Turbulent Flame Speed Research Articles

DNS of laboratory-scale turbulent premixed counterflow flames under elevated gravity conditions

In the present work, direct numerical simulations (DNSs) of laboratory-scale turbulent premixed counterflow flames are performed to understand the effect of elevated gravity on flame structures. In the DNS, a turbulent jet of CH4/N2/O2 mixture is injected in opposition to a stream of combustion product. Three cases were considered, i.e., one case with normal gravity and two cases with elevated gravity in different directions. The DNS results of the case under normal gravity conditions were found to be in good agreement with the experimental data. The turbulent flame speed is the highest in the case with elevated gravity aligned with the mean flame propagating direction, which is due to the increase in flame surface area and stretch factor. The gravity levels have a substantial influence on the flame extinction characteristics. It was observed that the flame front is pushed toward the direction of the gravity. When the elevated gravity antialigns with the mean flame propagating direction, the flame front travels across the stagnation plane toward the combustion product stream, leading to significant reactant dilution and flame extinction. The flame normal vector preferentially aligns with the most compressive strain rate for various gravity conditions. The tangential strain rate is the greatest (smallest) when elevated gravity antialigns (aligns) with the mean flame propagating direction. The impact of elevated gravity on the tangential strain rate at the flame fronts ultimately leads to a displacement in the flame front location relative to the stagnation plane, thereby influencing flame extinction.

Experimental Study of Turbulent Premixed Flames of Gas-to-Liquids (GTL) Fuel in a Fan-Stirred Combustion Bomb

This study experimentally investigates the turbulent flame speeds (St) of Gas-to-Liquids (GTL) fuel and a 50/50 blend with diesel across turbulence intensities (u') ranging from 0.5 to 3.0 m/s and equivalence ratios (Φ) from 0.7 to 1.3. Experiments were conducted using a cylindrical fan-stirred combustion bomb at an initial temperature (Ti) of 463 K and atmospheric pressure under near homogeneous and isotropic turbulence (HIT) conditions. A pressure transducer measured the St of the propagating GTL flame. High-speed imaging reveals that changes in flame brightness are associated with variations in flame temperature and soot incandescence. Compared to diesel, stoichiometric GTL and the 50/50 diesel-GTL blend showed peak combustion pressure reductions of 8.9% and 4.9%, respectively. Rich diesel fuel and lean GTL exhibited higher pressure rise rates, flame propagation, and St due to their Lewis numbers (Le) being less than one, enhancing flame-turbulence interaction. At u` of 0.5, 1.5, and 3.0 m/s, St for GTL increased by approximately 3.6%, 5.3%, and 2.8%, respectively, when compared to diesel, with these increases observed at Φ < 1.1. Transition from wrinkled to corrugated flamelets with increasing u` supports models predicting flame stability and potential industrial combustion efficiency improvements. Comparisons with numerical St results using the Zimont Turbulent Flame Speed Closure (Zimont TFC) model showed good agreement at lower (u' = 0.5 m/s) and mid-range (u' = 2.0 m/s) turbulence intensities but discrepancies at higher intensities (u' = 3.0 m/s), highlighting the need to refine numerical models for more accurate predictions across all turbulence levels.

Numerical investigation of premixed syngas/air flame evolution in a closed duct with tulip flame formation

In this study, the dynamic behavior of the premixed syngas/air flames in enclosed ducts with obstacles at the blockage ratio of 0.5, as well as without obstacles, was studied via large eddy simulation and the turbulent flame speed closure model. The hydrogen volume fraction (hydrogen content in fuel) was set to 50 %. Results indicate that the simulation reproduced the flame shapes, flame speed and explosion overpressure observed in the experiments. In cases without obstacles, a “tulip” flame with a tulip cusp can be seen, and secondary cusps emerge on tulip lips due to the high-pressure zone ahead of “tulip” flame lips, leading to in a distorted “tulip” flame. In cases with obstacles, the flame develops hemispherical and finger shapes upstream of the obstacles and evolves into a “tulip” flame downstream of the obstacles due to vortices and adverse pressure gradients. The alteration of the flow direction within the unburned and burned zones results in intricate flame shapes. Flame speed and overpressure show a close relationship and exhibit oscillations following the generation of the “tulip” flame. Typically, the velocity, pressure and flow fields adjacent to the flame front can affect dynamic behaviors of premixed syngas/air flames.

