Physics:G equation

In Combustion, G equation is a scalar [math]\displaystyle{ G(\mathbf{x},t) }[/math] field equation which describes the instantaneous flame position, introduced by Forman A. Williams in 1985 in the study of premixed turbulent combustion. The equation is derived based on the Level-set method. The equation was first studied by George H. Markstein, in a restrictive form for the burning velocity.

Mathematical description

Sources:

The G equation reads as

[math]\displaystyle{ \frac{\partial G}{\partial t} + \mathbf{v}\cdot\nabla G = U_L |\nabla G| }[/math]

where

• [math]\displaystyle{ \mathbf{v} }[/math] is the flow velocity field
• [math]\displaystyle{ U_L }[/math] is the local burning velocity

The flame location is given by [math]\displaystyle{ G(\mathbf{x},t)=G_o }[/math] which can be defined arbitrarily such that [math]\displaystyle{ G(\mathbf{x},t)\gt G_o }[/math] is the region of burnt gas and [math]\displaystyle{ G(\mathbf{x},t)\lt G_o }[/math] is the region of unburnt gas. The normal vector to the flame is [math]\displaystyle{ \mathbf{n}=-\nabla G /|\nabla G| }[/math].

Local burning velocity

The burning velocity of the stretched flame can be derived by subtracting suitable terms from the unstretched flame speed, for small curvature and small strain, as given by

[math]\displaystyle{ U_L = S_L - S_L \mathcal{L} \kappa - \mathcal{L} S }[/math]

where

• [math]\displaystyle{ S_L }[/math] is the burning velocity of unstretched flame
• [math]\displaystyle{ S=-\mathbf{n}\cdot\nabla\mathbf{v}\cdot\mathbf{n} }[/math] is the term corresponding to the imposed strain rate on the flame due to the flow field
• [math]\displaystyle{ \mathcal{L} }[/math] is the Markstein length, proportional to the laminar flame thickness [math]\displaystyle{ \delta_L }[/math], the constant of proportionality is Markstein number [math]\displaystyle{ \mathcal{M} }[/math]
• [math]\displaystyle{ \kappa = \nabla\cdot\mathbf{n} = -\frac{\nabla^2 G - \mathbf{n}\cdot\nabla(\mathbf{n}\cdot\nabla G)}{|\nabla G|} }[/math] is the flame curvature, which is positive if the flame front is convex with respect to the unburnt mixture and vice versa.

A simple example - Slot burner

The G equation has an exact expression for a simple slot burner. Consider a two-dimensional planar slot burner of slot width [math]\displaystyle{ b }[/math] with a premixed reactant mixture is fed through the slot with constant velocity [math]\displaystyle{ \mathbf{v}=(0,U) }[/math], where the coordinate [math]\displaystyle{ (x,y) }[/math] is chosen such that [math]\displaystyle{ x=0 }[/math] lies at the center of the slot and [math]\displaystyle{ y=0 }[/math] lies at the location of the mouth of the slot. When the mixture is ignited, a flame develops from the mouth of the slot to certain height [math]\displaystyle{ y=L }[/math] with a planar conical shape with cone angle [math]\displaystyle{ \alpha }[/math]. In the steady case, the G equation reduces to

[math]\displaystyle{ U\frac{\partial G}{\partial y} = U_L \sqrt{\left(\frac{\partial G}{\partial x}\right)^2+ \left(\frac{\partial G}{\partial y}\right)^2} }[/math]

If a separation of the form [math]\displaystyle{ G(x,y) = y + f(x) }[/math] is introduced, the equation becomes

[math]\displaystyle{ U = U_L\sqrt{1+ \left(\frac{\partial f}{\partial x}\right)^2}, \quad \text{or} \quad \frac{\partial f}{\partial x} = \frac{\sqrt{U^2-U_L^2}}{U_L} }[/math]

which upon integration gives

[math]\displaystyle{ f(x) = \frac{\left(U^2-U_L^2\right)^{1/2}}{U_L}|x| + C, \quad \Rightarrow \quad G(x,y) = \frac{\left(U^2-U_L^2\right)^{1/2}}{U_L}|x| + y+ C }[/math]

Without loss of generality choose the flame location to be at [math]\displaystyle{ G(x,y)=G_o=0 }[/math]. Since the flame is attached to the mouth of the slot [math]\displaystyle{ |x| = b/2, \ y=0 }[/math], the boundary condition is [math]\displaystyle{ G(b/2,0)=0 }[/math], which can be used to evaluate the constant [math]\displaystyle{ C }[/math]. Thus the scalar field is

[math]\displaystyle{ G(x,y) = \frac{\left(U^2-U_L^2\right)^{1/2}}{U_L}\left(|x|- \frac{b}{2}\right) + y }[/math]

At the flame tip, we have [math]\displaystyle{ x=0, \ y=L, \ G=0 }[/math], the flame height is easily determined as

[math]\displaystyle{ L = \frac{b\left(U^2-U_L^2\right)^{1/2}}{2U_L} }[/math]

and the flame angle [math]\displaystyle{ \alpha }[/math] is given by

[math]\displaystyle{ \tan \alpha = \frac{b/2}{L} = \frac{U_L}{\left(U^2-U_L^2\right)^{1/2}} }[/math]

Using the trigonometric identity [math]\displaystyle{ \tan^2\alpha = \sin^2\alpha/\left(1-\sin^2\alpha\right) }[/math], we have

[math]\displaystyle{ \sin\alpha = \frac{U_L}{U} }[/math]