Large-Eddy Simulation of spray combustion in aeronautical burners

Size: px
Start display at page:

Download "Large-Eddy Simulation of spray combustion in aeronautical burners"

Transcription

1 Large-Eddy Simulation of spray combustion in aeronautical burners Patricia Domingo Alvarez, Ghislain Lartigue and Vincent Moureau INSA Rouen, CNRS, CORIA, Rouen, France Large-Eddy Simulations (LES) of a novel lean-premixed (LP) injection system for aeronautical burners are performed and reported. In this LP injector, liquid kerosene is injected far upstream of the injector outlet in order to provide a gaseous fuel distribution as homogeneous as possible. In the LES, the full geometry of the injector is taken into account. A Lagrangian point-particle approach is chosen for the spray modeling. The two-phase combustion model relies on a tabulated chemistry approach. The flame/turbulence interactions are included via a presumed-pdf closure (PCM-FPI model). Numerical results are analyzed in order to bring more insights into the topology of such flames. I. Introduction The Advisory Council of Aeronautic Research in Europe (ACARE) gives mid and long-term objectives in order to limit pollutant emissions. As such, the 2020 objectives in the aviation transportation plan a reduction of CO 2 by 50% (i.e. a 50% reduction of the consumption) compared to the level in In the same time, an 80% reduction of the emission of NOx is required as well as a 50% reduction of the noise. Those requirements necessitate new innovative engine technologies using state-of-the-art injection technology such as LPP (Lean-Premixed-Prevaporized), LP (Lean Premixed), or variants. Those technologies enable to decrease dramatically pollution emissions by ensuring a better control of the combustion regime. The design of such injectors remains difficult and requires predictive models for their optimization. In order to achieve these requirements, industry is looking for new methods of combustion. The present work is lead by the industrial chair PERCEVAL (PowERing the future with Clean and Efficient Aero-engines) with the support of SAFRAN Tech, part of the SAFRAN group and the ANR (Agence Nationale de la Recherche). It is focused on advanced studies on the optimization of combustion efficiency and reduction of pollutants. This chair combines advanced time-resolved laser diagnostics and analysis techniques with highresolution time-resolved Large-Eddy simulation in a focused and highly-integrated program of scientific and technical research. Experiments are conducted in a new high-pressure high-temperature optical combustion facility fed with innovative lean fuel injection systems designed by the industrial partner. The chair is hosted by the CORIA joint laboratory organized between CNRS, University of Rouen and INSA-Rouen. Computational fluid dynamics (CFD) complements the experimental and theoretical approaches. It helps in interpreting and understanding theory and experiments and vice-versa. Numerical tools are very important when accompanied with theoretical knowledge in order to be able to study, analyze and understand complex combustion regimes at different conditions (as pressure, temperature) and configurations (1D, 2D, 3D). Several modeling frameworks exist: RANS (Reynolds Averaged Navier Stokes), DNS (Direct Numerical Simulation) or LES (Large Eddy Simulations). The latter is chosen in the present work, as it gives access to unsteady flow features that are present in such burners, at a reduced CPU cost. The paper is organized as follows: The theoretical framework of this study is presented in section II. The spray (sec. II.D), chemistry (sec. II.C) and combustion (sec. II.B) models are explained. The development of the experimental set-up is given in section III. The CFD code YALES2 which is used to accomplish the simulations is presented. Simulations corresponding to the aeronautical burner and the results obtained are exposed in section III.A. 1 of 10

