Mean Reaction Rate Research Articles

DNS of swirling hydrogen–air premixed flames

Direct numerical simulation is employed to investigate the turbulent flow characteristics and their effect on local flames for mean reaction rate modelling in turbulent swirling premixed flames. Two swirl numbers having significant effects on the formation of a central recirculation zone in the combustor are considered. The large velocity gradients in the higher swirl number case produce high turbulence intensity in a relatively upstream region compared to the lower swirl number case. The conditional Probability Density Functions (PDFs) of the reaction rate and dissipation rates of turbulent kinetic energy and scalar fluctuations are also examined. The PDFs show correlations between the turbulence energy dissipation and reaction rates and between the scalar dissipation and reaction rates, suggesting that the heat and radicals from the hot products trapped in the recirculation zones are mixed with the reactants, not only through scalar dissipation rate (i.e. scalar gradient) but also by small-scale processes of turbulence relevant to turbulent kinetic energy dissipation rate. Therefore, both scalar and velocity gradients have a strong influence on the chemical reactions through mixing of cold reactant and hot products. A conventional flamelet and EDC models are used to estimate the mean reaction rate, and to study the balance between these two mixing mechanisms. Although both models show a qualitative agreement with the DNS results, these models compensate their limitations each other, depending on the local turbulence and thermochemical conditions. A simple approach is proposed to exploit the advantages of these two models by considering the balance of two mixing mechanisms based on the chemical and turbulence time scales. The estimated mean reaction rate using the proposed model is significantly improved for the higher swirl number case, although the estimated value slightly shifts away from the DNS results for the lower swirl number case. The improved modelling estimate and the balance of turbulence and chemical time scales suggest that the locations of intense reaction zones are strongly related to the dissipation rates of both scalar and turbulent kinetic energy.

Effects of Lewis Number on Head on Quenching of Turbulent Premixed Flames: A Direct Numerical Simulation Analysis

The head on quenching of statistically planar turbulent premixed flames by an isothermal inert wall has been analysed using three-dimensional Direct Numerical Simulation (DNS) data for different values of global Lewis number Le(0.8, 1.0 and 1.2) and turbulent Reynolds number R e t . The statistics of head on quenching have been analysed in terms of the wall Peclet number P e (i.e. distance of the flame from the wall normalised by the Zel’dovich flame thickness) and the normalised wall heat flux Φ. It has been found that the maximum (minimum) value of Φ(P e) for the turbulent L e=0.8 cases are greater (smaller) than the corresponding laminar value, whereas both P e and Φ in turbulent cases remain comparable to the corresponding laminar values for L e=1.0 and 1.2. Detailed physical explanations are provided for the observed Le dependences of P e and Φ. The existing closure of mean reaction rate $\overline {\dot {\omega }}$ using the scalar dissipation rate (SDR) in the near wall region has been assessed based on a-priori analysis of DNS data and modifications to the existing closures of mean reaction rate and SDR have been suggested to account for the wall effects in such a manner that the modified closures perform well both near to and away from the wall.

Direct Numerical Simulation of Complex Fuel Combustion with Detailed Chemistry: Physical Insight and Mean Reaction Rate Modeling

Direct numerical simulation (DNS) of freely-propagating premixed flames of a multi-component fuel is performed using a skeletal chemical mechanism with 49 reactions and 15 species. The fuel consists of CO, H2, H2O, CH4, and CO2 in proportions akin to blast furnace gas or a low calorific value syngas. The simulations include low and high turbulence levels to elucidate the effect of turbulence on realistic chemistry flames. The multi-component fuel flame is found to have a more complex structure than most common flames, with individual species’ reaction zones not necessarily overlapping with each other and with a wide heat releasing zone. The species mass fractions and heat release rate show significant scatter, with their conditional average however remaining close to the laminar flame result. Probability density functions of displacement speed, stretch rate, and curvature are near-Gaussian. Five different mean reaction rate closures are evaluated in the RANS context using these DNS data, presenting perhaps the most stringent test to date of the combustion models. Significant quantitative differences are observed in the performance of the models tested, especially for the higher turbulence level case.

Transition from pulled to pushed fronts in premixed turbulent combustion: Theoretical and numerical study

