直喷柴油机火焰浮起长度的数值模拟研究

采用三维CFD模型模拟了直喷柴油机缸内喷雾燃烧过程,模拟缸压曲线得到了实验的验证.通过高温区与喷嘴之间的稳定距离来确定柴油机火焰浮起长度,研究在不同进气条件下火焰浮起长度的变化情况.该模型成功地预测了火焰浮起长度随着初始进气压力的增大而减少,随着进气温度的升高出现先增大后减少的趋势.同时模拟了在不同EGR率下柴油机缸内燃烧情况,发现火焰浮起长度和燃料着火延迟时间都随EGR率的增加而增大.     

Large Eddy Simulation of flame stabilisation dynamics and vortex control in a lifted H2/N2 jet flame

Flame stabilisation in (highly) preheated mixture is common in several industrial applications. When the reactants are injected separately in the device (usually at high-speed), the flame is lifted so that the fuel and oxidant first mix to give an ignitable mixture. If the temperature of the mixture is adequate, it auto-ignites stabilizing the flame. Here we focus on an academic lifted jet flame and Large Eddy Simulation (LES) is used to capture the flame and auto-ignition dynamics. Comparisons with experimental data show that LES simulates accurately high OH fluctuation levels at the stabilisation location. The vortex dynamics linked to these fluctuations is analyzed and it is found that small scale coherent structures play a vital role in the auto-ignition process. These structures are axial vorticity tubes (braids) and are located relatively far (in the radial direction) from the shear-layer. As a consequence, the lift-off height varies dramatically in time leading to OH fluctuations of the order of the mean OH concentration. This scenario is monitored in the compositional space highlighting the simultaneous evolution of OH, HO ₂ and temperature. Further, different strategies for open-loop control of the flame lift-off height are tested. In order to anchor the flame at different positions downstream of the nozzle, the vortex dynamics in the shear-layer was modified. Promoting successively vortex ring and braids, the auto-ignition region was moved significantly. In particular, modified nozzle geometries impacted the formation of braids and ensured a good premixing very close to the nozzle. As a consequence, it was possible to reduce significantly the lift-off height and stabilise the flame few diameters downstream of the nozzle.     

Large Eddy simulation of a pulverised coal jet flame

Large Eddy simulation (LES) has been applied to the pulverised coal jet flame studied at the Japanese Central Research Institute of Electric Power (CRIEPI). A working set of models to represent coal combustion, Lagrangian particle transport and radiative heat transfer in an LES framework has been implemented and tested. The simulation results of the flow field were compared to experimental data for both a reactive and non-reactive case, and an overall good agreement emerged. A simple method for replicating pyrometer measurements was developed for the LES and results obtained from the method were compared to the experimental data. Finally the species concentrations were compared to experimental results for CO₂, O₂ and N₂. The results show the potentials of using LES for pulverised coal combustion and open the way for further developments on the coal combustion models and the applications to more complex burners.     

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

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

Stabilized,flat iron flames on a hot counterflow burner

Metal powder combustion has traditionally been studied to mitigate the risk of industrial accidents and to determine the contributions of metals as additives to the performance of energetic materials. Recently, there has been growing interest in exploring the potential of metal powders as recyclable, zero-carbon energy carriers as an alternative to the hydrocarbons known to contribute to climate change. The present work introduces, for the first time, a stabilized flat iron flame. The counterflow burner used in this work is comprised of an inverted ceramic nozzle which sits above, and is aligned axially with, a lower nozzle producing a laminar flow of particles suspended in an oxidizing gas. A stabilized methane flame sits inside the top nozzle and the hot combustion products impinge upon the two-phase flow from the bottom nozzle, creating a stagnation plane. Spherical iron powder, with 90% of the particles less than 2.5 µm in size, is pre-loaded into a piston and dispersed using mixtures of 30% and 40% oxygen balanced in argon. Flame speeds are measured using particle image velocimetry (PIV), while flame temperatures are determined using multicolour pyrometry. It is found that flame speeds range between 30 cm/s and 45 cm/s for both oxidizing mixtures. Despite having fuel loadings below stoichiometric concentrations, the observed particle combustion temperatures are close to the adiabatic flame temperature of the stoichiometric mixture, indicating combustion in the diffusion-controlled regime for these small particles. Finally, the independence of the flame speeds with respect to oxygen concentration suggests flame propagation in the discrete regime.     

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

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

