Laminar flamelet model

Last updated May 16, 2024

The laminar flamelet model is a mathematical method for modelling turbulent combustion. The laminar flamelet model is formulated specifically as a model for non-premixed combustion ^[1]

Theory

The flamelet concept considers the turbulent flame as an aggregate of thin, laminar (Re < 2000), locally one-dimensional flamelet structures present within the turbulent flow field. Counterflow diffusion flame is a common laminar flame which is used to represent a flamelet in a turbulent flow. Its geometry consists of opposed and axi-symmetric fuel and oxidizer jets. As the distance between the jets is decreased and/or the velocity of the jets is increased, the flame is strained and departs from its chemical equilibrium until it eventually extinguishes. The mass fraction of species and temperature fields can be measured or calculated in laminar counterflow diffusion flame experiments. When calculated, a self-similar solution exists, and the governing equations can be simplified to only one dimension i.e. along the axis of the fuel and oxidizer jets. It is in this direction where complex chemistry calculations can be performed affordably.^[6]

Logic and formulae

To model a non-premixed combustion, governing equations for fluid elements are required. The conservation equation for the species mass fraction is as follows:-

{\partial Y_{k} \over \partial t}={\chi  \over 2}{1 \over Le_{k}}{\biggl (}{\partial ^{2}Y_{k} \over \partial Z^{2}}{\biggr )}+{\omega _{k} \over \rho }

(1)

Le_k → lewis number of kth species and the above formula was derived with keeping constant heat capacity. The energy equation with variable heat capacity:-

{c_{p}}{\partial T \over \partial t}={\chi  \over 2}{\biggl (}{\partial ^{2}h \over \partial Z^{2}}{\biggr )}-{\sum _{k=1}^{N}h_{k}}{\biggl \{}{\chi  \over 2}{\partial ^{2}Y_{k} \over \partial Z^{2}}+{\omega _{k} \over \rho }{\biggr \}}

(2)

As can be seen from above formulas that the mass fraction and temperature are dependent on

1. Mixture fraction Z

2. Scalar dissipation χ

3. Time

Many times we neglect the unsteady terms in above equation and assume the local flame structure having a balance between steady chemical equations and steady diffusion equation which result in Steady Laminar Flamelet Models (SLFM). For this, an average value of χ is computed known as Favre value ^[7]

{{\widetilde {\chi }}_{st}}={{\widetilde {\chi }} \over \int \limits _{0}^{1}{\frac {F(Z)}{F(Z_{st})}}\ P(Z)\,dZ}

(3)

F(Z)=\exp {\bigl (}-2{\bigl [}\operatorname {erfc} ^{-1}(2Z){\bigr ]}^{2}{\bigr )}

(3)

The basic assumption of a SLFM model is that a turbulent flame front behaves locally as a one dimensional, steady and laminar which proves to be a very useful while reducing the situation to a much simpler terms but it does create problems as few of the effects are not accounted for.

Advantages

The advantages of using this combustion model are as follows:

1. They have the advantage of showing strong coupling between chemical reactions and molecular transport.

2. The steady laminar flamelet model is also used to predict chemical non-equilibrium due to aerodynamic straining of the flame by the turbulence.

Disadvantages

The disadvantages of Steady Laminar Flamelet model due to above mentioned reason are:

1. It does not account for the curvature effects which can change the flame structure and is more detrimental while the structure hasn’t reached the quasi-steady state.

2. Such transient effects also arise in turbulent flow, the scalar dissipation experience a sudden change. As the flame structure take time to get stabilize.^[8]

To improve the above SLFM models, few more models has been proposed like Transient laminar flamelet model (TLFM) by Ferreira.

Related Research Articles

In physics, physical chemistry and engineering, fluid dynamics is a subdiscipline of fluid mechanics that describes the flow of fluids—liquids and gases. It has several subdisciplines, including aerodynamics and hydrodynamics. Fluid dynamics has a wide range of applications, including calculating forces and moments on aircraft, determining the mass flow rate of petroleum through pipelines, predicting weather patterns, understanding nebulae in interstellar space and modelling fission weapon detonation.

In combustion, a diffusion flame is a flame in which the oxidizer and fuel are separated before burning. Contrary to its name, a diffusion flame involves both diffusion and convection processes. The name diffusion flame was first suggested by S.P. Burke and T.E.W. Schumann in 1928, to differentiate from premixed flame where fuel and oxidizer are premixed prior to burning. The diffusion flame is also referred to as nonpremixed flame. The burning rate is however still limited by the rate of diffusion. Diffusion flames tend to burn slower and to produce more soot than premixed flames because there may not be sufficient oxidizer for the reaction to go to completion, although there are some exceptions to the rule. The soot typically produced in a diffusion flame becomes incandescent from the heat of the flame and lends the flame its readily identifiable orange-yellow color. Diffusion flames tend to have a less-localized flame front than premixed flames.