Turbulent Burning Velocity of High Hydrogen Flames

Abstract The turbulent burning velocity (ST) is one of the most important combustion properties controlling combustor operability limits, directly influencing blowoff, flashback, and combustion instabilities. Hydrogen has particularly significant influences on the turbulent flame speed. This paper presents new H2/CH4 data of high pressure, high hydrogen turbulent burning velocities. The datasets were designed to address fundamental questions as well as provide engineering/design relevant insights. This paper presents new scaling analysis of fuel composition, pressure, and preheat temperature effects on turbulent burning velocity. We also discuss the importance of considering what is being held constant (temperature, flame speed, Reynolds number, etc.) when one is analyzing these sensitivities. Data show that hydrogen fraction and pressure cause an increase in turbulent flame speed; whether quantified as raw ST,GC or normalized as ST,GC/SL,0 or ST,GC/SL,max (where SL,0 and SL,max are the unstretched and stretched laminar flame speed, respectively). We also propose that observed increases in ST,GC with pressure are due to increases in Reynolds number and not a kinetics/stretch sensitivity effect. With increasing preheat temperature, ST,GC increases while its normalized value (ST,GC/SL,0 and ST,GC/SL,max) can either increase or decrease, depending upon fuel composition. We also show how these sensitivities vary, depending on whether these comparisons are made at constant raw turbulence intensity urms or at constant normalized urms/SL,0 or urms/SL,max.

Turbulent flame propagation of C10 hydrocarbons/air expanding flames: Possible unified correlation based on the Markstein number

In this study, we experimentally analyzed the influence of turbulence and thermal-diffusional effect on the morphology and turbulent flame speed of premixed expanding flames. The turbulent flames of the n-decane/air and decalin/air mixtures were measured using a fan-stirred constant volume combustion bomb generating near-isotropic turbulence at elevated pressures (1, 2, and 5 bar), temperature (443 K), and a wide range of equivalence ratios (0.8–1.6). The results found that the wrinkling degree of the flame front was greater for Ma<0 compared to Ma>0. However, the distribution patterns of the curvature radius of the flame contours exhibited similarities across different equivalence ratios. Furthermore, an increase in turbulence intensity would aggravate the randomness of flame contour distribution. The turbulent expanding flames of n-decane/air and decalin/air mixtures were self-similar under different turbulence intensities and pressures, and this self-similar propagation followed a correlation between the normalized turbulent flame speed and the turbulent flame Reynolds number (ReT,f) to the one-half power. The similarity in normalized turbulent flame speeds was extended from normal alkanes (C4-C8) and isomeric alkanes (C8, C16) to longer straight-chain alkane (n-decane) and cycloalkane (decalin). As the increase of equivalence ratio, their normalized turbulent flame speeds increased nonlinearly due to the thermal-diffusional effect. Additionally, it was found that the classical, effective, and experimental Lewis numbers were confirmed as unsuitable parameters for characterizing the role of thermal-diffusional effect on the turbulent flame speed for hydrocarbon fuels with large molecular weight. Finally, two possible unified correlations were proposed based on the Markstein number, (d〈r〉/dt)/(σSL)=0.178e−0.231MaReT,f1/2 and (d〈r〉/dt)/(σSL)=(0.170−0.0447Ma+0.0067Ma2)ReT,f1/2, which could predict the turbulent flame speeds of hydrocarbon fuels with large molecular weight such as n-alkanes, iso-alkanes, and decalin.

Isolating effects of large and small scale turbulence on thermodiffusively unstable premixed hydrogen flames