2 II. Models II.A. LES framework In the present study, the Large-Eddy Simulation (LES) framework is considered. The objective of LES is to explicitly compute the large structures of the flow field, that are typically larger than the computational mesh size, whereas the effects of the smallest ones are modeled. LES is particularly well suited for turbulent combustion. 1 Even if reactions mostly occur at the sub-grid scale, LES can resolve the large-scale dynamics of the turbulent flame brush, which is important for burner performances such as pollutant emissions. In LES, variables are filtered in spectral space or in physical space and the filtered quantity f is defined as: f(x) = f(x )F (x x )dx (1) where F is the LES filter, where different type of filters are available. The filtered balance equations, coupled with the sub-grid models are numerically solved to determine the instantaneous filtered fields. The Navier-Stokes equations are written as follows: Mass conservation Momentum Energy ρ h s t + ρũ i h s x i ρũ i t = DP Dt + ρ t + (ρũ i ) = 0 (2) x i + (ρũ i ũ j ) + p = [τ ij ρ(ũ i u j ũ i ũ j )] (3) x i x i x i x i ( ) T λ th ρ(ũih s ũ i hs + ρ x i x i N V k,i Y k h s,k + ω T (4) k=1 Chemical species ρỹk t + x i ( ρũ i Ỹ k ) = x i [V k,i Y k ρ( u i Y k ũ i Ỹ k )] + ω k (5) II.B. Combustion model Two combustion regimes are usually defined: premixed and non-premixed flames. Moreover, intermediate regimes exists, as the partially premixed flames. A premixed flame occurs when the reactants are mixed before combustion. Fresh and burnt gases are separated by the flame front which propagates towards the fresh gases at a velocity of s L. The flame front thickness, δ L is composed by two different zones. The first one is called the pre-heat zone where diffusion phenomena occurs. The second one is the reaction zone characterized by a high heat release. The process evolution can be described by the normalized progress variable, c. This parameter is equal to 0 in the fresh gases and 1 in the burnt gases. c is generally built as a linear combination of species mass fractions. The mixture, perfectly homogeneous theoretically, is described by the equivalence ratio. Due to the temperature gradients, the flame propagates due to diffusion towards the fresh gases. The laminar flame speed as well as the flame thickness become fundamental parameters of the premixed combustion: s L = ω cdx and δ L,th = T b T u. (6) ρ u T max Non-premixed flames are characterized by the separated injection of fuel and oxidant. The flame develops in the diffusion zone and stabilizes at the stoichiometric iso-contour. Contrarily to premixed flames, nonpremixed flames are driven by the mixture. The parameter used to describe this type of flame is the mixture fraction, Z, equal to 0 for the oxidant and 1 for the fuel. Z is also generally built as a linear combination of species mass fractions. Another important parameter is the scalar dissipation rate. This parameter allows to know the characteristic time of diffusion. High dissipation rate implies fast species diffusion. 2 of 10

3 However, a third flame regime is encountered in most of the practical cases. In this regime of partially premixed combustion, 2 the mixture is partially mixed. The Takeno index 3 is used to determine the type of combustion regime in space. It is defined as: T = Y F Y O Y F Y O, (7) where Y F is the mass fraction of fuel and Y O the oxidant mass fraction. This index can be equal to 1, corresponding to premixed combustion zones and 1 in diffusion combustion zones. II.C. Tabulated chemistry: PCM-FPI In this type of model, thermodynamic data, i.e. species mass fractions and reactions rates, are tabulated as functions of a limited set of coordinates. Only this set of coordinates needs to be transported. For laminar adiabatic combustion with only two streams, the lookup table is a function of two coordinates: the progress variable, Y c, which evolves monotonically between fresh and burnt gases and the mixture fraction, Z, equal to zero in pure oxidizer and unity in pure fuel. Kerosene/air combustion involves a large number of species and reactions. Different strategies have been developed to reduce this complexity. Mass and Pope, 45 proposed the method of Intrinsic Low Dimensional Manifold (ILDM). It benefits from the large range of chemical time scales. The temporal evolution of a chemical system depends on its initial conditions, but after a short time, this system can be described by a reduced set of variables in the composition space. Then, on this low-dimensional manifold, species mass fractions and reaction rates are determined as functions of a reduced set of variables, which may contain a progress variable measuring the evolution of the reaction and a mixture fraction describing the mixing. A table is then created with this set of variables and searched by multi-linear interpolation. In order to describe correctly low temperature regions, Gicquel et al., 6 Oijen et al., 7 Fiorina et al. 8 and Goey et al. 9 proposed the Flame prolongation of ILDM (FPI) and Flamelet Generated Manifold (FGM) approaches. The basic idea is to generate look-up tables from simulations of one-dimensional laminar premixed flames using finite-rate chemical schemes. Reaction rates and species mass fractions are then tabulated as functions of a limited set of coordinates (progress variable, mixture fraction,...). Chemical tables may become large when additional parameters are introduced, leading to practical difficulties on massively parallel machines. On the one hand, ILDM is a method derived from the mathematical study of the chemical system. On the other hand, FPI assumes the flame front as a 1D laminar flame front. Hence, the flame structure does not depend on the local conditions. Tabulation can be applied then to premixed and non-premixed flames as well as partially premixed flames. In LES, FPI is combined with the Presumed Conditional Moments (PCM) method. It allows for taking into account the sub-grid flame-turbulence interaction. The look-up table is not parametrized by only two variables, as the mixture fraction and the progress variable, but also their second order moments, S Z and S c. This method belongs to the statistical approaches, where the source terms are modeled by a probability density function (PDF): ω k = ωp(z, c)dzdc, (8) where P(Z, c) is the PDF of Z and c. From Vervisch, 10 the PDF can be decomposed in: P(Z, c) = P(c Z)P(c)P(Z). (9) The source term must be written then as a function of two independent PDF, where the PDF shape must be presumed. The most popular PDF is the β-function, because of its ability to adopt the shapes encountered in mixing systems. Applied to the normalized variables ϕ, it reads P(ϕ) = ϕ (a 1) (1 ϕ) b (ϕa 1 )(1 ϕ) b 1 )dϕ, (10) 3 of 10

4 where a and b are parameters determined by the relationships : ( ) 1 a = ϕ 1 and b = ã a, (11) S ϕ ϕ where S ϕ is the segregation parameter of ϕ: S ϕ = ϕ 2 ϕ(1 ϕ). (12) The PDF is then a function of Z and c. The source term of any species and of the progress variable is a function of four parameters: ω k ( Z, c, Z 2, c 2 ) = ω k (Z, c)β Z,SZ (Z)β c,sc (c)dzdc, (13) and it is stored in a table with four inputs: Z, c, Z 2 and c 2. II.D. Spray model In many practical aeronautical applications, liquid fuel is injected into the combustion chamber resulting in fuel spray. Spray combustion involves many physical processes as vaporization, heat and mass transfer, droplet-air and vapor-air mixing, auto-ignition, flame extinction, turbulence, premixed and/or non-premixed flames and pollutant production. The spray is modeled with a point-particle Lagrangian approach. It consists on the individual treatment of each droplet in the computational domain. Mass and heat transfer, momentum and energy applies between the two phases. Particles have their own properties as temperature, diameter and velocity. Thermodynamic properties of the droplet are obtained from the conservation equations. Moreover, the mass transfer model of Spalding is adopted, 11 which makes the following assumptions: The droplet is perfectly spherical and isolated. Spherical liquid fuel droplet is made up of a single component. The liquid thermal conductivity is assumed infinitely large, leading to a uniform droplet temperature. Thermal diffusivity, D th = λ/(ρc p ), is small compared to the gas. Then, the gas phase thermal characteristic time is smaller than the characteristic time of the droplet. The thermal response is then assumed to be stationary in the gas phase and time derivatives are neglected. The droplet surface is assumed to be at an equilibrium state with the surrounding gas. Surrounding gas phase properties are assumed to be constant from the droplet surface to the far field. The droplet mass transfer is resulting from the integration of the mass conservation law from the droplet to the far field. dm p = πd p ρdsh log(1 + B M ), (14) dt where B M is the Spalding mass number, d p is the diameter of the particle, D is the diffusivity, ρ is the density and Sh is the Sherwood number which represents the ratio of the total rate mass transfer to the rate of diffusive mass transport. The diameter evolution is obtained from the mass transfer equation (14). Its temporal evolution is a function of an evaporation characteristic time τ m and the droplet initial diameter: dd 2 p dt = d2 p,0 2d p τ m. (15) The droplet temperature (T p ) evolution is obtained by integrating the energy conservation from the droplet s surface to the the far field, leading to: dt p = 1 ( ( T p T L )) vb T, (16) dt τ h C p,1/3 4 of 10

