The effect of fuel inlet turbulence intensity on H2/CH4 flame structure of MILD combustion using the LES method

Large eddy simulations (LES) are employed to investigate the effect of the inlet turbulence intensity on the H₂/CH₄ flame structure in a hot and diluted co-flow stream which emulates the (Moderate or Intense Low-oxygen Dilution) MILD combustion regime. In this regard, three fuel inlet turbulence intensity profiles with the values of 4%, 7% and 10% are superimposed on the annular mixing layer. The effects of these changes on the flame structure under the MILD condition are studied for two oxygen concentrations of 3% and 9% (by mass) in the oxidiser stream and three hot co-flow temperatures 1300, 1500 and 1750 K. The turbulence-chemistry interaction of the numerically unresolved scales is modelled using the (Partially Stirred Reactor) PaSR method, where the full mechanism of GRI-2.11 represents the chemical reactions. The influences of the turbulence intensity on the flame structure under the MILD condition are studied by using the profile of temperature, CO and OH mass fractions in both physical and mixture fraction spaces at two downstream locations. Also, the effects of this parameter are investigated by contours of OH, HCO and CH₂O radicals in an area near the nozzle exit zone. Results show that increasing the fuel inlet turbulence intensity has a profound effect on the flame structure particularly at low oxygen mass fraction. This increment weakens the combustion zone and results in a decrease in the peak values of the flame temperature and OH and CO mass fractions. Furthermore, increasing the inlet turbulence intensity decreases the flame thickness, and increases the MILD flame instability and diffusion of un-burnt fuel through the flame front. These effects are reduced by increasing the hot co-flow temperature which reinforces the reaction zone.     

Generalisation of the eddy-dissipation concept for jet flames with low turbulence and low Damköhler number

Moderate or intense low oxygen dilution (MILD) combustion has been the focus of a range of fundamental experimental and numerical studies. Reasonable agreement between experimental and numerical investigations, however, requires finite-rate chemistry models and, often, ad hoc model adjustment. To remedy this, an adaptive eddy dissipation concept (EDC) combustion model has previously been developed to target conditions encountered in MILD combustion; however, this model relies on a simplified, pre-defined assumption about the combustion chemistry. The present paper reports a generalised version of the modified EDC model without the need for an assumed, single-step chemical reaction or ad hoc coefficient tuning. The results show good agreement with experimental measurements of two CH₄/H₂ flames in hot coflows, showing improvements over the standard EDC model as well as the previously published modified EDC model. The updated version of the EDC model also demonstrates the capacity to reproduce the downstream transition in flame structure of a MILD jet flame seen experimentally, but which has previously proven challenging to capture computationally. Analyses of the previously identified dominant heat-release reactions provide insight into the structural differences between a conventional autoignitive flame and a flame in the MILD combustion regime, whilst highlighting the requirement for a generalised EDC combustion model.     

Additional criteria for MILD coal combustion

We proposed a theoretical basis for Moderate or Intense Low-oxygen Dilution (MILD) coal combustion based on the turbulent scalar energy spectra. This is motivated by the hypothesis that smallest scalar mixing length scales should be on the order of the particle size or smaller to ensure that mixing can occur to prevent formation of diffusion flames. Our proposed criterion is evaluated using several experimental datasets from the literature for coal combustion in both MILD and traditional combustion regimes. The experimental results confirm that the smallest mixing length scales should be of the order of or smaller than the particle diameter, η_mix?d_p, to breakup the heat and mass transfer boundary layers around particles in MILD coal combustion. Results indicate that poor mixing of species with small Schmidt numbers around small particles leads to the high luminous intensity in the reactor. The effects of inlet velocity and jet diameter on the mixing length scales are analyzed. Higher inlet velocity and smaller jet diameter are expected to reach MILD regime. The proposed criterion can be used to guide experimental design to achieve MILD conditions for coal combustion.     

Transported PDF simulation of turbulent CH4/H2 flames under MILD conditions with particle-level sensitivity analysis