Lean turbulent premixed hydrogen/air flames have substantially increased flame speeds, commonly attributed to differential diffusion effects. In this work, the effect of turbulence on lean hydrogen combustion is studied through Direct Numerical Simulation using detailed chemistry and detailed transport. Simulations are conducted at six Karlovitz numbers and four integral length scales. A general expression for the burning efficiency is proposed which depends on the conditional mean chemical source term and gradient of a progress variable, and the amount of superadiabatic burning. At a fixed Karlovitz number, the normalized turbulent flame speed and area both increase almost linearly with the integral length scale ratio. The effect on the mean source term profile is minimal, indicating that the increase in flame speed can solely be attributed to the increase in flame area. At a fixed integral length scale, both the flame speed and area first increase with Karlovitz number before decreasing. The qualitative observations and trends do not change when Soret effects are neglected. Specifically, neglecting Soret diffusion is shown to reduce the flame speed, area, and burning efficiency. At higher Karlovitz numbers, the diffusivity is enhanced due to penetration of turbulence into the reaction zone, significantly dampening differential diffusion effects.Novelty and SignificanceLean premixed hydrogen flames are subject to thermodiffusive instabilities, which can lead to system level instabilities such as flashback. A comprehensive study of the thermodiffusively unstable turbulent flames with detailed chemistry and detailed transport (including Soret diffusion) was conducted across a wide range of turbulent intensities and integral length scales. Varying these two parameters independently was necessary to isolate the effects of large- and small-scale turbulence. Using these results, we propose a general expression for the burning efficiency to explain the relationship between the turbulence intensity, flame speed, and flame area. This work is an important step in developing predictive models which can aid in the design of practical combustion devices.

Effect of pressure and H2 content on global consumption-based turbulent flame speed at gas turbine relevant conditions

This study focuses on the effect of pressure level and H2 content in the CH4–H2 fuel mixture on the global consumption-based turbulent flame speed (ST). The data considered in the analysis are collected for a variety of CH4–H2 fuel compositions (from 50% to 100% H2 vol.), pressure levels ranging from 0.5 to 1.5 MPa, and preheating temperatures up to 623.15 K. Different trends of ST/SL0 (SL0 is the adiabatic unstretched laminar flame speed) are observed as a function of the normalized turbulence intensity (u’/SL0), depending on the H2 content in CH4–H2 fuel blends. Yet, for all the CH4–H2 fuel mixtures considered, ST/SL0 consistently increases with H2 content. Rising pressure levels are responsible for larger ST/SL0, in direct proportion with the H2 content. These trends are mostly caused by the decrease in SL0, which causes changes to the characteristics of the flamelet that becomes more sensitive to hydrodynamic wrinkling through the high molecular diffusivity of H2.

Numerical Investigation of the Operating Characteristics of the Passive and Active Prechamber Jet Ignition in a Natural Gas Engine.

Prechamber jet ignition is a promising technology that enables stable ignition and fast combustion by combining thermal effects, chemical kinetics, and turbulent disturbance. The development and application of the prechamber ignition require a comprehensive and in-depth understanding of the operating characteristics of the prechamber ignition in the real engine working cycle. Therefore, numerical simulations are conducted to explore the operating performance of the prechamber ignition applied to a large-bore natural gas engine in this study. The differences between the passive prechamber (PPRE) and active prechamber (APRE) near the lean burn limit are compared. The results show that the jet ignition performance of the PPRE is hampered by the high residual gas coefficient and lean mixture in the prechamber under lean burn conditions. The ignition mode of the PPRE is similar to torch ignition, and the combustion process in the main chamber is mainly turbulent flame propagation. The ignition mechanism of the APRE is flame jet ignition. The main chamber combustion process presents a two-stage heat release characteristic, which can be subdivided into three phases: the initial flame development phase, the rapid combustion phase, and the late combustion phase. The heat release rate during the initial flame development phase depends on the physical and chemical properties of the prechamber jet and the mixture conditions in the main chamber. During the rapid combustion phase, the flame propagation along the radial direction of the jet largely depends on the time scale of the chemical reaction. The heat release rate depends on the coverage area of the jet, the jet residual momentum, and the turbulent flame speed. During the late combustion phase, the flame propagation is mainly affected by the turbulent flame speed. These results provide theoretical guidance for the subsequent application of prechamber ignition systems in various powertrains.

