Shear Layer Dynamics in a Reacting Jet in Crossflow

This study explores the effect of heat release on the growth of the shear layer vortical structures in a reacting jet in crossflow. Jets composed of mixtures of hydrogen, helium and nitrogen were used to independently vary the momentum flux ratio (J), jet to crossflow density ratio (S) and heat release. Velocity fields were obtained from 10?kHz high-speed stereoscopic particle image velocimetry (SPIV) and regions of elevated temperature/combustion products from simultaneous OH planar laser induced fluorescence (OH-PLIF). The shear layer vortices (SLV) originating from instabilities in the windward and leeward shear layers were identified using vortex identification indicator functions in order to track their spatial location and strength. The results show that the asymmetries in shear layer strength between the windward and leeward shear layers are dependent primarily on J, for both reacting and non-reacting flow-fields. The SLV growth rate dependencies on J and S is found to match trends noted by previous studies for non-reacting jets, where SLV growth rates increase with degree of global instability of the JICF. Heat release is also shown to suppress the SLV growth rates relative to non-reacting cases with the same jet parameters. Related to this point, the degree of lifting of the flame also has a significant impact on SLV growth. As flame lifting is directly related to autoignition times, this point shows strong coupling between kinetic rates and jet hydrodynamic stability.     

网格湍流对剪切层中预混火焰结构的影响

本文研究网格湍流对射流剪切层以及建立在其中的预混火焰的影响。利用热线风速仪测量射流的速度场,发现网格湍流使剪切层内湍流强度明显降低,抑制了低频速度脉动,同时增加了湍动能在小尺度脉动上的分配,使湍流更趋于各向同性,这表明网格湍流抑制了剪切层内的大涡和拟序结构。用细丝热电偶测量了火焰温度,结果显示网格湍流使火焰前峰的低频大幅摆动减少,小尺度皱褶增加,火焰区平均温度更高,说明网格湍流有利于剪切层中预混火焰的强化和稳定。     

Autoignition-controlled flame initiation and flame stabilization in a reacting jet in crossflow

This paper describes an analysis of the mechanisms of autoignition-controlled flame initiation and flame stabilization in a nonpremixed jet in crossflows, using simultaneous high-speed (10 kHz) tomographic particle image velocimetry, OH-PLIF and line-of-sight flame emissions. Measurements are conducted on a turbulent, transverse, reacting propane jet issued into a crossflow generated by combustion of natural gas at an equivalence ratio of 0.4 with the crossflow velocity of 10 m/s, the crossflow temperature of 1350 K and the jet momentum flux ratio of 41. While several prior studies have analyzed the lifted character of the flame in similar configurations, we show that several dynamic processes precede the leading edge of the lifted diffusion flame, including formation and evolution of “autoignition kernels”, “flame kernels” and “flame fragments”. “Autoignition kernels”, i.e., discrete compact reaction zones with the peak hydroxyl (OH) fluorescence intensity below that of the diffusion flame, initiate preferably at bulges along the jet periphery where the strain rates and the scalar dissipation rates are lower. The autoignition kernel grows in both size and the OH-fluorescence intensity as it convects downstream. An autoignition kernel transitions into a propagating flame kernel, which quickly gets distorted and elongated in the direction of the principal expansion strain rate to form a flame fragment. Neighboring flame fragments merge with each other and with the downstream diffusion flame via edge-flame propagation. Merging of upstream flame fragments with the downstream diffusion flame results in an upstream advancement of the diffusion-flame front. The diffusion flame front is intrinsically unsteady because of the rather random formation and evolution of autoignition kernels, flame kernels and flame fragments, presumably due to the stochastic velocity, the strain rate and mixture-fraction oscillations.     

A numerical study on the formation of diffusion flame islands in a turbulent hydrogen jet lifted flame