Transported probability density function (TPDF) simulation with sensitivity analysis has been conducted for turbulent non-premixed CH₄/H₂ flames of the jet-into-hot-coflow (JHC) burner, which is a typical model to emulate moderate or intense low oxygen dilution combustion (MILD). Specifically, two cases with different levels of oxygen in the coflow stream, namely HM1 and HM3, are simulated to reveal the differences between MILD and hot-temperature combustion. The TPDF simulation well predicts the temperature and species distributions including those of OH, CO and NO for both cases with a 25-species mechanism. The reduced reaction activity in HM1 as reflected in the peak OH concentration is well correlated to the reduced oxygen in the coflow stream. The particle-level local sensitivities with respect to mixing and chemical reaction further show dramatic differences in the flame characteristics. HM1 is less sensitive to mixing and reaction parameters than HM3 due to the suppressed combustion process. Specifically, for HM1 the sensitivities to mixing and chemical reactions have comparable magnitude, indicating that the combustion progress is controlled by both mixing and reaction in MILD combustion. For HM3, there is however a change in the combustion mode: during the flame initialization, the combustion progress is more sensitive to chemical reactions, indicating that finite-rate chemistry is the controlling process during the autoignition process for flame stabilization; at further downstream where the flame has established, the combustion progress is controlled by mixing, which is characteristic of nonpremixed flames. An examination of the particles with the largest sensitivities reveals the difference in the controlling mixtures for flame stabilization, namely, the stoichiometric mixtures are important for HM1, whereas, fuel-lean mixtures are controlling for HM3. The study demonstrates the potential of TPDF simulations with sensitivity analysis to investigate the effects of finite-rate chemistry on the flame characteristics and emissions, and reveal the controlling physio-chemical processes in MILD combustion.     

Numerical study of influence of molecular diffusion in the Mild combustion regime

In this paper, the importance of molecular diffusion versus turbulent transport in the moderate or intense low-oxygen dilution (Mild) combustion mode has been numerically studied. The experimental conditions of Dally et al. [Proc. Combust. Inst. 29 (2002) 1147–1154] were used for modelling. The EDC model was used to describe the turbulence–chemistry interaction. The DRM-22 reduced mechanism and the GRI 2.11 full mechanism were used to represent the chemical reactions of an H₂/methane jet flame. The importance of molecular diffusion for various O₂ levels, jet Reynolds numbers and H₂ fuel contents was investigated. Results show that the molecular diffusion in Mild combustion cannot be ignored in comparison with the turbulent transport. Also, the method of inclusion of molecular diffusion in combustion modelling has a considerable effect on the accuracy of numerical modelling of Mild combustion. By decreasing the jet Reynolds number, decreasing the oxygen concentration in the airflow or increasing H₂ in the fuel mixture, the influence of molecular diffusion on Mild combustion increases.     

Simulations of laminar non-premixed flames of kerosene with hot combustion products as oxidiser

The effects of hot combustion product dilution in a pressurised kerosene-burning system at gas turbine conditions were investigated with laminar counterflow flame simulations. Hot combustion products from a lean (φ = 0.6) premixed flame were used as an oxidiser with kerosene surrogate as fuel in a non-premixed counterflow flame at 5, 7, 9 and 11 bar. Kerosene-hot product flames, referred to as ‘MILD’, exhibit a flame structure similar to that of kerosene–air flames, referred to as ‘conventional’, at low strain rates. The Heat Release Rate (HRR) of both conventional and MILD flames reflects the pyrolysis of the primary and intermediate fuels on the rich side of the reaction zone. Positive HRR and OH regions in mixture fraction space are of similar width to conventional kerosene flames, suggesting that MILD flames are thin fronts. MILD flames do not exhibit typical extinction behaviour, but gradually transition to a mixing solution at very high rates of strain (above A = 160, 000 s^?1 for all pressures). This is in agreement with literature that suggests heavily preheated and diluted flames have a monotonic S-shaped curve. Despite these differences in comparison with kerosene–air flames, MILD flames follow typical trends as a function of both strain and pressure. Further still, the peak locations of the overlap of OH and CH₂O mass fractions in comparison with the peak HRR indicate that the pixel-by-pixel product of OH- and CH₂O-PLIF signals is a valid experimental marker for non-premixed kerosene MILD and conventional flames.     