5 where τ h is a thermal characteristic time, T is the temperature of the surrounding gas, B T the Spalding thermal number, L v the latent heat and C p,1/3 the heat capacity at a preponderated temperature T 1/3 used to characterize the intermediate state of the gas around the droplet. This approach is validated and confronted to the experimental results of Ghassemi 12 in Fig. 1. An isolated droplet of kerosene with an initial diameter of d 0 = 1.52 mm is placed in an atmosphere at rest at a temperature of T = 973 K and a pressure of P = 1 MPa Ghassemi et al. (2006) Drop Diam (Y2) Drop T (Y2) (d/d 0 ) 2 [-] Temperature [K] t/d 0 2 [s mm -2 ] Figure 1. Isolated droplet evaporating in a stagnant atmosphere II.E. Sub-grid fluxes modeling The sub-grid Reynolds stresses, scalar transport fluxes and enthalpy transport fluxes are closed using models developed for non-reactive flows. In the Smagorinsky subgrid-scale model, 13 unresolved momentum fluxes are expressed according to the Boussinesq assumption. Boussinesq 14, 15 postulates that the momentum transfer caused by turbulent eddies can be modeled with an eddy viscosity and it states that the Reynolds stress tensor, T ij, is proportional to the trace-less mean strain rate tensor, Sij : T ij δ ij 3 T kk = ν t where ν t is a subgrid scale viscosity. It is modeled as: ( ui + u ) j = 2ν j S ij, (17) x j x i ν t = C 2 S 4/3 l 2/3 t S = C 2 S 4/3 l 2/3 t (2S ij S ij ) 1/2, (18) where l t is the turbulence integral length scale, C S a model constant, and S the resolved shear stress. Equation 18 is simplified by assuming that the integral length scale l t is of the order of the grid size l t : ν t = (C S ) 2 S = (C S ) 2 (2S ij S ij ) 1/2. (19) The isotropic contribution τ kk in Eq. 17, corresponding to twice the sub-grid scale turbulent kinetic energy, is unknown and is usually absorbed into the filtered pressure P. In the case of homogeneous isotropic turbulence, the model constant is estimated as C S However, C S depends on the flow configuration. Germano 16 introduced a modification in the Smagorinsky model. The constant C S is automatically computed in space and time. A second filter bigger than is introduced. The sub-grid tensor τ ij = ρ( u i u j ũ i ũ j ) and the sub-grid tensor based on the two times filtered velocity is τ ij = ρ( u i u j ˆũ i ˆũ j ) can be written from the Smagorinsky model as: τ ij = 2ρ(C S ) 2 S S ij and τ ij = 2ˆρ(C S ) 2 ˆ S ˆ Sij. (20) 5 of 10

6 Finally, the Germano identity, which relates both filtering levels, can be computed as: ( ) L ij = τ ij τ ij = ρ ˆũ i ˆũ j ũ i ũ j. (21) From Eqs. 20 and 21, the Smagorinsky constant can be obtained from the filtered velocity fields. III. Numerical set-up LES is performed with the YALES2 code. It is a LES and DNS semi-industrial code developed at CORIA laboratory. It is a low-mach number solver based on the Finite Volume Method (FVM). It is capable to solve the filtered Navier-Stokes equations on unstructured meshes of billions of elements. 17 The innovative lean fuel injection system is designed to operate at high-pressure and high-temperature conditions. As a first step, it is tested at atmospheric pressure in a dedicated combustion chamber. The numerical domain is equivalent to the experimental facility at CORIA laboratory. 18 The geometry is first meshed with 7.1 millions tetrahedral elements (mesh M1). One refinement stage is then applied homogeneously, leading to a 59.8 million elements mesh (mesh M2). The smallest element size is 20 µm in M2. The combustion chamber domain is cm 3. One operating point is computed. The air mass flow rate is ṁ air = 30 g/s at atmospheric pressure and a temperature of T = 473 K. The liquid kerosene fuel is injected by the LP (Lean Premixed) swirled injector resulting in an equivalence ratio of φ = 0.7 in the combustion chamber. In this study, atomization is not considered. Injection is modeled as a fog of droplets directly injected with a prescribed diameter PDF and velocity profile. Experimentally, it has been observed that the fuel hits the injector wall creating a film of fuel. This film is afterwards atomized and re-injected into the combustion chamber. In the modeling of the injection, two injectors are configured, see Fig. 2. This configuration is chosen according to the experimental data. The first one, which injects 3/4 of the fuel mass flow rate, that is ṁ 1 = 1.2 g/s with a mean diameter of 10 µm and the second injector, which injects the remaining 1/4 of the fuel mass flow rate, ṁ 2 = g/s with a mean diameter of 20 µm. In both cases, the Rosin-Rammler distribution is used to model the droplet size distribution, see Eq. 22, F v (D) = 1 exp [ ( D D m ) γ ], (22) where D m is the mean diameter and γ is a constant describing the material uniformity and hence called the uniformity constant. 19 This distribution fits well to various fuel injection types. The atomization here is a consequence of the high velocity of the air, that is, air blast. Results are sensitive to the initial droplet size distribution, which makes it an important step. However, there is no granulometry data available for the moment, which makes the determination of the droplet size a challenge. Figure 2. Fuel injection modeling 6 of 10

7 III.A. Results and discussion The instantaneous and mean velocity fields are presented in Figs.3 and 3(a). Two recirculation zones characteristic of swirled flows, appear. The Outer Recirculation Zone (ORZ), where fresh and burnt gases mix, and the Central Recirculation Zones (CRZ) far from the injection point because of the low swirl number (lower than 0.6). These recirculation zones are marked by the iso-contour at zero mean axial velocity. The spray flame topology is represented by the spray fog, the temperature, and the progress variable in Figs. 4 and 5. The position of the flame can be observed by the temperature field and the progress variable. Flame is anchored to the external wall of the injector outlet. Its V-shape is typical of swirl burners. In Figure 6, the kerosene mass fraction, the Takeno index and the mixture fraction are represented. The evaporated kerosene mass fraction is not homogeneous upstream of the flame front but it exhibits some richer pockets corresponding to remaining evaporating droplets. These pockets are illustrated in Fig. 6(c) with a red iso-contour at Z stoichiometric = These pockets interact with the fresh gases, evaporating and giving the mixture air/kerosene that feeds the flame. The mixture fraction source term with the particles colored by their diameter are represented. Particles are abundant on the jet. They are not evaporated immediately. They are evaporated in the flame front and in the burnt gases. The flame therefore burns in a partially premixed regime. To confirm this observation, the Takeno index is analyzed. The Takeno index gives information about the combustion regime. The zones where Takeno index is equal to 1 are the zones of premixed combustion while in the zones where the Takeno index is equal to 1 combustion is done by diffusion. The mixture fraction allows to know the local equivalence ratio in the chamber. OH mass fraction is a useful marker to locate the transient reaction zones in highly turbulent flows. In Fig. 7, instantaneous and mean OH mass fraction are presented. The V-shape of the flame is clearly visible. Moreover, zones with strong OH concentration gradients are present. The computation time for the simulations has been approximately hours using the resources of the computing center CRIANN 20 with the cluster ANTARES. Despite this computed time, the results are not fully converged, which prevents a further analysis of mean quantities. (a) (b) Figure 3. Midplane instantaneous field of the velocity (a) and mean axial velocity field with iso-lines of velocity null representing recirculation zones (b) IV. Conclusion Industry seeks new combustion technologies to reduce pollution. Experimental research is being carried out in CORIA laboratory with a novel LP injection system. The aim of this work was to perform numerical simulations of this configuration with the LES solver YALES2 to help interpret such physical experiments, and even to ascertain basic phenomenological aspects of the experiments, which are not apparent from the laboratory data. The topology, the aerodynamic and pollutants of this configuration have been studied. Results show a 7 of 10