Direct numerical simulations of pure and partially cracked ammonia/air turbulent premixed jet flames

Ammonia has been identified as a promising fuel to diminish greenhouse gas emission. However, ammonia combustion presents certain challenges including low reactivity and high NO emission. In the present study, three-dimensional direct numerical simulations (DNS) of ammonia/air premixed slot jet flames with varying Karlovitz numbers (Ka) and cracking ratios were performed. Three cases were considered, including two pure ammonia/air flames with different turbulence intensities and one partially cracked ammonia/air flame with high turbulence intensity. The effects of turbulence intensity and partial ammonia cracking on turbulence–flame interactions and NO emission characteristics of the flames were investigated. It was shown that the turbulent flame speed is higher for the flames with high turbulence intensity. In general, the flame displacement speed is negatively correlated with curvature in negative curvature regions, while the correlation is weak in the positive curvature regions for highly turbulent flames. Most flame area is consumed in negatively curved regions and produced in positively curved regions. It was found that the NO mass fraction is higher in the flame with partial ammonia cracking compared to the pure ammonia/air flames. The NO pathway analysis shows that the NH → NO pathway is enhanced, while the NO consumption pathway is suppressed in the partially cracked ammonia/air flame. The NO mass fraction is higher in regions of negative curvature than positive curvature. Interestingly, the NO mass fraction is found to be negatively correlated with the local equivalence ratio, which is consistent in both the DNS and the corresponding laminar premixed flames.

Turbulent partially cracked ammonia/air premixed spherical flames

The combustion of ammonia requires, for most energy conversion systems, a combustion promoter such as hydrogen to guarantee the start-up, stability and combustion efficiency. Partially cracked ammonia (PCA) can provide sufficient hydrogen concentrations to enhance the burning velocity in comparison with pure ammonia. However, little work exists on the use of PCA blends operating under relevant turbulent conditions. To that end the outwardly propagating spherical flame configuration was employed to determine the laminar and turbulent flame propagation characteristics of PCA (NH3/(H2+N2)) and corresponding binary (NH3/H2) mixtures across various turbulent combustion regimes. First, PCA and ammonia-hydrogen blends exhibit similar flame propagation rates under various turbulent intensities, even for the laminar case. The highest turbulent burning velocity was observed at leanest conditions, as opposed to laminar flames which exhibited highest flame speed at conditions above stoichiometry. Under rich conditions, no substantial flame enhancement due to turbulence was measured irrespective of the hydrogen content. This lack of flame enhancement under turbulent conditions is attributed to the effect of preferential diffusion with good agreement observed with trends in measured Markstein numbers. The normalized turbulent flame speed is dominated by the enhanced molecular diffusivity afforded by the presence of hydrogen up to 15 % enrichment, prior to decreasing upon further hydrogen addition under lean and stoichiometric conditions. This ‘bending’ phenomenon may be the contribution of several factors including; the transitioning between combustion regimes associated with low Damköhler numbers (Da) and flame thickening; merging of flamelets due to the presence of ammonia enhancing wrinkling; and combined changes in laminar burning velocity and preferential diffusional behavior. Furthermore, good agreement for turbulent flame speed is observed with a correlation that includes the influence of turbulent stretch (Ka) and non-equidiffusion (Le), with the agreement reducing with decreasing chemical to turbulent time scale ratios (Da << 1).

Investigation of differential diffusion in lean, premixed, hydrogen-enriched swirl flames