This paper presents a numerical study on the formation of diffusion flame islands in a hydrogen jet lifted flame. A real size hydrogen jet lifted flame is numerically simulated by the DNS approach over a period of about 0.5 ms. The diameter of hydrogen injector is 2 mm, and the injection velocity is 680 m/s. The lifted flame is composed of a stable leading edge flame, a vigorously turbulent inner rich premixed flame, and a number of outer diffusion flame islands. The relatively long-term observation makes it possible to understand in detail the time-dependent flame behavior in rather large time scales, which are as large as the time scale of the leading edge flame unsteadiness. From the observation, the following three findings are obtained concerning the formation of diffusion flame islands. (1) A thin oxygen diffusion layer is developed along the outer boundary of the lifted flame, where the diffusion flame islands burn in a rather flat shape. (2) When a diffusion flame island comes into contact with the fluctuating inner rich premixed flame, combustion is intensified due to an increase in the hydrogen supply by molecular diffusion. This process also works for the production of the diffusion flame islands in the oxygen diffusion layer. (3) When a large unburned gas volume penetrates into the leading edge flame, the structure of the leading edge flame changes. In this transformation process, a diffusion flame island comes near the leading edge flame. The local deficiency of oxygen plays an important role in this production process.     

Modelling of turbulent lifted jet flames using flamelets: a priori assessment and a posteriori validation

This study focuses on the modelling of turbulent lifted jet flames using flamelets and a presumed Probability Density Function (PDF) approach with interest in both flame lift-off height and flame brush structure. First, flamelet models used to capture contributions from premixed and non-premixed modes of the partially premixed combustion in the lifted jet flame are assessed using a Direct Numerical Simulation (DNS) data for a turbulent lifted hydrogen jet flame. The joint PDFs of mixture fraction Z and progress variable c, including their statistical correlation, are obtained using a copula method, which is also validated using the DNS data. The statistically independent PDFs are found to be generally inadequate to represent the joint PDFs from the DNS data. The effects of Z–c correlation and the contribution from the non-premixed combustion mode on the flame lift-off height are studied systematically by including one effect at a time in the simulations used for a posteriori validation. A simple model including the effects of chemical kinetics and scalar dissipation rate is suggested and used for non-premixed combustion contributions. The results clearly show that both Z–c correlation and non-premixed combustion effects are required in the premixed flamelets approach to get good agreement with the measured flame lift-off heights as a function of jet velocity. The flame brush structure reported in earlier experimental studies is also captured reasonably well for various axial positions. It seems that flame stabilisation is influenced by both premixed and non-premixed combustion modes, and their mutual influences.     

Stratified turbulent Bunsen flames: flame surface analysis and flame surface density modelling

In this paper it is investigated whether the Flame Surface Density (FSD) model, developed for turbulent premixed combustion, is also applicable to stratified flames. Direct Numerical Simulations (DNS) of turbulent stratified Bunsen flames have been carried out, using the Flamelet Generated Manifold (FGM) reduction method for reaction kinetics. Before examining the suitability of the FSD model, flame surfaces are characterized in terms of thickness, curvature and stratification.

All flames are in the Thin Reaction Zones regime, and the maximum equivalence ratio range covers 0.1?φ?1.3. For all flames, local flame thicknesses correspond very well to those observed in stretchless, steady premixed flamelets. Extracted curvature radii and mixing length scales are significantly larger than the flame thickness, implying that the stratified flames all burn in a premixed mode. The remaining challenge is accounting for the large variation in (subfilter) mass burning rate.

In this contribution, the FSD model is proven to be applicable for Large Eddy Simulations (LES) of stratified flames for the equivalence ratio range 0.1?φ?1.3. Subfilter mass burning rate variations are taken into account by a subfilter Probability Density Function (PDF) for the mixture fraction, on which the mass burning rate directly depends. A priori analysis point out that for small stratifications (0.4?φ?1.0), the replacement of the subfilter PDF (obtained from DNS data) by the corresponding Dirac function is appropriate. Integration of the Dirac function with the mass burning rate m=m(φ), can then adequately model the filtered mass burning rate obtained from filtered DNS data. For a larger stratification (0.1?φ?1.3), and filter widths up to ten flame thicknesses, a β-function for the subfilter PDF yields substantially better predictions than a Dirac function. Finally, inclusion of a simple algebraic model for the FSD resulted only in small additional deviations from DNS data, thereby rendering this approach promising for application in LES.     