Large-eddy simulation of MILD combustion using partially stirred reactor approach

Subgrid-scale (SGS) parameterization and method for calculating filtered reaction rate are critical components of an accurate large-eddy simulation (LES) of turbulent flames. In this study, we integrate gradient-type structural SGS models with a partially stirred reactor approach by using detailed chemical kinetics to simulate a turbulent methane/hydrogen jet flame under moderate or intense low-oxygen dilution (MILD) conditions. The study examines two oxygen dilution levels. The framework is assessed through a systematic and comprehensive comparison of temperature, and mass fractions of major and minor species with experimental data and other reference simulation results. Overall, the statistics of the combustion field show excellent agreement with measurements at different axial locations, and a significant improvement compared to some previous simulations. It suggests that the proposed nonlinear LES framework is able to accurately model MILD combustion with reasonable computational cost.     

In situ Raman characterization of nanoparticle aerosols during flame synthesis

Raman spectroscopy is applied to diagnose nanoparticle presence and characteristics in a gaseous flow field. Specifically, in situ monitoring of the Raman-active modes of TiO₂ and Al₂O₃ nanoparticles in aerosol form is demonstrated in high-temperature flame environments. This technique serves as a sensitive and reliable way to characterize particle composition and crystallinity (e.g. anatase versus rutile) and delineate the phase conversion of nanoparticles as they evolve in the flow field. The effect of temperature on the solid-particle Raman spectra is investigated by seeding nanoparticles into a co-flow jet diffusion flame, where local gas-phase temperatures are correlated by shape-fitting the N₂ vibrational Stokes Q-branch Raman spectra. Applying the technique to a flame synthesis environment, the results demonstrate that in situ Raman of as-formed nanoparticles can be readily applied to other gas-phase synthesis systems, especially as an on-line diagnostic.     

Validation of flow simulation and gas combustion sub-models for the CFD-based prediction of NOx formation in biomass grate furnaces

While reasonably accurate in simulating gas phase combustion in biomass grate furnaces, CFD tools based on simple turbulence–chemistry interaction models and global reaction mechanisms have been shown to lack in reliability regarding the prediction of NO_x formation. Coupling detailed NO_x reaction kinetics with advanced turbulence–chemistry interaction models is a promising alternative, yet computationally inefficient for engineering purposes. In the present work, a model is proposed to overcome these difficulties. The model is based on the Realizable k–? model for turbulence, Eddy Dissipation Concept for turbulence–chemistry interaction and the HK97 reaction mechanism. The assessment of the sub-models in terms of accuracy and computational effort was carried out on three laboratory-scale turbulent jet flames in comparison with the experimental data. Without taking NO_x formation into account, the accuracy of turbulence modelling and turbulence–chemistry interaction modelling was systematically examined on Sandia Flame D and Sandia CO/H₂/N₂ Flame B to support the choice of the associated models. As revealed by the Large Eddy Simulations of the former flame, the shortcomings of turbulence modelling by the Reynolds averaged Navier–Stokes (RANS) approach considerably influence the prediction of the mixing-dominated combustion process. This reduced the sensitivity of the RANS results to the variations of turbulence–chemistry interaction models and combustion kinetics. Issues related to the NO_x formation with a focus on fuel bound nitrogen sources were investigated on a NH₃-doped syngas flame. The experimentally observed trend in NO_x yield from NH₃ was correctly reproduced by HK97, whereas the replacement of its combustion subset by that of a detailed reaction scheme led to a more accurate agreement, but at increased computational costs. Moreover, based on results of simulations with HK97, the main features of the local course of the NO_x formation processes were identified by a detailed analysis of the interactions between the nitrogen chemistry and the underlying flow field.     

介质阻挡体放电对甲烷扩散火焰特性的影响