8 Figure 4. Particles injected in the combustion chamber (a) (b) Figure 5. Midplane instantaneous temperature (a) and progress variable (b) V-shape topology of the flame with a partially-premixed combustion regime. To continue this work, a PhD thesis has started. The same injector is being simulated in a new aeronautical burner configuration. A high-pressure high-temperature combustion chamber is being investigated. First, a reference point of 8.33 bar has been established. The choice of this point is because of the good characterization of the flame experimentally for a wide range of FAR (Fuel-Air Ratio) before the limit of extinction of the flame. In addition, modeling of the optical diagnostics will be also performed in simple and LES configurations. The objective is to help experimentalists designing their optical diagnostics. 8 of 10

9 (a) (b) (c) (d) Figure 6. Intantaneous kerosene mass fraction (a), Takeno index (b) mixture fraction with iso-contour at Zs toi (c) and mixture fraction source term with particles colored by their diameter [µm] (d) (a) (b) Figure 7. Midplane instantaneous OH mass fraction(a) and mean OH mass fraction (b) 9 of 10

10 References 1 Poinsot, T. and Veynante, D., Theoretical and numerical combustion, T. Poinsot, S.l., Peters, N., Turbulent Combustion, Cambridge University Press, Aug Yamashita, H., Shimada, M., and Takeno, T., A numerical study on flame stability at the transition point of jet diffusion flames, Symposium (International) on Combustion, Vol. 26, Elsevier, 1996, pp Maas, U. and Pope, S. B., Simplifying Chemical Kinetics: Intrinsic Low-Dimensional Manifolds in Composition Space, Combustion and Flame, Vol. 88, No. 3-4, March 1992, pp Pope, S. B., PDF methods for turbulent reactive flows, Progress in Energy and Combustion Science, Vol. 11, No. 2, 1985, pp Gicquel, O., Darabiha, N., and Thévenin, D., Laminar premixed hydrogen/air counterflow flame simulations using Flame Prolongation of ILDM with differential diffusion, Proceedings of the Combustion Institute, Vol. 28, No. 2, Jan. 2000, pp Oijen, J. a. V. and Goey, L. P. H. D., Modelling of Premixed Laminar Flames using Flamelet-Generated Manifolds, Combustion Science and Technology, Vol. 161, No. 1, Dec. 2000, pp Fiorina, B., Gicquel, O., Vervisch, L., Carpentier, S., and Darabiha, N., Premixed turbulent combustion modeling using tabulated detailed chemistry and PDF, Proceedings of the Combustion Institute, Vol. 30, No. 1, Jan. 2005, pp Goey, L. P. H. D., Oijen, J. A. V., Bongers, H., and Groot, G. R. A., New flamelet based reduction methods: the bridge between chemical reduction techniques and flamelet methods, European Combustion Meeting. 10 Vervisch, L., Hauguel, R., Domingo, P., and Rullaud, M., Three facets of turbulent combustion modelling: DNS of premixed V-flame, LES of lifted nonpremixed flame and RANS of jet-flame, Journal of Turbulence, Vol. 5, Jan Spalding, D. B., The combustion of liquid fuels, Symposium (International) on Combustion, Vol. 4, No. 1, Jan. 1953, pp Ghassemi, H., Baek, S. W., and Khan, Q. S., Experimental Study on Evaporation of Kerosene Droplets at Elevated Pressures and Temperatures, Combustion Science and Technology, Vol. 178, No. 9, Sept. 2006, pp Smagorinsky, J., General circulation experiments with the primitive equations, Monthly Weather Review, Vol. 91, No. 3, March 1963, pp Boussinesq, J..-. A. d. t., Essai sur la théorie des eaux courantes / par J. Boussinesq, Impr. nationale, Paris, Schmitt, F. G., About Boussinesq s turbulent viscosity hypothesis: historical remarks and a direct evaluation of its validity, Comptes Rendus Mécanique, Vol. 335, No. 910, Sept. 2007, pp A dynamic subgridscale eddy viscosity model, Physics of Fluids A: Fluid Dynamics, Vol. 3, No. 7, July 1991, pp Moureau, V., Domingo, P., and Vervisch, L., Design of a massively parallel CFD code for complex geometries, Comptes rendus mecanique, Vol. 339, No. 2, Feb. 2011, pp Salaün, E., Malbois, P., Vandel, A., Godard, G., Grisch, F., Renou, B., Cabot, G., and Boukhalfa, A. M., Experimental investigation of a spray swirled flame in gas turbine model combustor, 18th International Symposium on the Application of Laser and Imaging Techniques to Fluid Mechanics, Lisbon, Portugal, Vesilind, P. A., The Rosin-Rammler particle size distribution, Resource Recovery and Conservation, Vol. 5, No. 3, Sept. 1980, pp Documentation IBM cluster idataplex - ANTARES CRIANN,. 10 of 10

DNS and LES of Turbulent Combustion

DNS and LES of Turbulent Combustion Computational Fluid Dynamics In Chemical Reaction Engineering IV June 19-24, 2005 Barga, Italy DNS and LES of Turbulent Combustion Luc Vervisch INSA de Rouen, IUF, CORIA-CNRS Pascale Domingo, Julien Réveillon

More information

Predicting NO Formation with Flamelet Generated Manifolds

Predicting NO Formation with Flamelet Generated Manifolds Predicting NO Formation with Flamelet Generated Manifolds J. A. van Oijen and L. P. H. de Goey Dept. Mechanical Engineering, Technische Universiteit Eindhoven P.O. Box, 6 MB Eindhoven, The Netherlands

More information

Coupling tabulated chemistry with large-eddy simulation of turbulent reactive flows

Coupling 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 information

Self-similar behavior of chemistry tabulation in laminar and turbulent multi-fuel injection combustion systems