This paper extends a previous theoretical study (Sabelnikov and Lipatnikov, 2013) of the influence of countergradient transport (CGT) on the speed of a statistically stationary, planar, 1D premixed flame that passively propagates in homogenous turbulence in the form of a traveling wave, i.e. retains its mean thickness and structure. While two particular models of the mean rate of product creation were addressed in the previous article, with the shape of the rate as a function of the Favre-averaged combustion progress variable being concave in both cases, the present paper deals with a more general model that subsumes both concave functions and functions with an inflection point, i.e. a point where the function changes from being concave to convex or vice versa. In this more general case, transition from pulled (flame speed is controlled by processes localized to the flame leading edge) to pushed (flame speed is controlled by processes within the entire flame brush) flames can occur both due to interplay of the nonlinear reaction term and a nonlinear convection term associated with CGT and due to the change of the shape of the reaction term in the absence of CGT. Explicit pushed traveling wave solutions to the studied problem are theoretically derived and conditions under that developing flames approach either pushed or pulled traveling wave solution are obtained by analyzing the governing equations at the flame leading edge and invoking the steepness selection criterion which highlights traveling wave with the steepest profile at the leading edge. Other analytical results include conditions for transition from pulled to pushed premixed turbulent flames, dependence of flame speed on the magnitude of the CGT term and the shape of the mean reaction rate, analytical expressions for the mean thickness of the pushed flames and turbulent scalar flux within the pushed flames. All these theoretical findings are validated by results of unsteady numerical simulations of the initial boundary value problem with steep initial wave profiles.

DNS assessment of relation between mean reaction and scalar dissipation rates in the flamelet regime of premixed turbulent combustion

The linear relation between the mean rate of product creation and the mean scalar dissipation rate, derived in the seminal paper by K.N.C. Bray [‘The interaction between turbulence and combustion’, Proceedings of the Combustion Institute, Vol. 17 (1979), pp. 223–233], is the cornerstone for models of premixed turbulent combustion that deal with the dissipation rate in order to close the reaction rate. In the present work, this linear relation is straightforwardly validated by analysing data computed earlier in the 3D Direct Numerical Simulation (DNS) of three statistically stationary, 1D, planar turbulent flames associated with the flamelet regime of premixed combustion. Although the linear relation does not hold at the leading and trailing edges of the mean flame brush, such a result is expected within the framework of Bray's theory. However, the present DNS yields substantially larger (smaller) values of an input parameter cm (or K2 = 1/(2cm − 1)), involved by the studied linear relation, when compared to the commonly used value of cm = 0.7 (or K2 = 2.5). To gain further insight into the issue and into the eventual dependence of cm on mixture composition, the DNS data are combined with the results of numerical simulations of stationary, 1D, planar laminar methane–air flames with complex chemistry, with the results being reported in terms of differently defined combustion progress variables c, i.e. the normalised temperature, density, or mole fraction of CH4, O2, CO2 or H2O. Such a study indicates the dependence of cm both on the definition of c and on the equivalence ratio. Nevertheless, K2 and cm can be estimated by processing the results of simulations of counterpart laminar premixed flames. Similar conclusions were also drawn by skipping the DNS data, but invoking a presumed beta probability density function in order to evaluate cm for the differently defined c's and various equivalence ratios.

Conditional statistics in a planar liquid jet with a second-order chemical reaction

Conditional statistics, conditioned on mixture fraction, are experimentally investigated in a planar liquid jet with a second-order isothermal chemical reaction A+B→R. Reactants A and B are contained in the jet and ambient flows, respectively, and are supplied into the test section under the non-premixed condition. Streamwise velocity, mixture fraction and concentrations of all reactive species are simultaneously measured by using an I-type hot-film probe and an optical fibre probe based on light absorption spectrometry. The cross-streamwise profiles of conditional mean concentrations and conditional mean reaction rate show that the conditional mean concentrations and the conditional mean reaction rate near the jet exit change with the cross-streamwise position, whereas they are independent of the cross-streamwise position in the downstream regions. On the jet centreline, the conditional mean reaction rate has a peak value on the condition that the conditional mean concentration of product R has a peak value. A comparison between the conditional mean reaction rate and the covariance of conditional concentration fluctuations shows that the covariance of conditional concentration fluctuations is so small that a first moment closure model for the conditional mean reaction rate is valid. The conditional mean streamwise velocity is almost linear to the mixture fraction for small mixture fraction fluctuations. The conditional scalar dissipation rate is calculated from the budget of the conditional moment closure equation. The results show that the conditional scalar dissipation rate has a single peak value in the upstream region, whereas in the downstream region, it has two peaks for the large and small mixture fraction values. The conventional mean scalar dissipation rate is calculated from the probability density function of the mixture fraction and the conditional scalar dissipation rate. The conventional mean scalar dissipation rate on the jet centreline decreases in the downstream direction as (x∗)-2.9 (where x∗ is the distance from the virtual origin), which is almost the same as that expected from scaling arguments.

