The Matalon–Matkowsky–Clavin–Joulin theory or MMCJ theory refers to a longwave hydrodynamic model of a premixed flame with a large-amplitude flame wrinkling, developed independently by Moshe Matalon & Bernard J. Matkowsky and Paul Clavin & Guy Joulin, following the pioneering study by Paul Clavin and Forman A. WilliamsClavin, P., & Williams, F. A. (1982). Effects of molecular diffusion and of thermal expansion on the structure and dynamics of premixed flames in turbulent flows of large scale and low intensity. Journal of fluid mechanics, 116, 251-282., Pierre Pelcé and Paul ClavinPelce, P., & Clavin, P. (1988). Influence of hydrodynamics and diffusion upon the stability limits of laminar premixed flames. In Dynamics of curved fronts (pp. 425-443). Academic Press. and by Gregory Sivashinsky. The theory, for the first time, calculated the burning rate of the curved flame that differs from the burning rate of the planar flame due to flame stretch, associated with the flame curvature and the strain imposed on the flame by the flow field. Specifically, the theory unraveled two important results
The usual Rankine–Hugoniot conditions, applicable across the flame, gets a correction due to the inner structure of the premixed flame. The burning-rate of the flame is influenced by the intrinsic curvature (seen by an observer moving with the flame) and tangential straining of the flame due to flow non-uniformities.
Corrections to the interfacial conditions
Kinematic condition and the burning rate formula
The kinematic condition, i.e., G equation, at the flame interface is given by
where is the unit normal to the flame surface (pointing towards the burnt gas side), is the flow velocity field evaluated at the flame surface. According to Matalon–Matkowsky–Clavin–Joulin theory, if and are the laminar burning speed and thickness of a planar flame (and be the corresponding flame residence time), then the burning speed for the curved flame with respect to the unburnt gas is given by
where is the surface divergence of the tangential velocity ; with flow being incompressible outside the flame, . Moreover, and are the two Markstein numbers, associated with the curvature and tangential straining, which are given, according to Clavin–Williams formula by
where
Here, is the ratio of non-dimensional temperature, is the density scaled by its unburnt gas value, is the ratio of density-thermal conductivity product to its value in the unburnt gas, is the Zel'dovich number and || is the effective Lewis number of the deficient reactant (either fuel or oxidizer or a combination of both).
Normal mass flux
Due to the finite thickness of the flame, the normal mass flux condition acquires a correction. Let be the flame-front speed along (towards the unburnt gas) in the laboratory frame. Then the usual Rankine–Hugoniot conditions for normal mass flux across the interface is , where . However, we now arrive at
where
See also G equation Markstein number
References