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.     

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.     

小火焰模型在贫燃预混火焰中的研究

由层流小火焰库引入详细化学反应机理,通过简化的PDF方法计算组分浓度、平均温度和密度等变量,以钝体火焰稳定燃烧室和某燃气轮机上的燃烧室为例,模拟甲烷/空气贫燃条件下预混燃烧的平均火焰位置和火焰厚度,计算结果与实验结果吻合良好,这表明此方法能够较好计算出平均湍流火焰的主要特征。     

Parametric study on a compound-drop spray flame

The introduction of compound-drop spray in a combustion system is a new concept. These droplets bear two gasification stages to cause an integral positive or negative effect on a premixed flame to raise or lower the local temperature of the gasification region. In this paper, we adopt a compound drop which contains a water core encased by a layer of shell fuel. A one-dimensional homogeneous lean or rich premixed flame with the dilute compound-drop spray was investigated by using large activation energy asymptotic analysis. The compound-drop spray burning mode was defined and divided into completely pre-vaporised burning (CPB), shell pre-vaporised burning (SPB) and shell partially pre-vaporised (SPP) burning modes by way of the gasification zones of the shell fuel and the core water relative to the flame position. The influences of the initial droplet radius, the shell-fuel mass fraction and the liquid loading of the compound-drop spray on the lean and rich flames were analysed. By means of the normalisation parameter of flame propagation mass flux (), enhancement, suppression or extinction of the compound-drop spray flame can be represented clearly. Furthermore, from the observation of extinction, the necessary conditions of extinction of a lean spray flame by the internal heat transfer are that the spray is a negative effect and causes a sufficient heat loss rate at flame sheet downstream side. For a rich spray flame, three extinction patterns were observed; they occur in SPP, SPB or at the critical SPB mode, but do not in CPB. The extinction maps of the compound-drop spray demarcate the patterns and also indicate the limitations and corresponding conditions of the flame extinction.     

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.     

Laminar premixed spray flame analysis with thermally sensitive intermediate kinetics

A new thermo-diffusive analysis of one-dimensional laminar lean or rich off-stoichiometric premixed spray flames has been performed using a chain branching/chain breaking chemical kinetic scheme and under the assumption that the fuel droplets evaporate in a sharp front. The sensitivity of the flame structure, speed and the location of the evaporation front to the initial droplet load have been demonstrated. A linear stability analysis reveals the way in which the spray's presence modifies the neutral stability curves.     

Prediction of soot and thermal radiation in a model gas turbine combustor burning kerosene fuel spray at different swirl levels

Combustion of kerosene fuel spray has been numerically simulated in a laboratory scale combustor geometry to predict soot and the effects of thermal radiation at different swirl levels of primary air flow. The two-phase motion in the combustor is simulated using an Eulerian–Lagragian formulation considering the stochastic separated flow model. The Favre-averaged governing equations are solved for the gas phase with the turbulent quantities simulated by realisable k–? model. The injection of the fuel is considered through a pressure swirl atomiser and the combustion is simulated by a laminar flamelet model with detailed kinetics of kerosene combustion. Soot formation in the flame is predicted using an empirical model with the model parameters adjusted for kerosene fuel. Contributions of gas phase and soot towards thermal radiation have been considered to predict the incident heat flux on the combustor wall and fuel injector. Swirl in the primary flow significantly influences the flow and flame structures in the combustor. The stronger recirculation at high swirl draws more air into the flame region, reduces the flame length and peak flame temperature and also brings the soot laden zone closer to the inlet plane. As a result, the radiative heat flux on the peripheral wall decreases at high swirl and also shifts closer to the inlet plane. However, increased swirl increases the combustor wall temperature due to radial spreading of the flame. The high incident radiative heat flux and the high surface temperature make the fuel injector a critical item in the combustor. The injector peak temperature increases with the increase in swirl flow mainly because the flame is located closer to the inlet plane. On the other hand, a more uniform temperature distribution in the exhaust gas can be attained at the combustor exit at high swirl condition.     

Problems of predicting turbulent burning rates