Large-eddy simulation of evaporating spray in a coaxial combustor

Large-eddy simulation of an evaporating isopropyl alcohol spray in a coaxial combustor is performed. The Favre-averaged, variable density, low-Mach number Navier-Stokes equations are solved on unstructured grids with dynamic subgrid scale model to compute the turbulent gas-phase. The original incompressible flow algorithm for LES on unstructured grids by [Mahesh et al., J. Comp. Phys. 197 (2004) 215–240] is extended to include density variations and droplet evaporation. An efficient particle-tracking scheme on unstructured meshes is developed to compute the dispersed phase. Experimentally measured droplet size distribution and size-velocity correlation near the nozzle exit are used as the inlet conditions for the spray. The predictive capability of the LES approach on unstructured grids together with Lagrangian droplet dynamics models to capture the droplet dispersion characteristics, size distributions, and the spray evolution is examined in detail. The mean and turbulent quantities for the gas and particle phases are compared to experimental data to show good agreement. It is shown that for low evaporation rates considered in the present study, a well resolved large-eddy simulation together with simple subgrid models for droplet evaporation and motion provides good agreement of the mean and turbulent quantities for the gas and droplet phases compared to the experimental data. This work represents an important first step to assess the predictive capability of the unstructured grid LES approach applied to spray vaporization. The novelty of the results presented is that they establish a baseline fidelity in the ability to simulate complex flows on unstructured grids at conditions representative of gas-turbine combustors.     

Diesel flame lift-off stabilization in the presence of laser-ignition: a numerical study

Diesel flame lift-off and stabilization in the presence of laser-ignition were numerically investigated with the method of Eulerian stochastic fields. The aim was to scrutinise the interaction between the lifted diesel flame and an ignition kernel upstream of the lifted flame. The numerical simulation was carried out in a constant-volume combustion vessel with n-heptane as fuel. The process was studied previously in an experiment employing Diesel #2 as the fuel in the same combustion vessel. In the experiment a lifted flame was first established at a position downstream of the nozzle. An ignition kernel was then initiated using a high-energy pulse laser at a position upstream of the natural lift-off position of the diesel flame. The laser-ignition kernel was modelled using a high-temperature (～2000 K) hot spot. In both experiment and simulations the upstream front of the ignition kernel was shown to remain around the initial laser ignition site for a substantially long period of time, while the downstream front of the ignition kernel propagates rapidly towards the natural lift-off position downstream of the laser ignition site. The lift-off position oscillated before the final stabilization at the natural lift-off position. The structures and the propagation speed of the reaction fronts in the laser-ignition kernel and the main flame were analysed. Two different stabilization mechanisms, the auto-ignition mechanism and the flame propagation mechanism, were identified for the naturally lifted flame and the laser-induced reaction front, respectively. A mechanism was proposed to explain the oscillation of the lift-off position.     

LES studies of the flow in a swirl gas combustor

Environmental and other practical concerns have led to the development of compact gas turbine combustors burning lean mixtures leading to potentially low CO and NO_x emissions. The compact design requires efficient atomization and mixing together with a compact premixed flame. Associated with these requirements are higher temperatures, increased heat transfer, and thermal load, thus increasing the danger of combustion instabilities (causing performance deterioration and excessive mechanical loads), and possible off-design operation. Numerical simulations of reacting flows are well suited to address these issues. To this end, large eddy simulation (LES) is particularly promising. The philosophy behind LES is to explicitly simulate the large scales of the flow and the thermochemistry, affected by boundary conditions whilst modeling only the small scales, including the interaction between the flow and the combustion processes. Here, we examine the flow and the flame in a model gas turbine combustor (General Electric’s lean premixed dry low NO_x LM6000) to evaluate the potential of LES for design studies of engineering applications and to study the effects of the combustor confinement geometry on the flow and on the flame dynamics. Two LES models, a Monotone Integrated LES model with 1 and 2 step Ahrrenius chemistry, and a fractal flame-wrinkling LES model coupled to a conventional one-equation eddy-viscosity subgrid model, are used. Reasonable agreement is found when comparing predictions with experimental data and with other LES computations of the same case. Furthermore, the combustor confinement geometry is found to strongly affect the vortical flow, and hence also the flame and its dynamics.     