Self-similar behavior of chemistry tabulation in laminar and turbulent multi-fuel injection combustion systems Center for Turbulence Research Proceedings of the Summer Program 26 39 Self-similar behavior of chemistry tabulation in laminar and turbulent multi-fuel injection combustion systems By A. Naudin, B. Fiorina,

More information

Analysis of dynamic models for turbulent premixed combustion

Analysis 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 information

A validation study of the flamelet approach s ability to predict flame structure when fluid mechanics are fully resolved

A validation study of the flamelet approach s ability to predict flame structure when fluid mechanics are fully resolved Center for Turbulence Research Annual Research Briefs 2009 185 A validation study of the flamelet approach s ability to predict flame structure when fluid mechanics are fully resolved By E. Knudsen AND

More information

LES Approaches to Combustion

LES 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 information

D. VEYNANTE. Introduction à la Combustion Turbulente. Dimanche 30 Mai 2010, 09h00 10h30

D. 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 information

LARGE-EDDY SIMULATION OF PARTIALLY PREMIXED TURBULENT COMBUSTION

LARGE-EDDY SIMULATION OF PARTIALLY PREMIXED TURBULENT COMBUSTION LARGE-EDDY SIMULATION OF PARTIALLY PREMIXED TURBULENT COMBUSTION Heinz Pitsch Mechanical Engineering Department Stanford University Stanford, CA 94305, USA h.pitsch@stanford.edu ABSTRACT The development

More information

Towards regime identification and appropriate chemistry tabulation for computation of autoigniting turbulent reacting flows

Towards 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 information

Large-eddy simulation of an evaporating and reacting spray

Large-eddy simulation of an evaporating and reacting spray Center for Turbulence Research Annual Research Briefs 2008 479 Large-eddy simulation of an evaporating and reacting spray By T. Lederlin AND H. Pitsch 1. Motivation and objectives 1.1. Evaporative spray

More information

Flame / wall interaction and maximum wall heat fluxes in diffusion burners

Flame / 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 information

Flow Structure Investigations in a "Tornado" Combustor

Flow 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 information

FLAME WRINKLING FACTOR DYNAMIC MODELING FOR LARGE EDDY SIMULATIONS OF TURBULENT PREMIXED COMBUSTION

FLAME 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 information

Lecture 9 Laminar Diffusion Flame Configurations

Lecture 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 information

LEAN PREMIXED TURBULENT COMBUSTION MODELING USING FLAME TABULATED CHEMISTRY AND A PRESUMED PDF APPROACH

LEAN 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 information

Flow and added small-scale topologies in a turbulent premixed flame

Flow and added small-scale topologies in a turbulent premixed flame Flow and added small-scale topologies in a turbulent premixed flame L. Cifuentes*, A. Kempf* and C. Dopazo** luis.cifuentes@uni-due.de *University of Duisburg-Essen, Chair of Fluid Dynamics, Duisburg -

More information

A comparison between two different Flamelet reduced order manifolds for non-premixed turbulent flames

A comparison between two different Flamelet reduced order manifolds for non-premixed turbulent flames 8 th U. S. National Combustion Meeting Organized by the Western States Section of the Combustion Institute and hosted by the University of Utah May 19-22, 2013 A comparison between two different Flamelet

More information

Simulation of Turbulent Lifted Flames and their Transient Propagation

Simulation 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 information

Flamelet Analysis of Turbulent Combustion

Flamelet 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 information

ANSYS Advanced Solutions for Gas Turbine Combustion. Gilles Eggenspieler 2011 ANSYS, Inc.

ANSYS 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 information

A filtered tabulated chemistry model for LES of premixed combustion

A filtered tabulated chemistry model for LES of premixed combustion A filtered tabulated chemistry model for LES of premixed combustion Benoit Fiorina, Ronan Vicquelin, Pierre Auzillon, Nasser Darabiha, Olivier Gicquel, Denis Veynante To cite this version: Benoit Fiorina,

More information

Dynamics of Lean Premixed Systems: Measurements for Large Eddy Simulation

Dynamics 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 information

arxiv: v1 [physics.flu-dyn] 25 Nov 2018

arxiv: 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 information

Lecture 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. 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 information

Overview of Turbulent Reacting Flows

Overview of Turbulent Reacting Flows Overview of Turbulent Reacting Flows Outline Various Applications Overview of available reacting flow models LES Latest additions Example Cases Summary Reacting Flows Applications in STAR-CCM+ Ever-Expanding

More information

Large-eddy simulation of an industrial furnace with a cross-flow-jet combustion system

Large-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 information

A first investigation on using a species reaction mechanism for flame propagation and soot emissions in CFD of SI engines

A first investigation on using a species reaction mechanism for flame propagation and soot emissions in CFD of SI engines A first investigation on using a 1000+ species reaction mechanism for flame propagation and soot emissions in CFD of SI engines F.A. Tap *, D. Goryntsev, C. Meijer, A. Starikov Dacolt International BV

More information

A priori Tabulation of Turbulent Flame Speeds via a Combination of a Stochastic Mixing Model and Flamelet Generated Manifolds 5

A priori Tabulation of Turbulent Flame Speeds via a Combination of a Stochastic Mixing Model and Flamelet Generated Manifolds 5 Konrad-Zuse-Zentrum für Informationstechnik Berlin Takustraße 7 D-14195 Berlin-Dahlem Germany HEIKO SCHMIDT 1 MICHAEL OEVERMANN 2 ROB J.M. BASTIAANS 3 ALAN R. KERSTEIN 4 A priori Tabulation of Turbulent

More information

Representing Turbulence/Chemistry interaction with strained planar premixed

Representing Turbulence/Chemistry interaction with strained planar premixed Representing Turbulence/Chemistry interaction with strained planar premixed flames CES Seminar Andreas Pesch Matr. Nr. 300884 Institut für technische Verbrennung Prof. H. Pitsch Rheinisch-Westfälische

More information

Application of FGM to DNS of premixed turbulent spherical flames

Application of FGM to DNS of premixed turbulent spherical flames Application of FGM to DNS of premixed turbulent spherical flames R.J.M. Bastiaans, G.R.A Groot, J.A. van Oijen and L.P.H. de Goey, Section Combustion Technology, Department of Mechanical Engineering, Eindhoven

More information

A Priori Model for the Effective Lewis Numbers in Premixed Turbulent Flames

A 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 information

NUMERICAL ANALYSIS OF TURBULENT FLAME IN AN ENCLOSED CHAMBER

NUMERICAL 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 information

Advanced Turbulence Models for Emission Modeling in Gas Combustion