Characterization of multi-regime reaction zones in a piloted inhomogeneous jet flame with local extinction

Gradient free regime identification (GFRI) is applied to 1D Raman/Rayleigh/LIF measurements of temperature and major species from the intermediate velocity case of the Sydney piloted inhomogeneous jet flame series to better understand the structure of reaction zones and the downstream evolution of multi-regime characteristics. The GFRI approach allows local reaction zones to be detected and characterized as premixed, dominantly premixed, multi-regime, dominantly non-premixed, or non-premixed flame structures, based on flame markers (mixture fraction, chemical mode, and heat release rate) derived from the experimental data. The statistics of chemical mode zero-crossings, which mark premixed reaction zones, and the relative populations of flame structures are shown to be sensitive to the state of mixing in the near field of the flame and to the level of local extinction farther downstream. Multi-regime structures, where premixed and non-premixed reaction zones occur in close proximity and both contribute to overall heat release, account for nearly half the total population at streamwise locations within the first several jet diameters. There is a rapid transition within the near field whereby the relative population of non-premixed and dominantly non-premixed structures grows from 0.05 to nearly 0.5, and the population of premixed and dominantly premixed structures decreases correspondingly as fluid entering the reaction zone becomes progressively fuel-rich. Local extinction and re-ignition bring a resurgence in premixed-type structures, many of which occur at fuel-lean conditions. There are also modest populations of multi-regime structures, having chemical mode zero-crossings at lean conditions, which would not exist in a fully burning jet flame.     

Turbulence/flame/wall interactions in non-premixed inclined slot-jet flames impinging at a wall using direct numerical simulation

In the present work, three-dimensional turbulent non-premixed oblique slot-jet flames impinging at a wall were investigated using direct numerical simulation (DNS). Two cases are considered with the Damköhler number (Da) of case A being twice that of case B. A 17 species and 73-step mechanism for methane combustion was employed in the simulations. It was found that flame extinction in case B is more prominent compared to case A. Reignition in the lower branch of combustion for case A occurs when the scalar dissipation rate relaxes, while no reignition occurs in the lower branch for case B due to excessive scalar dissipation rate. A method was proposed to identify the flame quenching edges of turbulent non-premixed flames in wall-bounded flows based on the intersections of mixture fraction and OH mass fraction iso-surfaces. The flame/wall interactions were examined in terms of the quenching distance and the wall heat flux along the quenching edges. There is essentially no flame/wall interaction in case B due to the extinction caused by excessive turbulent mixing. In contrast, significant interactions between flames and the wall are observed in case A. The quenching distance is found to be negatively correlated with wall heat flux as previously reported in turbulent premixed flames. The influence of chemical reactions and wall on flow topologies was identified. The FS/U and FC/U topologies are found near flame edges, and the NNN/U topology appears when reignition occurs. The vortex-dominant topologies, FC/U and FS/S, play an increasingly important role as the jet turbulence develops.     

Temperature evolution,flame behavior and extinction inside compartment with a ceiling vent under the ambient wind