Flame characterisation of gas-assisted pulverised coal combustion using FPV-LES

A multiphase flamelet/progress variable (FPV) model for the large eddy simulation (LES) of gas-assisted pulverised coal combustion (PCC) is developed. The target of the simulation is the Darmstadt turbulent gas-assisted swirling solid fuel combustion chamber. The coal particles are treated as Lagrangian point particles, the position, momentum and energy of which are tracked. The gas phase is described by the low-Mach Navier-Stokes equations alongside the Eulerian transport equations of the governing variables for the FPV model. The set of chemical states of the PCC flame is pre-tabulated in a six-dimensional flamelet table and determined by the mixing of the primary fuel stream, volatiles and char off-gases with the oxidising air, the progress of chemical reactions, the interphase heat transfer, as well as sub-grid scale variations. A presumed $\beta$-PDF approach for the total mixture fraction is applied to capture sub-grid scale effects. The discrete ordinate method (DOM) with the weighted sum of grey gases model (WSGGM) is employed to model radiation. The FPV-LES results are validated against the experimental evidence and a good agreement of the predicted mean and RMS velocities, as well as the mean gas temperature between experiments and simulations is obtained. The contributions of the pilot, volatile and char off-gas fuel streams to the coal flame are analysed. It is found that most regions of the furnace are dominated by either pilot or volatile combustion, while char conversion only occurs in the far downstream and outer furnace regions. The pilot gas dominates the near-wall region inside the quarl, whereas the volatile gas mainly released from small particles dominates a first volatile combustion zone in the interior of the internal recirculation zone. Larger particles heat up more slowly and release their volatile content further downstream, leading to a secondary volatile combustion zone.     

Modeling of flame lift-off length in diesel low-temperature combustion with multi-dimensional CFD based on the flame surface density and extinction concept

Low-Temperature Combustion (LTC) is becoming a promising technology for simultaneously reducing soot and NOx emissions from diesel engines. LTC regimes are evaluated by the flame lift-off length – the distance from the injector orifice to the location of hydroxyl luminescence closest to the injector in the flame jet. Various works have been dedicated to successful simulations of lifted flames of a diesel jet by use of various combustion modeling approaches. In this work, flame surface density and flamelet concepts were used to model the diesel lift-off length under LTC conditions. Numerical studies have been performed with the ECFM3Z model, n-Heptane and diesel fuels to determine the flame lift-off length and its correlation with soot formation under quiescent conditions. The numerical results showed good agreement with experimental data, which were obtained from an optically accessible constant volume chamber and presented at the Engine Combustion Network (ECN) of Sandia National Laboratories. It was shown that at a certain distance downstream from the injector orifice, stoichiometric scalar dissipation rate matched the extinction scalar dissipation rate. This computed extinction scalar dissipation rate correlated well with the flame lift-off length. For the range of conditions investigated, adequate quantitative agreement was obtained with the experimental measurements of lift-off length under various ambient gas O₂ concentrations, ambient gas temperatures, ambient gas densities and fuel injection pressures. The results showed that the computed lift-off length values for most of the conditions lay in a reasonable range within the quasi-steady lift-off length values obtained from experiments. However, at ambient temperatures lower than 1000 K, the lift-off length values were under-predicted by the numerical analysis. This may be due to the use of the droplet evaporation model as it is believed that evaporation has a strong effect on the lift-off length.     

LES/CMC modelling of ignition and flame propagation in a non-premixed methane jet

The Large Eddy Simulation (LES) / Conditional Moment Closure (CMC) model with detailed chemistry is used for modelling spark ignition and flame propagation in a turbulent methane jet in ambient air. Two centerline and one off-axis ignition locations are simulated. We focus on predicting the flame kernel formation, flame edge propagation and stabilization. The current LES/CMC computations capture the three stages reasonably well compared to available experimental data. Regarding the formation of flame kernel, it is found that the convection dominates the propagation of its downstream edge. The simulated initial downstream and radial flame propagation compare well with OH-PLIF images from the experiment. Additionally, when the spark is deposited at off-centerline locations, the flame first propagates downstream and then back upstream from the other side of the stoichiometric iso-surface. At the leading edge location, the chemical source term is larger than others in magnitude, indicating its role in the flame propagation. The time evolution of flame edge position and the final lift-off height are compared with measurements and generally good agreement is observed. The conditional quantities at the stabilization point reflect a balance between chemistry and micro-mixing. This investigation, which focused on model validation for various stages of spark ignition of a turbulent lifted jet flame through comparison with measurements, demonstrates that turbulent edge flame propagation in non-premixed systems can be reasonably well captured by LES/CMC.     