Pore to Core Scale Simulation of the Mass Transfer with Mineral Reaction in Porous Media

Pore Network Model (PNM) is used to simulate mass transfer with mineral reaction in a single phase flow through porous medium which is here a sandstone sample from the reservoir formation of the Pakoslaw gas field. The void space of the porous medium is represented by an idealized geometry of pore-bodies joined by pore-throats. Parameters defining the pore-bodies and the pore-throats distribution are determined by an optimization process aiming to match the experimental Mercury Intrusion Capillary Pressure (MICP) curve and petrophysical properties of the rock such as intrinsic permeability and formation factor. The generated network is used first to simulate the multiphase flow by solving Kirchhoff’s laws. The capillary pressure and relative permeability curves are derived. Then, reactive transport is addressed under asymptotic regime where the solute concentration undergoes an exponential evolution with time. The porosity/permeability relationship and the three phenomenological coefficients of transport, namely the solute velocity, the dispersion and the mean reaction rate are determined as functions of Peclet and Peclet-Damkohler dimensionless numbers. Finally, the role of the dimensionless numbers on the reactive flow properties is highlighted.

Direct Monomer-By-Monomer Observation of Gelsolin-Mediated Actin Filament Nucleation

Actin nucleation is a fundamental biological process involved in several cellular events. Actin filament nucleation requires several hundred nano-molars to micro-molar amounts of polymerizable actin monomers. This requirement makes the impossible use of standard single molecule techniques. Here, we use zero-mode waveguides to reduce the observation volume and monitor, with single molecule resolution, the initial steps in actin filament formation. We immobilized gelsolin, a capping and nucleating protein at the bottom of the waveguides and added up to several micro-molar amounts of fully labeled actin monomers in polymerization conditions. We observe the first steps during filament formation as step-wise increases in detected fluorescence intensity. Mean reaction rates are then estimated for each oligomer size. Arrival time distribution of individual steps reveals the presence of kinetic intermediates. Furthermore, analysis of oligomer size distribution as a function of actin concentration reveals transitions at around 600 and 900 nM where the system dynamics change concomitant with the appearance of filament elongation. The change in dynamics at the transition is associated with an increase in dissociation rates. Such transition is blocked by Latrunculin A, suggesting a functional role for ATP hydrolysis and/or monomer flattening. Our results point towards a gelsolin-mediated actin filament nucleation mechanism where protein-protein interactions trigger ‘nuclei’ instability necessary for filament elongation.

244 スワール予混合火炎における乱流混合の影響

Direct numerical simulation is employed to investigate the turbulent flow characteristics and their effect on local flames for mean reaction rate modelling of turbulent swirling premixed flames. Two swirl number conditions having significant effects on the formation of a central recirculation zone in the combustor are considered. The mean turbulent and heat release rate fields suggest that the heat and radicals from the hot products trapped in the recirculation zones are mixed with the reactants, not only through scalar dissipation rate (i.e.scalar gradient) but also by small-scale processes of turbulence relevant to turbulent kinetic energy dissipation rate. Therefore, both scalar and velocity gradients have a strong influence on the chemical reactions through mixing of cold reactant and hot products. A conventional flamelet and non-flamelet models are used to estimate the mean reaction rate, and to study the balance between these two mixing mechanisms. Although both models show a qualitative agreement with the DNS results, these models compensate their limitations each other, depending on the local turbulence and thermochemical conditions. A simple approach is proposed to exploit the advantages of these two models by considering the balance of two mixing mechanisms based on the chemical and turbulence time scales.

Reaction and internal energy relaxation rates in viscous thermochemically non-equilibrium gas flows

