Toward the study of combustion instabilities in aircraft engines
|
|
- Allen Freeman
- 5 years ago
- Views:
Transcription
1 Toward the study of combustion instabilities in aircraft engines A. Chatelier The prediction of combustion instabilities is crucial for aircraft engines manufacturers to ensure the reliability and the life span of their engines. To study these phenomena, numerical simulations are more and more often performed. They offer advantages: a wide range of parameters is available, they are safe and inexpensive compared to experiments. This work presents a first step in the general study of combustion instabilities. Firstly, several tools for acoustics and combustion are validated. The tests are performed on simple cases which have reference solutions from analytical resolution or measurements, such as one dimensional flames for combustion and rectangular duct for acoustics. The parameters of these cases are close to the more complex ones. Then, the dynamics of a laminar flame are studied. Finally, a laboratory scale configuration is explored. Kinetics for a methane-air mixture are validated as well as the TFLES model for combustion. The AVSP code for acoustics correctly determines the eigenmodes of a simple configuration. The TFLES model on a pulsed flame has a significant impact on its dynamics, depending on the thickening factor and the model used. An unstable mode is found for reacting RANS computations of a laboratory scale configuration. An unstable mode is also predicted by AVSP computations with an active flame model. Even if further work is needed to develop these tools, the first results indicate that they can quickly yield predictions about combustion instabilities. I. Introduction The environmental legislation becoming more and more strict, aircraft engine manufacturers developed new concepts of less polluting combustion chambers with a focus on NOx emissions. However, these new chambers, designed to operate at low equivalence ratio, may experience combustion instabilities that can damage or, even worse, destroy them. Trying to predict numerically these instabilities is thus a crucial challenge for the developers of CFD codes. The first scientist to study combustion instabilities, i.e. the interaction between acoustics and combustion, was Lord Rayleigh with his singing flame. 1 He defined what is called nowadays the Rayleigh criterion in combustion. Another major researcher to focus on combustion instabilities was Crocco, 2 during the 1950s. In the 1990s, McManus et al. 3 made a review of the progress on the control of these instabilities. Recently, researchers focused on ways to model flames for acoustic purposes, for instance Giauque 4 on the Flame Transfer Function (FTF). With the currently available computation resources, unsteady numerical simulations of reacting flows, on technical scale configurations can not resolve the entirety of flame characteristic lengths and their interactions with turbulence, due to computational costs. For the conditions usually encountered inside combustion chambers, the heat thickness of the flame front is very small (usually below one millimetre) and the reaction zone even smaller (one tenth of the heat thickness). Standard LES (Large Eddy Simulation) meshes, designed initially to resolve the most energetic turbulent structures, have mostly a characteristic mesh size larger than the heat thickness. Thus they are unable to resolve properly the instantaneous flame front. This is why numerous models, such as TFLES (Thickened Flame for LES) have been developed to take into account the limits of standard LES meshes in combustion. However, these models have been designed to recover integral properties of the flame, for example its consumption speed, but that can be done to the detriment of its dynamics. Unfortunately, this dynamics has a key role in the combustion instabilities, in which the heat Intern at DEFA (Department of Fundamental and Applied Energetics), ONERA, 29 avenue de la Division Leclerc, 92322, Chatillon, France and Double-Degree Student at KTH, Stockolm, Sweden and ENSTA ParisTech, Palaiseau, France 1 of 17
2 release coming from the combustion is coupled with the acoustics of the chamber: a poor resolution of the flame dynamics modifies this coupling and leads to an incorrect prediction of instabilities, if they exist. The structure of this article is the following: the numerical tools are presented, then the combustion and acoustic models are validated. The interaction between acoustics and combustion is studied on a pulsed laminar flame. Finally, the combustion instabilities are explored on a laboratory scale configuration, using the tools validated previously. II. Numerical Tools II.A. II.A.1. CEDRE Main characteristics of CEDRE CEDRE is a research code developed at ONERA and designed to solve multi-physics problems for energetics and propulsion, with several solvers. The solver used in this work is CHARME, which solves the Navier- Stokes equations for multispecies reacting flows. Using Favre s decomposition (for more details, see the appendix A on the filtering of Navier-Stokes equations) and the Einstein s notations, the filtered or averaged Navier-Stokes equations are recast in RANS and LES models as Eqs (1), (2), (3) and (4): ρẽ t t ρũ i t ρỹk t + ρũ iỹk x i (ρũ i ) = 0, (1) x i + J k,i ω k = ( ) x i x ρũ i Y k = ( ) J t/sgs k,i, (2) i x i ) = ( ) τ t/sgs ij, (3) x i ρ t + + (ρũ i ũ j ) + p τ ij = ( x j x i x j x ρũ i u j j + (ρũ i ẽ t ) [(τ pδ ij )u i ] + J e t,i = [ x i x j x i x ρũ i e t + ( τ ij ) ] δ p ij u i j (4) = ( ) Je t/sgs x t + τ t/sgs ij ũ i. (5) j The terms on the right hand side (RHS) need to be closed, i.e. modelled, in order to solve these equations. The closure uses classical analogies (e.g. Chaouat 5 ). The first term that has to be closed is the RHS of the momentum equation, Eq. (3). In RANS, an analogy is made between the Reynolds tensor τij t and the viscous strain tensor τ ij, as shown in Eq. (6) in BVM (Boussinesq Viscosity Model) approach: ( τij t ρũ i u j = µ ũi t + ũ j 2 ) ũ k δ ij 2 x j x i 3 x k 3 ρkδ ij. (6) In the present work, turbulent viscosity µ t is modelled with the two-equation k L model: In LES, the approach is similar and uses a classical Smagorinsky model: ( τ sgs ij ρũ i u j = µ ũi sgs + ũ j 2 ũ k x j x i 3 µ t = ρl k. (7) x k δ ij ). (8) The subgrid scale (sgs) viscosity µ sgs is defined in Eq. (9) with C s being a constant of the model (usually 0.1 < C s < 0.2), the size of the spatial filter and S the resolved velocity gradient tensor: µ sgs = ρ(c s ) 2 2 S : S. (9) The closure of the specie equation (Eq. (2)) and total energy equation (Eq. (4)) is made on the assumption that the turbulent flux follows a gradient law as an analogy with the Fick and Fourier laws, as in Eqs (10) and (11). It is valid in both RANS and LES. This model requires the introduction of two non-dimensional numbers, the turbulent Schmidt number Sc t/sgs,k and the turbulent Prandtl number P r t/sgs,k : J t/sgs k,i = D t/sgs Ỹk k x i with D t/sgs k = µ t/sgs /(ρsc t/sgs,k ), (10) 2 of 17
3 II.A.2. Je t/sgs t = q t/sgs n e i + k=1 h k J t/sgs k,i = λ t/sgs k Numerical methods in CHARME T n e + x i k=1 h k J t/sgs k,i with λ t/sgs k = µ t/sgs c p /P r t/sgs,k. (11) CHARME is based on a cell-centered finite volume methodology for unstructured meshes of general polyhedra. The interpolation is based on MUSCL (Monotone Upstream Schemes for Conservation Laws) schemes. Both explicit and implicit methods of time integration are implemented in CHARME. The method mainly used in this work is the backward Euler implicit method. II.B. II.B.1. AVSP Main characteristics of AVSP The AVSP code, from CERFACS, is designed to solve linear acoustic problems on 2D or 3D configurations. It solves the Helmholtz equations assuming harmonic waves (Eqs (12), (13), (14) and (15) from Poinsot and Veynante 6 for non-reacting flows) where p ω is the harmonic pressure variation, u ω the harmonic velocity variation, ρ ω the harmonic density variation, k = ω/c 0 the wave number, ω the pulsation and c 0 the speed of sound: iωp ω + ρ 0 c 2 0. u ω = 0, (12) ρ 0 iω u ω + p ω = 0, (13) which leads to: The quantities p 1, u 1 and ρ 1 are defined in Eq. (16): p ω = c 2 0ρ ω, (14) 2 p ω + k 2 p ω = 0. (15) p 1 = R(p ω e iωt ), u 1 = R( u ω e iωt ), ρ 1 = R(ρ ω e iωt ). (16) These quantities are used to define the acoustic impedance Z and the reflection coefficient R in Eq. (17): II.B.2. Numerical methods in AVSP Z = 1 ρ 0 c 0 p 1 u 1, R = Z + 1 Z 1. (17) AVSP deals with triangles in 2D and tetrahedra in 3D. Two options are available: multiple modes, where modes are computed with increasing frequencies or single mode, where the target frequency is specified as an input. In multiple modes option, two boundary conditions are available. Either the acoustic pressure variation p 1 is null, which corresponds to an open end or the acoustic velocity variation u 1 is zero, which is equivalent to a wall. In single mode option, more boundary conditions are available for specific purposes, for instance a given impedance can be set as a constant or from experimental data. Also, an active flame model is available for the study of combustion instabilities. III. Validation of combustion models This section aims to compare computational and experimental results on one-dimensional laminar premixed flames. It is possible to express the laminar flame speed S L as Eq. (18). The derivation of this formula in shown in the appendix C. This expression is valid for the fuel, as well as for the oxidizer: S L = ω k ρ u (Y kb Y ku ). (18) 3 of 17
4 Values of ρ u, Y kb and Y ku are computed once, knowing the composition of both the burned and unburned gases. ω k is the integral on the domain of the source term ω k. The finer the grid along the x direction, the better the term ω k is approximated. The experimental results for this part of the study come from the experiments of Vagelopoulos 7 for the methane-air premixed flame. The computations are done using CEDRE. The computational domain is a two-dimensional rectangle, m-wide in the x direction and 1 m-high in the y direction. Each cell is m-wide in the x direction and 1 m-high in the y direction. The domain is therefore composed of 1000 cells. The initial position of the flame front is at x = m with the origin on the left boundary. On the left side of the flame, the domain is initialised with fresh gases while the right side is initialised with burnt gases. The equivalence ratio of the computation is φ = The kinetics used in this section is CM2-CH4-2R. 8 It involves six species and two reactions, shown in Eq. (19). Table [1] sums up the composition and state in fresh and burned gases for this equivalence ratio. { CH O 2 CO + 2H 2 O CO O (19) 2 CO 2 Quantity Value in fresh gases Value in burned gases Y CH4 [ ] Y O2 [ ] Y H2O [ ] Y CO2 [ ] Y CO [ ] Y N2 [ ] T [K] P [Pa] Table 1: Specie composition, temperature and pressure in fresh and burned gases III.A. Methane-air chemical scheme Flame speed (m/s) Figure 1: 0.15 Experimental data GRI-Mech 3.0 2S_CH4_BFER with PEA correction CM2_CH4_2R Equivalence ratio (-) Comparison of the laminar flame speed between experiments and computations. Fig. 1 shows a comparison between three chemical schemes (GRI-Mech 3.0, 9 2S CH4 BFER with PEA correction described by Franzelli et al. 10 and CM2 CH4 2R) and experimental results from Vagelopoulos 4 of 17
5 and Egolfopoulos 7 for an equivalence ratio between 0.7 and 1.4. GRI-Mech is the reference detailed chemical scheme for natural gases, including methane. One can notice that the experimental data and the GRI-Mech 3.0 do not give exactly the same flame speed. The discrepancy may have two origins: firstly, GRI-Mech 3.0 is a detailed scheme but does not describe perfectly all the elementary reactions that occur in reality; secondly, the experimental results have uncertainties of measurement. Despite this, the evolution is well reproduced. 2S CH4 BFER and CM2 CH4 2R are two reduced chemistry schemes that are very similar. CM2 CH4 2R leads to poor results for rich mixture (i.e. φ > 1.1). Franzelli et al. proposed a PEA (Pre Exponential factor Adjustment) correction to fit 2S CH4 BFER with GRI-Mech 3.0 for rich mixtures. This is confirmed by the results shown in Fig. 1. The reduced schemes reproduce well the flame speed for lean mixture (φ < 1). The scheme that is used in this work is CM2 CH4 2R, dedicated to lean combustion. III.B. III.B.1. TFLES model Presentation of the model TFLES stands for Thickened Flame for LES. For common LES simulations, the size of the flame front is usually smaller than the smallest cell size. It means that either the mesh can be refined or the flame front can be thickened. The first option is most of the time not realistic due to computational costs. The second option is the principle of TFLES. The laminar speed S L and thickness δ L of a flame can be expressed as shown in Eq. (20): S L D u th A, δ L Du th S L. (20) Assuming that the heat diffusion coefficient in Eq. (20) can be considered as D th regardless the position in space and that Le k = 1, so that the diffusion coefficients can be written as a unique coefficient D, the previous equation becomes Eq. (21): S L DA, δ L D S L. (21) In order to thicken the flame front δ L, the coefficient diffusion D is multiplied by a factor F. Meanwhile, the laminar flame speed has to remain constant, so the preexponential factor A is divided by F, as in Eq. (22): S L F D A F = DA S L = S L, δ L F D S δ L = F δ L. (22) L The next step in the use of the TLFES model is the application of the thickening factor. One has two possibilities: either the TFLES model is applied everywhere in the domain, or the model is applied specifically in the region where the reaction occurs (i.e. the flame zone). The first case is the simplest one and is a reasonable first step when studying the effect of the TFLES model on a flame. It will be called here standard TFLES. The second case requires a criterion that determines where to apply the model. This one will be called dynamic TFLES. The idea behind dynamic TFLES is that the thickening factor is computed in each cell of the domain. In order to apply the model only in the flame, one requires a flame sensor. Eq. (23) is an example of a pseudo reaction rate that has its maximum in the flame. Subscripts F and O stand for Fuel and Oxidizer. Γ is a constant and is usually set to 0.5. T a is the activation temperature for the considered reaction. Eq. (24) is the expression of the sensor α. Ω 0 is the pseudo-reaction rate computed in the laminar 1D flame. C α is a parameter that determines the thickness of the thickened zone. Finally, Eq. (25) gives the expression of the thickening factor as a function of the maximum thickening factor F max and the sensor α: ( Ω = Y F Y O exp Γ T ) a, (23) T ( ) Ω α = tanh C α, (24) Ω 0 F (α) = 1 + (F max 1)α. (25) III.B.2. Validation of the TFLES model on a 1D flame In order to confirm that the TFLES model gives a correct flame speed, it is tested on the same case as in section III.A with a coarser mesh (the cells are four times larger). The results are shown in Fig. 2: 5 of 17
6 Laminar flame speed (m/s) Flame thickness (m) F = 1 F = 2 F = 3 F = Iteration F = 1 F = 2 F = 3 F = Iteration Figure 2: Convergence of the computation of S L and δ L for several thickening factors F Fig. 2 shows that the convergence of the flame speed depends on the thickening factor. Computations with F = 1 and F = 2 present oscillations which show that F is not large enough to make the model converge on this given mesh. Also, the flame thickness is not correct for these two thickening factors. The results for F = 3 are much closer to what was expected. Finally for F = 4, the flame speed S L0 is 40.4 cm/s while the flame thickness δ L0 is approximately m. Therefore the flame speed is underestimated (by less than 3%) and the flame thickness is overestimated (by 4%) compared to the results without the TFLES model. This margin of error remains acceptable. IV. Validation of acoustic models IV.A. Boundary conditions and acoustics The study of the interaction between acoustics and combustion implies that the acoustic waves are properly propagated in the computational domain. The boundary conditions must be able to reproduce the correct acoustic impedance of the configuration. Physical domain length Physical inlet Inlet Outlet Physical outlet Computational domain length Figure 3: Acoustic impedance modelling for an ideal one-dimensional duct Fig. 3 presents an ideal case for acoustic modelling. Two boundary conditions must be set to reproduce the correct acoustic behaviour: the inlet and the outlet. The computational domain is usually different from the experimental domain in order to decrease the computational cost. If the pressure is fixed at a boundary, then p 1 = 0 which leads to Z = 0 and R = 1. If the velocity is fixed, then u 1 = 0 so that Z = and R = 1. These two conditions will be called standard boundary conditions. If the physical boundaries can be considered very far from the computational ones, non-reflecting boundary conditions can be used. This means that acoustic waves travel without reflection, so that R = 0 or R = (and Z = 1 or Z = 1). The previous conditions are ideal. For some configurations, it is not possible to use them because the actual boundary does not fit any of the previous cases. In order to set properly the boundary conditions, the acoustic impedance must be measured during experiments. 6 of 17
7 IV.B. Behaviour of CEDRE boundary conditions Pulsed pressure at the outlet boundary In this case, the pressure of the outlet boundary of a one dimensional domain is pulsed with a single 4 Pa-amplitude pulse compared to the initial value ( Pa). Symmetry Inlet 2 cm Probe position Outlet Symmetry Figure 4: Computational domain for the test of boundary conditions Pressure (Pa) Pressure (Pa) Time (s) Figure 5: Pressure over time at a given position for the standard boundary condition at 20% of the total length of the domain Time (s) Figure 6: Pressure over time at a given position for the non-reflecting boundary condition at 20% of the total length of the domain. Fig. 5 shows the behaviour of the standard boundary condition on the pressure. The position is roughly at 20% of the total length of the domain. The first peak in pressure corresponds to the normal propagation from the right of the pressure wave. However the second peak short after shows that the pressure peak was reflected entirely by the boundary condition. The last two pressure drops come from the reflection of the pulse on the outlet boundary. They are inverted because the pressure is fixed at the outlet boundary at Pa. Fig. 6 shows what happens with non-reflecting boundary conditions. The position is the same as in Fig. 5. In this case, the peak reaches the limit but is not reflected and leaves the domain. This validates the behaviour of the non-reflecting boundary condition. IV.C. AVSP validation In order to validate the results of AVSP, the code is tested on simple configurations that can be solved analytically. The configuration considered here is a 2D rectangular duct closed on the left and opened on the right and enclosed by walls above and below, as shown in Fig. 7. The resolution of the Helmholtz equation for harmonic waves gives the eigenfrequencies for this configuration in Eq. (26), where m and n are integers: (2m ) ( n ) 2. f m,n = c 0 + (26) 4L 2h For L = 1 m and h = 0.5 m, the results of AVSP and the analytical solution are shown in Table 2. 7 of 17
8 Wall, u 1 = 0 Closed end, u 1 = 0 h Opened end, p 1 = 0 Figure 7: Wall, u 1 = 0 L Reference configuration for the validation of the AVSP code Mode (m,n) Analytical solution (Hz) AVSP solution (Hz) Error (in %) (0,0) (1,0) (2,0) (3,0) (0,1) (1,1) (2,1) (3,1) (0,2) (1,2) Table 2: Comparison of eigenfrequencies between analytical and AVSP solutions The results given by AVSP on this configuration are globally close from the analytical results given by linear acoustics. One can notice that the highest the frequency, the largest the error tends to be. Also, when the eigenfrequency is double (for example the modes (2,0) and (1,1) have the same eigenfrequency), AVSP finds two different frequencies. Figure 8: (2,1)) Pressure variations in arbitrary units for f = Hz (mode (1,0)) and f = Hz (mode Fig. 8 shows two examples of results from AVSP on the absolute pressure variation (the scale is arbitrary). The first one corresponds to the so-called three-quarter wave mode (because f = 3/4 c 0 /L) and is purely longitudinal. The second one shows a combination of the second longitudinal and the first transversal modes (therefore begin the (2,1) mode). 8 of 17
9 Figure 9: Pressure variations in arbitrary units for f = Hz and f = Hz (modes (2,0) and (1,1)) Fig. 9 shows what happens when an eigenfrequency is double. The figure on the left should show a purely longitudinal mode (similar to the Fig. 8 on the right without the transversal component). The second one should show a combination of the first longitudinal and first transversal modes. This result can be explained with the method of resolution of AVSP. The code solves the Helmholtz equation for eigenfrequencies. When an eigenfrequency is double, the subspace associated to this frequency is of dimension two. The code is not able to distinguish between the two physical eigenmodes that exist for this frequency. Therefore the result is a linear combination of the two pressure modes associated with the same frequency. V.A. V. Interaction between acoustics and combustion n τ model in AVSP Combustion instabilities exist because of the interaction between acoustics and combustion. A first step in their analysis is the study of eigenmodes of the configuration. To do so, the code AVSP from CERFACS has been used. Firstly, the eigenmodes of a configuration are computed with no combustion. The assumptions in this computation are those of linear acoustics. Secondly, the eigenmodes are computed again with an active flame model. The n τ model introduced by Crocco 11 is used to perform this. As its name suggests, it requires two parameters, n and τ. n is called the interaction index, and describes the intensity of a change in velocity on the heat release. τ represents the time delay between variations in velocity and heat release. This model links the linearised heat release term Ω 1 T with the linearised acoustic velocity u 1 (Eq. (16)) at a certain position a, as shown in Eq. (27): Ω 1 T = S 1ρ 1 c 2 1 γ 1 nu 1(a, t τ). (27) One has several ways to determine these parameters. The first one is the analysis of experimental data. The second one is to use values within an order of magnitude. The complete description is available in the technical report of CERFACS written by Ghani. 12 The introduction of the n τ model modifies the eigenfrequencies. Their real part is shifted by the presence of the flame and their imaginary part can become positive, which can lead to a combustion instability. V.B. Pulsed flame The computations of this section are derived from the EM2C experiment 13 for the original experimental setup and results and from the work of Auzillon 14 for the input data and computation results. The problem is considered two-dimensional and axisymmetrical. The goal of this section is to study the dynamics of a laminar flame when pulsed by a sinusoidal velocity (Eq. (28)) for several TFLES factors: v(t) = v + 2 v cos(2πft). (28) The computational domain is 400 mm-wide and 600 mm-high and is composed of two physical domains. The first one represents the burner, with 30mm height and 11mm radius. The second one is the atmosphere 9 of 17
10 with a 1mm-wide and 30mm-high wall to separate it from the burner. The two domains are connected at the exit of the burner. Outlet 400 mm Axis of symmetry 11 mm CH4/air mixture 30 mm 600 mm Inlet (a) Computational domain and boundary conditions. (b) Close up view of the mesh Figure 10: Pulsed flame configuration Fig. 10a shows the computational domain and the boundary conditions and Fig. 10b shows a part of the mesh. The data used to perform the computations are shown in Tables [1] and [3]. The kinetics used in this section is CM2-CH4-2R (Eq. (19)). Input data Symbol Value(s) Equivalence ratio φ 1.05 Diameter of the burner (mm) d 22 Mean injection velocity (m/s) v 0.97 Variable injection velocity (m/s) v 0.19 Pulsation frequency (Hz) f 62.5 Pulsation period (s) T Thickening factors for standard TFLES F Mesh sizes (number of cells) for each F 240k 180k 150k Characteristic length of the smallest cells l c for each F (in the flame region) (m) Table 3: Data for the pulsed flame. 10 of 17
11 F = 2.5 F = 4 F = 5.7 DNS data Y (m) Y (m) Y (m) Y (m) X (m) X (m) X (m) X (m) t = t i t = t i + T 0 /4 t = t i + T 0 /2 t = t i + 3T 0 /4 Figure 11: Flame front (c = 0.5) for three thickening factors over a cycle of velocity excitation Fig. 11 shows the comparison between the DNS data from Auzillon 14 and the standard TFLES model for three thickening factors. The thickening factor has a crucial impact on the flame front, as it was expected. The curvature of the flame front is smoothed when F increases. Over a cycle of excitation, the thickening factor also has an effect on the phase shift between the thickened flame and the DNS results: the flame front is higher along the y axis when F increases for the first two instants when the speed increases while it is lower or at the same level when the speed decreases. VI. Application to a laboratory scale configuration: PRECCINSTA VI.A. The PRECCINSTA configuration The PRECCINSTA configuration, shown in Fig. 12a, is an experiment designed by DLR in Stuttgart with a Turbomeca helicopter injector (presented for instance by Lartigue and Meier 15 ). It is composed of a plenum (where the air is injected), an injector (where the fuel is injected) and a chamber (where the combustion occurs). The plenum is used as a tranquillisation chamber. The injector is a swirler and was initially designed by Turbomeca for helicopter combustors. This injector was then modified in order to inject methane. It was chosen because it was likely to create combustion instabilities in certain regimes. In this paper, the studied case is simplified. The fuel and the air are considered premixed and injected in the plenum. 11 of 17
12 (a) The PRECCINSTA configuration. (b) A 2D view of the computational domain for the PRECCINSTA configuration. Figure 12 The computational domain used in this section is shown in Figure 12b. The large cone is the atmosphere. Input data Symbol Value(s) Equivalence ratio φ 0.75 Diameter of the inlet (mm) d 25 Flow rate (kg/m 2 /s) ρu Mesh size for RANS (number of cells) 970,000 Characteristic length of the smallest cells (m) l c Turbulence model for RANS k l Turbulence rate in RANS T u Characteristic length of turbulence (m) in RANS l 0.02 Table 4: Data for the PRECCINSTA case Quantity Value in fresh gases Value in burned gases Y CH4 [ ] Y O2 [ ] Y H2O [ ] Y CO2 [ ] Y CO [ ] Y N2 [ ] T [K] P [Pa] Table 5: Specie composition, temperature and pressure in fresh and burned gases VI.B. RANS computations The data used to initialize the RANS computations are shown in Table of 17
13 Figure 13: Axial velocity field and streamlines at z = 0 from RANS computations. Fig. 13 shows the axial velocity field at z = 0 for RANS computations. The effect of the swirler on the velocity is seen. A large recirculation zone appears downstream of the injector, and two smaller are located at the bottom of the chamber. VI.C. RANS computations with combustion The data used to perform the RANS computations with combustion are shown in Tables 4 and 5. The kinetics used are the same as previously. The initial state used is the result of the RANS computations without combustion. Figure 14: Temperature field at z = 0 from RANS computations with combustion before the instability appears Fig. 14 shows the temperature field given by the RANS computations with combustion. The flame was initialised just after the exit of the injector, by changing the composition from fresh gases to burnt gases. 13 of 17
14 Figure 15: Pressure field at z = 0 for two instants, showing the coupling Fig. 15 shows what happens during the computations: an acoustic coupling appears between the plenum and the chamber Plenum pressure Chamber pressure 1.8e+06 FFT result e e e+06 Pressure (Pa) Amplitude (a.u.) 1e Time (s) Frequency (Hz) Figure 16: Pressure in the chamber and the plenum over time and Fast Fourier Transform of one signal Fig. 16 shows the pressure in the plenum and the chamber over time and the Fast Fourier Transform of one signal. One can note that they are almost in opposition of phase. When the pressure reaches the maximum in the plenum, it is close to the minimum in the chamber. VI.D. AVSP computations The geometry of the 3D PRECCINSTA configuration in the plane z = 0 has been simplified as a 2D domain in order to compute the eigenmodes in AVSP. The modes that are found are therefore only the ones in this plane. This approach is consistent with the mode observed in the RANS computations with combustion. Figure 17: Absolute pressure variations for the first mode f = Hz without active flame The first mode observed in AVSP without active flame is shown in Fig. 17 and is mainly longitudinal. 14 of 17
15 The active flame model is then enabled. It is initialized as a compact flame (a field of n and τ), 0.01 m wide along x and as large as the domain along y. It is located directly after the injector. The resolution of the eigenproblem with an active flame is made on a single point, i.e. a target frequency is set before the computations are run. In this section, the computations are focused on the first mode (initially f 0 = Hz). The period associated to this frequency is T 0 = s. The n parameter is set at 10 8 J/m. This value is chosen based on the typical order of magnitude of n. 12 Two values of τ have been chosen: τ 1 = s and τ 2 = s. The first value is below T 0 /2 while the second one is between T 0 /2 and T 0. According to Poinsot and Veynante, 6 one might expect that computations with τ 1 converge on a stable mode, while the ones with τ 2 on an unstable mode. An effect that has to be studied is the choice of the target frequency (f t ). Two values have been studied: f t,1 = 200 Hz and f t,2 = 400 Hz. Case τ (s) f t (Hz) Resulting eigenfrequency (Hz) Table 6: 1 τ 1 f t, i τ 1 f t, i τ 2 f t, i τ 2 f t, i Results of the AVSP computations with active flame Table 6 presents the results from the AVSP computations with active flame. One can notice that in all cases, the real part of the eigenfrequency is identical. However the imaginary part differs. In cases 1 and 3, the target frequency is set far from the original frequency f 0. The algorithm in these cases correctly predicts the real part, but the imaginary part is not converged with the default parameters of the algorithm. When setting the target frequency much closer to f 0, the imaginary part is different. In case 1, it is negative which denotes a stable mode while in case 4 it is positive thus unstable. The FFT of the signals of Fig. 16 gives a maximum amplitude around 400 Hz. However the signals have a length of 0.01 seconds, which is not sufficient to have a good frequency resolution, given that the phenomenon occurs at arout 400 Hz. Further computations must be done to improve the accuracy of the results. VII. Conclusion The present work shows that several tools are available to predict combustion instabilities. The use of an acoustic code on a simple geometry can quickly highlight the unstable behaviour of a configuration. Meanwhile, this unstable mode can be found with computations on the complete configuration. However, one must be careful with the model used for combustion. For instance, the TFLES model has a crucial impact on the dynamics of the flame. A. Filtering of Navier-Stokes equations Appendix The Navier-Stokes equations are shown in Eqs (29), (30), (31) and (32): ρ DY k Dt = ρy k t ρ DV Dt = ρv t ρ + (ρv ) = 0, (29) t + (ρy k V ) = J k + ω k, (30) + (ρv V ) = (τ pi) + S m, (31) ρ De t Dt = ρe t + (ρe t V ) = J et + [(τ pi)v ] + S et. (32) t In Eq. (30), ω k is the production (if positive) or consumption (if negative) rate of specie k. J k is the diffusion flux of specie k and can be written as J k = ρv k Y k with V k being the diffusion velocity of specie 15 of 17
16 k. n e k=1 J k = 0 and n e k=1 ω k = 0 have to be verified in order to ensure the global mass conservation. J k is based on the phenomenologic diffusion law of Fick, which is a common approximation, with a diffusion coefficient D k, as shown in Eq. (33): J k = ρd k Y k. (33) In Eq. (31), τ is the viscous strain tensor and S m is the source term that can be expressed as S m = ρ n e k=1 Y kf k +..., where f k are the volume forces. In Eq. (32), J et = q + n e k=1 h kj k is the energy flux, where h k is the enthalpy of each specie and q is the heat flux, S et = S m (V + V k ) + S rad +... is the source term, where S rad is the source term from radiation. It is assumed that any quantity ϕ can be decomposed as the sum of a main part ϕ and a fluctuation ϕ. In RANS computations, ϕ is the mean part of the quantity, while it is the resolved part of the quantity in LES. In RANS, Eq. (34) are verified, but not in LES, where new correlations appear. The Favre decomposition (or filtering) is an extension of the first decomposition, where ϕ = ϕ+ϕ, with ϕ and ϕ defined in Eq. (35): ϕ = ϕ, ϕ = 0, (34) ϕ = ρϕ ρ and ρϕ = 0, (35) For example: ρu i u j = ρũ i ũ j + ρũ i u j. (36) It is designed to take into account the variability of the density in compressible flows. Using this filtering, there will be no correlated part in the equation for conversation of mass. B. Classical quantities in combustion The Prandtl number P r k, which compares the kinematic viscosity with the heat conductivity, is defined for each specie k in Eq. (37): P r k = µ kc pk λ k = ν D th. (37) The Schmidt number Sc k, which compares the kinematic viscosity with the specie diffusion, is defined for each specie k in Eq. (38): Sc k = µ k ρ k D k = ν D k. (38) The equivalence ratio φ is defined for a generic reaction between a fuel F and an oxidant O, resulting in a product P, as shown in Eq. (39): F + so (1 + s)p. (39) The parameter s depends on the reaction and is called the mass stoichiometric coefficient. The general definition of φ is given in Eq. (40). For a premixed laminar flame, the fuel and the oxidant are mixed and injected at the same mass rate ṁ. Therefore the equivalence ratio has a simplified expression in Eq. (41): φ = s Y F ṁ Y O ṁ = s Y ( F YF = Y O Y O φ = s ṁf, (40) ṁ O ) ( ) YF / Y O. (41) st The subscript st denotes the stoichiometry. The progress variable c is introduced to be able to locate a position in the flame in an unambiguous way. For instance, c for premixed methane-air flame can be defined as in Eq. (42): Y CO + Y CO2 c =. (42) (Y CO + Y CO2 ) eq(φ) The composition of the mixture when the reaction reaches equilibrium (subscript eq) depends on the equivalence ratio φ defined in Eq. (41). The definition of c is based on the necessity to have an increasing and monotonous variable in the flame to describe its structure. 16 of 17
17 C. Derivation of the laminar flame speed for a simplified 1D-flame For a simplified 1D-flame (incompressible, stationary and using the Fick law of diffusion), Eq. (30) becomes Eq. (43). Under these assumptions, the flow rate ṁ is constant and can be expressed as ρu or ρ u S L, where u is the subscript for unburned gases (while b corresponds to burned gases): ρu dy k dx = d dx ( ρd k Y k dx ) + ω k. (43) Eq. (43) can be recast as Eq. (44) and integrated in the x direction as in Eq. (45): + ρ u S L dy k = ρ u S L dy k dx = dj k dx + ω k, (44) + + dj k + ω k dx. (45) The diffusion flow rate J k is zero far from the flame (since it is a premixed flame). Therefore the first term of the right hand side is equal to zero. By writing the second term as ω k, Eq. (18) is obtained. Acknowledgements The author would like to thank Dr. Nicolas Bertier and Aymeric Boucher for their patience, support and invaluable scientific input, the DEFA department and ONERA for their welcome and Dr. David Eller for his supervision. References 1 John William Strutt 3rd Baron Rayleigh. The Theory of Sound Luigi Crocco and Sin-I Cheng. Theory of combustion instability in liquid propellant rocket motors. Cambridge Univ Press, K. R. McManus, T. Poinsot, and S. M. Candel. A review of active control of combustion instabilities. Progress in Energy and Combustion Science, 19(1):1 29, Alexis Giauque. Flame transfer function and disturbance energies in gaseous reacting flows. PhD thesis, March Bruno Chaouat. Modelisation et simulation des ecoulements turbulents dans les propulseurs a propergol solide [Modelling and simulation of turbulent flows in solid propergol propulsors]. PhD thesis, Thierry Poinsot and Denis Veynante. Theoretical and Numerical Combustion. R.T. Edwards Inc., Christine M. Vagelopoulos and Fokion N. Egolfopoulos. Direct experimental determination of laminar flame speeds. Symposium (International) on Combustion, 27(1): , Twenty-Seventh Sysposium (International) on Combustion Volume One. 8 G. Boudier. Methane/air flame with 2-step chemistry: 2s-ch4-cm2. Technical report, CERFACS technical report, F. Frenklach, H. Wang, C.-L. Yu, M. Goldenberg, C. Bowman, R. Hanson, D. Davidson, E. Chang, G. Smith, D. Golden, W. Gardiner, and V. Lissianski. GRI-Mech [Online; accessed 27-August- 2014]. 10 Benedetta Franzelli, Eleonore Riber, Laurent Y. M. Gicquel, and Thierry Poinsot. Large eddy simulation of combustion instabilities in a lean partially premixed swirled flame. Combustion and Flame, 159(2): , February L. Crocco. Aspects of combustion stability in liquid propellant rocket motors part i: Fundamentals. low frequency instability with monopropellants. Journal of the American Rocket Society, 21(6): , A. Ghani. Avsp active flame: Tool descriptions. Technical report, CERFACS internal technical report. 13 Sébastien Ducruix, Daniel Durox, and Sébastien Candel. Theoretical and experimental determinations of the transfer function of a laminar premixed flame. Proceedings of the Combustion Institute, 28(1): , Pierre Auzillon. Modeling chemical flame structure and combustion dynamics in large eddy simulation. PhD thesis, G. Lartigue and U. Meier. Experimental and numerical investigation of self-excited combustion oscillations in a scaled gas turbine combustor. Applied Thermal Engineering, pages , of 17
XXXVIII Meeting of the Italian Section of the Combustion Institute
Coupling a Helmholtz solver with a Distributed Flame Transfer Function (DFTF) to study combustion instability of a longitudinal combustor equipped with a full-scale burner D. Laera*, S.M. Camporeale* davide.laera@poliba.it
More informationExperimental analysis and large eddy simulation to determine the response of non premixed flame submitted to acoustic forcing
Experimental analysis and large eddy simulation to determine the response of non premixed flame submitted to acoustic forcing B. Varoquié, J.P. Légier, F. Lacas, D. Veynante and T. Poinsot Laboratoire
More informationADVANCED DES SIMULATIONS OF OXY-GAS BURNER LOCATED INTO MODEL OF REAL MELTING CHAMBER
ADVANCED DES SIMULATIONS OF OXY-GAS BURNER LOCATED INTO MODEL OF REAL MELTING CHAMBER Ing. Vojtech Betak Ph.D. Aerospace Research and Test Establishment Department of Engines Prague, Czech Republic Abstract
More informationThermoacoustic Instabilities Research
Chapter 3 Thermoacoustic Instabilities Research In this chapter, relevant literature survey of thermoacoustic instabilities research is included. An introduction to the phenomena of thermoacoustic instability
More informationA Novel FEM Method for Predicting Thermoacoustic Combustion Instability
Excerpt from the Proceedings of the COMSOL Conference 009 Milan A Novel FEM Method for Predicting Thermoacoustic Combustion Instability G. Campa *, S.M. Camporeale Politecnico di Bari * campa@imedado.poliba.it,
More informationBest Practice Guidelines for Combustion Modeling. Raphael David A. Bacchi, ESSS
Best Practice Guidelines for Combustion Modeling Raphael David A. Bacchi, ESSS PRESENTATION TOPICS Introduction; Combustion Phenomenology; Combustion Modeling; Reaction Mechanism; Radiation; Case Studies;
More informationStudies of mean and unsteady flow in a swirled combustor using experiments, acoustic analysis and Large Eddy Simulations
Studies of mean and unsteady flow in a swirled combustor using experiments, acoustic analysis and Large Eddy Simulations S. Roux a,, G. Lartigue a, T. Poinsot a,b, U. Meier c and C. Bérat d a CERFACS,
More informationHEAT TRANSFER IN A RECIRCULATION ZONE AT STEADY-STATE AND OSCILLATING CONDITIONS - THE BACK FACING STEP TEST CASE
HEAT TRANSFER IN A RECIRCULATION ZONE AT STEADY-STATE AND OSCILLATING CONDITIONS - THE BACK FACING STEP TEST CASE A.K. Pozarlik 1, D. Panara, J.B.W. Kok 1, T.H. van der Meer 1 1 Laboratory of Thermal Engineering,
More informationFlame / wall interaction and maximum wall heat fluxes in diffusion burners
Flame / wall interaction and maximum wall heat fluxes in diffusion burners de Lataillade A. 1, Dabireau F. 1, Cuenot B. 1 and Poinsot T. 1 2 June 5, 2002 1 CERFACS 42 Avenue Coriolis 31057 TOULOUSE CEDEX
More informationA Ghost-fluid method for large-eddy simulations of premixed combustion in complex geometries
Center for Turbulence Research Annual Research Briefs 2005 269 A Ghost-fluid method for large-eddy simulations of premixed combustion in complex geometries By V. Moureau, P. Minot, C. Bérat AND H. Pitsch
More informationPremixed MILD Combustion of Propane in a Cylindrical. Furnace with a Single Jet Burner: Combustion and. Emission Characteristics
Premixed MILD Combustion of Propane in a Cylindrical Furnace with a Single Jet Burner: Combustion and Emission Characteristics Kin-Pang Cheong a, c, Guochang Wang a, Jianchun Mi a*, Bo Wang a, Rong Zhu
More informationAPPLICATION OF HYBRID CFD/CAA TECHNIQUE FOR MODELING PRESSURE FLUCTUATIONS IN TRANSONIC FLOWS
TASK QUARTERLY Vol. 17, Nos 3 4, 2013, pp. 145 154 APPLICATION OF HYBRID CFD/CAA TECHNIQUE FOR MODELING PRESSURE FLUCTUATIONS IN TRANSONIC FLOWS SŁAWOMIR DYKAS, WŁODZIMIERZ WRÓBLEWSKI AND SEBASTIAN RULIK
More informationNumerical simulations of a massively separated reactive flow using a DDES approach for turbulence modelling
Numerical simulations of a massively separated reactive flow using a DDES approach for turbulence modelling Bruno Sainte-Rose, Nicolas Bertier, Sébastien Deck and Francis Dupoirieux Abstract Computations
More informationS. Kadowaki, S.H. Kim AND H. Pitsch. 1. Motivation and objectives
Center for Turbulence Research Annual Research Briefs 2005 325 The dynamics of premixed flames propagating in non-uniform velocity fields: Assessment of the significance of intrinsic instabilities in turbulent
More informationLEAN PREMIXED TURBULENT COMBUSTION MODELING USING FLAME TABULATED CHEMISTRY AND A PRESUMED PDF APPROACH
Seoul, Korea, 22-24 June 29 LEAN PREMIXED TURBULENT COMBUSTION MODELING USING FLAME TABULATED CHEMISTRY AND A PRESUMED PDF APPROACH Julien Savre, Nicolas Bertier, Daniel Gaffie Department of Fundamental
More informationNumerical Investigation of Ignition Delay in Methane-Air Mixtures using Conditional Moment Closure
21 st ICDERS July 23-27, 27 Poitiers, France Numerical Investigation of Ignition Delay in Methane-Air Mixtures using Conditional Moment Closure Ahmad S. El Sayed, Cécile B. Devaud Department of Mechanical
More informationDynamics of Lean Premixed Systems: Measurements for Large Eddy Simulation
Dynamics of Lean Premixed Systems: Measurements for Large Eddy Simulation D. Galley 1,2, A. Pubill Melsió 2, S. Ducruix 2, F. Lacas 2 and D. Veynante 2 Y. Sommerer 3 and T. Poinsot 3 1 SNECMA Moteurs,
More informationA G-equation formulation for large-eddy simulation of premixed turbulent combustion
Center for Turbulence Research Annual Research Briefs 2002 3 A G-equation formulation for large-eddy simulation of premixed turbulent combustion By H. Pitsch 1. Motivation and objectives Premixed turbulent
More informationTOPICAL PROBLEMS OF FLUID MECHANICS 97
TOPICAL PROBLEMS OF FLUID MECHANICS 97 DOI: http://dx.doi.org/10.14311/tpfm.2016.014 DESIGN OF COMBUSTION CHAMBER FOR FLAME FRONT VISUALISATION AND FIRST NUMERICAL SIMULATION J. Kouba, J. Novotný, J. Nožička
More informationCombustion basics... We are discussing gaseous combustion in a mixture of perfect gases containing N species indexed with k=1 to N:
Combustion basics... T. Poinsot poinsot@imft.fr Only the things you should know to understand the following courses Mainly elements of laminar flame theory for premixed and diffusion flames 1 Copyright
More informationTurbulent Boundary Layers & Turbulence Models. Lecture 09
Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects
More informationAnalysis of dynamic models for turbulent premixed combustion
Center for Turbulence Research Proceedings of the Summer Program 2012 387 Analysis of dynamic models for turbulent premixed combustion By D. Veynante, T. Schmitt, M. Boileau AND V. Moureau Very few attempts
More informationLecture 14. Turbulent Combustion. We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing.
Lecture 14 Turbulent Combustion 1 We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing. In a fluid flow, turbulence is characterized by fluctuations of
More informationThermal NO Predictions in Glass Furnaces: A Subgrid Scale Validation Study
Feb 12 th 2004 Thermal NO Predictions in Glass Furnaces: A Subgrid Scale Validation Study Padmabhushana R. Desam & Prof. Philip J. Smith CRSIM, University of Utah Salt lake city, UT-84112 18 th Annual
More informationSimulation of Turbulent Lifted Flames and their Transient Propagation
25 th ICDERS August 2-7th, 2015 Leeds, UK Simulation of Turbulent Lifted Flames and their Transient Propagation S. Ruan, Z. Chen, N. Swaminathan University of Cambridge Cambridge, UK 1 Introduction Turbulent
More informationLarge Eddy Simulations for the Flame Describing Function of a premixed turbulent swirling flame
Large Eddy Simulations for the Flame Describing Function of a premixed turbulent swirling flame Davide LAERA, and Aimee S. MORGANS Department of Mechanical Engineering, Imperial College London, London,
More informationTurbulence Modeling. Cuong Nguyen November 05, The incompressible Navier-Stokes equations in conservation form are u i x i
Turbulence Modeling Cuong Nguyen November 05, 2005 1 Incompressible Case 1.1 Reynolds-averaged Navier-Stokes equations The incompressible Navier-Stokes equations in conservation form are u i x i = 0 (1)
More informationRouen LBV 2012 ACCURACY OF TWO METHODS TO MEASURE LAMINAR FLAME SPEEDS: (1) STEADY BUNSEN BURNER FLAMES AND (2) SPHERICAL FLAMES IN BOMBS.
Rouen LBV 2012 ACCURACY OF TWO METHODS TO MEASURE LAMINAR FLAME SPEEDS: (1) STEADY BUNSEN BURNER FLAMES AND (2) SPHERICAL FLAMES IN BOMBS. A. Bonhomme, T. Boushaki*, L. Selle, B. Ferret and T. Poinsot
More informationA mixed acoustic-entropy combustion instability in a realistic gas turbine
Center for Turbulence Research Proceedings of the Summer Program 2012 449 A mixed acoustic-entropy combustion instability in a realistic gas turbine By E. otheau, L. Selle, Y. ery, T. Poinsot AND F. Nicoud
More informationFor a long time, engine noise has been dominated by fan and jet noise. With their
Challenges in Combustion for Aerospace Propulsion M. Huet, F. Vuillot, N. Bertier (ONERA) M. Mazur, N. Kings, W. Tao, P. Scouflaire, F. Richecoeur, S. Ducruix (Laboratoire EM2C, CNRS, Centrale Supélec,
More informationFlow Structure Investigations in a "Tornado" Combustor
Flow Structure Investigations in a "Tornado" Combustor Igor Matveev Applied Plasma Technologies, Falls Church, Virginia, 46 Serhiy Serbin National University of Shipbuilding, Mikolayiv, Ukraine, 545 Thomas
More informationLOW TEMPERATURE MODEL FOR PREMIXED METHANE FLAME COMBUSTION
ISTP-16, 2005, PRAGUE 16TH INTERNATIONAL SYMPOSIUM ON TRANSPORT PHENOMENA LOW TEMPERATURE MODEL FOR PREMIXED METHANE FLAME MBUSTION M. Forman, J.B.W.Kok,M. Jicha Department of Thermodynamics and Environmental
More informationA REDUCED-ORDER METHANE-AIR COMBUSTION MECHANISM THAT SATISFIES THE DIFFERENTIAL ENTROPY INEQUALITY
THE PUBLISHING HOUSE PROCEEDINGS OF THE ROMANIAN ACADEMY, Series A, OF THE ROMANIAN ACADEMY Special Issue/2018, pp. 285 290 A REDUCED-ORDER METHANE-AIR COMBUSTION MECHANISM THAT SATISFIES THE DIFFERENTIAL
More informationUniversità degli Studi di Firenze Dipartimento di Energetica Sergio Stecco
Università degli Studi di Firenze Dipartimento di Energetica Sergio Stecco Thermo-Acoustic Analysis of an Advanced Lean Injection System in a Tubular Combustor Configuration A. Andreini 1, B. Facchini
More informationA NUMERICAL ANALYSIS OF COMBUSTION PROCESS IN AN AXISYMMETRIC COMBUSTION CHAMBER
SCIENTIFIC RESEARCH AND EDUCATION IN THE AIR FORCE-AFASES 2016 A NUMERICAL ANALYSIS OF COMBUSTION PROCESS IN AN AXISYMMETRIC COMBUSTION CHAMBER Alexandru DUMITRACHE*, Florin FRUNZULICA ** *Institute of
More informationWall-Functions and Boundary Layer Response to Pulsating and Oscillating Turbulent Channel Flows
K. Hanjalić, Y. Nagano and S. Jakirlić (Editors) Wall-Functions and Boundary Layer Response to Pulsating and Oscillating Turbulent Channel Flows D. Panara 1, M. Porta 2,R. Dannecker 1, and B. Noll 1 1
More informationNUMERICAL ANALYSIS OF TURBULENT FLAME IN AN ENCLOSED CHAMBER
NUMERICAL ANALYSIS OF TURBULENT FLAME IN AN ENCLOSED CHAMBER Naveen Kumar D 1*, Pradeep R 2 and Bhaktavatsala H R 3 1 Assistant Professor Department of Mechanical Engineering, M S Engineering College,
More informationCoupling tabulated chemistry with large-eddy simulation of turbulent reactive flows
Center for Turbulence Research Proceedings of the Summer Program 2008 237 Coupling tabulated chemistry with large-eddy simulation of turbulent reactive flows By R. Vicquelin, B. Fiorina, N. Darabiha, D.
More informationCFD Simulation of Internal Flowfield of Dual-mode Scramjet
CFD Simulation of Internal Flowfield of Dual-mode Scramjet C. Butcher, K. Yu Department of Aerospace Engineering, University of Maryland, College Park, MD, USA Abstract: The internal flowfield of a hypersonic
More informationHELMHOLTZ RESONATORS FOR DAMPING COMBUSTOR THERMOACOUSTICS
HELMHOLTZ RESONATORS FOR DAMPING COMBUSTOR THERMOACOUSTICS Dong Yang and Aimee S. Morgans Department of Aeronautics, Imperial College London, London, UK, SW7 AZ email: d.yang13@imperial.ac.uk Helmholtz
More informationFLAME WRINKLING FACTOR DYNAMIC MODELING FOR LARGE EDDY SIMULATIONS OF TURBULENT PREMIXED COMBUSTION
August 8 -, Poitiers, France FLAME WRINKLING FACTOR DYNAMIC MODELING FOR LARGE EDDY SIMULATIONS OF TURBULENT PREMIXED COMBUSTION Thomas Schmitt, Matthieu Boileau, Denis Veynante Laboratoire EMC CNRS -
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationMass Transfer in Turbulent Flow
Mass Transfer in Turbulent Flow ChEn 6603 References: S.. Pope. Turbulent Flows. Cambridge University Press, New York, 2000. D. C. Wilcox. Turbulence Modeling for CFD. DCW Industries, La Caada CA, 2000.
More informationA Priori Model for the Effective Lewis Numbers in Premixed Turbulent Flames
Paper # 070LT-0267 Topic: Turbulent Flames 8 th US National Combustion Meeting Organized by the Western States Section of the Combustion Institute and hosted by the University of Utah May 19-22, 2013.
More informationThermoacoustic Instabilities in a Gas Turbine Combustor
Thermoacoustic Instabilities in a Gas Turbine Combustor The Royal Institute of Technology School of Technological Sciences Department of Vehicle and Aeronautical Engineering The Marcus Wallenberg Laboratory
More informationLarge Eddy Simulation of a turbulent premixed flame stabilized by a backward facing step.
Large Eddy Simulation of a turbulent premixed flame stabilized by a backward facing step. Angelo Murrone and Dominique Scherrer ONERA, 29 Avenue de la Division Leclerc, BP. 72, 92322 Châtillon, France
More informationCFD Analysis of Vented Lean Hydrogen Deflagrations in an ISO Container
35 th UKELG Meeting, Spadeadam, 10-12 Oct. 2017 CFD Analysis of Vented Lean Hydrogen Deflagrations in an ISO Container Vendra C. Madhav Rao & Jennifer X. Wen Warwick FIRE, School of Engineering University
More informationANSYS Advanced Solutions for Gas Turbine Combustion. Gilles Eggenspieler 2011 ANSYS, Inc.
ANSYS Advanced Solutions for Gas Turbine Combustion Gilles Eggenspieler ANSYS, Inc. 1 Agenda Steady State: New and Existing Capabilities Reduced Order Combustion Models Finite-Rate Chemistry Models Chemistry
More informationValida&on of the FDF methodology with a mul&dimensional Helmholtz solver
Valida&on of the FDF methodology with a mul&dimensional Helmholtz solver Alexis Cuquel, Frédéric Boudy, Daniel Durox, Thierry Schuller and Sébas;en Candel Laboratoire EM2C CNRS Ecole Centrale Paris. KIAI
More informationAME 513. " Lecture 8 Premixed flames I: Propagation rates
AME 53 Principles of Combustion " Lecture 8 Premixed flames I: Propagation rates Outline" Rankine-Hugoniot relations Hugoniot curves Rayleigh lines Families of solutions Detonations Chapman-Jouget Others
More informationCombustion Behind Shock Waves
Paper 3F-29 Fall 23 Western States Section/Combustion Institute 1 Abstract Combustion Behind Shock Waves Sandeep Singh, Daniel Lieberman, and Joseph E. Shepherd 1 Graduate Aeronautical Laboratories, California
More informationSpontaneous Oscillations in LNGT Combustors: CFD Simulation
Spontaneous Oscillations in LNGT Combustors: CFD Simulation V. Di Sarli, A. Di Benedetto and F. S. Marra Istituto di Ricerche sulla Combustione - C.N.R., Naples - ITALY INTRODUCTION The development of
More informationEVALUATION OF FOUR TURBULENCE MODELS IN THE INTERACTION OF MULTI BURNERS SWIRLING FLOWS
EVALUATION OF FOUR TURBULENCE MODELS IN THE INTERACTION OF MULTI BURNERS SWIRLING FLOWS A Aroussi, S Kucukgokoglan, S.J.Pickering, M.Menacer School of Mechanical, Materials, Manufacturing Engineering and
More informationModeling instabilities in lean premixed turbulent combustors using detailed chemical kinetics
Accepted for publication in Combustion Science and Technology Modeling instabilities in lean premixed turbulent combustors using detailed chemical kinetics Bjørn Lilleberg, Ivar S. Ertesvåg and Kjell Erik
More informationTowards the prediction of soot in aero-engine combustors with large eddy simulation
Center for Turbulence Research Proceedings of the Summer Program 2014 117 Towards the prediction of soot in aero-engine combustors with large eddy simulation By B. Cuenot, E. Riber AND B. Franzelli A simplified
More informationTurbulence: Basic Physics and Engineering Modeling
DEPARTMENT OF ENERGETICS Turbulence: Basic Physics and Engineering Modeling Numerical Heat Transfer Pietro Asinari, PhD Spring 2007, TOP UIC Program: The Master of Science Degree of the University of Illinois
More informationLarge-Eddy Simulation of spray combustion in aeronautical burners
Large-Eddy Simulation of spray combustion in aeronautical burners Patricia Domingo Alvarez, Ghislain Lartigue and Vincent Moureau INSA Rouen, CNRS, CORIA, 76000 Rouen, France Large-Eddy Simulations (LES)
More informationLES Approaches to Combustion
LES Approaches to combustion LES Approaches to combustion LES Approaches to Combustion W P Jones Department of Mechanical Engineering Imperial College London Exhibition Road London SW7 2AZ SIG on Combustion
More informationHybrid RANS/LES simulations of a cavitating flow in Venturi
Hybrid RANS/LES simulations of a cavitating flow in Venturi TSFP-8, 28-31 August 2013 - Poitiers - France Jean Decaix jean.decaix@hevs.ch University of Applied Sciences, Western Switzerland Eric Goncalves
More informationD. VEYNANTE. Introduction à la Combustion Turbulente. Dimanche 30 Mai 2010, 09h00 10h30
D. VEYNANTE Introduction à la Combustion Turbulente Dimanche 30 Mai 2010, 09h00 10h30 Introduction to turbulent combustion D. Veynante Laboratoire E.M2.C. CNRS - Ecole Centrale Paris Châtenay-Malabry France
More informationKinetic study of combustion behavior in a gas turbine -Influence from varying natural gas composition
Kinetic study of combustion behavior in a gas turbine -Influence from varying natural gas composition Catharina Tillmark April 18, 2006 Lund University Dept. of Energy Sciences P.O.Box 118, SE-221 00 Lund
More informationLES modeling of heat and mass transfer in turbulent recirculated flows E. Baake 1, B. Nacke 1, A. Umbrashko 2, A. Jakovics 2
MAGNETOHYDRODYNAMICS Vol. 00 (1964), No. 00, pp. 1 5 LES modeling of heat and mass transfer in turbulent recirculated flows E. Baake 1, B. Nacke 1, A. Umbrashko 2, A. Jakovics 2 1 Institute for Electrothermal
More informationModelling thermo-acoustic instabilities in an oxy-fuel premixed burner
Modelling thermo-acoustic instabilities in an oy-fuel premied burner Nils Erland L. Haugen Øyvind Langørgen Sigurd Sannan Acknowledgements This presentation forms a part of the BIGCO project, performed
More informationImpact of numerical method on auto-ignition in a temporally evolving mixing layer at various initial conditions
Journal of Physics: Conference Series PAPER OPEN ACCESS Impact of numerical method on auto-ignition in a temporally evolving mixing layer at various initial conditions To cite this article: A Rosiak and
More informationPublished in: Proceedings of the Fluent Benelux User Group Meeting, 6-7 October 2005, Wavre, Belgium
The application of Flamelet Generated Manifolds in partailly-premixed flames Ramaekers, W.J.S.; Albrecht, B.A.; van Oijen, J.A.; de Goey, L.P.H.; Eggels, R.L.G.M. Published in: Proceedings of the Fluent
More informationCombustion Instability Modelling Using Different Flame Models
Combustion Instability Modelling Using Different Flame Models Somayeh Nosrati Shoar, Abbas Fakhrtabatabaei MAPNA Turbine Engineering and Manufacturing Company (TUGA MAPNA Building, No., Mirdamad Ave, Tehran,
More informationTowards regime identification and appropriate chemistry tabulation for computation of autoigniting turbulent reacting flows
Center for Turbulence Research Annual Research Briefs 009 199 Towards regime identification and appropriate chemistry tabulation for computation of autoigniting turbulent reacting flows By M. Kostka, E.
More informationLarge eddy simulation of spark ignition in a turbulent methane jet
Large eddy simulation of spark ignition in a turbulent methane jet Guilhem Lacaze, Edward Richardson, Thierry Poinsot To cite this version: Guilhem Lacaze, Edward Richardson, Thierry Poinsot. Large eddy
More informationLecture 9 Laminar Diffusion Flame Configurations
Lecture 9 Laminar Diffusion Flame Configurations 9.-1 Different Flame Geometries and Single Droplet Burning Solutions for the velocities and the mixture fraction fields for some typical laminar flame configurations.
More informationENTROPY PRODUCTION MODELING IN CFD OF TURBULENT COMBUSTION FLOW
ENTROPY PRODUCTION MODELING IN CFD OF TURBULENT COMBUSTION FLOW Ivar S. Ertesvåg, Jostein Kolbu Department of Energy and Process Engineering Norwegian University of Science and Technology NO-7491 Trondheim,
More informationLES AND ACOUSTIC ANALYSIS OF COMBUSTION INSTABILITIES IN GASTURBINES
LES AND ACOUSTIC ANALYSIS OF COMBUSTION INSTABILITIES IN GASTURBINES T. Poinsot and L. Selle IMF and CERFACS, Toulouse, France Abstract Keywords RANS (Reynolds Averaged) techniques for combustion constitute
More informationDNS of Reacting H 2 /Air Laminar Vortex Rings
46th AIAA Aerospace Sciences Meeting and Exhibit 7-10 January 2008, Reno, Nevada AIAA 2008-508 DNS of Reacting H 2 /Air Laminar Vortex Rings Jeff Doom and Krishnan Mahesh University of Minnesota, Minneapolis,
More informationLarge Eddy Simulation of Flame Flashback by Combustion Induced Vortex Breakdown
June 30 - July 3, 2015 Melbourne, Australia 9 1C-5 Large Eddy Simulation of Flame Flashback by Combustion Induced Vortex Breakdown Eike Tangermann Institut für Mathematik und Rechneranwendung Universität
More informationLarge-eddy simulation of an industrial furnace with a cross-flow-jet combustion system
Center for Turbulence Research Annual Research Briefs 2007 231 Large-eddy simulation of an industrial furnace with a cross-flow-jet combustion system By L. Wang AND H. Pitsch 1. Motivation and objectives
More informationEffects of Variation of the Flame Area and Natural Damping on Primary Acoustic Instability of Downward Propagating Flames in a Tube
5 th ICDERS August 7, 015 Leeds, UK Effects of Variation of the Flame Area and Natural Damping on Primary Acoustic Instability of Downward Propagating Flames in a Tube Sung Hwan Yoon and Osamu Fujita Division
More informationDirect numerical simulation of a turbulent reacting jet
Center for Turbulence Research Annual Research Briefs 999 59 Direct numerical simulation of a turbulent reacting jet By B. J. Boersma. Motivation and objectives Turbulent reacting jets are important in
More informationModeling chemical flame structure and combustion dynamics in LES
Author manuscript, published in "Proceeding of the combustion Institute (2011) In Press" Modeling chemical flame structure and combustion dynamics in LES P. Auzillon a,, B. Fiorina a, R. Vicquelin a,b,
More informationAdvanced Numerical Methods for non-premixed Flames
0 Advanced Numerical Methods for non-premixed Flames Annarita Viggiano Department of Environmental Engineering and Physics University of Basilicata Italy 1. Introduction Engine designers are under increasing
More informationIMPROVED POLLUTANT PREDICTIONS IN LARGE-EDDY SIMULATIONS OF TURBULENT NON-PREMIXED COMBUSTION BY CONSIDERING SCALAR DISSIPATION RATE FLUCTUATIONS
Proceedings of the Combustion Institute, Volume 9, 00/pp. 1971 1978 IMPROVED POLLUTANT PREDICTIONS IN LARGE-EDDY SIMULATIONS OF TURBULENT NON-PREMIXED COMBUSTION BY CONSIDERING SCALAR DISSIPATION RATE
More informationarxiv: v1 [physics.flu-dyn] 25 Nov 2018
Combustion regimes in sequential combustors: Flame propagation and autoignition at elevated temperature and pressure O. Schulz,a, N. Noiray,a a CAPS Laboratory, Department of Mechanical and Process Engineering,
More informationMUSCLES. Presented by: Frank Wetze University of Karlsruhe (TH) - EBI / VB month review, 21 September 2004, Karlsruhe
MUSCLES Modelling of UnSteady Combustion in Low Emission Systems G4RD-CT-2002-00644 R&T project within the 5 th Framework program of the European Union: 1 Numerical computations of reacting flow field
More informationLarge eddy simulations of gaseous flames in gas turbine combustion chambers
Large eddy simulations of gaseous flames in gas turbine combustion chambers L.Y.M. Gicquel,G.Sta elbachandt.poinsot CERFACS, 42 Avenue G. Coriolis, 31057 Toulouse Cedex 1, France Abstract Recent developments
More informationAn Unsteady/Flamelet Progress Variable Method for LES of Nonpremixed Turbulent Combustion
43rd AIAA Aerospace Sciences Meeting and Exhibit, -3 Jan 25, Reno, NV An Unsteady/Flamelet Progress Variable Method for LES of Nonpremixed Turbulent Combustion Heinz Pitsch and Matthias Ihme Stanford University,
More informationComputational Fluid Dynamics 2
Seite 1 Introduction Computational Fluid Dynamics 11.07.2016 Computational Fluid Dynamics 2 Turbulence effects and Particle transport Martin Pietsch Computational Biomechanics Summer Term 2016 Seite 2
More informationDEVELOPMENT OF CFD MODEL FOR A SWIRL STABILIZED SPRAY COMBUSTOR
DRAFT Proceedings of ASME IMECE: International Mechanical Engineering Conference & Exposition Chicago, Illinois Nov. 5-10, 2006 IMECE2006-14867 DEVELOPMENT OF CFD MODEL FOR A SWIRL STABILIZED SPRAY COMBUSTOR
More information(Received?? and in revised form??)
Under consideration for publication in J. Fluid Mech. Large-eddy simulation and experimental study of heat transfer, nitric oxide emissions and combustion instability in a swirled turbulent high pressure
More informationFlamelet Analysis of Turbulent Combustion
Flamelet Analysis of Turbulent Combustion R.J.M. Bastiaans,2, S.M. Martin, H. Pitsch,J.A.vanOijen 2, and L.P.H. de Goey 2 Center for Turbulence Research, Stanford University, CA 9435, USA 2 Eindhoven University
More informationNumerical Modelling of Aerodynamic Noise in Compressible Flows
Open Journal of Fluid Dynamics 22 2 65-69 http://dx.doi.org/.4236/ofd.22.237 Published Online September 22 (http://www.scirp.org/ournal/ofd) Numerical Modelling of Aerodynamic Noise in Compressible Flows
More informationAcoustic and Large Eddy Simulation studies of azimuthal modes in annular combustion chambers
Acoustic and Large Eddy Simulation studies of azimuthal modes in annular combustion chambers Pierre Wolf, Gabriel Staffelbach, Laurent Y.M. Gicquel, Jens-Dominik Müller, Thierry Poinsot To cite this version:
More informationExperimental Study on the Non-reacting Flowfield of a Low Swirl Burner
Experimental Study on the Non-reacting Flowfield of a Low Swirl Burner Hang Yin & Ren Dai School of Energy and Powering Engineering, University of Shanghai for Science and Technology Box 25, 516# Jungong
More informationLarge-Eddy Simulation of Turbulent Flames in Syn-Gas Fuel-Air Mixtures
45th AIAA Aerospace Sciences Meeting and Exhibit 8-11 January 2007, Reno, Nevada AIAA 2007-1435 Large-Eddy Simulation of Turbulent Flames in Syn-Gas Fuel-Air Mixtures Baris A. Sen and Suresh Menon School
More informationInvestigation of an implicit solver for the simulation of bubble oscillations using Basilisk
Investigation of an implicit solver for the simulation of bubble oscillations using Basilisk D. Fuster, and S. Popinet Sorbonne Universités, UPMC Univ Paris 6, CNRS, UMR 79 Institut Jean Le Rond d Alembert,
More informationReaction Rate Closure for Turbulent Detonation Propagation through CLEM-LES
5 th ICDERS August 7, 05 Leeds, UK through CLEM-LES Brian Maxwell, Matei Radulescu Department of Mechanical Engineering, University of Ottawa 6 Louis Pasteur, Ottawa, KN 6N5, Canada Sam Falle School of
More informationOpen Archive Toulouse Archive Ouverte (OATAO)
Open Archive Toulouse Archive Ouverte (OATAO) OATAO is an open access repository that collects the work of Toulouse researchers and makes it freely available over the web where possible. This is an author-deposited
More informationLiquid-Rocket Transverse Triggered Combustion Instability: Deterministic and Stochastic Analyses
Liquid-Rocket Transverse Triggered Combustion Instability: Deterministic and Stochastic Analyses by W. A. Sirignano Mechanical and Aerospace Engineering University of California, Irvine Collaborators:
More informationDispersion of Entropy Perturbations Transporting through an Industrial Gas Turbine Combustor
Flow Turbulence Combust (218) 1:481 52 DOI 1.17/s1494-17-9854-6 Dispersion of Entropy Perturbations Transporting through an Industrial Gas Turbine Combustor Yu Xia 1 Ignacio Duran 2 Aimee S. Morgans 1
More informationNumerical Simulation of Entropy Generation in Hydrogen Enriched Swirl Stabilized Combustion
Saqr & Wahid CFD Letters Vol. 5(1) 13 www.cfdl.issres.net Vol. 5 (1) March 13 Numerical Simulation of Entropy Generation in Hydrogen Enriched Swirl Stabilized Combustion Khalid M. Saqr 1,* and Mazlan A.
More informationIntroduction Flares: safe burning of waste hydrocarbons Oilfields, refinery, LNG Pollutants: NO x, CO 2, CO, unburned hydrocarbons, greenhouse gases G
School of Process, Environmental and Materials Engineering Computational study of combustion in flares: structure and emission of a jet flame in turbulent cross-flow GERG Academic Network Event Brussels
More informationLARGE EDDY SIMULATION AND FLOW CONTROL OVER A 25 RAMP MODEL
LARGE EDDY SIMULATION AND FLOW CONTROL OVER A 25 RAMP MODEL 09/11/2017 Paolo Casco Stephie Edwige Philippe Gilotte Iraj Mortazavi LES and flow control over a 25 ramp model : context 2 Context Validation
More informationPhysics of turbulent flow
ECL-MOD 3A & MSc. Physics of turbulent flow Christophe Bailly Université de Lyon, Ecole Centrale de Lyon & LMFA - UMR CNRS 5509 http://acoustique.ec-lyon.fr Outline of the course A short introduction to
More information