Spectral kinetic energy transfer in turbulent premixed reacting flows

Simple criterion of importance of laminar flame instabilities in premixed turbulent combustion of mixtures characterized by low Lewis numbers

By (i) highlighting the mitigation effect of strain rates on laminar flame instabilities and (ii) comparing peak growth rates of laminar flame instabilities with strain rates generated by small-scale turbulent eddies, a simple criterion of importance of the influence of the instabilities on an increase in premixed flame surface area in turbulent flows is suggested. The criterion implies that, even in lean hydrogen-air mixtures, laminar flame instabilities can significantly affect the flame area only in weak or moderate turbulence (the Karlovitz number defined using laminar flame speed, thermal flame thickness, and Kolmogorov time scale is on the order of 10 or less under room conditions).

A unified mechanism for unconfined deflagration-to-detonation transition in terrestrial chemical systems and type Ia supernovae

The nature of type Ia supernovae (SNIa)—thermonuclear explosions of white dwarf stars—is an open question in astrophysics. Virtually all existing theoretical models of normal, bright SNIa require the explosion to produce a detonation in order to consume all of stellar material, but the mechanism for the deflagration-to-detonation transition (DDT) remains unclear. We present a unified theory of turbulence-induced DDT that describes the mechanism and conditions for initiating detonation both in unconfined chemical and thermonuclear explosions. The model is validated by using experiments with chemical flames and numerical simulations of thermonuclear flames. We use the developed theory to determine criteria for detonation initiation in the single-degenerate Chandrasekhar-mass SNIa model and show that DDT is almost inevitable at densities of 107 to 108 grams per cubic centimeter.

Surface-averaged quantities in turbulent reacting flows and relevant evolution equations

While quantities conditioned to an isosurface of reaction progress variable c, which characterizes fluid state in a turbulent reacting flow, have been attracting rapidly growing interest in the recent literature, a mathematical and physical framework required for research into such quantities has not yet been elaborated properly. This paper aims at filling two fundamental gaps in this area, i.e., (i) ambiguities associated with a definition of a surface-averaged quantity and (ii) the lack of rigorous equations that describe evolutions of such quantities. In the first (theoretical) part of the paper, (a) analytical relations between differently defined (area-weighted and unweighted) surface-averaged quantities are obtained and differences between them (quantities) are discussed, (b) a unified method for deriving an evolution equation for bulk area-weighted surface-averaged value of a local characteristic ϕ of a turbulent reacting flow is developed, and (c) the method is applied for deriving evolution equations for the bulk area-weighted surface-averaged reaction-surface density |∇c|, local reaction-wave thickness 1/|∇c|, and local displacement speed S_{d}, i.e., the speed of an isosurface of the c(x,t) field with respect to the local flow. In the second (numerical) part of the paper, direct numerical simulation data obtained recently from a highly turbulent reaction wave are analyzed in order to (1) highlight substantial differences between area-weighted and unweighted surface-averaged quantities and (2) show that various terms in the derived evolution equations are amenable to accurate numerical evaluation in spite of appearance of the so-called zero-gradient points [C. H. Gibson, Phys. Fluids 11, 2305 (1968)PFLDAS0031-917110.1063/1.1691820] in a highly turbulent medium. Finally, the obtained analytical and numerical results are used to shed light on the paradox of local flame thinning and broadening which is widely discussed in the turbulent combustion literature.

Direct numerical simulation study of statistically stationary propagation of a reaction wave in homogeneous turbulence

A three-dimensional (3D) direct numerical simulation (DNS) study of the propagation of a reaction wave in forced, constant-density, statistically stationary, homogeneous, isotropic turbulence is performed by solving Navier-Stokes and reaction-diffusion equations at various (from 0.5 to 10) ratios of the rms turbulent velocity U^{'} to the laminar wave speed, various (from 2.1 to 12.5) ratios of an integral length scale of the turbulence to the laminar wave thickness, and two Zeldovich numbers Ze=6.0 and 17.1. Accordingly, the Damköhler and Karlovitz numbers are varied from 0.2 to 25.1 and from 0.4 to 36.2, respectively. Contrary to an earlier DNS study of self-propagation of an infinitely thin front in statistically the same turbulence, the bending of dependencies of the mean wave speed on U^{'} is simulated in the case of a nonzero thickness of the local reaction wave. The bending effect is argued to be controlled by inefficiency of the smallest scale turbulent eddies in wrinkling the reaction-zone surface, because such small-scale wrinkles are rapidly smoothed out by molecular transport within the local reaction wave.