Hydrogen combustion is a promising alternative to fossil fuel combustion in an effort to reduce our carbon footprint. However, hydrogen combustion is prone to thermodiffusive instabilities largely dependent on differential diffusion, a phenomenon that can lead to higher probabilities of flashback in industrial burners, given hydrogen’s high reactivity and diffusivity. This paper evaluates low-swirl flames of methane and air enriched with hydrogen to highlight the onset of differential diffusion. Testing was conducted in a fully controllable swirl burner, where bulk velocity Uav = 13 m/s and swirl number S = 0.6 were kept constant for each hydrogen–methane blend to isolate increases in flame surface area from increases in turbulence intensity. Furthermore, each fuel blend of hydrogen and methane is evaluated at the same laminar flame speed of SLo = 0.267 m/s to isolate flame stretch effects on the turbulent burning rate. Combined hydroxyl (OH) PLIF and stereoscopic PIV at the National Research Council of Canada were used to analyze the OH fluorescence in a 2D-3C velocity field for each flame condition. High-speed PIV at McGill University was used to resolve local flame phenomena, such as local flame displacement velocity and flame stretch rate. Using these techniques, it can be observed that the flame displaces axially in response to turbulent flame speed while exhibiting increases in flamefront wrinkling. This increased corrugation due to flame stretch is highlighted in the PDFs of local curvature and κSf and is further evidenced by a shift towards positive curvatures (κ> 0) for increasing H2 volume fraction. This trend suggests that there is a strong correlation with increases in turbulent burning rate and positive curvature as a result of differential diffusion, but it is not necessarily a control mechanism of the most forward propagating points proposed by the theory of leading points.

Machine learning models for the prediction of turbulent combustion speed for hydrogen-natural gas spark ignition engines

The work carried out in this paper focused on “Machine learning models for the prediction of turbulent combustion speed for hydrogen-natural gas spark ignition engines”. The aim of this work is to develop and verify the ability of machine learning models to solve the problem of estimating the turbulent flame speed for a spark-ignition internal combustion engine operating with a hydrogen-natural gas mixture, then evaluate the relevance of these models in relation to the usual approaches. The novelty of this work is the possibility of a direct calculation of turbulent combustion speed with a good precision, using only machine learning model. The obtained models are also compared to each other by considering in turn as a comparison criterion: the precision of the result, calculation time, and the ability to assimilate original data (which has not undergone preprocessing). An important particularity of this work is that the input variables of the machine learning models were chosen among the variables directly measurable experimentally, based on the opinion of experts in combustion in internal combustion engines and not on the usual approaches to dimensionality reduction on a dataset. The data used for this work was taken from a MINSEL 380, a 380-cc single-cylinder engine. The results show that all the machine learning models obtained are significantly faster than the usual approach and Random Forest (R2: R-squared = 0.9939 and RMSE: Root Mean Square Error = 0.4274) gives the best results. With a forecasting accuracy of over 90 %, both approaches can make reasonable predictions for most industrial applications such as designing engine monitoring and control systems, firefighting systems, simulation, and prototyping tools.

Numerical modeling of combustion in gas engines with prechamber ignition

A numerical model is developed to predict the combustion processes in large-bore gas engines featuring prechamber ignition systems and being coupled with Computational Fluid Dynamics (CFD). The proposed combustion model is based on the simultaneous modeling of premixed and partially-premixed flames by an extended progress variable approach. By introducing an additional fuel scalar to track the externally injected fuel of the scavenged prechamber and depending on a gradient criteria of the fuel scalar, regions of locally premixed or partially-premixed state can be differentiated. The reaction rate for premixed combustion is described through a turbulent flame speed closure approach, whereas the partially-premixed combustion is described by a pre-tabulated flamelet chemistry approach. For the validation of the combustion model and its performance in the different flame propagation phases, that is, prechamber flame ignition and main-chamber flame propagation, experimental data of two large-bore gas engines with different operating conditions, prechamber configurations and engine geometries are taken into account. A good agreement of the simulations with the experimental results is shown for the variety of operating conditions and engine configurations. The developed combustion model is able to predict the combustion process in the prechamber as well as the ignition of the main chamber charge by means of the protruding flame jets through the prechamber nozzles. The prechamber ignition system accelerates the early flame phase and hence shortens the burning duration due to the deep-penetrating and turbulence-inducing flame jets in comparison to a conventional spark plug engine.