Different approaches to the modelling of turbulent combustion first are reviewed briefly. A unified, stretched flamelet approach then is presented. With Reynolds stress modelling and a generalized probability density function (PDF) of strain rate, it enables a source term, in the form of a probability of burning function, P_b, to be expressed as a function of Markstein numbers and the Karlovitz stretch factor. When P_b is combined with some turbulent flame fractal considerations, an expression is obtained for the turbulent burning velocity. When it is combined with the profile of the unstretched laminar flame volumetric heat release rate plotted against the reaction progress variable and the PDF of the latter, an expression is obtained for the mean volumetric turbulent heat release rate. Through these relationships experimental values of turbulent burning velocity might be used to evaluate P_b and hence the CFD source term, the mean volumetric heat release rate.

Different theoretical expressions for the turbulent burning velocity, including the present one, are compared with experimental measurements. The differences between these are discussed and this is followed by a review of CFD applications of these flamelet concepts to premixed and non-premixed combustion. The various assumptions made in the course of the analyses are scrutinized in the light of recent direct numerical simulations of turbulent flames and the applications to the flames of laser diagnostics. Remaining problem areas include a sufficiently general combination of strain rate and flame curvature PDFs to give a single PDF of flame stretch rate, the nature of flame quenching under positive and negative stretch rates, flame responses to changing stretch rates and the effects of flame instabilities.     

Propagation and extinction of an unsteady spherical spray flame front

Experimental evidence seems to indicate that the life of a laminar spherical flame front propagating through a fresh mixture of air and liquid fuel droplets can be roughly split into three stages: (1) ignition, (2) radial propagation with a smooth flame front and (3) propagation with flame front cellularization and/or pulsation. In this work, the second stage is analysed using the slowly varying flame approach, for a fuel rich flame. The droplets are presumed to vaporize in a sharp front ahead of the reaction front. Evolution equations for the flame and evaporation fronts are derived. For the former the combined effect of heat loss due to droplet vaporization and radiation plays a dominant explicit role. In addition, the structure of the evaporation front is deduced using asymptotics based on a large parameter associated with spray vaporization. Numerical calculations based on the analysis point to the way in which the spray modifies conditions for flame front extinction. Within the framework of the present simplified model the main relevant parameters turn out to be the initial liquid fuel load in the fresh mixture and/or the latent heat of vaporization of the fuel.     

Combustion of partially premixed spray jets

Gas turbines, liquid rocket motors, and oil-fired furnaces utilize the spray combustion of continuously injected liquid fuels. In most cases, the liquid spray is mixed with an oxidizer prior to combustion, and further oxidizer is supplied from the outside of the spray to complete diffusion combustion. This rich premixed spray is called “partially premixed spray.” Partially premixed sprays have not been studied systematically although they are of practical importance. In the present study, the burning behavior of partially premixed sprays was experimentally studied with a newly developed spray burner. A fuel spray and an oxidizer, diluted with nitrogen, was injected into the air. The overall equivalence ratio of the spray jet was set larger than unity to establish partially premixed spray combustion. In the present burner, the mean droplet diameter of the atomized liquid fuel could be varied without varying the overall equivalence ratio of the spray jet. Two combustion modes with and without an internal flame were observed. As the mean droplet diameter was increased or the overall equivalence ratio of the spray jet was decreased, the transition from spray combustion only with an external group flame to that with the internal premixed flame occurred. The results suggest that the internal flame was supported by flammable mixture through the vaporization of fine droplets, and the passage of droplet clusters deformed the internal flame and caused internal flame oscillation. The existence of the internal premixed flame enhanced the vaporization of droplets in the post-premixed-flame zone within the external diffusion flame.     

A comparison of different approaches to integrate flamelet tables with presumed-shape PDF in flamelet models for turbulent flames

Flamelet models for turbulent combustion modelling make use of presumed-shape probability density functions (PDFs) for integrating laminar flamelet solutions to obtain an integrated flamelet table that can readily be used for turbulent flame calculations. The existence of non-unique approaches for such an integration has rarely been investigated before. For the first time, this work studies systematically the non-uniqueness of the flamelet table integration approaches. A flamelet model called the flamelet/progress variable model is used in the study, although the issue exists generally in many other flamelet models. Two classes of table integration approaches are investigated, one preserving the laminar flamelet structures during integration and the other not. Three different table integration approaches are examined and compared in detail to provide a thorough understanding of the different approaches. A partially stirred reactor is used as a test case for examining the different approaches. A method based on the transported PDF method is also employed to provide a reference for the assessment of the different flamelet table integration approaches. It is found in general that the flamelet preserving integration approach yields a more reasonable joint PDF of the mixture fraction and the progress variable, and the prediction results are closer to the referenced transported PDF results.     

