Liñán's diffusion flame theory explained

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.[1] [2] [3] 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,

Mathematical description

The theory is well explained in the simplest possible model. Thus, assuming a one-step irreversible Arrhenius law for the combustion chemistry with constant density and transport properties and with unity Lewis number reactants, the governing equation for the non-dimensional temperature field

T(y)

in the stagnation point flow reduces to
d2T
dy2

+y

dT
dy

=-DayFyO

-Ta/T
e

,Z=

1erfc\left(
2
y
\sqrt2

\right)

where

Z

is the mixture fraction,

Da

is the Damköhler number,

Ta=E/R

is the activation temperature and the fuel mass fraction and oxidizer mass fraction are scaled with their respective feed stream values, given by

\begin{align} yF&=Z+To-T\\ yO&=(1-Z)/S+To-T \end{align}

with boundary conditions

T(-infty)=T(infty)=To

. Here,

To

is the unburnt temperature profile (frozen solution) and

S

is the stoichiometric parameter (mass of oxidizer stream required to burn unit mass of fuel stream). The four regime are analyzed by trying to solve above equations using activation energy asymptotics and Damköhler number asymptotics. The solution to above problem is multi-valued. Treating mixture fraction

Z

as independent variable reduces the equation to
d2T
dZ2

=-2\pi

y2
e

DayFyO

-Ta/T
e

with boundary conditions

T(0)=T(1)=To

and

y=\sqrt{2}erfc-1(2Z)

.

Extinction Damköhler number

The reduced Damköhler number is defined as follows

\delta=8\pi

2
Z
s
2
y
s
e\left(
2
T
s
Ta

\right)3

-Ta/T
Dae

where

ys=\sqrt{2}erfc-1(2Zs),Zs=1/(S+1)

and

Ts=To+Zs

. The theory predicted an expression for the reduced Damköhler number at which the flame will extinguish, given by

\deltaE=e\left[(1-\gamma)-(1-\gamma)2+0.26(1-\gamma)3+0.055(1-\gamma)4\right]

where

\gamma=1-2(1-\alpha)(1-Zs)

.

See also

Notes and References

  1. Linan, A. (1974). The asymptotic structure of counterflow diffusion flames for large activation energies. Acta Astronautica, 1(7-8), 1007-1039.
  2. Williams, F. A. (1985). Combustion Theory, (1985). Cummings Publ. Co.
  3. 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.