Darrieus–Landau instability

The Darrieus–Landau instability or density fingering refers to an instability of chemical fronts propagating into a denser medium, named after Georges Jean Marie Darrieus and Lev Landau.^[1][2] This instability is one of the key 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. 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.^[3][4]

The instability analysis behind the Darrieus–Landau instability considers a planar, premixed flame front subjected to very small perturbations.^[5] It is useful to think of this arrangement as one in which the unperturbed flame is stationary, with the reactants (fuel and oxidizer) directed towards the flame and perpendicular to it with a velocity u1, and the burnt gases leaving the flame also in a perpendicular way but with velocity u2. The analysis assumes that the flow is an incompressible flow, and that the perturbations are governed by the linearized Euler equations and, thus, are inviscid. With these considerations, the main result of this analysis is that, if the density of the burnt gases is less than that of the reactants, which is the case in practice due to the thermal expansion of the gas produced by the combustion process, the flame front is unstable to perturbations of any wavelength. Another result is that the rate of growth of the perturbations is inversely proportional to their wavelength; thus small flame wrinkles (but larger than the characteristic flame thickness) grow faster than larger ones. In practice, however, diffusive and buoyancy effects that are not taken into account by the analysis of Darrieus and Landau may have a stabilizing effect.^[6][7][8][9]

Dispersion relation

If the disturbances to the steady planar flame sheet are of the form ${\displaystyle e^{i\mathbf {k} \cdot \mathbf {x} _{\bot }+\sigma t}}$, where ${\displaystyle \mathbf {x} _{\bot }}$ is the transverse coordinate system that lies on the undisturbed stationary flame sheet, ${\displaystyle t}$ is the time, ${\displaystyle \mathbf {k} }$ is the wavevector of the disturbance and ${\displaystyle \sigma }$ is the temporal growth rate of the disturbance, then the dispersion relation is given by^[10]

${\displaystyle {\frac {\sigma }{S_{L}k}}={\frac {r}{r+1}}\left({\sqrt {1+{\frac {r^{2}-1}{r}}}}-1\right)}$

where ${\displaystyle S_{L}}$ is the laminar burning velocity (or, the flow velocity far upstream of the flame in a frame that is fixed to the flame), ${\displaystyle k=|\mathbf {k} |}$ and ${\displaystyle r=\rho _{u}/\rho _{b}}$ is the ratio of burnt to unburnt gas density. In combustion ${\displaystyle r>1}$ always and therefore the growth rate ${\displaystyle \sigma >0}$ for all wavenumbers. This implies that a plane sheet of flame with a burning velocity ${\displaystyle S_{L}}$ is unstable for all wavenumbers. In fact, Amable Liñán and Forman A. Williams quote in their book^[11][12] that in view of laboratory observations of stable, planar, laminar flames, publication of their theoretical predictions required courage on the part of Darrieus and Landau.

If the buoyancy forces are taken into account (in others words, accounts of Rayleigh–Taylor instability are considered) for planar flames that are perpendicular to the gravity vector, then some level of stability can be anticipated for flames propagating vertically downwards (or flames that held stationary by a vertically upward flow) since in these cases, the denser unburnt gas lies beneath the lighter burnt gas mixture. Of course, flames that are propagating vertically upwards or those that are held stationary by a vertically downward flow, both the Darrieus–Landau mechanism and the Rayleigh–Taylor mechanism contributes to the destabilizing effect. The dispersion relation when buoyance forces are included becomes

${\displaystyle {\frac {\sigma }{S_{L}k}}={\frac {r}{r+1}}\left[{\sqrt {1+\left({\frac {r^{2}-1}{r}}\right)\left(1-{\frac {g}{S_{L}^{2}rk}}\right)}}-1\right]}$