Advanced Turbulence Models for Emission Modeling in Gas Combustion 1 Advanced Turbulence Models for Emission Modeling in Gas Combustion Ville Tossavainen, Satu Palonen & Antti Oksanen Tampere University of Technology Funding: Tekes, Metso Power Oy, Andritz Oy, Vattenfall

More information

XXXVIII Meeting of the Italian Section of the Combustion Institute

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 information

EVALUATION 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 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 information

DARS overview, IISc Bangalore 18/03/2014

DARS overview, IISc Bangalore 18/03/2014 www.cd-adapco.com CH2O Temperatur e Air C2H4 Air DARS overview, IISc Bangalore 18/03/2014 Outline Introduction Modeling reactions in CFD CFD to DARS Introduction to DARS DARS capabilities and applications

More information

Large eddy simulations of gaseous flames in gas turbine combustion chambers

Large 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 information

Modeling chemical flame structure and combustion dynamics in LES

Modeling 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 information

LES of the Sandia Flame D Using an FPV Combustion Model

LES of the Sandia Flame D Using an FPV Combustion Model Available online at www.sciencedirect.com ScienceDirect Energy Procedia 82 (2015 ) 402 409 ATI 2015-70th Conference of the ATI Engineering Association LES of the Sandia Flame D Using an FPV Combustion

More information

Experimental 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 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 information

ADVANCED 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 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 information

Lecture 8 Laminar Diffusion Flames: Diffusion Flamelet Theory

Lecture 8 Laminar Diffusion Flames: Diffusion Flamelet Theory Lecture 8 Laminar Diffusion Flames: Diffusion Flamelet Theory 8.-1 Systems, where fuel and oxidizer enter separately into the combustion chamber. Mixing takes place by convection and diffusion. Only where

More information

Turbulent Boundary Layers & Turbulence Models. Lecture 09

Turbulent 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 information

Best Practice Guidelines for Combustion Modeling. Raphael David A. Bacchi, ESSS

Best 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 information

Large eddy simulation of hydrogen-air propagating flames

Large eddy simulation of hydrogen-air propagating flames Loughborough University Institutional Repository Large eddy simulation of hydrogen-air propagating flames This item was submitted to Loughborough University's Institutional Repository by the/an author.

More information

Investigation of ignition dynamics in a H2/air mixing layer with an embedded vortex

Investigation of ignition dynamics in a H2/air mixing layer with an embedded vortex Paper # 070LT-0211 The 8th US National Meeting of the Combustion Institute, Park City, UT, May 19-22, 2013 Investigation of ignition dynamics in a H2/air mixing layer with an embedded vortex S.K. Menon

More information

A Priori Testing of Flamelet Generated Manifolds for Turbulent Partially Premixed Methane/Air Flames

A Priori Testing of Flamelet Generated Manifolds for Turbulent Partially Premixed Methane/Air Flames DOI 10.1007/s10494-009-9223-1 A Priori Testing of Flamelet Generated Manifolds for Turbulent Partially Premixed Methane/Air Flames W. J. S. Ramaekers J. A. van Oijen L. P. H. de Goey Received: 7 September

More information

A G-equation formulation for large-eddy simulation of premixed turbulent combustion

A 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 information

A Ghost-fluid method for large-eddy simulations of premixed combustion in complex geometries

A 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 information

Evaluation of Liquid Fuel Spray Models for Hybrid RANS/LES and DLES Prediction of Turbulent Reactive Flows

Evaluation of Liquid Fuel Spray Models for Hybrid RANS/LES and DLES Prediction of Turbulent Reactive Flows Evaluation of Liquid Fuel Spray Models for Hybrid RANS/LES and DLES Prediction of Turbulent Reactive Flows by Ali Afshar A thesis submitted in conformity with the requirements for the degree of Masters

More information

LES modeling for lifted turbulent jet flames

LES modeling for lifted turbulent jet flames Center for Turbulence Research Proceedings of the Summer Program 1998 83 LES modeling for lifted turbulent jet flames By Luc Vervisch 1 AND Arnaud Trouvé 2 The LES method is an attractive approach for

More information

Simulation of a lean direct injection combustor for the next high speed civil transport (HSCT) vehicle combustion systems

Simulation of a lean direct injection combustor for the next high speed civil transport (HSCT) vehicle combustion systems Center for Turbulence Research Annual Research Briefs 27 241 Simulation of a lean direct injection combustor for the next high speed civil transport (HSCT) vehicle combustion systems By H. El-Asrag, F.

More information

Turbulent Premixed Combustion

Turbulent Premixed Combustion Turbulent Premixed Combustion Combustion Summer School 2018 Prof. Dr.-Ing. Heinz Pitsch Example: LES of a stationary gas turbine velocity field flame 2 Course Overview Part II: Turbulent Combustion Turbulence

More information

Published in: Proceedings of the Fluent Benelux User Group Meeting, 6-7 October 2005, Wavre, Belgium

Published 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 information

DNS of droplet evaporation and combustion in a swirling combustor

DNS of droplet evaporation and combustion in a swirling combustor Center for Turbulence Research Annual Research Briefs 28 DNS of droplet evaporation and combustion in a swirling combustor By K. Luo, O. Desjardins AND H. Pitsch. Motivation and objective Turbulent multi-phase

More information

Thermal NO Predictions in Glass Furnaces: A Subgrid Scale Validation Study

Thermal 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 information

Fluid Dynamics and Balance Equations for Reacting Flows

Fluid Dynamics and Balance Equations for Reacting Flows Fluid Dynamics and Balance Equations for Reacting Flows Combustion Summer School 2018 Prof. Dr.-Ing. Heinz Pitsch Balance Equations Basics: equations of continuum mechanics balance equations for mass and

More information

Survey of Turbulent Combustion Models for Large Eddy Simulations of Propulsive Flowfields

Survey of Turbulent Combustion Models for Large Eddy Simulations of Propulsive Flowfields AIAA SciTech 5-9 January 2015, Kissimmee, Florida 53rd AIAA Aerospace Sciences Meeting AIAA 2015-1379 Survey of Turbulent Combustion Models for Large Eddy Simulations of Propulsive Flowfields Downloaded

More information

Numerical Methods in Aerodynamics. Turbulence Modeling. Lecture 5: Turbulence modeling

Numerical Methods in Aerodynamics. Turbulence Modeling. Lecture 5: Turbulence modeling Turbulence Modeling Niels N. Sørensen Professor MSO, Ph.D. Department of Civil Engineering, Alborg University & Wind Energy Department, Risø National Laboratory Technical University of Denmark 1 Outline