Analysis of a spray diffusion flame and its extinction in a shear layer

A new semi-analytical solution for a laminar spray diffusion flame in the shear layer between fuel and oxidant streams is developed. The Stokes number is identified as a small spray droplet-related parameter to be used in a perturbation analysis of the liquid phase governing equations. Appropriate specification of an additional parameter ensures that similarity is achieved so that the concentration field of the liquid in the spray can be readily found. The coupled gas-phase equations are treated using the usual inverse of the large Zeldovitch number for the asymptotic analysis. Numerical results demonstrate how the distribution of the liquid phase in the developing shear layer between two unidirectional gas streams flowing over one another with (the possibility of) dissimilar velocities in their respective free-streams influences the flame shape, location, fuel vapour and temperature fields. An extinction analysis enables a parametric mapping of conditions for extinguishment of the spray diffusion flames to be drawn.     

Polydisperse spray diffusion flames in oscillating flow

The phenomenon of droplet clustering or grouping found when a spray of droplets is moving in an oscillating host flow field is investigated for the case of a polydisperse spray that fuels a laminar co-flow diffusion flame. A mathematical solution is developed for the liquid phase based on use of small Stokes numbers for size sections into which the polydisperse spray size distribution is divided. Droplet clustering in the oscillatory flow field is accounted for by constructing a special model for the sectional vaporization Damkohler numbers in accordance with droplet size. Combining this with a formal solution for a gas phase Schvab-Zel'dovich variable yields the means whereby flame dynamics can be described. Results calculated from this solution demonstrate that preferential droplet size behaviour (with smaller droplets tending to cluster to a greater extent and reduce the vaporization Damkohler number more than larger ones) can have a major impact on the flame dynamics through local droplet enrichment with attendant consequences on the production of fuel vapour. The dynamics of the sort of flame (over- or under-ventilated) and the occurrence of flame pinching leading to multiple flame sheets are altered under these circumstances. However, potential control of the actual initial spray polydispersity may reduce the intensity of such effects.     

PDF mixing time scales for premixed combustion in the laminar flame limit

A new mixing time scale for PDF calculations of premixed combustion in the laminar flame limit is introduced. It is based on the characteristics of the “random walk” diffusion process and accurately captures scalar micro-mixing such that physical features of the premixed flames are preserved. The speed of a freely-propagating laminar flame can be captured accurately. Extreme stochastic events may, however, lead to flame acceleration, but flame stability can be recovered with the introduction of a conditioning variable that relaxes towards an Eulerian reference field. With conditioning the predicted flame speed is quite insensitive towards the exact value of the only modelling parameter. This modelling parameter then allows to control the minor scalar fluctuations and deviations from a flamelet structure which is thought to be one of the key features of conditioning methods.     

甲醇喷雾湍流反应流的JPDF数值模拟

喷雾湍流燃烧过程中,液滴、湍流和化学反应之间强烈耦合,物理化学机理非常复杂。本文将速度-标量-频率联合概率密度函数JPDF输运方程方法应用于两相喷雾湍流反应流问题,利用火焰面模型解耦流动和化学反应动力学的耦合关系,建立起相应的数值计算模型。采用Monte-Carlo数值计算方法,针对澳大利亚悉尼大学Masri等人以甲醇为燃料所进行的湍流喷雾燃烧值班火焰这一试验进行了数值模拟,通过与Fluent下的计算结果及试验结果的对比分析,验证了本文所建模型的准确性。     

钝体驻定湍流扩散火焰的数值研究——燃烧模型比较

分别采用标量联合的概率密度函数方法、稳态火焰面模型、Euler非稳态火焰面模型和基于有限体积/Monte Carlo混合算法的完备PDF模型对钝体驻定的Sydney湍流扩散火焰HM1进行数值模拟,以比较不同燃烧模型的性能,并比较标量联合的概率密度函数方法和Euler非稳态火焰面模型对氮氧化物排放预测的差异.计算结果和实验数据的比较表明,采用概率密度函数方法计算化学反应可以得到更好的结果但计算量较大,而用火焰面模型求解计算量较小,在接近完全燃烧的情形下,其计算结果比较合理.