A premixed flame is a flame formed under certain conditions during the combustion of a premixed charge of fuel and oxidiser. Since the fuel and oxidiser—the key chemical reactants of combustion—are available throughout a homogeneous stoichiometric premixed charge, the combustion process once initiated sustains itself by way of its own heat release. The majority of the chemical transformation in such a combustion process occurs primarily in a thin interfacial region which separates the unburned and the burned gases. The premixed flame interface propagates through the mixture until the entire charge is depleted. The propagation speed of a premixed flame is known as the flame speed which depends on the convection-diffusion-reaction balance within the flame, i.e. on its inner chemical structure. The premixed flame is characterised as laminar or turbulent depending on the velocity distribution in the unburned pre-mixture.

In combustion engineering and explosion studies, the Markstein number characterizes the effect of local heat release of a propagating flame on variations in the surface topology along the flame and the associated local flame front curvature. The dimensionless Markstein number is defined as:

The Sugden Award is an annual award for contributions to combustion research. The prize is awarded by the British Section of The Combustion Institute for the published paper with at least one British Section member as author, which makes the most significant contribution to combustion research. The prize is named after Sir Morris Sugden.

The eddy break-up model (EBU) is used in combustion engineering. Combustion modeling has a wide range of applications. In most of the combustion systems, fuel and oxygen are separately supplied in the combustion chamber. Due to this, chemical reaction and combustion occur simultaneously in the combustion chamber. However, the rate of the chemical reaction is faster than the rate of mixing fuel and oxygen. Therefore, that rate of combustion is controlled by rate of mixing. Such cases, where formation of pre-mixture is difficult, are called diffusion combustion or diffusion flames.

Combustion models for CFD refers to combustion models for computational fluid dynamics. Combustion is defined as a chemical reaction in which a hydrocarbon fuel reacts with an oxidant to form products, accompanied with the release of energy in the form of heat. Being the integral part of various engineering applications like: internal combustion engines, aircraft engines, rocket engines, furnaces, and power station combustors, combustion manifests itself as a wide domain during the design, analysis and performance characteristics stages of the above-mentioned applications. With the added complexity of chemical kinetics and achieving reacting flow mixture environment, proper modeling physics has to be incorporated during computational fluid dynamic (CFD) simulations of combustion. Hence the following discussion presents a general outline of the various adequate models incorporated with the Computational fluid dynamic code to model the process of combustion.

Flamelet-Generated Manifold (FGM) is a combustion chemistry reduction technique. The approach of FGM is based on the idea that the most important aspects of the internal structure of the flame front should be taken into account. In this view, a low-dimensional chemical manifold is created on the basis of one-dimensional flame structures, including nearly all of the transport and chemical phenomena as observed in three-dimensional flames. In addition, the progress of the flame is generally described by transport equations for a limited number of control variables.

Chemical reaction models transform physical knowledge into a mathematical formulation that can be utilized in computational simulation of practical problems in chemical engineering. Computer simulation provides the flexibility to study chemical processes under a wide range of conditions. Modeling of a chemical reaction involves solving conservation equations describing convection, diffusion, and reaction source for each component species.

The Darrieus–Landau instability or hydrodynamic instability is an instrinsic flame instability that occurs in premixed flames, caused by the density variation due to the thermal expansion of the gas produced by the combustion process. In simple terms, the stability inquires whether a steadily propagating plane sheet with a discontinuous jump in density is stable or not. It was predicted independently by Georges Jean Marie Darrieus and Lev Landau. Yakov Zeldovich notes that Lev Landau generously suggested this problem to him to investigate and Zeldovich however made error in calculations which led Landau himself to complete the work.

In combustion, the Karlovitz number is defined as the ratio of chemical time scale $to Kolmogorov time scale, named after Béla Karlovitz. The number reads as$

Paul Andrews Libby was a professor of mechanical and aerospace engineering at the University of California, San Diego, a specialist in the field of combustion and aerospace engineering.

Liñán diffusion flame theory is a theory developed by Amable Liñán in 1974 to explain the diffusion flame structure using activation energy asymptotics and Damköhler number asymptotics. Liñán used counterflowing jets of fuel and oxidizer to study the diffusion flame structure, analyzing for the entire range of Damköhler number. His theory predicted four different types of flame structure as follows,