本文在高频交流激励模式下,采用同轴圆柱构型激励器,开展了介质阻挡体放电对空气/甲烷同轴剪切扩散火焰燃烧特性影响实验研究。激励器敷设在外喷嘴环缝以电离空气,采用纹影系统和B型热电偶分别获取流场形态和火焰温度,激励频率为8 kHz,通过改变气体流量和放电电压,分析了不同工况下射流流场、火焰结构和火焰温度在等离子体作用下的变化规律。结果表明:等离子体气动效应能有效增强射流湍流强度,强化空气/甲烷掺混,增大射流角,并随激励电压提高作用效果逐渐增强,实验中未形成明显扩张流动的初始射流在放电电压30 kV时其射流角最大为23.5°。贫燃条件下等离子体激励会改善火焰形态,增强燃烧稳定性,并在流量较低时缩短火焰长度。此外,富燃火焰下游温度会随着激励强度增大不断升高,而贫燃火焰下游温度变化受上游燃烧强度影响存在升高和降低两种情况。     

A transported probability density function/photon Monte Carlo method for high-temperature oxy–natural gas combustion with spectral gas and wall radiation

A computational fluid dynamics model for high-temperature oxy–natural gas combustion is developed and exercised. The model features detailed gas-phase chemistry and radiation treatments (a photon Monte Carlo method with line-by-line spectral resolution for gas and wall radiation – PMC/LBL) and a transported probability density function (PDF) method to account for turbulent fluctuations in composition and temperature. The model is first validated for a 0.8 MW oxy–natural gas furnace, and the level of agreement between model and experiment is found to be at least as good as any that has been published earlier. Next, simulations are performed with systematic model variations to provide insight into the roles of individual physical processes and their interplay in high-temperature oxy–fuel combustion. This includes variations in the chemical mechanism and the radiation model, and comparisons of results obtained with versus without the PDF method to isolate and quantify the effects of turbulence–chemistry interactions and turbulence–radiation interactions. In this combustion environment, it is found to be important to account for the interconversion of CO and CO₂, and radiation plays a dominant role. The PMC/LBL model allows the effects of molecular gas radiation and wall radiation to be clearly separated and quantified. Radiation and chemistry are tightly coupled through the temperature, and correct temperature prediction is required for correct prediction of the CO/CO₂ ratio. Turbulence–chemistry interactions influence the computed flame structure and mean CO levels. Strong local effects of turbulence–radiation interactions are found in the flame, but the net influence of TRI on computed mean temperature and species profiles is small. The ultimate goal of this research is to simulate high-temperature oxy–coal combustion, where accurate treatments of chemistry, radiation and turbulence–chemistry–particle–radiation interactions will be even more important.     

A comparative study of H2-air premixed flame in micro combustors with different physical and boundary conditions

A numerical study of H₂-air premixed combustion in the micro channels with a detailed chemical reaction mechanism is performed by solving the two-dimensional fully elliptic governing equations of continuity, momentum, energy and species, coupled with the energy equation in the solid wall. A reference case is defined as the combustion in a cylindrical tube with 0.8 mm inner diameter and 8 mm length with a non-slip wall and a uniform velocity profile at the inlet plane. Different physical and boundary conditions have been applied in order to investigate their respective effects on the flame temperature. The conditions studied in the current paper include the combustor size and geometry, inlet velocity profile, axial heat conduction in the solid wall and slip-wall and temperature jump at the gas–solid interface. It is noted that effects of Knudsen number (slip-wall and temperature jump) on the thermal and fluid field are not very significant in a d = 0.4 mm micro combustor. Furthermore, the qualitative effects of Knudsen number on the flame temperature are analysed. The results of this paper indicate that these various boundary and physical conditions have effects on the flame temperature to different extent and should be carefully monitored when applied for different applications.     

Stability and emissions control using air injection and H2 addition in premixed combustion