This article presents an experimental and CFD simulation investigation and analysis on temperature evolution and flame behavior inside compartment with a ceiling vent under ambient wind. Experiments were conducted employing a reduced-scale model containing a cubic fire compartment with a ceiling vent under the external wind generated by a wind tunnel. The temperature of windward and leeward inside compartment as well as the flame behavior were recorded for various vent dimensions, heat release rates and wind speeds for a total of 720 test conditions. Results show that there are two types of fire behavior regimes inside the compartment for relatively small- (Bernoulli flow regime) or large vents (oscillatory exchange flow regime). With a relatively small vent, the flame is located at the center of compartment and the temperatures of windward and leeward are almost the same with or without the ambient wind. Then the flame turns to extinct due to lack of oxygen with increasing of fuel supply. However, with a relatively large vent, there is a flame transition from windward to leeward accompanied by the transition of temperature distribution with increasing of fuel supply when subject to ambient wind. Such extinction and transition mechanisms are interpreted by the aid of CFD simulation of global equivalence ratio and flow/oxygen field inside compartment. The critical heat release rate for the occurrence of flame transition decreases with raising wind speed while increases with vent size. Their complex dependence is found to be well represented in terms of a non-dimensional heat release rate as a function of the wind Froude number, employing the vent area-equivalent characteristic diameter as length scale. These new findings facilitate the understanding of the compartment fire evolution with a ceiling vent subject to ambient wind.     

Combustion characteristics of a dual-mode scramjet combustor with cavity flameholder

Combustion characteristics of a laboratory dual-mode ramjet/scramjet combustor were studied experimentally. The combustor consists of a sonic fuel jet injected into a supersonic crossflow upstream of a wall cavity pilot flame. These fundamental components are contained in many dual-mode combustor designs. Experiments were performed with an isolator entrance Mach number of 2.2. Air stagnation temperatures were varied from 1040 to 1490 K, which correspond to flight Mach numbers of 4.3–5.4. Both pure hydrogen and a mixture of hydrogen and ethylene fuels were used. High speed imaging of the flame luminosity was performed along with measurements of the isolator and combustor wall pressures. For ramjet mode operation, two distinct combustion stabilization locations were found for fuel injection a sufficient distance upstream of the cavity. At low T₀, the combustion was anchored at the leading edge of the cavity by heat release in the cavity shear layer. At high T₀, the combustion was stabilized a short distance downstream of the fuel injection jet in the jet-wake. For an intermediate range of T₀, the reaction zone oscillated between the jet-wake and cavity stabilization locations. Wall pressure measurements showed that cavity stabilized combustion was the steadiest, followed by jet-wake stabilized, and the oscillatory case. For fuel injection close to the cavity, a hybrid stabilization mode was found in which the reaction zone locations for the two stabilization modes overlapped. For this hybrid stabilization, cavity fueling rate was an important factor in the steadiness of the flow field. Scramjet mode combustion was found to only exist in the cavity stabilized location for the conditions studied.     

Experimental and numerical study of a conical turbulent partially premixed flame

The structure and dynamics of a turbulent partially premixed methane/air flame in a conical burner were investigated using laser diagnostics and large-eddy simulations (LES). The flame structure inside the cone was characterized in detail using LES based on a two-scalar flamelet model, with the mixture fraction for the mixing field and level-set G-function for the partially premixed flame front propagation. In addition, planar laser induced florescence (PLIF) of CH and chemiluminescence imaging with high speed video were performed through a glass cone. CH and CH₂O PLIF were also used to examine the flame structures above the cone. It is shown that in the entire flame the CH layer remains very thin, whereas the CH₂O layer is rather thick. The flame is stabilized inside the cone a short distance above the nozzle. The stabilization of the flame can be simulated by the triple-flame model but not the flamelet-quenching model. The results show that flame stabilization in the cone is a result of premixed flame front propagation and flow reversal near the wall of the cone which is deemed to be dependent on the cone angle. Flamelet based LES is shown to capture the measured CH structures whereas the predicted CH₂O structure is somewhat thinner than the experiments.     

Critical conditions at liftoff limit of a laminar n-butane-air jet diffusion flame