Activation energy asymptotics (AEA), also known as large activation energy asymptotics, is an asymptotic analysis used in the combustion field utilizing the fact that the reaction rate is extremely sensitive to temperature changes due to the large activation energy of the chemical reaction.

In combustion, flame stretch is a quantity which measures the amount of stretch of the flame surface due to curvature and due to the outer velocity field strain. The early concept of flame stretch was introduced by Karlovitz in 1953, although the correct definition was introduced two decades later by Forman A. Williams in 1975. George H. Markstein studied flame stretch by treating the flame surface as a hydrodynamic discontinuity. The flame stretch is also discussed by Bernard Lewis and Guenther von Elbe in their book. All these discussions treated flame stretch as an effect of flow velocity gradients. The stretch can be found even if there is no velocity gradient, but due to the flame curvature. So, the definition required a more general formulation and its precise definition was first introduced by Forman A. Williams in 1975 as the ratio of rate of change of flame surface area to the area itself

In the study of diffusion flame, Liñán's equation is a second-order nonlinear ordinary differential equation which describes the inner structure of the diffusion flame, first derived by Amable Liñán in 1974. The equation reads as

Thierry Poinsot, is a French researcher, research director at the CNRS, researcher at the Institute of Fluid Mechanics in Toulouse, scientific advisor at CERFACS and senior research fellow at Stanford University. He has been a member of the French Academy of sciences since 2019.

Mixture fraction is a quantity used in combustion studies that measures the mass fraction of one stream of a mixture formed by two feed streams, one the fuel stream and the other the oxidizer stream. Both the feed streams are allowed to have inert gases. The mixture fraction definition is usually normalized such that it approaches unity in the fuel stream and zero in the oxidizer stream. The mixture-fraction variable is commonly used as a replacement for the physical coordinate normal to the flame surface, in nonpremixed combustion.

David S-K Ting is a Canadian academic, author and researcher. He is a professor of mechanical, automotive & materials engineering at the University of Windsor. He is the founder of the Turbulence & Energy Laboratory.

Diffusive–thermal instability or thermo–diffusive instability is an intrinsic flame instability that occurs both in premixed flames and in diffusion flames and arises because of the difference in the diffusion coefficient values for the fuel and heat transport, characterized by non-unity values of Lewis numbers. The instability mechanism that arises here is the same as in Turing instability explaining chemical morphogenesis, although the mechanism was first discovered in the context of combustion by Yakov Zeldovich in 1944 to explain the cellular structures appearing in lean hydrogen flames. Quantitative stability theory for premixed flames were developed by Gregory Sivashinsky (1977), Guy Joulin and Paul Clavin (1979) and for diffusion flames by Jong S. Kim and Forman A. Williams (1996,1997).

References

↑ Cuenot, B. (2011). The Flamelet Model for Non-Premixed Combustion. In: Echekki, T., Mastorakos, E. (eds) Turbulent Combustion Modeling. Fluid Mechanics and Its Applications, vol 95. Springer, Dordrecht. https://doi.org/10.1007/978-94-007-0412-1_3
↑ Williams, F. A. (1975). Recent advances in theoretical descriptions of turbulent diffusion flames. In Turbulent mixing in nonreactive and reactive flows (pp. 189–208). Springer, Boston, MA.
↑ Peters, N. (1983). Local quenching due to flame stretch and non-premixed turbulent combustion. Combustion Science and Technology, 30(1–6), 1–17.
↑ Peters, N., & Williams, F. A. (1983). Liftoff characteristics of turbulent jet diffusion flames. AIAA journal, 21(3), 423–429.
↑ Peters, N. (1984). Laminar diffusion flamelet models in non-premixed turbulent combustion. Progress in energy and combustion science, 10(3), 319–339.
↑ ANSYS, Fluent. "Laminar flamelet models theory". FLUENT. ANSYS. Archived from the original on 6 November 2014. Retrieved 6 November 2014.
↑ Pfuderer, D.G.; Neuber, A.A.; Fruchtel, G.; Hassel, E.P.; Janicka, J. (1996). "Combustion flame". 106: 301–317.{{cite journal}}: Cite journal requires |journal= (help)
↑ Pitsch, H.; Peters, N. (1998). "Unsteady Flamelet Modelling Of Turbulent Hydrogen-Air Diffusion Flames". Twenty-Seventh Symposium (International) on Combustion/The Combustion Institute. pp. 1057–1064.