We experimentally study lean premixed combustion stabilized behind a backward-facing step. For a propane–air mixture, the lean blowout limit is associated with strong pressure fluctuation arising simultaneously with strong flame–vortex interactions, which have been shown to constitute the mechanism of heat release dynamics in this flow. A high-speed air jet, issuing from a small slot and injected perpendicular to the main flow near the step, is used to disrupt this mechanism. For momentum ratio of jet to main flow below unity, the jet dilutes the mixture, further destabilizing the flame or leading to complete blowout. Above unity, the flame becomes more stable, and the pressure oscillations are suppressed. Flow visualization and OH*/CH* chemiluminescence measurements show that a strong jet produces a more compact flame that is less driven by the wake vortex, anchored closer to the step, and deflected upwards away from the lower wall of the channel. This renders the flame less vulnerable to heat loss and strong strains, which improves its stability and extends the flammability limit. Adding hydrogen to the main fuel improves the flame stability over the entire range of the air jet mass flow, with better results for momentum ratio larger than 1; H₂ pulls the flame further upstream, away from the shear zone and the unsteady vortex. NO_x emission benefits from the air jet, while, with H₂ addition, NO_x concentration is higher in the products as the overall burning temperature rises. However, hydrogen addition enables extending the flammability limit further by increasing air supply in the primary stream, hence achieving lower NO_x. The study suggests a simpler, almost passive, multi-objective combustion control technique and indicates that hydrogen addition can be a successful in situ approach for NO_x reduction.     

Stochastic low-order modelling of hydrogen autoignition in a turbulent non-premixed flow

Autoignition risk in initially non-premixed flowing systems, such as premixing ducts, must be assessed to help the development of low-NO_x systems and hydrogen combustors. Such situations may involve randomly fluctuating inlet conditions that are challenging to model in conventional mixture-fraction-based approaches. A Computational Fluid Dynamics (CFD)-based surrogate modelling strategy is presented here for fast and accurate predictions of the stochastic autoignition behaviour of a hydrogen flow in a hot air turbulent co-flow. The variability of three input parameters, i.e., inlet fuel and air temperatures and average wall temperature, is first sampled via a space-filling design. For each sampled set of conditions, the CFD modelling of the flame is performed via the Incompletely Stirred Reactor Network (ISRN) approach, which solves the reacting flow governing equations in post-processing on top of a Large Eddy Simulation (LES) of the inert hydrogen plume. An accurate surrogate model, namely a Gaussian Process, is then trained on the ISRN simulations of the burner, and the final quantification of the variability of autoignition locations is achieved by querying the surrogate model via Monte Carlo sampling of the random input quantities. The results are in agreement with the observed statistics of the autoignition locations. The methodology adopted in this work can be used effectively to quantify the impact of fluctuations and assist the design of practical combustion systems.     

Large eddy simulation of a pulverized coal jet flame ignited by a preheated gas flow

Large eddy simulation (LES) is applied to a pulverized coal jet flame ignited by a preheated gas flow. The simulation results are compared to experimental data obtained for the inlet stoichiometric ratios of 0.14, 0.22, and 0.36. An accurate and computationally inexpensive devolatilization model suitable for combustion simulation in LES is proposed and incorporated into the LES. The numerical results of gas temperature and coal burnout on the centerline show good agreement with the experimental data. Two kinds of lift-off heights are introduced to verify the combustion simulation. One is the height from the primary nozzle exit to the starting point of the growing flame region. The other is the height from the primary nozzle exit to the starting point of the continuous flame region. The calculated results of the two lift-off heights show good agreement with the experimental data. In contrast to LES, the standard k–ε model overestimates the lift-off heights because it calculates time-averaged temperature which does not contain information about local flame structure. The stoichiometric ratio in the gas phase at the starting point of the growing flame region is found to be independent of the inlet stoichiometric ratio in the range from 0.14 to 0.36.     

Role of radicals on MILD combustion inception

The criterion used to define MILD combustion in non-premixed condition is analysed using Direct Numerical Simulation (DNS) of MILD combustion of methane-diluted air established with internal exhaust gas recirculation. The simulations reveal multiple interacting reaction zones in MILD combustion which are extremely different from conventional combustion. Furthermore, DNS deduced S-curves highlight the role of chemically active species. Specifically, the temperature rise is accompanied with an increase in the scalar dissipation rate of mixture fraction, which is quite contrasting to the classical S-curve from the classical flame theories. This observation is explained on a physical basis.     

Ethanol turbulent spray flame response to gas velocity modulation