RANS predictions of turbulent diffusion flames: comparison of a reactor and a flamelet combustion model to the well stirred approach

The flame stabilisation process in turbulent non-premixed flames is not fully understood and several models have been developed to describe the turbulence–chemistry interaction. This work compares the performance of the multiple Representative Interactive Flamelet (mRIF) model, the Volume Reactor Fraction Model (VRFM), and the Well-Stirred reactor (WS) model in describing such flames. The predicted ignition delay and flame lift-off length of n-heptane sprays are compared to experimental results published within the Engine Combustion Network (ECN). All of the models predict the trend of ignition delay reasonably well. At a low gas pressure (42 bar) the ignition delay is overpredicted compared to the experimental data, but the difference between the models is not significant. However, the predicted lift-off lengths differ. At high pressure (87 bar) the difference between the models is small. All models slightly underpredict the lift-off length compared to the experimental data. At low gas pressure (42 bar) the mRIF model gives the best results. The VRFM and WS models predict excessively short lift-off lengths, but the VRFM model gives better results than the WS model. The flame structures of the models are also compared. The WS model and the VRFM model yield a well defined flame stabilisation point whereas the mRIF model does not. The flame of the mRIF model is more diffuse and the model is not able to predict flame propagation. All models were able to predict the experimental trends in lift-off and ignition delay, but certain differences between them are demonstrated.     

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.     

A comparison of computational and experimental lift-off heights of coflow laminar diffusion flames

As a sensitive marker of changes in flame structure, the number densities of excited-state CH (denoted CH^*), and excited-state OH (denoted OH^*) are imaged in coflow laminar diffusion flames. Measurements are made both in normal gravity and on the NASA KC-135 reduced-gravity aircraft. The spatial distribution of these radicals provides information about flame structure and lift-off heights that can be directly compared with computational predictions. Measurements and computations are compared over a range of buoyancy and fuel dilution levels. Results indicate that the lift-off heights and flame shapes predicted by the computations are in excellent agreement with measurement for both normal gravity (1g) and reduced gravity flames at low dilution levels. As the fuel mixture is increasingly diluted, however, the 1g lift-off heights become underpredicted. This trend continues until the computations predict stable flames at highly dilute fuel mixtures beyond the 1g experimental blow-off limit. To better understand this behavior, an analysis was performed, which indicates that the lift-off height is sensitive to the laminar flame speed of the corresponding premixed mixture at the flame edge. By varying the rates of two key “flame speed” controlling reactions, we were able to modify the predicted lift-off heights so as to be in closer agreement with the experiments. The results indicate that reaction sets that work well in low dilution systems may need to be modified to accommodate high dilution flames.     

喷孔内流对射流初始扰动影响的大涡模拟研究

本文利用大涡模拟方法对喷孔内流和近喷孔区域的燃油雾化过程进行数值模拟研究。为了能更加准确地模拟喷孔内流,本文根据实际喷油器结构设定了计算域,能够考虑针阀、油囊和喷孔结构对喷孔内流的影响。结果表明:在考虑了具体的喷油器结构参数情况下,采用气液两相流大涡模拟能够获得与实验非常接近的喷束结构;喷孔内流的非轴对称分布是引起射流非轴对称型扰动的原因;这种非轴对称扰动使喷束变为非轴对称结构。结果还显示这种非轴对称扰动具有远大于表面波的振幅和波长。     

基于LES-PDF方法的双旋流模型燃烧室数值模拟