In the present paper, reaction and energy relaxation rates as well as the normal stress are studied for viscous gas flows with vibrational and chemical non-equilibrium. Using the modified Chapman-Enskog method, multi-temperature models based on the Treanor and Boltzmann vibrational distributions are developed for the general case taking into account all kinds of vibrational energy transitions, exchange reactions, dissociation, and recombination. Integral equations specifying the first-order corrections to the normal mean stress and reaction rates are derived, as well as approximate systems of linear equations for their numerical computation. Generalized thermodynamic driving forces associated with all non-equilibrium processes are introduced. It is shown that normal stresses and rates of non-equilibrium processes can be expressed in terms of the same driving forces; the symmetry of kinetic coefficients in these expressions is proven. The developed general model is applied to a particular case of a pure N2 viscous flow with slow VT relaxation. Normal stress and rates of vibrational relaxation are studied for various ratios of vibrational and translational temperatures. The cross effects between different vibrational transitions in viscous flows are evaluated, along with the influence of anharmonicity and flow compressibility on the first-order corrections to the relaxation rate. Limits of validity for the widely used Landau–Teller model of vibrational relaxation are indicated.

The balances of mixing ratios and segregation intensity: a case study from the field (ECHO 2003)

Abstract. An inhomogeneous mixing of reactants causes a reduction of their chemical removal compared to the homogeneously mixed case in turbulent atmospheric flows. This can be described by the intensity of segregation IS being the covariance of the mixing ratios of two species divided by the product of their means. Both terms appear in the balance equation of the mixing ratio and are discussed for the reaction between isoprene and OH for data of the field study ECHO 2003 above a deciduous forest. For most of these data, IS is negatively correlated with the fraction of mean OH mixing ratio reacting with isoprene. IS is also negatively correlated with the isoprene standard deviation. Both findings agree with model results discussed by Patton et al. (2001) and others. The correlation coefficient between OH and isoprene and, therefore, IS increases with increasing mean reaction rate. In addition, the balance equation of the covariance between isoprene and OH is applied as the theoretical framework for the analysis of the same field data. The storage term is small, and, therefore, a diagnostic equation for this covariance can be derived. The chemical reaction term Rij is dominated by the variance of isoprene times the quotient of mixing ratios of OH and isoprene. Based on these findings a new diagnostic equation for IS is formulated. Comparing different terms of this equation, IS and Rij show a relation also to the normalised isoprene standard deviation. It is shown that not only chemistry but also turbulent and convective mixing and advection – considered in a residual term – influence IS. Despite this finding, a detection of the influence of coherent eddy transport above the forest according to Katul et al. (1997) on IS fails, but a relation to the turbulent and advective transport of isoprene variance is determined. The largest values of IS are found for most unstable conditions with increasing buoyant production, confirming qualitatively model predictions by Ouwersloot et al. (2011).

Experimental study on the reaction rate of a second‐order chemical reaction in a planar liquid jet

Instantaneous concentrations of reactive species are simultaneously measured in a planar liquid jet with a second‐order chemical reaction to investigate the statistical properties of the chemical reaction rate and the validity of models which have been proposed for concentration correlation. The jet flow contains the reactant A, and the ambient flow contains the reactant B. The results show that the concentration correlation of the reactants makes a negative contribution to the mean reaction rate, and this contribution is important in the downstream direction. The concentration correlation changes owing to the chemical reaction. The effects of the chemical reaction on the concentration correlation change with the flow location and the Damköhler number. The concentration correlation predicted by the Toor's model and the three‐environment model are compared with the experimental results. The results show that these models fail to accurately estimate the concentration correlation. © 2014 American Institute of Chemical Engineers AIChE J, 60: 3969–3988, 2014

Evaluation of a reduced mechanism for turbulent premixed combustion

In this study, 3D direct numerical simulations of a multi-component fuel consisting of CO,H2,H2O,CO2 and CH4 reacting with air are performed. A freely propagating turbulent premixed stoichiometric flame is simulated for both low and high turbulence conditions i.e., the rms values of turbulent velocity fluctuations normalised by the laminar flame speed are of order 1 and 10. A skeletal mechanism involving 49 reactions and 15 species, and a 5-step reduced mechanism with 9 species, are used in order to evaluate the performance of the reduced mechanism under turbulent conditions. The 5-step mechanism incurs significantly lower computational expenses compared to the skeletal mechanism. The majority of species mean mass fractions and mean reaction rates computed using these two mechanisms are in good agreement with one another. The mean progress variable and heat release rate variations across the flame brush are also recovered by the reduced mechanism. No major differences are observed in flame response to curvature or strain effects induced by turbulence, although some differences are observed in instantaneous flame structure. These differences are studied using a correlation coefficient and detailed analysis suggests that this comes from the fluctuating heat release induced effects in the case with higher turbulence level. Further considerations based on instantaneous reaction rate and local displacement speed are discussed to evaluate the suitability of the reduced mechanism.

