G equation

Last updated

In Combustion, G equation is a scalar field equation which describes the instantaneous flame position, introduced by Forman A. Williams in 1985 [1] [2] 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. [3] [4] [5]

Contents

Mathematical description

Sources: [6] [7]

The G equation reads as

where

The flame location is given by which can be defined arbitrarily such that is the region of burnt gas and is the region of unburnt gas. The normal vector to the flame is .

Local burning velocity

According to Matalon–Matkowsky–Clavin–Joulin theory, the burning velocity of the stretched flame, for small curvature and small strain, is given by

where

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 with a premixed reactant mixture is fed through the slot with constant velocity , where the coordinate is chosen such that lies at the center of the slot and 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 with a planar conical shape with cone angle . In the steady case, the G equation reduces to

If a separation of the form is introduced, the equation becomes

which upon integration gives

Without loss of generality choose the flame location to be at . Since the flame is attached to the mouth of the slot , the boundary condition is , which can be used to evaluate the constant . Thus the scalar field is

At the flame tip, we have , the flame height is easily determined as

and the flame angle is given by

Using the trigonometric identity , we have

Related Research Articles

<span class="mw-page-title-main">Lorentz force</span> Force acting on charged particles in electric and magnetic fields

In physics, specifically in electromagnetism, the Lorentz force is the combination of electric and magnetic force on a point charge due to electromagnetic fields. A particle of charge q moving with a velocity v in an electric field E and a magnetic field B experiences a force of It says that the electromagnetic force on a charge q is a combination of (1) a force in the direction of the electric field E, and (2) a force at right angles to both the magnetic field B and the velocity v of the charge.

<span class="mw-page-title-main">Laplace's equation</span> Second-order partial differential equation

In mathematics and physics, Laplace's equation is a second-order partial differential equation named after Pierre-Simon Laplace, who first studied its properties. This is often written as or where is the Laplace operator, is the divergence operator, is the gradient operator, and is a twice-differentiable real-valued function. The Laplace operator therefore maps a scalar function to another scalar function.

<span class="mw-page-title-main">Navier–Stokes equations</span> Equations describing the motion of viscous fluid substances

The Navier–Stokes equations are partial differential equations which describe the motion of viscous fluid substances. They were named after French engineer and physicist Claude-Louis Navier and the Irish physicist and mathematician George Gabriel Stokes. They were developed over several decades of progressively building the theories, from 1822 (Navier) to 1842–1850 (Stokes).

<span class="mw-page-title-main">Potential flow</span> Velocity field as the gradient of a scalar function

In fluid dynamics, potential flow or irrotational flow refers to a description of a fluid flow with no vorticity in it. Such a description typically arises in the limit of vanishing viscosity, i.e., for an inviscid fluid and with no vorticity present in the flow.

In mathematics, the Laplace operator or Laplacian is a differential operator given by the divergence of the gradient of a scalar function on Euclidean space. It is usually denoted by the symbols , (where is the nabla operator), or . In a Cartesian coordinate system, the Laplacian is given by the sum of second partial derivatives of the function with respect to each independent variable. In other coordinate systems, such as cylindrical and spherical coordinates, the Laplacian also has a useful form. Informally, the Laplacian Δf (p) of a function f at a point p measures by how much the average value of f over small spheres or balls centered at p deviates from f (p).

In the calculus of variations, a field of mathematical analysis, the functional derivative relates a change in a functional to a change in a function on which the functional depends.

<span class="mw-page-title-main">Green's function</span> Impulse response of an inhomogeneous linear differential operator

In mathematics, a Green's function is the impulse response of an inhomogeneous linear differential operator defined on a domain with specified initial conditions or boundary conditions.

A Newtonian fluid is a fluid in which the viscous stresses arising from its flow are at every point linearly correlated to the local strain rate — the rate of change of its deformation over time. Stresses are proportional to the rate of change of the fluid's velocity vector.

<span class="mw-page-title-main">Euler equations (fluid dynamics)</span> Set of quasilinear hyperbolic equations governing adiabatic and inviscid flow

In fluid dynamics, the Euler equations are a set of partial differential equations governing adiabatic and inviscid flow. They are named after Leonhard Euler. In particular, they correspond to the Navier–Stokes equations with zero viscosity and zero thermal conductivity.

In mathematics, the covariant derivative is a way of specifying a derivative along tangent vectors of a manifold. Alternatively, the covariant derivative is a way of introducing and working with a connection on a manifold by means of a differential operator, to be contrasted with the approach given by a principal connection on the frame bundle – see affine connection. In the special case of a manifold isometrically embedded into a higher-dimensional Euclidean space, the covariant derivative can be viewed as the orthogonal projection of the Euclidean directional derivative onto the manifold's tangent space. In this case the Euclidean derivative is broken into two parts, the extrinsic normal component and the intrinsic covariant derivative component.