The stability mechanism of laminar coflow jet diffusion flames in normal gravity has been studied computationally and experimentally. N-butane, the heaviest alkane in a gaseous state at ambient temperature and pressure, is used as the fuel since the reaction mechanism is similar to that of higher (liquid) hydrocarbons. The critical mean n-butane jet and coflowing air velocities at flame stability limits are measured using a small fuel tube burner (0.8 mm inner diameter). The time-dependent, axisymmetric numerical code with a detailed reaction mechanism (58 species and 540 reactions), molecular diffusive transport, and a radiation model, reveals a flame structure. A fuel-lean peak reactivity spot (i.e., reaction kernel), possessing the hybrid nature of diffusion-premixed flame structure at a constant temperature of ≈1560 K, is formed at the flame base and controls the flame stability. In a near-quiescent environment, the flame base resides below the fuel tube exit plane and thereby premixing is limited. As the coflowing air velocity is increased incrementally under a fixed fuel jet velocity, the flame base moves slightly above (≈1 mm) the burner exit and vigorous premixed combustion becomes prevailing. The local heat-release rate at the reaction kernel nearly doubles due to the increased convective oxygen flux (i.e., a blowing effect). The local Damköhler number, newly defined as a ratio of the square root of the local heat-release rate and the local velocity, decreases gradually first and drops abruptly at a critical threshold value and the flame base lifts off from the burner rim. The calculated coflow air velocity at liftoff is ≈0.38 m/s at the fuel jet velocity of 2 m/s, which is consistent with an extrapolated measured value of 0.41 m/s. This work has determined the critical Damköhler number at the stability limit quantitatively, for the first time, for laminar jet diffusion flames.     

湍流分层燃烧的直接数值模拟和小火焰模型研究

湍流分层燃烧广泛应用于工业燃烧装置,但是目前还比较缺乏适用于湍流分层燃烧的高精度数值模型。本文利用直接数值模拟数据库,对高Karlovitz数分层射流火焰的小火焰模型表现进行了先验性评估。考虑了两种小火焰模型,一种是基于自由传播层流预混火焰的小火焰模型M1,另一种是基于分层对冲小火焰的小火焰模型M2。研究发现M1和M2在c-Z空间的结果与直接数值模拟在定性上是一致的。在物理空间,M2对过程变量反应速率脉动值的预测结果要优于M1.     

A generalized flame surface density modelling approach for the auto-ignition of a turbulent non-premixed system

Auto-ignition of turbulent non-premixed systems is encountered in practical devices such as diesel internal combustion engines. It remains a challenge for modellers, as it exhibits specific features such as unsteadiness, flame propagation and combustion far from stoichiometric conditions. In this paper, a two-dimensional DNS database of an igniting H₂/O₂/N₂ mixing layer, including detailed chemistry and transport, is extensively post-processed in order to gain physical insight into the flame structure and dynamics during auto-ignition. The results are used as a framework for the development of a generalized flame surface density modelling approach by integrating the equations over all possible mixture fraction values. The mean reaction rate is split into two contributions: a generalized flame surface density and a mean reaction rate per unit generalized flame surface density. The unsteadiness of the ignition phenomenon is accounted for via a generalized progress variable. Closures for the generalized surface average of the reaction rate and for the generalized progress variable are proposed, and the modelling approach is tested a priori versus the DNS data. The use of a laminar database for the chemistry coupled to the mean turbulent field via the generalized progress variable shows very promising results, capturing the correct ignition delay and the premixed peak in the turbulent mean heat release rate evolution. This allows confidence in future inclusion and validation of this approach in a RANS-CFD code.     

Auto-ignition and flame stabilization of pulsed methane jets in a hot vitiated coflow studied with high-speed laser and imaging techniques