Experimentally Closing the Balance of Progress of Reaction in Premixed Turbulent Combustion in the Thin Flame Regime

We investigate the possibility of determining the local turbulent flame speed by measuring the individual terms in the balance of a mean progress of reaction variable for the case of a low turbulence methane-air Bunsen flame in the thin flame regime. Velocity distributions and flame edge positions were measured by particle image velocimetry techniques at 3 kHz for a flame stabilized by a surrounding pilot of the same stoichiometry, for a turbulent Reynolds number around 66 and Karlovitz numbers of the order of 4. The conservation equation for mean progress variable was analyzed along different streamlines as a balance of terms expressed as velocities, including terms for convection, turbulent diffusion, mean reaction, and turbulent and molecular diffusion. Each term was estimated from local velocities and flame locations using a thin flame approximation, and their uncertainty was evaluated based on propagation of experimentally measured statistical correlations. The largest terms were the convective and reaction terms, as expected, with smaller roles for turbulent and molecular diffusion across the flame brush. Countergradient diffusion and transition to gradient diffusion were observed across the flame brush. Closure of the balance of terms in the conservation equations using independently measured terms was not consistently achieved across the flame brush within the reckoned uncertainties, arriving at a balance within 20–30% of the absolute value. Testable hypotheses are offered for the possible reasons for the mismatch, including the role of spatial filtering and 3D effects on the reaction rate term. Finally, the experiments identify the inaccuracies in measuring a true local turbulent flame speed, and suggest a consistent methodology to reduce errors in such estimations. This is the first time such a detailed experimental closure is attempted for any configuration. The results suggest that the significant improvements in spatial resolution are necessary for a full closure.

Lean ammonia-fueled engine operation enabled by hydrogen-assisted turbulent jet ignition

Anhydrous ammonia (NH3) use in internal combustion engines represents a zero-carbon energy solution that is fully sustainable if NH3 is generated renewably. An active hydrogen-fueled pre-chamber to induce turbulent jet ignition is investigated in this work as a means to enhance ignition energy and turbulent flame speed in an NH3 fueled engine. The strength of the turbulent jets, and thus their effectiveness in igniting the main-chamber and enhancing combustion, is highly dependent on pre-chamber equivalence ratio and hydrogen fraction. Local pre-chamber mixtures are varied in the present study by investigating a range of pre-mixed intake NH3-air equivalence ratios (ϕ = 0.5–1) under a consistent hydrogen direct injection strategy in the pre-chamber. Additionally, given the knock-resistance of NH3, multiple compression ratios were studied to investigate the impact on efficiency, emissions, and the combustion process. Results show a clear trade-off where leaner intake equivalence ratios enhance the reactivity of the pre-chamber (greater local hydrogen fraction and closer to stoichiometry) but reduce the reactivity of the main-chamber (lean and slow flame speed). Spark timing optimizes the trade-off under a fixed injection strategy; advancing spark provides more time for combustion to occur in the main-chamber but inhibits pre-chamber reactivity for a less energetic ignition of the main chamber. Optimal indicated thermal efficiency and minimum unburned NH3 and N2O emissions occur around 0.7–0.8 equivalence ratio for all compression ratios. Conversely, NOx is highest at these equivalence ratios but could theoretically be eliminated using selective catalytic reduction aftertreatment using the NH3 present in the exhaust.

Heat release and flame scale effects on turbulence dynamics in confined premixed flows