A directional derivative is a concept in multivariable calculus that measures the rate at which a function changes in a particular direction at a given point.

In Hamiltonian mechanics, a canonical transformation is a change of canonical coordinates (q, p) → that preserves the form of Hamilton's equations. This is sometimes known as form invariance. Although Hamilton's equations are preserved, it need not preserve the explicit form of the Hamiltonian itself. Canonical transformations are useful in their own right, and also form the basis for the Hamilton–Jacobi equations and Liouville's theorem.

In physics, the Hamilton–Jacobi equation, named after William Rowan Hamilton and Carl Gustav Jacob Jacobi, is an alternative formulation of classical mechanics, equivalent to other formulations such as Newton's laws of motion, Lagrangian mechanics and Hamiltonian mechanics.

In differential geometry, the four-gradient is the four-vector analogue of the gradient from vector calculus.

The Vlasov equation is a differential equation describing time evolution of the distribution function of plasma consisting of charged particles with long-range interaction, such as the Coulomb interaction. The equation was first suggested for the description of plasma by Anatoly Vlasov in 1938 and later discussed by him in detail in a monograph.

In calculus, the Leibniz integral rule for differentiation under the integral sign, named after Gottfried Wilhelm Leibniz, states that for an integral of the form where and the integrands are functions dependent on the derivative of this integral is expressible as where the partial derivative indicates that inside the integral, only the variation of with is considered in taking the derivative.

<span class="mw-page-title-main">Premixed flame</span>

A premixed flame is a flame formed under certain conditions during the combustion of a premixed charge of fuel and oxidiser. Since the fuel and oxidiser—the key chemical reactants of combustion—are available throughout a homogeneous stoichiometric premixed charge, the combustion process once initiated sustains itself by way of its own heat release. The majority of the chemical transformation in such a combustion process occurs primarily in a thin interfacial region which separates the unburned and the burned gases. The premixed flame interface propagates through the mixture until the entire charge is depleted. The propagation speed of a premixed flame is known as the flame speed which depends on the convection-diffusion-reaction balance within the flame, i.e. on its inner chemical structure. The premixed flame is characterised as laminar or turbulent depending on the velocity distribution in the unburned pre-mixture.

<span class="mw-page-title-main">Covariant formulation of classical electromagnetism</span> Ways of writing certain laws of physics

The covariant formulation of classical electromagnetism refers to ways of writing the laws of classical electromagnetism in a form that is manifestly invariant under Lorentz transformations, in the formalism of special relativity using rectilinear inertial coordinate systems. These expressions both make it simple to prove that the laws of classical electromagnetism take the same form in any inertial coordinate system, and also provide a way to translate the fields and forces from one frame to another. However, this is not as general as Maxwell's equations in curved spacetime or non-rectilinear coordinate systems.

In fluid dynamics, the Oseen equations describe the flow of a viscous and incompressible fluid at small Reynolds numbers, as formulated by Carl Wilhelm Oseen in 1910. Oseen flow is an improved description of these flows, as compared to Stokes flow, with the (partial) inclusion of convective acceleration.

In continuum mechanics, a compatible deformation tensor field in a body is that unique tensor field that is obtained when the body is subjected to a continuous, single-valued, displacement field. Compatibility is the study of the conditions under which such a displacement field can be guaranteed. Compatibility conditions are particular cases of integrability conditions and were first derived for linear elasticity by Barré de Saint-Venant in 1864 and proved rigorously by Beltrami in 1886.

References

  1. Williams, F. A. (1985). Turbulent combustion. In The mathematics of combustion (pp. 97-131). Society for Industrial and Applied Mathematics.
  2. Kerstein, Alan R., William T. Ashurst, and Forman A. Williams. "Field equation for interface propagation in an unsteady homogeneous flow field." Physical Review A 37.7 (1988): 2728.
  3. GH Markstein. (1951). Interaction of flow pulsations and flame propagation. Journal of the Aeronautical Sciences, 18(6), 428-429.
  4. Markstein, G. H. (Ed.). (2014). Nonsteady flame propagation: AGARDograph (Vol. 75). Elsevier.
  5. Markstein, G. H., & Squire, W. (1955). On the stability of a plane flame front in oscillating flow. The Journal of the Acoustical Society of America, 27(3), 416-424.
  6. Peters, Norbert. Turbulent combustion. Cambridge university press, 2000.
  7. Williams, Forman A. "Combustion theory." (1985).