Large eddy simulations of partially premixed ethanol dilute spray flames using the flamelet generated manifold model

The paper presents Large Eddy Simulations (LESs) for the Sydney ethanol piloted turbulent dilute spray flames ETF2, ETF6, and ETF7. The Flamelet Generated Manifold (FGM) approach is employed to predict mixing and burning of the evaporating fuel droplets. A methodology to match the experimental inflow spray profiles is presented. The spray statistical time-averaged results show reasonable agreement with mean and RMS data. The Particle Size Distribution (PSD) shows a good match downstream of the nozzle exit and up to x/D = 10. At x/D = 20 and 30 the PSD is under-predicted for droplets with mean diameter D10 > 20μm and over-predicted for the smaller size droplets. The simulations reasonably predict the reported mean flame structure and length. The effect of increasing the carrier velocity (ETF2–ETF7) or decreasing the liquid fuel injection mass flow rate (ETF2–ETF6) is found to result in a leaner, shorter flame and stronger spray–flow interactions. Higher tendency to local extinction is observed for ETF7 which is closer to blow-off compared to ETF2 and has higher scalar dissipation rates, higher range of Stokes number, and faster droplet response. The possible sources of LES-FGM deviations from the measurements are discussed and highlighted. In particular, the spray time-averaged statistical error contribution is quantified and the impact of the inflow uncertainty is studied. Sensitivity analysis to the pre-vaporized nozzle fuel mass fraction show that such small inflow perturbations (by ±?2% for the ETF2 flame) have a strong impact on the flame structure, and the droplets’ dynamics. Conditional scatter plots show that the flame exhibits wide range of mixing conditions and bimodal mixing lines particularly at upstream locations (x/D?相似文献   

Partially premixed flamelet in LES of acetone spray flames

An effective partially premixed flamelet model for large eddy simulation (LES) of turbulent spray combustion is formulated. Different flame regimes are identified with a flame index defined by budget terms in a 2-D multi-phase flamelet formulation, and the application in LES of partially pre-vaporized spray flames shows a favorable agreement with experiments. Simulations demonstrate that, compared to the conventional single-regime flamelets, the present partially premixed flamelet formulation shows its ability in capturing the subgrid regime transitions, yielding a well prediction of peak gas temperature and the downstream flame spreading. A propagating premixed flame front is found coupled with a trailing diffusion burning through the spray evaporation, and the spray effect on regime discrimination is manifested with transport budget analysis. A two-phase regime indicator is then proposed, by which the evaporation-dictated regime is properly described. Its intended use will rely on both gas and spray flamelet structures.     

Effects of fine fuel droplets on a laminar flame stabilized in a partially prevaporized spray stream

A partially prevaporized spray burner was developed to investigate the interaction between fuel droplets and a flame. Monodispersed partially prevaporized ethanol sprays with narrow diameter distribution were generated by the condensation method using rapid pressure reduction of a saturated ethanol vapor–air mixture. A tilted flat flame was stabilized at the nozzle exit using a hot wire. Particle tracking velocimetry (PTV) was applied to measurements of the droplet velocity; the laminar burning velocity was obtained from gas velocity derived from the droplet velocity. Observations were made of flames in partially prevaporized spray streams with mean droplet diameters of 7 μm and the liquid equivalence ratios of 0.2; the total equivalence ratio was varied. In all cases, a sharp vaporization plane was observed in front of the blue flame. Flame oscillation was observed on the fuel-rich side. At strain rates under 50 s⁻¹, the change in the burning velocity with the strain rate is small in fuel-lean spray streams. In spray streams of 0.7 and 0.8 in the total equivalence ratio, burning velocity increases with strain rates of greater than 50 s⁻¹. However, in spray streams with 0.9 and 1.0 in the total equivalence ratio, burning velocity decreases as the strain rate increases. At strain rates greater than 80 s⁻¹, burning velocity decreases with an increased gas equivalence ratio. The effect of mean droplet diameter, and the entry length of droplets into a flame on the laminar burning velocity, were also investigated to interpret the effect of the strain rate on the laminar burning velocity of partially prevaporized sprays.     

An experimental study of partially premixed flame sound

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

Experimental and numerical study of a conical turbulent partially premixed flame

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

Structure of partially premixed n-heptane–air counterflow flames

To avoid the complexities associated with the droplet/vapor transport and nonuniform evaporation processes, a fundamental investigation of liquid fuel combustion in idealized configurations is very useful. An experimental–computational investigation of prevaporized n-heptane nonpremixed and partially premixed flames established in a counterflow burner is described. There is a general agreement between various facets of our nonpremixed flame measurements and the literature data. The partially premixed flames are characterized by a double flame structure. This becomes more distinct as the strain rate decreases and partial premixing increases, which also increases the separation distance between the two reaction zones. The peak partially premixed flame temperature increases with increasing premixing of the fuel stream. The peak CO₂ and H₂O concentrations are relatively insensitive to partial premixing. The CO and H₂ peak concentrations on the premixed side increase as the fuel-side equivalence ratio decreases. These species are transported to the nonpremixed reaction zone where they oxidize. The C₂ species have peaks in the premixed reaction zone. The concentrations of olefins are ten times larger than those of the corresponding paraffins. The oxidizer is present in partially premixed flames throughout the combustion system and there are no regions characterized by simultaneous high temperature and high fuel concentration. As a result, pyrolysis reactions leading to soot formation are greatly diminished.     

Measurements of flame orientation and scalar dissipation in turbulent partially premixed methane flames

Three turbulent flames were studied using a new experimental facility developed at Sandia National Laboratories. Line imaging of Raman and Rayleigh scattering and CO laser-induced fluorescence (LIF) yielded information on all major species, temperature, mixture fraction, and a 1D surrogate measure of scalar dissipation. Simultaneously, crossed planar OH LIF imaging provided information on the instantaneous flame orientation, allowing estimation of the full 3D (flame-normal) scalar dissipation rate. The three flames studied were methane–air piloted jet flames (Sandia flames C, D, and E), which cover a range in Reynolds number from 13,400 to 33,600. The statistics of the instantaneous flame orientation are examined in the different flames, with the purpose of studying the prevailing kinematics of isoscalar contours. The 1D and 3D results for scalar dissipation rate are examined in detail, both in the form of conditional averages and in the form of probability density functions. The effect of overall strain and Reynolds number on flame suppression and eventual extinction is also investigated, by examining the doubly conditional statistics of temperature in the form of S-shaped curves. This latter analysis reveals that double conditioning of temperature on both mixture fraction and scalar dissipation does not collapse the data from these flames onto the same curve at low scalar dissipation rates, as might be expected from simple flamelet concepts.     

Identification of droplet burning modes in lean, partially prevaporized swirl-stabilized spray flames

Experimental and numerical investigations of single droplet burning modes in a lean, partially prevaporized swirl-stabilized spray flame are reported. In the experiment single droplet flames have been visualized by CH-PLIF and simultaneous recording of the Mie signal. Two single droplet burning modes were identified: the envelope flame is a spherical diffusion flame burning at near-stoichiometric conditions. The wake flame is a potentially lean, partially premixed flame located downstream of the droplet. The droplet burning mode is of practical relevance, since it has significant impact on NO formation due to incomplete prevaporization.The droplet burning mode is determined by the ratio of chemical and convective time scales. The convective time scale is related to the droplet slip velocity. The impact of turbulent gas phase velocity fluctuations on droplet mechanics and droplet burning is discussed, based on a previous numerical investigation. In the present study the droplet slip velocity was measured with the 3D Phase Doppler (3D-PD) technique. For the measured slip velocities and ambient conditions in the hot gas region of the spray flame, simulations of single droplet burning were performed utilizing detailed models for chemical reaction, diffusive transport and vaporization. An agreement between the droplet burning modes predicted by the simulation and the droplet burning modes observed in the experiments was found.     

Influence of micro-explosions on the stability of laminar premixed water-in-fuel emulsion spray flames

A model is presented for a one-dimensional laminar premixed flame, propagating into a rich, off-stoichiometric, fresh homogenous mixture of water-in-fuel emulsion spray, air and inert gas. Due to its relatively large latent heat of vaporisation, the water vapour acts to cool the flame that is sustained by the prior release of fuel vapour. To simplify the inherent complexity that characterises the analytic solution of multi-phase combustion processes, the analysis is restricted to fuel-rich laminar premixed water-in-fuel flames, and assumes a single-step global chemical reaction mechanism. The main purpose is to investigate the steady-state burning velocity and burnt temperature as functions of parameters such as initial water content in the emulsified droplet and total liquid droplet loading. In particular, the influence of micro-explosion of the spray’s droplets on the flame’s characteristics is highlighted for the first time. Steady-state analytical solutions are obtained and the sensitivity of the flame temperature and the flame propagating velocity to the initial water content of the micro-exploding emulsion droplets is established. A linear stability analysis is also performed and reveals the manner in which the micro-explosions influence the neutral stability boundaries of both cellular and pulsating instabilities.     

Flame structure and flame spread rate over a solid fuel in partially premixed atmospheres

We have investigated the downward flame spread over a thin solid fuel. Hydrogen, methane, or propane, included in the gaseous product of pyrolysis reaction, is added in the ambient air. The fuel concentration is kept below the lean flammability limit to observe the partially premixing effect. Both experimental and numerical studies have been conducted. Results show that, in partially premixed atmospheres, both blue flame and luminous flame regions are enlarged, and the flame spread rate is increased. Based on the flame index, a so-called triple flame is observed. The heat release rate ahead of the original diffusion flame is increased by adding the fuel, and its profile is moved upstream. Here, we focus on the heat input by adding the fuel in the opposed air, which could be a direct factor to intensify the combustion reaction. The dependence of the flame spread rate on the heat input is almost the same for methane and propane/air mixtures, but larger effect is observed for hydrogen/air mixture. Since the deficient reactant in lean mixture is fuel, the larger effect of hydrogen could be explained based on the Lewis number consideration. That is, the combustion is surely intensified for all cases, but this effect is larger for lean hydrogen/air mixture (Le < 1), because more fuel diffuses toward the lean premixed flame ahead of the original diffusion flame. Resultantly, the pyrolysis reaction is promoted to support the higher flame spread rate.     

A comparative study on partially premixed combustion (PPC) and reactivity controlled compression ignition (RCCI) in an optical engine

Partially premixed combustion (PPC) and reactivity controlled compression ignition (RCCI) are two new combustion modes in compression-ignition (CI) engines. However, the detailed in-cylinder ignition and flame development process in these two CI modes were not clearly understood. In the present study, firstly, the fuel stratification, ignition and flame development in PPC and RCCI were comparatively studied on a light-duty optical engine using multiple optical diagnostic techniques. The overall fuel reactivity (PRF number) and concentration (fuel-air equivalence ratio) were kept at 70 and 0.77 for both modes, respectively. Iso-octane and n-heptane were separately used in the port-injection (PI) and direct-injection (DI) for RCCI, while PRF70 fuel was introduced through direct-injection (DI) for PPC. The DI timing for both modes was fixed at –25°CA ATDC. Secondly, the combustion characteristics of PPC and RCCI with more premixed charge were explored by increasing the PI mass fraction for RCCI and using the split DI strategy for PPC. In the first part, results show that RCCI has shorter ignition delay than PPC due to the fuel reactivity stratification. The natural flame luminosity, formaldehyde and OH PLIF images prove that the flame front propagation in the early stage of PPC can be seen, while there is no distinct flame front propagation in RCCI. In the second part, the higher premixed ratio results in more auto-ignition sites and faster combustion rate for PPC. However, the higher premixed ratio reduces the combustion rate in RCCI mode and the flame front propagation can be clearly seen, the flame speed of which is similar to that in spark ignition engines but lower than that in PPC. It can be concluded that the ratio of flame front propagation and auto-ignition in RCCI and PPC can be modulated by the control over the fuel stratification degree through different fuel-injection strategies.     

Direct numerical simulations of rich premixed turbulent n-dodecane/air flames at diesel engine conditions

Rich premixed turbulent n-dodecane/air flames at diesel engine conditions are analyzed using direct numerical simulations. The conditions correspond to a parametric variation of the Engine Combustion Network Spray A (pressure 60 atm; oxidizer oxygen level and temperature 21% and 900 K, respectively; fuel temperature 363 K). Three simulations with equivalence ratios of 3, 5, and 7 are performed with a Karlovitz number (Ka, based on flame time) of order 100 to match the estimated Ka of the rich premixed combustion region in Spray A. At these conditions, the reference laminar flames exhibit a complex structure which involves both low-temperature chemistry (LTC) and high-temperature chemistry over a wide range of length scales. In the presence of turbulence, the flame structure is strongly affected in physical space and the reaction zone exhibits a very complex structure in which broken, distributed, and thin regions co-exist, especially for the leanest case. However, the contribution of the LTC pathway is only weakly affected by turbulence. In progress variable space, the mean flame structure, including the chemical source terms, is found to match remarkably well that of the corresponding unity Lewis number laminar flame, particularly for the $?=$ 3 and 5 cases. This behavior is attributed to the strong turbulent mixing occurring throughout the flames/reaction zones, which suppresses differential diffusion effects. Nevertheless, large conditional fluctuations around the mean chemical source terms are identified. These are found to correlate very well with radical species mass fractions such as OH. In addition, a similar functional dependence is obtained from counterflow laminar flames. As such, it appears from these results that laminar flame models have a potential to be used to represent the thermochemical state of rich premixed turbulent flames under diesel engine conditions.     

A theoretical study of premixed turbulent flame development

Flame development in a statistically stationary and uniform, planar, one-dimensional turbulent flow is theoretically studied. A generalized balance equation for the mean combustion progress variable, which includes turbulent diffusion and pressure-driven transport terms, as well as the mean rate of product creation, is introduced and analyzed by invoking the sole assumption of a self-similar flame structure, well-supported by numerous experiments. The assumption offers the opportunity to simplify the problem by splitting the aforementioned partial differential equation into two ordinary differential equations, which separately model spatial variations of the progress variable and time variations of flame speed and thickness. The self-similar profile of the progress variable, obtained in numerous experiments, is theoretically predicted. Closures of the normalized pressure-driven transport term and mean rate of product creation are obtained. The closed balance equation shows that turbulent diffusion dominates during the initial stage of flame development, followed by the transition to counter-gradient transport in a sufficiently developed flame. A criterion of the transition is derived. The transition is promoted by the heat release and pressure-driven transport. Fully developed mean flame brush thickness and speed are shown to decrease when either density ratio or pressure-driven transport increases. Solutions for the development of the thickness are obtained. The development is accelerated by the pressure-driven transport and heat release.     

燃料敏感性对内燃机部分预混燃烧(PPC)光谱特性的影响研究

在一台光学发动机上,利用火焰高速成像技术和自发光光谱分析法,研究了燃料敏感性(S)为0和6时对发动机缸内火焰发展和燃烧发光光谱的影响。试验过程中,通过改变喷油时刻(SOI=-25,-15和-5°CA ATDC)使燃烧模式从部分预混燃烧过渡到传统柴油燃烧模式。通过使用正庚烷、异辛烷、乙醇混合燃料来改变燃料敏感性。结果表明,在PPC模式下(-25°CA ATDC),火焰发展过程是从近壁面区域开始着火,而后向燃烧室中心发展,即存在类似火焰传播过程,同时在燃烧室下部未燃区域也形成新的着火自燃点。敏感性对燃烧相位影响较大,对缸内燃烧火焰发展历程影响较小;高敏感性燃料OH和CH带状光谱出现的时刻推迟,表明高敏感性燃料高温反应过程推迟,且光谱强度更低,表明碳烟辐射强度减弱。在PPC到CDC之间的过渡区域(-15°CA ATDC),燃烧火焰发光更亮,燃烧反应速率比-25°CA ATDC时刻的反应速率更快。高、低敏感性燃料对缸压放热率的影响规律与-25°CA ATDC相近,此时的燃烧反应更剧烈,放热率更高,碳烟出现时刻更早。该喷油时刻下的光谱强度高于PPC模式下的光谱强度,说明此时的CO氧化反应与碳烟辐...     

A flamelet analysis of the burning velocity of premixed turbulent expanding flames

Direct numerical simulation is a very powerful tool to evaluate the validity of new models and theories for turbulent combustion. In this paper, direct numerical simulations of spherically expanding premixed turbulent flames in the corrugated flamelet regime are performed. The flamelet-generated manifold method is used to deal with detailed reaction kinetics. The numerical method is validated for both laminar and turbulent expanding flames. The computational results are analyzed by using an extended flame stretch theory. It is investigated whether this theory is able to describe the influence of flame stretch and curvature on the local burning velocity of the flame. If the full profiles of flame stretch and curvature through the flame front are included in the theory, the local mass burning rate is predicted accurately. The influence of several approximations, which are used in other existing theories, is studied. When flame stretch is assumed to be constant through the flame front or when curvature of the flame front is neglected, the theory fails to predict the local mass burning rate.     

Large scale effects in weakly turbulent premixed flames

In this study we numerically investigate large scale premixed flames in weakly turbulent flow fields. A large scale flame is classified as such based on a reference hydrodynamic lengthscale being larger than a neutral (cutoff) lengthscale for which the hydrodynamic or Darrieus–Landau (DL) instability is balanced by stabilizing diffusive effects. As a result, DL instability can develop for large scale flames and is inhibited otherwise. Direct numerical simulations of both large scale and small scale three-dimensional, weakly turbulent flames are performed at constant Karlovitz and turbulent Reynolds number, using two paradigmatic configurations, namely a statistically planar flame and a slot Bunsen flame. As expected from linear stability analysis, DL instability induces its characteristic cusp-like corrugation only on large scale flames. We therefore observe significant morphological and topological differences as well as DL-enhanced turbulent flame speeds in large scale flames. Furthermore, we investigate issues related to reaction rate modeling in the context of flame surface density closure. Thicker flame brushes are observed for large scale flames resulting in smaller flame surface densities and overall larger wrinkling factors.     

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.     

Investigation of flame structure and burning intensity of partially premixed methane enrichment of syngas using OH-PLIF and kinetic simulation

Various experiments were conducted to study the combustion characteristics of partially premixed methane enrichment of syngas by using the OH-PLIF technique. Experiments were conducted on a co-flow burner, and the methane concentration (X_CH4 = CH₄/(H₂+CO+CH₄)) was varied from 0 to 20%, the overall equivalence ratio was varied from 0.4 to 1.2 and the inner equivalence ratio was varied from 1.5 to 3.5. Kinetic simulation was conducted by using OPPDIF module of CHEMKIN-Pro software. Results show that an increase in X_CH4 and ?_overall weakens the OH signal intensity. Adding methane into the fuel greatly increases the height of the inner flame front, and the increase of methane concentration has a negative effect on flame propagation speed. Meanwhile, simulation results remain consistent with the experiments. The main OH radical production reaction changes from R46: H+HO₂ = 2OH to R38: H+O₂ = O+OH when methane concentration contained in the fuel mixture increases. Sensitivity analysis also indicates that reaction which plays a dominant effect on temperature changes with the increase of methane concentration.     

Numerical study of transient evolution of lifted jet flames: partially premixed flame propagation and influence of physical dimensions

Three-dimensional (3D) unsteady Reynolds-averaged Navier–Stokes simulations of a spark-ignited turbulent methane/air jet flame evolving from ignition to stabilisation are conducted for different jet velocities. A partially premixed combustion model is used involving a correlated joint probability density function and both premixed and non-premixed combustion mode contributions. The 3D simulation results for the temporal evolution of the flame's leading edge are compared with previous two-dimensional (2D) results and experimental data. The comparison shows that the final stabilised flame lift-off height is well predicted by both 2D and 3D computations. However, the transient evolution of the flame's leading edge computed from 3D simulation agrees reasonably well with experiment, whereas evident discrepancies were found in the previous 2D study. This difference suggests that the third physical dimension plays an important role during the flame transient evolution process. The flame brush's leading edge displacement speed resulting from reaction, normal and tangential diffusion processes are studied at different typical stages after ignition in order to understand the effect of the third physical dimension further. Substantial differences are found for the reaction and normal diffusion components between 2D and 3D simulations especially in the initial propagation stage. The evolution of reaction progress variable scalar gradients and its interaction with the flow and mixing field in the 3D physical space have an important effect on the flame's leading edge propagation.