文章使用基于LES-PDF方法的AECSC两相程序,对模型燃烧室GTMC进行了数值模拟,以此验证AECSC程序对燃烧室模拟的可行性和可信度,并对旋流燃烧室的流动和燃烧特性进行分析.首先,分别用商用软件Fluent 18.1和AECSC程序的LES方法对GTMC的冷态工况进行了模拟.与实验结果相比,程序计算的轴向、径向、切向速度峰值的相对误差在大多数统计点上在20%以内,3个方向速度峰值位置的相对误差基本都在10%以内.同时发现,AECSC程序的计算结果在一定程度上比Fluent的计算结果更接近实验值.文章进一步使用AECSC程序对GTMC的热态工况进行了模拟.在计算结果中,整体的时间平均温度在数值和分布上和实验结果很接近,并很好地再现了实验的"V形"火焰和内外低温回流区,但是高温区出现的位置相比实验结果有些靠前.总体来说,热态结果和实验结果比较吻合,尤其在液雾的模拟方面,计算结果和实验结果相比有着很好的一致性.这说明将LES和PDF相结合的方法在模拟湍流燃烧方面有较强的优势,可以成为未来的研究和发展方向.      

Simulation of spray combustion in a lean-direct injection combustor

Large-eddy simulation (LES) of a liquid-fueled lean-direct injection (LDI) combustor is carried out by resolving the entire inlet flow path through the swirl vanes and the combustor. A localized dynamic subgrid closure is combined with a subgrid mixing and combustion model so that no adjustable parameters are required for both non-reacting and reacting LES. Time-averaged velocity predictions compare well with the measured data. The unsteady flow features that play a major role in spray dispersion, fuel–air mixing and flame stabilization are identified from the simulation data. It is shown that the vortex breakdown bubble (VBB) is smaller with more intense reverse flow when there is heat release. The swirling shear layer plays a major role in spray dispersion and the VBB provides an efficient flameholding mechanism to stabilize the flame.     

Experimental and numerical characterization of a turbulent spray flame

A turbulent ethanol spray flame is characterized through quantitative experiments using laser-based imaging techniques. The data set is used to validate a numerical code for the simulation of spray combustion. The spray burner has been designed to generate a stable flame without the use of a bluff body or a pilot flame facilitating numerical simulations. The experiments include spatially-resolved measurements of droplet sizes (Mie/LIF-dropsizing and PDA), droplet velocity (PDA), liquid-phase temperature (2-color LIF temperature imaging with Rhodamine B) and gas-phase temperature (multi-line NO-LIF temperature imaging). The measurements close to the nozzle exit are used to determine the initial conditions for numerical simulations. An Eulerian–Lagrangian model including spray flamelet modeling is applied to calculate the development of the spray. Good agreement with the experimental data is found. The experimental data set and the numerical results will be published on a website to allow other groups to evaluate their experimental and/or numerical data.     

A spray flamelet/progress variable approach combined with a transported joint PDF model for turbulent spray flames

A spray flamelet/progress variable approach is developed for use in spray combustion with partly pre-vaporised liquid fuel, where a laminar spray flamelet library accounts for evaporation within the laminar flame structures. For this purpose, the standard spray flamelet formulation for pure evaporating liquid fuel and oxidiser is extended by a chemical reaction progress variable in both the turbulent spray flame model and the laminar spray flame structures, in order to account for the effect of pre-vaporised liquid fuel for instance through use of a pilot flame. This new approach is combined with a transported joint probability density function (PDF) method for the simulation of a turbulent piloted ethanol/air spray flame, and the extension requires the formulation of a joint three-variate PDF depending on the gas phase mixture fraction, the chemical reaction progress variable, and gas enthalpy. The molecular mixing is modelled with the extended interaction-by-exchange-with-the-mean (IEM) model, where source terms account for spray evaporation and heat exchange due to evaporation as well as the chemical reaction rate for the chemical reaction progress variable. This is the first formulation using a spray flamelet model considering both evaporation and partly pre-vaporised liquid fuel within the laminar spray flamelets. Results with this new formulation show good agreement with the experimental data provided by A.R. Masri, Sydney, Australia. The analysis of the Lagrangian statistics of the gas temperature and the OH mass fraction indicates that partially premixed combustion prevails near the nozzle exit of the spray, whereas further downstream, the non-premixed flame is promoted towards the inner rich-side of the spray jet since the pilot flame heats up the premixed inner spray zone. In summary, the simulation with the new formulation considering the reaction progress variable shows good performance, greatly improving the standard formulation, and it provides new insight into the local structure of this complex spray flame.