Turbulence–chemistry interactions in a heavy-duty compression–ignition engine

The influences of unresolved turbulent fluctuations in composition and temperature (turbulence–chemistry interactions – TCI) on heat release, flame structure, and emissions are explored in unsteady Reynolds-averaged simulations for a heavy-duty diesel engine at four operating conditions. TCI are isolated and quantified by comparing results from a transported probability density function (PDF) method with those from a model that neglects the influence of fluctuations on local mean reaction rates (a well-stirred-reactor – WSR-model), with all other aspects of the modeling being the same (e.g., spray model, gas-phase chemical mechanism, and soot model). The simulations feature standard fuel-spray and turbulence models, skeletal-level gas-phase chemistry, and a semi-empirical two-equation soot model. Computed pressure and heat-release traces, turbulent flame structure, and emissions from the WSR and PDF models show marked differences, with the PDF-model results being in closer agreement with experiment in most cases. The soot results are especially striking. Computed soot levels from the PDF model are within a factor of five of the measured engine-out particulate matter, and computed soot levels from the WSR and PDF models differ by up to several orders of magnitude, with the PDF-model results being in much closer agreement with experiment. The results show that TCI are important in compression–ignition engines, and that accurate accounting for turbulent fluctuations is at least as important as accurate kinetic rate parameters in the gas-phase chemistry and soot models.

Direct numerical simulation and reaction rate modelling of premixed turbulent flames

Direct numerical simulation (DNS) is employed to investigate the flame characteristics and the mean chemical source term ω˙k˜ of hydrogen/air premixed turbulent combustion in the thin reaction zones regime. Three flames with different initial radius and equivalence ratio are considered. The probability density functions (PDFs) of the flame curvature, shape factor and turbulent flame velocity are plotted and compared for different flames. The distributions of the flame curvature and shape factor indicate that the local geometries for both of the flame kernel and planar flame are similar, irrespective of their difference in the global geometry. The modified laminar flamelet PDF (MLF-PDF) shows improvements over the widely used β-PDF for the progress variable distribution. Several closures for the conditional mean reaction rate are integrated with the presumed PDFs and compared with the DNS data. It is demonstrated that the predictions of the first-order reaction rate integrated with the MLF-PDF are satisfactory for H2 but not for the intermediate species OH. The second-order closure integrated with the MLF-PDF is proved to be promising in predicting both of the major and intermediate species.

Reaction Zones and Their Structure in MILD Combustion

Three-dimensional direct numerical simulation (DNS) of turbulent combustion under moderate and intense low-oxygen dilution (MILD) conditions has been carried out inside a cuboid with inflow and outflow boundaries on the upstream and downstream faces, respectively. The initial and inflowing mixture and turbulence fields are constructed carefully to be representative of MILD conditions involving partially mixed pockets of unburned and burned gases. The combustion kinetics are modeled using a skeletal mechanism for methane-air combustion, including non-unity Lewis numbers for species and temperature-dependent transport properties. The DNS data is analyzed to study the MILD reaction zone structure and its behavior. The results show that the instantaneous reaction zones are convoluted and the degree of convolution increases with dilution and turbulence levels. Interactions of reaction zones occur frequently and are spread out in a large portion of the computational domain due to the mixture non-uniformity and high turbulence level. These interactions lead to local thickening of reaction zones yielding an appearance of distributed combustion despite the presence of local thin reaction zones. A canonical MILD flame element, called MIFE, is proposed, which represents the averaged mass fraction variation for major species reasonably well, although a fully representative canonical element needs to include the effect of reaction zone interactions and associated thickening effects on the mean reaction rate.

Simulation of turbulent explosion of hydrogen–air mixtures