The auto-ignition of a pulsed methane jet issuing into a laminar coflow of hot exhaust products of a lean premixed hydrogen/air flat flame was examined using high-speed laser and optical measurement techniques with frame rates of 5 kHz or more. OH* chemiluminescence was used to determine the downstream location of the first auto-ignition kernel as well as the stabilization height of the steady-state lifted jet flame. OH planar laser-induced fluorescence (PLIF) was used to determine further details of the auto-ignition with a higher spatial resolution. Simultaneous imaging of broadband luminosity from a viewing angle perpendicular to the OH* chemiluminescence was applied, to three-dimensionally reconstruct the ignition kernel location in space and to determine whether the first occurrence of the kernel was within or beyond the PLIF laser sheet. The development and expansion of the jet was characterized by high-speed Schlieren imaging. Statistics have been compiled for both the ignition time as well as the downstream location of the first auto-ignition kernel and the stabilization height of the steady-state lifted jet flame. From the PLIF images it was found that auto-ignition tended to occur at the interface between bulges of the inflowing jet and the coflow. For steady-state conditions, auto-ignition kernels were observed frequently below the flame base, emphasizing that the lifted jet flame is stabilized by auto-ignition.     

Turbulent flame visualization using direct numerical simulation

Combustion phenomena are of high scientific and technological interest, in particular for energy generation and transportation systems. Direct Numerical Simulations (DNS) have become an essential and well established research tool to investigate the structure of turbulent flames, since they do not rely on any approximate turbulence models. In this work two complementary DNS codes are employed to investigate different types of fuels and flame configurations. The code is π³ is a 3-dimensional DNS code using a low-Mach number approximation. Chemistry is described through a tabulation, using two coordinates to enter a database constructed for example with 29 species and 141 reactions for methane combustion. It is used here to investigate the growth of a turbulent premixed flame in a methane-air mixture (Case 1). The second code,Sider is an explicit three-dimensional DNS code solving the fully compressible reactive Navier-Stokes equations, where the chemical processes are computed using a complete reaction scheme, taking into account accurate diffusion properties. It is used here to compute a hydrogen/air turbulent diffusion flame (Case 2), considering 9 chemical species and 38 chemical reactions.     

An experimental study of partially premixed flame sound

The occurrence of oscillating combustion and combustion instability has led to resurgence of interest in the causes, mechanisms, suppression, and control of combustion noise. Noise generated by enclosed flames is of greater practical interest but is more complicated than that by open flames, which itself is not clearly understood. Studies have shown that different modes of combustion, premixed and non-premixed, differ in their sound generation characteristics. However, there is lack of understanding of the region bridging these two combustion modes. This study investigates sound generation by partially premixed flames. Starting from a non-premixed flame, air was gradually added to achieve partial premixing while maintaining the fuel flow rate constant. Methane, ethylene, and ethane partially premixed flames were studied with hydrogen added for flame stabilization. The sound pressure generated by methane partially premixed flames scales with M⁵ compared to M³ for turbulent non-premixed methane flames. Also, the sound pressure generated by partially premixed flames of ethane and ethylene scales as M^4.5. With progressive partial premixing, spectra level increases at all frequencies with a greater increase in the high-frequency region compared to the low-frequency region; flames develop a peak and later a constant level plateau in the low frequency region. The partially premixed flames of methane, ethylene, and ethane generate a similar SPL as a function of equivalence ratio when the fuel volume flow rate is matched. However, when fuel mass flow rate is matched, the ethane and ethylene flames produce a similar SPL, which is lower than that produced by the methane flame.     

Large Eddy Simulations of a piloted lean premix jet flame using finite-rate chemistry

