Markstein number

In combustion engineering and explosion studies, the Markstein number characterizes the effect of local heat release of a propagating flame on variations in the surface topology along the flame and the associated local flame front curvature. The dimensionless Markstein number is defined as:


 * $$\mathcal{M} = \frac{\mathcal{L}}{\delta_L}$$

where $$\mathcal{L}$$ is the Markstein length, and $$\delta_L$$ is the characteristic laminar flame thickness. The larger the Markstein length, the greater the effect of curvature on localised burning velocity. It is named after George H. Markstein (1911—2011), who showed that thermal diffusion stabilized the curved flame front and proposed a relation between the critical wavelength for stability of the flame front, called the Markstein length, and the thermal thickness of the flame. Phenomenological Markstein numbers with respect to the combustion products are obtained by means of the comparison between the measurements of the flame radii as a function of time and the results of the analytical integration of the linear relation between the flame speed and either flame stretch rate or flame curvature. The burning velocity is obtained at zero stretch, and the effect of the flame stretch acting upon it is expressed by a Markstein length. Because both flame curvature and aerodynamic strain contribute to the flame stretch rate, there is a Markstein number associated with each of these components.

Clavin–Williams formula
The Markstein number with respect to the unburnt gas mixture was derived by Paul Clavin and Forman A. Williams in 1982, using activation energy asymptotics. The Clavin–Williams formula is given by


 * $$\mathcal{M} = \frac{r}{r-1} \mathcal{J} + \frac{\beta(Le_{\mathrm{eff}}-1)}{2(r-1)} \mathcal{I},$$

where


 * $$\mathcal{J} = \int_1^r \frac{\lambda(\theta)}{\theta}d\theta, \quad \mathcal{I} = \int_1^r \frac{\lambda(\theta)}{\theta} \ln\frac{r-1}{\theta-1}\,d\theta.$$

Here

The function $$\lambda(\theta)$$, in most cases, is simply given by $$\lambda =\theta^m$$, where $$m=0.7$$, in which case, we have


 * $$\mathcal{J} = \frac{1}{m}(r^m-1), \quad \mathcal{I} = \frac{1}{m}(r^m-1)\ln(r-1)-\int_{1}^r \theta^{m-1}\ln(\theta-1) d\theta.$$

In the constant transport coefficient assumption, $$\lambda=1$$, in which case, we have


 * $$\mathcal{J} =\ln r, \quad \mathcal I = -\mathrm{Li_2}[-(r-1)]$$

where $$\mathrm{Li_2}$$ is the dilogarithm function.

Second Markstein number
In general, Markstein number for the curvature effects $$\mathcal{M}_c$$ and strain effects $$\mathcal{M}_s$$ are not same in real flames. In that case, one defines a second Markstein number as


 * $$\mathcal{M}_2=\mathcal{M}_c-\mathcal{M}_s.$$