A numerical investigation of the interaction between a spray flame and an acoustic forcing of the velocity field is presented in this paper. In combustion systems, a thermoacoustic instability is the result of a process of coupling between oscillations in heat released and acoustic waves. When liquid fuels are used, the atomisation and the evaporation process also undergo the effects of such instabilities, and the computational fluid dynamics of these complex phenomena becomes a challenging task. In this paper, an acoustic perturbation is applied to the mass flow of the gas phase at the inlet and its effect on the evaporating fuel spray and on the flame front is investigated with unsteady Reynolds averaged Navier-Stokes numerical simulations. Two flames are simulated: a partially premixed ethanol/air spray flame and a premixed pre-vaporised ethanol/air flame, with and without acoustic forcing. The frequencies used to perturb the flames are 200 and 2500 Hz, which are representative for two different regimes. Those regimes are classified based on the Strouhal number St = (D/U)f_f: at 200 Hz, St = 0.07, and at 2500 Hz, St = 0.8. The exposure of the flame to a 200 Hz signal results in a stretching of the flame which causes gas field fluctuations, a delay of the evaporation and an increase of the reaction rate. The coupling between the flame and the flow excitation is such that the flame breaks up periodically. At 2500 Hz, the evaporation rate increases but the response of the gas field is weak and the flame is more stable. The presence of droplets does not play a crucial role at 2500 Hz, as shown by a comparison of the discrete flame function in the case of spray and pre-vaporised flame. At low Strouhal number, the forced response of the pre-vaporised flame is much higher compared to that of the spray flame.     

高焓横向射流尾迹区掺混燃烧的数值特性

尾迹区作为横向射流流场的重要结构受到广泛关注,其掺混和燃烧特性对近壁面区域的流场特性有重要影响.文章在对仿真充分验证的基础上,采用Reynolds平均模拟方法对Ma=8飞行条件下高焓横向射流尾迹区中的掺混和燃烧特性进行了数值研究.探究了冷热流场尾迹区中的氢气掺混特性,冷态流场尾迹区中的激波结构对氢气分布产生一定影响,热态流场尾迹区中存在多种氢气掺混路径.V形回流区中的高浓度氢气对燃烧产生了一定的阻碍作用.定量测量了尾迹区中的火焰结构,尾迹区火焰的顶点位置随高度增加向下游线性移动,受射流主流影响,尾迹区火焰的展向宽度在距离壁面一定高度后开始增大.对冷热流场中的主要参数进行了对比,燃烧消耗了氢气使温度升高,但是尾迹区中的流动速度没有明显增加,燃烧放出的热量没有完全转化为流体的动能.      

The response of turbulent premixed flames to normal acoustic excitation

To model the thermo-acoustic excitation of flames in practical combustion systems, it is necessary to know how a turbulent flame front responds to an incident acoustic wave. This will depend partly on the way in which the burning velocity responds to the wave. In this investigation, the response of CH₄/air and CH₄/H₂/air mixtures has been observed in a novel flame stabilisation configuration, in which the premixture of fuel and air is made to decelerate under controlled conditions in a wide-angle diffuser. Control is provided by an annular wall-jet of air and by turbulence generators at the inlet. Ignition from the outlet of the diffuser allows an approximately flat flame to propagate downwards and stabilise at a height that depends on the turbulent burning velocity. When the flow is excited acoustically, the ensemble-averaged height oscillates. The fluctuations in flow velocity and flame height are monitored by phase-locked particle image velocimetry and OH-planar laser induced fluorescence, respectively. The flame stabilised against a lower incident velocity as the acoustic amplitude increased. In addition, at the lowest frequency of 52 Hz, the fluctuations in turbulent burning velocity (as represented by the displacement speed) were out-of-phase with the acoustic velocity. Thus, the rate of displacement of the flame front relative to the flow slowed as the flow accelerated, and so the flame movement was bigger than it would have been if the burning velocity had not responded to the acoustic fluctuation. With an increase in frequency to 119 Hz, the relative flame movement became even larger, although the phase-difference was reduced, so the effect on burning velocity was less dramatic. The addition of hydrogen to the methane, so as to maintain the laminar burning velocity at a lower equivalence ratio, suppressed the response at low amplitude, but at a higher amplitude, the effect was reversed.