Liñán's equation

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[1]. The equation reads as

${\displaystyle {\frac {d^{2}y}{d\zeta ^{2}}}=(y^{2}-\zeta ^{2})e^{-\delta ^{-1/3}(y+\gamma \zeta )}}$

subjected to the boundary conditions

${\displaystyle {\begin{aligned}\zeta \rightarrow -\infty$ :&\quad {\frac {dy}{d\zeta }}=-1,\\\zeta \rightarrow \infty :&\quad {\frac {dy}{d\zeta }}=1\end{aligned}}}

where ${\displaystyle \delta }$ is the reduced or rescaled Damköhler number and ${\displaystyle \gamma }$ is the ratio of excess heat conducted to one side of the reaction sheet to the total heat generated in the reaction zone. If ${\displaystyle \gamma >0}$, more heat is transported to the oxidizer side, thereby reducing the reaction rate on the oxidizer side (since reaction rate depends on the temperature) and consequently greater amount of fuel will be leaked into the oxidizer side. Whereas, if ${\displaystyle \gamma <0}$, more heat is transported to the fuel side of the diffusion flame, thereby reducing the reaction rate on the fuel side of the flame and increasing the oxidizer leakage into the fuel side. When ${\displaystyle \gamma \rightarrow 1}$ ${\displaystyle (\gamma \rightarrow -1)}$, all the heat is transported to the oxidizer (fuel) side and therefore the flame sustains extremely large amount of fuel (oxidizer) leakage.[2]

The equation is, in some aspects, universal (also called as the canonical equation of the diffusion flame) since although Liñán derived the equation for stagnation point flow, assuming unity Lewis numbers for the reactants, the same equation is found to represent the inner structure for general laminar flamelets[3][4][5], having arbitrary Lewis numbers[6][7][8].

Existence of solutions

Near the extinction of the diffusion flame, ${\displaystyle \delta }$ is order unity. The equation has no solution for ${\displaystyle \delta <\delta _{E}}$, where ${\displaystyle \delta _{E}}$ is the extinction Damköhler number. For ${\displaystyle \delta >\delta _{E}}$ with ${\displaystyle |\gamma |<1}$, the equation possess two solutions, of which one is an unstable solution. Unique solution exist if ${\displaystyle |\gamma |>1}$ and ${\displaystyle \delta >\delta _{E}}$. The solution is unique for ${\displaystyle \delta >\delta _{I}}$, where ${\displaystyle \delta _{I}}$ is the ignition Damköhler number.

Liñán also gave a correlation formula for the extinction Damköhler number, which is increasingly accurate for ${\displaystyle 1-\gamma \ll 1}$,

${\displaystyle \delta _{E}=e[(1-\gamma )-(1-\gamma )^{2}+0.26(1-\gamma )^{3}+0.055(1-\gamma )^{4}].}$

Generalized Liñán's equation

The generalized Liñán's equation is given by

${\displaystyle {\frac {d^{2}y}{d\zeta ^{2}}}=(y-\zeta )^{m}(y+\zeta )^{n}e^{-\delta ^{-1/3}(y+\gamma \zeta )}}$

where ${\displaystyle m}$ and ${\displaystyle n}$ are constant reaction orders of fuel and oxidizer, respectively.

Large Damköhler number limit

In the Burke–Schumann limit, ${\displaystyle \delta \rightarrow \infty }$. Then the equation reduces to

${\displaystyle {\frac {d^{2}y}{d\zeta ^{2}}}=(y-\zeta )^{m}(y+\zeta )^{n},\quad \zeta \rightarrow \pm \infty$ :\,{\frac {dy}{d\zeta }}=\pm 1.}

An approximate solution to this equation was developed by Liñán himself using integral method in 1963 for his thesis[9],

${\displaystyle y(\zeta )=y_{m}+(\zeta -\zeta _{m})\operatorname {erf} [k(\zeta -\zeta _{m})]-{\frac {1}{{\sqrt {\pi }}k}}\left[1-e^{-k^{2}(\zeta -\zeta _{m})^{2}}\right],}$

where ${\displaystyle \mathrm {erf} }$ is the error function and

${\displaystyle {\begin{aligned}\zeta _{m}&={\frac {n-m}{n+m}}y_{m},\\y_{m}&={\frac {m+n}{2}}\left[{\frac {2}{\pi m^{2}n^{2}}}\left({\sqrt {1+{\frac {\pi (m+n)}{2mn}}}}-1\right)\right]^{\frac {1}{m+n+1}},\\k&={\frac {\sqrt {\pi }}{2}}m^{m}n^{n}\left({\frac {2y_{m}}{m+n}}\right)^{m+n}.\end{aligned}}}$

Here ${\displaystyle \zeta =\zeta _{m}}$ is the location where ${\displaystyle y(\zeta )}$ reaches its minimum value ${\displaystyle y(\zeta _{m})=y_{m}}$. When ${\displaystyle m=n=1}$, ${\displaystyle \zeta _{m}=0}$, ${\displaystyle y_{m}=0.8702}$ and ${\displaystyle k=0.6711}$.

See also

• Liñán's diffusion flame theory

References

1. Linan, A. (1974). The asymptotic structure of counterflow diffusion flames for large activation energies. Acta Astronautica, 1(7-8), 1007-1039.
2. Gubernov, V., & Kim, J. S. (2006). On the fast-time oscillatory instabilities of Linan's diffusion-flame regime. Combustion Theory and Modelling, 10(5), 749-770.
3. Peters, N., & Williams, F. A. (1983). Liftoff characteristics of turbulent jet diffusion flames. AIAA journal, 21(3), 423-429.
4. Peters, N. (1983). Local quenching due to flame stretch and non-premixed turbulent combustion. Combustion Science and Technology, 30(1–6), 1–17.
5. Peters, N. (1986). Laminar Flamelet concept in turbulent combustion Twenty-First Symposium (International) on Combustion-The Combustion Institute 1231.
6. Seshadri, K., & Trevino, C. (1989). The influence of the Lewis numbers of the reactants on the asymptotic structure of counterflow and stagnant diffusion flames. Combustion science and technology, 64(4-6), 243-261.
7. Cheatham, S., & Matalon, M. (2000). A general asymptotic theory of diffusion flames with application to cellular instability. Journal of Fluid Mechanics, 414, 105–144.
8. Liñán, A., Martínez-Ruiz, D., Vera, M., & Sánchez, A. L. (2017). The large-activation-energy analysis of extinction of counterflow diffusion flames with non-unity Lewis numbers of the fuel. Combustion and Flame, 175, 91–106.
9. Linan, A. (1963). On the Structure of Laminar Diffusion Flames. Caltech. https://thesis.library.caltech.edu/4833/1/Linan_a_1963.pdf http://oa.upm.es/2108/2/LINAN_MONO_1963_01a.pdf