As industry transitions to a net-zero carbon future, turbulent premixed combustion will remain an integral process for power generating gas turbines, aviation engines, and high-speed propulsion due to their ability to minimize pollutant emissions. However, accurately predicting the behavior of a turbulent reacting flow field remains a challenge. To better understand the dynamics of premixed reacting flows, this study experimentally investigates the effects of combustion heat release and flame scales on the evolution of turbulence in a high-speed, confined bluff-body combustor. The combustor is operated across a range of equivalence ratios from 0.7 to 1 to isolate the role of chemical heat release, flame speed, and flame thickness on the evolution of turbulence as the flow progresses from reactants to products. High-speed particle image velocimetry and CH* chemiluminescence imaging systems are simultaneously employed to quantify turbulent flame and flow dynamics. The results notably demonstrate that the flame augments turbulence fluctuations as the flow evolves from reactants to products for all cases, which opposes most simulations of premixed turbulent reactions. Notably, turbulence fluctuations increase monotonically with the heat of combustion and corresponding turbulent flame speed. Spatial profiles of turbulence statistics are conditioned on the mean flame front, and nondimensionalizing the turbulence profiles using laminar flame properties is shown to collapse all conditions onto a single curve. The resulting nondimensional profile confirms that turbulence dynamics scales with the heat of combustion and was used to develop a novel correlation to predict the increase in turbulent fluctuations across the premixed flame. A Reynolds averaged Navier Stokes decomposition is also explored to further characterize the effects of combustion heat release on the dominant mechanisms of turbulent energy transport. The cumulative results can guide modeling capabilities to better predict flame and flow dynamics and accelerate design strategies for premixed turbines with carbon-free fuels.

Modeling of high-speed, methane–air, turbulent combustion, Part I: One-dimensional turbulence modeling with comparison to DNS

The ability of the one-dimensional turbulence (ODT) model to serve as a surrogate direct numerical simulation (DNS) is assessed for highly turbulent flames. The ODT model is applied to freely propagating premixed methane–air flames at Karlovitz numbers 10, 102, 103, and 104, and results are compared with DNS. The ODT model solves the conservation equations for momentum, energy, and species on a one-dimensional domain, which corresponds to a streamwise line of sight spanning the DNS domain. The effects of turbulent advection are modeled via a stochastic process, in which the Kolmogorov and reactive length and time scales are explicitly resolved. Molecular transport and chemical kinetics are concurrently advanced in time. Both the ODT and DNS simulations use a 21-species skeletal chemical model for methane combustion. The accuracy of the ODT model is assessed by comparing its predictions of several key characteristics of the flames for each Karlovitz number tested, including the turbulent flame speed and width and the joint probability density functions (jPDFs) of major and selected minor species as well as the heat release rate conditioned on temperature with the results of DNS under comparable conditions. The ODT model is shown to yield qualitative and quantitative agreement with the DNS data for most of the above flame characteristics. Discrepancies are observed primarily for the jPDFs of several minor species examined. Overall, the ODT approach is shown to be an effective surrogate of DNS, potentially useful for guiding chemical reaction model reduction and for assessing the sensitivities of the flame structure and the burning rate to chemistry under highly turbulent conditions.Novelty and Significance:The direct numerical simulations (DNS) of premixed turbulent methane–air flames presented in this work span a uniquely wide range of turbulent intensities, from relatively modest corresponding to Karlovitz number Ka=10 to ultra-high intensities at Ka=104. This represents virtually the entire range of turbulent intensities that could be encountered in any realistic situation. This is also the first time that such a wide range of conditions is probed for methane in high-fidelity, fully resolved simulations, which use a fully compressible set of flow equations. The one-dimensional turbulence (ODT) model utilizes the same forcing that is present in the DNS enabling a direct comparison between the ODT and DNS. The results show that ODT captures the key features of the DNS results. ODT is shown to be an effective surrogate for DNS and may be useful in guiding chemical reaction model reduction, where many simulations are required.

Development of Machine Learning Models for Studying the Premixed Turbulent Combustion of Gas-To-Liquids (GTL) Fuel Blends