More information

Process Chemistry Toolbox - Mixing

Process Chemistry Toolbox - Mixing Process Chemistry Toolbox - Mixing Industrial diffusion flames are turbulent Laminar Turbulent 3 T s of combustion Time Temperature Turbulence Visualization of Laminar and Turbulent flow http://www.youtube.com/watch?v=kqqtob30jws

More information

TURBINE BURNERS: Engine Performance Improvements; Mixing, Ignition, and Flame-Holding in High Acceleration Flows

TURBINE BURNERS: Engine Performance Improvements; Mixing, Ignition, and Flame-Holding in High Acceleration Flows TURBINE BURNERS: Engine Performance Improvements; Mixing, Ignition, and Flame-Holding in High Acceleration Flows Presented by William A. Sirignano Mechanical and Aerospace Engineering University of California

More information

Spray evaporation model sensitivities

Spray evaporation model sensitivities Center for Turbulence Research Annual Research Briefs 20 23 Spray evaporation model sensitivities By Shashank, E. Knudsen AND H. Pitsch. Motivation and objective The energy density of solid- and liquid-phase

More information

Introduction Flares: safe burning of waste hydrocarbons Oilfields, refinery, LNG Pollutants: NO x, CO 2, CO, unburned hydrocarbons, greenhouse gases G

Introduction 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 information

HEAT 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 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 information

Soot formation in turbulent non premixed flames

Soot formation in turbulent non premixed flames Soot formation in turbulent non premixed flames A. Cuoci 1, A. Frassoldati 1, D. Patriarca 1, T. Faravelli 1, E. Ranzi 1, H. Bockhorn 2 1 Dipartimento di Chimica, Materiali e Ing. Chimica, Politecnico

More information

Introduction to Turbulence and Turbulence Modeling

Introduction to Turbulence and Turbulence Modeling Introduction to Turbulence and Turbulence Modeling Part I Venkat Raman The University of Texas at Austin Lecture notes based on the book Turbulent Flows by S. B. Pope Turbulent Flows Turbulent flows Commonly

More information

New sequential combustion technologies for heavy-duty gas turbines

New sequential combustion technologies for heavy-duty gas turbines New sequential combustion technologies for heavy-duty gas turbines Conference on Combustion in Switzerland 07.09.2017 ETH Zurich Nicolas Noiray, Oliver Schulz CAPS Lab D-MAVT ETH Nicolas Noiray 07/09/17

More information

New Developments in Large Eddy Simulation of Complex Flows

New Developments in Large Eddy Simulation of Complex Flows 14th Australasian Fluid Mechanics Conference Adelaide University, Adelaide, Australia 1-14 December 21 New Developments in Large Eddy Simulation of Complex Flows P. Moin Center for Turbulence Research

More information

DEVELOPMENT OF CFD MODEL FOR A SWIRL STABILIZED SPRAY COMBUSTOR

DEVELOPMENT 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

Modeling of Wall Heat Transfer and Flame/Wall Interaction A Flamelet Model with Heat-Loss Effects

Modeling of Wall Heat Transfer and Flame/Wall Interaction A Flamelet Model with Heat-Loss Effects 9 th U. S. National Combustion Meeting Organized by the Central States Section of the Combustion Institute May 17-20, 2015 Cincinnati, Ohio Modeling of Wall Heat Transfer and Flame/Wall Interaction A Flamelet

More information

CFD Analysis of Vented Lean Hydrogen Deflagrations in an ISO Container

CFD 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 information

An evaluation of a conservative fourth order DNS code in turbulent channel flow

An evaluation of a conservative fourth order DNS code in turbulent channel flow Center for Turbulence Research Annual Research Briefs 2 2 An evaluation of a conservative fourth order DNS code in turbulent channel flow By Jessica Gullbrand. Motivation and objectives Direct numerical

More information

DNS of Reacting H 2 /Air Laminar Vortex Rings

DNS 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 information

Lecture 10 Turbulent Combustion: The State of the Art

Lecture 10 Turbulent Combustion: The State of the Art Lecture 10 Turbulent Combustion: The State of the Art 10.-1 Engineering applications are typically turbulent turbulence models These models use systematic mathematical derivations based on the Navier-

More information

LES of an auto-igniting C 2 H 4 flame DNS

LES of an auto-igniting C 2 H 4 flame DNS Center for Turbulence Research Annual Research Briefs 2011 237 LES of an auto-igniting C 2 H 4 flame DNS By E. Knudsen, E. S. Richardson, J. H. Chen AND H. Pitsch 1. Motivation and objectives Large eddy

More information

An Unsteady/Flamelet Progress Variable Method for LES of Nonpremixed Turbulent Combustion

An 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 information

A dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries

A dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries Center for Turbulence Research Annual Research Briefs 2006 41 A dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries By D. You AND P. Moin 1. Motivation

More information

Impact of numerical method on auto-ignition in a temporally evolving mixing layer at various initial conditions

Impact 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 information

Construction of Libraries for Non-Premixed Tabulated Chemistry Combustion Models including Non-Adiabatic Behaviour due to Wall Heat Losses

Construction of Libraries for Non-Premixed Tabulated Chemistry Combustion Models including Non-Adiabatic Behaviour due to Wall Heat Losses Sonderforschungsbereich/Transregio 40 Annual Report 2016 193 Construction of Libraries for Non-Premixed Tabulated Chemistry Combustion Models including Non-Adiabatic Behaviour due to Wall Heat Losses By

More information

TOPICAL PROBLEMS OF FLUID MECHANICS 97

TOPICAL 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 information

Modelling Detailed-Chemistry Effects on Turbulent Diffusion Flames using a Parallel Solution-Adaptive Scheme

Modelling Detailed-Chemistry Effects on Turbulent Diffusion Flames using a Parallel Solution-Adaptive Scheme Modelling Detailed-Chemistry Effects on Turbulent Diffusion Flames using a Parallel Solution-Adaptive Scheme by Pradeep Kumar Jha A thesis submitted in conformity with the requirements for the degree of

More information

LOW TEMPERATURE MODEL FOR PREMIXED METHANE FLAME COMBUSTION

LOW 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 information

A DEDICATED LES EXPERIMENTAL DATABASE FOR THE ASSESSMENT OF LES SGS MODELS: THE PULSATILE JET IMPINGEMENT IN TURBULENT CROSS FLOW