Spherically expanding turbulent premixed hydrogen–air flames are computed using the Reynolds-Averaged Navier Stokes (RANS) approach. The mean reaction rate is modelled using unstrained and strained flamelets, and an algebraic model. Since the temperature and mass fraction evolve differently in hydrogen flames because of non-unity Lewis numbers, two reaction progress variables are used in the calculations. The computed turbulent burning velocity is compared to measured values to validate the computational models. The strained flamelets model captured the experimental variations quite well while the other two models overestimated the burning velocity of stoichiometric hydrogen–air flames. All of these models underestimated the burning velocity for a lean flame, which is thermo-diffusively unstable, indicating a need to include these instability effects in turbulent combustion modelling. A comparative analyses of stoichiometric hydrogen– and methane–air flames are also performed to understand their relative behaviour for a given turbulent combustion condition signified by u′/sL0 and Λ/δ. Although the gross behaviour is observed to be similar, there is a significant difference in the normalised turbulent flame speeds, st/sL0, of these flames. A detailed analysis showed that the mean reaction rate increases very sharply, at least by two orders of magnitude, near the leading edge for the hydrogen flames compared to the methane flames. Since this behaviour is controlled not only by turbulence and its interaction with flame but also by chemical kinetics, one must also consider chemical kinetics parameters such as activation temperature, in addition to u′/sL0 and Λ/δ to characterise st/sL0.

Local Strain Rate and Curvature Dependences of Scalar Dissipation Rate Transport in Turbulent Premixed Flames: A Direct Numerical Simulation Analysis

The statistical behaviours of the instantaneous scalar dissipation rateNcof reaction progress variablecin turbulent premixed flames have been analysed based on three-dimensional direct numerical simulation data of freely propagating statistically planar flame and V-flame configurations with different turbulent Reynolds numberRet. The statistical behaviours ofNcand different terms of its transport equation for planar and V-flames are found to be qualitatively similar. The mean contribution of the density-variation termT1is positive, whereas the molecular dissipation term(-D2)acts as a leading order sink. The mean contribution of the strain rate termT2is predominantly negative for the cases considered here. The mean reaction rate contributionT3is positive (negative) towards the unburned (burned) gas side of the flame, whereas the mean contribution of the diffusivity gradient term(D)assumes negative (positive) values towards the unburned (burned) gas side. The local statistical behaviours ofNc,T1,T2,T3,(-D2), andf(D)have been analysed in terms of their marginal probability density functions (pdfs) and their joint pdfs with local tangential strain rateaTand curvaturekm. Detailed physical explanations have been provided for the observed behaviour.

Computations of turbulent lean premixed combustion using conditional moment closure

Conditional Moment Closure (CMC) is a suitable method for predicting scalars such as carbon monoxide with slow chemical time scales in turbulent combustion. Although this method has been successfully applied to non-premixed combustion, its application to lean premixed combustion is rare. In this study the CMC method is used to compute piloted lean premixed combustion in a distributed combustion regime. The conditional scalar dissipation rate of the conditioning scalar, the progress variable, is closed using an algebraic model and turbulence is modelled using the standard k–ɛ model. The conditional mean reaction rate is closed using a first order CMC closure with the GRI-3.0 chemical mechanism to represent the chemical kinetics of methane oxidation. The PDF of the progress variable is obtained using a presumed shape with the Beta function. The computed results are compared with the experimental measurements and earlier computations using the transported PDF approach. The results show reasonable agreement with the experimental measurements and are consistent with the transported PDF computations. When the compounded effects of shear-turbulence and flame are strong, second order closures may be required for the CMC.

Numerical estimation of blowout, flashback, and flame position in MIT micro gas-turbine chamber

Combustion of hydrogen–air mixture has been simulated numerically inside the MIT (Massachusetts Institute of Technology) micro gas-turbine chamber. Blowout, flashback, and flame position have been studied for different equivalence ratios. Some of the considerations in this simulation are applying a 9-species, 19-step hydrogen–air reaction mechanism, thermal coupling of reacting flow and solid structure of the combustor, considering radiation and convection heat loss from the outer surface of the combustor, and exerting physical boundary conditions on 3D geometry of the combustion chamber. To solve the simulating equations for 3D computational fluid dynamics model, finite volume method has been implemented, and parallel processing has been performed on 6 compute nodes. To validate employed simulating models, the simulation results have been compared with experiment results reported from MIT laboratory and also with simulation results obtained by another research team. The comparison shows that using eddy dissipation concept model (EDC) with disabled turbulence productions and turbulent viscosity terms in k and ε transport equations and solving equations with remaining terms can predict range of mass flow for stable combustion much closer to experimental results (more than 200% improvement in simulation results), which implies that it can be considered as a relatively reliable method for modeling mean reaction rate of micro-combustion.