Liñán's equation

Short description: Type of ordinary differential 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. The equation reads as

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

subjected to the boundary conditions

[math]\displaystyle{ \begin{align} \zeta\rightarrow -\infty : &\quad \frac{dy}{d\zeta}=-1,\\ \zeta\rightarrow \infty : &\quad \frac{dy}{d\zeta}=1 \end{align} }[/math]

where [math]\displaystyle{ \delta }[/math] is the reduced or rescaled Damköhler number and [math]\displaystyle{ \gamma }[/math] is the ratio of excess heat conducted to one side of the reaction sheet to the total heat generated in the reaction zone. If [math]\displaystyle{ \gamma\gt 0 }[/math], 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 [math]\displaystyle{ \gamma\lt 0 }[/math], 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 [math]\displaystyle{ \gamma\rightarrow 1 }[/math] [math]\displaystyle{ (\gamma\rightarrow -1) }[/math], all the heat is transported to the oxidizer (fuel) side and therefore the flame sustains extremely large amount of fuel (oxidizer) leakage.

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, having arbitrary Lewis numbers.

Existence of solutions

Near the extinction of the diffusion flame, [math]\displaystyle{ \delta }[/math] is order unity. The equation has no solution for [math]\displaystyle{ \delta\lt \delta_E }[/math], where [math]\displaystyle{ \delta_E }[/math] is the extinction Damköhler number. For [math]\displaystyle{ \delta\gt \delta_E }[/math] with [math]\displaystyle{ |\gamma|\lt 1 }[/math], the equation possess two solutions, of which one is an unstable solution. Unique solution exist if [math]\displaystyle{ |\gamma|\gt 1 }[/math] and [math]\displaystyle{ \delta\gt \delta_E }[/math]. The solution is unique for [math]\displaystyle{ \delta\gt \delta_I }[/math], where [math]\displaystyle{ \delta_I }[/math] is the ignition Damköhler number.

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

[math]\displaystyle{ \delta_E = e[(1-\gamma)-(1-\gamma)^2+0.26(1-\gamma)^3 + 0.055(1-\gamma)^4]. }[/math]

Generalized Liñán's equation

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

[math]\displaystyle{ \frac{d^2y}{d\zeta^2} =(y-\zeta)^m (y+\zeta)^ne^{-\delta^{-1/3}(y+\gamma \zeta)} }[/math]

where [math]\displaystyle{ m }[/math] and [math]\displaystyle{ n }[/math] are constant reaction orders of fuel and oxidizer, respectively.

Large Damköhler number limit

In the Burke–Schumann limit, [math]\displaystyle{ \delta\rightarrow\infty }[/math]. Then the equation reduces to

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

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

[math]\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], }[/math]

where [math]\displaystyle{ \mathrm{erf} }[/math] is the error function and

[math]\displaystyle{ \begin{align} \zeta_m &= \frac{n-m}{n+m} y_m,\\ y_m &= \frac{m+n}{2}\left[\frac{2}{\pi m^2n^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{align} }[/math]

Here [math]\displaystyle{ \zeta=\zeta_m }[/math] is the location where [math]\displaystyle{ y(\zeta) }[/math] reaches its minimum value [math]\displaystyle{ y(\zeta_m)=y_m }[/math]. When [math]\displaystyle{ m=n=1 }[/math], [math]\displaystyle{ \zeta_m=0 }[/math], [math]\displaystyle{ y_m=0.8702 }[/math] and [math]\displaystyle{ k=0.6711 }[/math].

See also

• Liñán's diffusion flame theory