A DEDICATED LES EXPERIMENTAL DATABASE FOR THE ASSESSMENT OF LES SGS MODELS: THE PULSATILE JET IMPINGEMENT IN TURBULENT CROSS FLOW A DEDICATED LES EXPERIMENTAL DATABASE FOR THE ASSESSMENT OF LES SGS MODELS: THE PULSATILE JET IMPINGEMENT IN TURBULENT CROSS FLOW Hubert Baya Toda Energy Applications Techniques IFP Energie Nouvelles Rueil

More information

Numerical Investigation of Ignition Delay in Methane-Air Mixtures using Conditional Moment Closure

Numerical 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 information

An Introduction to Theories of Turbulence. James Glimm Stony Brook University

An Introduction to Theories of Turbulence. James Glimm Stony Brook University An Introduction to Theories of Turbulence James Glimm Stony Brook University Topics not included (recent papers/theses, open for discussion during this visit) 1. Turbulent combustion 2. Turbulent mixing

More information

PDF Modeling and Simulation of Premixed Turbulent Combustion

PDF Modeling and Simulation of Premixed Turbulent Combustion Monte Carlo Methods Appl. Vol. No. (), pp. 43 DOI 5 / MCMA.7. c de Gruyter PDF Modeling and Simulation of Premixed Turbulent Combustion Michael Stöllinger and Stefan Heinz Abstract. The use of probability

More information

Modeling flame brush thickness in premixed turbulent combustion

Modeling flame brush thickness in premixed turbulent combustion Center for Turbulence Research Proceedings of the Summer Program 2006 299 Modeling flame brush thickness in premixed turbulent combustion By E. Knudsen, O. Kurenkov, S. Kim, M. Oberlack AND H. Pitsch Turbulent

More information

Accounting for spray vaporization in turbulent combustion modeling

Accounting for spray vaporization in turbulent combustion modeling Center for Turbulence Research Proceedings of the Summer Program 1998 25 Accounting for spray vaporization in turbulent combustion modeling By J. Réveillon 1 AND L. Vervisch 1 Three dimensional Direct

More information

Numerical Simulation of Hydrogen Gas Turbines using Flamelet Generated Manifolds technique on Open FOAM

Numerical Simulation of Hydrogen Gas Turbines using Flamelet Generated Manifolds technique on Open FOAM Numerical Simulation of Hydrogen Gas Turbines using Flamelet Generated Manifolds technique on Open FOAM Alessio Fancello (M.Sc.) Department of Mechanical Engineering Combustion Technology Technische Universiteit

More information

Insights into Model Assumptions and Road to Model Validation for Turbulent Combustion

Insights into Model Assumptions and Road to Model Validation for Turbulent Combustion Insights into Model Assumptions and Road to Model Validation for Turbulent Combustion Venke Sankaran AFRL/RQR 2015 AFRL/RQR Basic Research Review UCLA Jan 20, 2015 AFTC PA Release# 15011, 16 Jan 2015 Goals

More information

Large-eddy simulation of supersonic reacting flows

Large-eddy simulation of supersonic reacting flows Sonderforschungsbereich/Transregio 4 Proceedings of the Summer Program 23 2 Large-eddy simulation of supersonic reacting flows By G. Ribert, L. Bouheraoua AND P. Domingo CORIA - UMR664 CNRS, INSA and Université

More information

Topics in Other Lectures Droplet Groups and Array Instability of Injected Liquid Liquid Fuel-Films

Topics in Other Lectures Droplet Groups and Array Instability of Injected Liquid Liquid Fuel-Films Lecture Topics Transient Droplet Vaporization Convective Vaporization Liquid Circulation Transcritical Thermodynamics Droplet Drag and Motion Spray Computations Turbulence Effects Topics in Other Lectures

More information

hydrogen auto ignition in a turbulent co flow of heated air with LES and CMC approach.

hydrogen auto ignition in a turbulent co flow of heated air with LES and CMC approach. biblio.ugent.be The UGent Institutional Repository is the electronic archiving and dissemination platform for all UGent research publications. Ghent University has implemented a mandate stipulating that

More information

Budget analysis and model-assessment of the flamelet-formulation: Application to a reacting jet-in-cross-flow

Budget analysis and model-assessment of the flamelet-formulation: Application to a reacting jet-in-cross-flow Center for Turbulence Research Proceedings of the Summer Program 212 397 Budget analysis and model-assessment of the flamelet-formulation: Application to a reacting jet-in-cross-flow By W. L. Chan, Y.

More information

Investigation of ignition dynamics in a H2/air mixing layer with an embedded vortex

Investigation of ignition dynamics in a H2/air mixing layer with an embedded vortex Paper # 070LT-0211 The 8th US National Meeting of the Combustion Institute, Park City, UT, May 19-22, 2013 Investigation of ignition dynamics in a H2/air mixing layer with an embedded vortex S.K. Menon

More information

Numerical Simulations of Hydrogen Auto-ignition in a Turbulent Co-flow of Heated Air with a Conditional Moment Closure

Numerical Simulations of Hydrogen Auto-ignition in a Turbulent Co-flow of Heated Air with a Conditional Moment Closure Numerical Simulations of Hydrogen Auto-ignition in a Turbulent Co-flow of Heated Air with a Conditional Moment Closure I. Stanković*, 1, A. Triantafyllidis, E. Mastorakos, C. Lacor 3 and B. Merci 1, 4

More information

AER1310: TURBULENCE MODELLING 1. Introduction to Turbulent Flows C. P. T. Groth c Oxford Dictionary: disturbance, commotion, varying irregularly

AER1310: TURBULENCE MODELLING 1. Introduction to Turbulent Flows C. P. T. Groth c Oxford Dictionary: disturbance, commotion, varying irregularly 1. Introduction to Turbulent Flows Coverage of this section: Definition of Turbulence Features of Turbulent Flows Numerical Modelling Challenges History of Turbulence Modelling 1 1.1 Definition of Turbulence

More information

S. Kadowaki, S.H. Kim AND H. Pitsch. 1. Motivation and objectives

S. 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 information

FLAME AND EDDY STRUCTURES IN HYDROGEN AIR TURBULENT JET PREMIXED FLAME

FLAME AND EDDY STRUCTURES IN HYDROGEN AIR TURBULENT JET PREMIXED FLAME FLAME AND EDDY STRUCTURES IN HYDROGEN AIR TURBULENT JET PREMIXED FLAME M. Shimura, K. Yamawaki, Y.-S. Shim, M. Tanahashi and T. Miyauchi Department of Mechanical and Aerospace Engineering Tokyo Institute

More information