AbstractStudying the spatial and temporal evolution in turbulent flames represents one of the most challenging problems in the combustion community. Based on previous 3D numerical analyses, this study aims to develop data-driven machine learning (ML) models for predicting the flame radius evolution and turbulent flame speeds for diesel, gas-to-liquids (GTL), and their 50/50 blend (by volumetric composition) under different thermodynamic and turbulence operating conditions. Two ML models were developed in this study. Model 1 predicts the variations of the flame radius with time, equivalence ratio, and turbulence intensity, whereas model 2 predicts the variations of the turbulence flame speed with the operating parameters. The k-fold cross-validation technique is used for model training, and the developed neural network-based model is used to investigate the effects of operating parameters on the premixed turbulent flames. In addition, the possible minimum and maximum values of responses at the corresponding operating parameters are found using a genetic algorithm (GA) approach. Model 1 could capture the computational fluid dynamics (CFD) outputs with high precision at different flame radiuses and time instants with a maximum absolute error percentage of 5.46%. For model 2, the maximum absolute error percentage was 6.58%. Overall, this study demonstrates the applicability and promising performance of the proposed ML models, which will be used in subsequent research to analyze turbulent flames a posteriori.

Experimental study of the influence of Lewis number, laminar flame thickness, temperature, and pressure on turbulent flame speed using hydrogen and methane fuels

To experimentally explore the influence of Lewis number Le, laminar flame thickness δL, pressure P, and unburned gas temperature Tu on turbulent flame speed ST, a set of conditions is designed by adjusting nitrogen mole fraction in lean H2/O2/N2 (Le=0.35) and stoichiometric CH4/O2/N2 (Le≈1) mixtures. The adjustment is performed by simulating complex-chemistry laminar flames to obtain the same laminar flame speeds SL not only for different fuels, but also for different pressures (1, 2, and 5 atm). The mixtures are characterized by significantly different SL at Tu=300 K and 400 K, whereas variations in δL with the temperature are sufficiently weak. Moreover, laminar flame thicknesses are approximately equal for H2-based and CH4-based mixtures at the same P, but are significantly decreased with increasing pressure. For this set of conditions, ST is measured by applying schlieren imaging techniques to film expansion of centrally ignited, statistically spherical flames in homogeneous isotropic turbulence generated by a dual-chamber, constant-pressure, fan-stirred explosion facility. Analyses of the measured data show the following trends. First, turbulent flame speed is increased by both P and Tu, whereas ST/SL is decreased with increasing Tu. Second, turbulent flame speed measured at different P and Tu can be predicted by allowing for SL(P,Tu) and δL(P,Tu). Thus, the present data do not call for explicitly substituting normalized pressure or temperature into a turbulent flame speed approximation. Third, ST is increased with decreasing laminar flame thickness. Fourth, speeds of the lean H2/O2/N2 flames are higher when compared to the stoichiometric CH4/O2/N2 flames, with this difference is increased (reduced) by P (Tu, respectively). Fifth, all measured data on ST can quantitatively be described by substituting SL and δL with the counterpart characteristics of highly strained twin laminar flames. The latter finding supports leading point concept of premixed turbulent combustion.

Propagation characteristics of lean turbulent premixed ammonia–hydrogen flames

Direct numerical simulations (DNS) of fuel-lean turbulent premixed NH3-H2-N2-air flames are analyzed to investigate propagation and flame structural characteristics under fixed velocity and length ratios. To comprehensively assess the impact of diffusive-thermal imbalances on hydrogen–enriched ammonia flames, additional solutions with unity-Lewis-number transport were analyzed and compared with those obtained using the mixture-averaged transport model. The increase of H2 fraction in the fuel leads to elevated mean turbulent flame speed and stretch factor, indicating the impact of thermal-diffusive instability. The turbulent flame speed of the 60%NH3-25%H2-15%N2-air flame displays pronounced oscillations, a phenomenon absent in other mixtures considered in the current study. This behavior is attributed to the preferential diffusion of H2 mixed with the low-reactive NH3 in moderate quantities, resulting in higher generation of flame elements extending into the product side and dynamic evolution of H2. The flame structure analysis, in terms of conditional averages, revealed a distinctive variation in H2 and H atom distributions. The flames with a higher H2 fraction (40%NH3-45%H2-15%N2-air) produced a second peak of H2O2 in the trailing edge region, indicating additional production in the intense reaction zone. Additionally, in the 60%NH3-25%H2-15%N2-air flame, the reaction rate of H2 exhibited a unique behavior, with H2 being produced in the intermediate flame zone and rapidly consumed in the reaction zone, differing from other cases.