where ${\displaystyle g>0}$ corresponds to gravitational acceleration for flames propagating downwards and ${\displaystyle g<0}$ corresponds to gravitational acceleration for flames propagating upwards. The above dispersion implies that gravity introduces stability for downward propagating flames when ${\displaystyle k^{-1}>l_{b}=S_{L}^{2}r/g}$, where ${\displaystyle l_{b}}$ is a characteristic buoyancy length scale. For small values of ${\displaystyle r-1}$, the growth rate becomes

${\displaystyle {\frac {\sigma }{S_{L}k}}={\frac {1}{2}}(r-1)+\cdots \quad k^{-1}\ll l_{b},}$

${\displaystyle {\frac {\sigma }{S_{L}k}}={\frac {1}{2}}(r-1)\left(1-{\frac {g}{S_{L}^{2}k}}\right)+\cdots \quad k^{-1}\sim l_{b},}$

${\displaystyle {\frac {\sigma }{g/S_{L}}}=-{\frac {1}{2}}(r-1)+\cdots \quad k^{-1}\gg l_{b}.}$

Limitations

Darrieus and Landau's analysis treats the flame as a plane sheet to investigate its stability with the neglect of diffusion effects, whereas in reality, the flame has a definite thickness, say the laminar flame thickness ${\displaystyle k^{-1}\sim \delta _{L}=D_{T}/S_{L}}$, where ${\displaystyle D_{T}}$ is the thermal diffusivity, wherein diffusion effects cannot be neglected. Accounting for the flame structure, as first envisioned by George H. Markstein, are found to stabilize the flames for small wavelengths ${\displaystyle k^{-1}\sim \delta _{L}}$, except when fuel diffusion coefficient and thermal diffusivity differ from each other significantly leading to the so-called (Turing) diffusive-thermal instability.

Darrieus–Landau instability manifests in the range ${\displaystyle \delta _{L}\ll k^{-1}\ll l_{b}}$ for downward propagating flames and ${\displaystyle \delta _{L}\ll k^{-1}}$ for upward propagating flames.

Dispersion relation under Darcy's law

The classical dispersion relation was based on the assumption that the hydrodynamics is governed by Euler equations. In strongly confinement system such as a Hele-Shaw cell or in porous media, the hydrodynamics is however governed by Darcy's law. The dispersion relation based on Darcy's law was derived by J. Daou and P. Rajamanickam.^[13][14] The dispersion relation under Darcy's law reads

${\displaystyle {\frac {\sigma }{S_{L}k}}={\frac {r-1}{1+m}}+{\frac {1-m}{1+m}}{\frac {V}{S_{L}}}-{\frac {r-1}{1+m}}{\frac {\rho _{b}g\kappa _{b}}{S_{L}\mu _{b}}}}$

where ${\displaystyle r=\rho _{u}/\rho _{b}>1}$ is the density ratio, ${\displaystyle m=(\mu _{u}/\kappa _{u})/(\mu _{b}/\kappa _{b})<1}$ is the ratio of friction factor which involves the viscosity ${\displaystyle \mu }$ and the permeability ${\displaystyle \kappa }$ (in Hele-Shaw cells, ${\displaystyle \kappa _{u}=\kappa _{b}=h^{2}/12}$, where ${\displaystyle h}$ is the cell width, so that ${\displaystyle m=\mu _{u}/\mu _{b}}$ is simply the viscosity ratio) and ${\displaystyle V}$ is the speed of a uniform imposed flow. When ${\displaystyle V>0}$, the imposed flow opposes flame propagation and when ${\displaystyle V<0}$, it aids flame propagation. As before, ${\displaystyle g>0}$ corresponds to downward flame propagation and ${\displaystyle g<0}$ to upward flame propagation. The three terms in the above formula, respectively, corresponds to Darrieus–Landau instability (density fingering), Saffman–Taylor instability (viscous fingering) and Rayleigh–Taylor instability (gravity fingering), in the context of Darcy's law. The Saffman–Taylor instability is specific to confined flames and does not exist in unconfined flames.

See also

• Michelson–Sivashinsky equation
• Clavin–Garcia equation