A Large Eddy Simulation (LES) model capable of accurately representing finite-rate chemistry effects in turbulent premixed combustion is presented. The LES computations use finite-rate chemistry and implicit LES combustion modelling to simulate an experimentally well-documented lean-premixed jet flame stabilized by a stoichiometric pilot. The validity of the implicit LES assumption is discussed and criteria are expressed in terms of subgrid scale Damköhler and Karlovitz numbers. Simulation results are compared to experimental data for velocity, temperature and species mass fractions of CH₄, CO and OH. The simulation results highlight the validity and capability of the present approach for the flame and in general the combustion regime examined. A sensitivity analysis to the choice of the finite-rate chemistry mechanism is reported, this analysis indicates that the one and two-step global reaction mechanisms evaluated fail to capture the reaction layer with sufficient accuracy, while a 20-species skeletal mechanism reproduces the experimental observations accurately including the key finite-rate chemistry indicators CO and OH. The LES results are shown to be grid insensitive and that the grid resolution within the bounds examined is far less important compared to the sensitivity of the finite-rate chemistry representation. The results are analyzed in terms of the flame dynamics and it is shown that intense small scale mixing (high Karlovitz number) between the pilot and the jet is an important mechanism for the stabilization of the flame.     

Vortex dynamics in different combustion regions of a cavity-based scramjet

A hybrid RANS/LES study of a cavity-based scramjet was performed and reasonable agreements were found between simulation results and experimental measurements. In the current case, the flame was stabilized by the subsonic cavity shear layer and propagated downstream into the supersonic flow. The vortex dynamic in the flow, mixing, and combustion regions was comparatively investigated. The averaged vorticity in the combustion regions was lower by 55% compared to the mixing region, primarily due to dilatation as a result of the heat release. Furthermore, the combustion zone was decomposed into four regions based on premixed/diffusion flame and subsonic/supersonic combustion. Then the vorticity and its transport in the four regions were compared. The averaged vorticity in the premixed combustion regions was only slightly larger than that in the diffusion combustion regions. However, the averaged heat release rate was nearly 3 times larger in the premixed regions, leading to higher contributions of dilatation and baroclinic torque in the premixed regions, with an overall weak positive impact on the vorticity generation. In the subsonic combustion regions, the vorticity was three times larger than that in the supersonic combustion regions, despite similar heat release rates on average. It could be explained by the relatively large magnitude of dilatation and baroclinic torque in the supersonic flow. Vortex stretching and dilatation were comparable in the supersonic flame but the former became two times larger than the latter in the subsonic flame. Moreover, the baroclinic torque had larger contributions than diffusion in the supersonic flame whereas the opposite trend was found in the subsonic flame. The results highlight that the subsonic combustion regions in the cavity shear layer and near the walls significantly contribute to the vortex dynamics and mixing process, in addition to flame stabilization.     

A multi-zone chemistry mapping approach for direct numerical simulation of auto-ignition and flame propagation in a constant volume enclosure

A direct numerical simulation (DNS) coupling with multi-zone chemistry mapping (MZCM) is presented to simulate flame propagation and auto-ignition in premixed fuel/air mixtures. In the MZCM approach, the physical domain is mapped into a low-dimensional phase domain with a few thermodynamic variables as the independent variables. The approach is based on the fractional step method, in which the flow and transport are solved in the flow time steps whereas the integration of the chemical reaction rates and heat release rate is performed in much finer time steps to accommodate the small time scales in the chemical reactions. It is shown that for premixed mixtures, two independent variables can be sufficient to construct the phase space to achieve a satisfactory mapping. The two variables can be the temperature of the mixture and the specific element mass ratio of H atom for fuels containing hydrogen atoms. An aliasing error in the MZCM is investigated. It is shown that if the element mass ratio is based on the element involved in the most diffusive molecules, the aliasing error of the model can approach zero when the grid in the phase space is refined. The results of DNS coupled with MZCM (DNS-MZCM) are compared with full DNS that integrates the chemical reaction rates and heat release rate directly in physical space. Application of the MZCM to different mixtures of fuel and air is presented to demonstrate the performance of the method for combustion processes with different complexity in the chemical kinetics, transport and flame–turbulence interaction. Good agreement between the results from DNS and DNS-MZCM is obtained for different fuel/air mixtures, including H₂/air, CO/H₂/air and methane/air, while the computational time is reduced by nearly 70%. It is shown that the MZCM model can properly address important phenomena such as differential diffusion, local extinction and re-ignition in premixed combustion.