Extended flamelet model and improved interaction of chemistry and turbulence for modeling partially premixed combustion with a joint PDF method

Turbulent combustion is commonly categorized into premixed, non-premixed and partially premixed combustion. For nonpremixed combustion simulations the laminar flamelet concept proved to be very valuable while for the more complex case of partially premixed combustion this model shows considerable deficiencies. Here, the classical laminar flamelet approach is extended to the partially premixed combustion regime. For that, the joint statistics of mixture fraction, scalar dissipation rate and a progress variable, calculated with a joint probability density function (PDF) method, is used to get the statistics of the compositions and of the chemical energy source term from pre-processed flame tables. This approach can be compared with the unsteady flamelet concept; the main differences consists of the way the progress variable evolution is computed and in the pre-computed flame tables. The progress variable describes the point of time a fluid parcel is consumed by a flame front. The fluid parcels are represented by computational particles, which are used for PDF methods. The pre-computed flame tables are computed from steady solutions 2D stabilized flames propagating into an unburnt mixture with varying mixture fraction. The corresponding position of a fluid particle in such a 2D laminar flame is determined by its mixture fraction and a burning time; both to be modeled for each computational particle in the PDF simulation. Numerical experiments of turbulent diffusion jet flames demonstrate that this approach can be employed for challenging test cases. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

An Approach to Model Partially Premixed Turbulent Combustion withProbability Density Function (PDF) Methods

In turbulent combustion one distinguishes between premixed, non-premixed and partially premixed combustion. While laminar flamelet models proved to be extremely valuable for a wide range of non-premixed flame simulations, similar approaches are more problematic in the partially premixed regime. Here the laminar flamelet concept for non-premixed turbulent combustion simulations is generalized for the partially premixed regime. Similar as in the unsteady flamelet approach, the joint statistics of a progress variable, mixture fraction and scalar dissipation rate is used to obtain the joint statistics of the compositions from pre-computed flame tables. The required distribution is computed with a joint PDF method and the main differences between the new approach and previous ones, are the pre-computed tables and the way the evolution of the progress variable is calculated. Instead of evolving 1D flamelets, steady 2D solutions of burning flamelets propagating into unburned mixtures with varying mixture fraction are considered. The location of a fluid particle in this 2D laminar flame is defined by its mixture fraction and a burning time, which are modeled for each computational particle used in the PDF method. Numerical experiments of a turbulent lifted diffusion flame and a premixed Bunsen flame demonstrate that this approach can be employed for a wide range of applications. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Lagrangian Joint PDF Approach for Turbulent Premixed Combustion

A model for premixed turbulent combustion based on a joint velocity probability density function (PDF) method and a progress variable is presented. Compared with other methods employing progress variables, the advantage here is that turbulent mixing of the progress variable requires no modeling. Moreover, by applying scale separation, the Lagrangian framework allows to account for the embedded, quasi laminar flame structure in a very natural way. The numerical results presented here are based on a simple closure of the progress variable source term and it is demonstrated that the new modeling approach is robust and shows the correct qualitative behavior. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Mixing modeling for stochastic simulations of turbulent premixed flames

The use of probability density function (PDF) methods for turbulent combustion simulations is very attractive because arbitrary finite-rate chemistry can be exactly taken into account. However, many real flames involve a variety of mixing regimes (non-premixed, partially-premixed and premixed turbulent combustion), and the development of PDF methods for partiallypremixed and premixed turbulent combustion turned out to be a very challenging task. A promising way to extend the range of applicability of PDF methods to the fast flamelet chemistry of turbulent premixed flames is described here. Simulation results of three turbulent premixed flames demonstrate the suitability of this approach. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Examination of errors of table integration in flamelet/progress variable modeling of a turbulent non-premixed jet flame

A flamelet model implementation consists of two steps: the generation of a set of laminar flamelet solutions and the integration of the laminar flamelet solutions with presumed shape probability density functions (PDFs) to produce a flamelet table for turbulent flame simulations. Many studies have been done in the past to examine the effect of different flamelet modeling strategies including the effect of employing different laminar flamelet solutions for the modeling. However, little work has been done to examine the effect of different presumed PDF table integration approaches on different flamelet model predictions. This work aims at investigating the source of errors arising from the flamelet table integration. The flamelet/progress variable model is chosen as a representative flamelet model, and three different presumed PDF table integration approaches are compared to examine the effect of table integration on flamelet model predictions. A laboratory-scale turbulent non-premixed jet flame (Sandia flame D) is chosen as a test case for the examination. In general, some evident sensitivity of the modeling results to the different flamelet table integration approaches is observed. The underlying reasons for the performance difference of different approaches are explored, and it is found that a model that preserves the one-dimensional laminar flamelet structure during the presumed-PDF table integration can improve the model prediction accuracy. Different sources of errors involved in flamelet model implementation are investigated, including numerical integration errors, flamelet table errors, and the errors in the predictions of the flamelet independent variables.     

Modelling Intermediate Species in Partially Premixed Turbulent Combustion Based on the LES-FGM Method北大核心CSCD

The LES of partially premixed turbulent flame MRB in TU Darmstadt was conducted based on the flamelet-tabulated combustion model FGM, and effects of premixed and partially premixed tabulations on the modelling results were studied. The results show that, different methods of tabulation exhibit limited influences on the predictions of the flame structure, velocity, and major species, but using a partially premixed tabulation largely improves the reliability of modelling intermediate minor species CO and H2. The underlying reason lies in a better inclusion of the fuel-air mixing effects through the partially premixed tabulation, which is built based on laminar counter-flow flames. Adding extra transport equations for the intermediate species improves the predictions of intermediate species, especially given a premixed tabulation adopted; meanwhile, the stretch effects in this turbulent flame are ignorable. The results are significant to guide the high-fidelity simulation of partially premixed turbulent flames based on the flamelet-tabulated combustion model. © 2023 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

RANS-simulation of premixed turbulent combustion using the level set approach

A model for premixed turbulent combustion is investigated using a RANS-approach. The evolution of the flame front is described with the help of the level set approach [1] which is used for tracking of propagating interfaces in free-surface flows, geodesics, grid generation and combustion. The fluid properties are conditioned on the flame front position using a burntunburnt probability function across the flame front. Computations are performed using the code FASTEST-3D which is a flow solver for a non-orthogonal, block-structured grid. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Re.ned Modeling of the pressure redistribution/scrambling under consistent consideration of scalar turbulence effects in turbulent combustion

In practical turbulent flow problems of engineering importance the coupling between velocity and scalar turbulence along with the variable density plays a non negligible role. For computations using second moment closure approach, the pressure redistribution/scrambling is the most critical term to be modeled as well known. Almost all existing models consist in rescating models derived on a constant density basis in a density weighted form. With regard to turbulent premixed combustion in fact, the application of such models to a range of transient one‐dimensional and two‐dimensional premixed flames in the flamelet regime has been found to yield unsatisfactory results, see [1]. As pointed out by Sadiki [2], the use of the Favre method must be consistently considered as far as open thermodynamic systems are concerned. Furthermore, the need for maintaining certain invariance properties, physical and mathematical realizability conditions in formulating turbulence models is well accepted. Because turbulent processes are irreversible, these efforts demand a carefull consideration of thermodynamic concepts. Based on the results in [1] and following [2], this work aims to derive a physically consistent formulation of the pressure redistribution/scrambling term under consideration of the variable density. Considering the case of premixed flames, the thermochemistry is included by means of a single reactive scalar ‐ the reaction progress variable. The accuracy of the model extensions proposed is demonstrated by comparing the numerical results with experimental data in opposed jet premixed flame configuration.     

A self-similar premixed turbulent flame model

This paper is devoted to premixed combustion modelling in turbulent flow. First, we derive a model for the turbulent flame velocity based on the observed self-similarity of the turbulent flame. The model uses the local flame brush width as a fundamental parameter and, therefore, we show how it can be retrieved for numerical implementation. The diffusive property of the brush width is treated in such a way as to theoretically let the brush have a clearly defined boundary propagating at finite velocity. The model, implemented in Star-CD CFD software through user programming, is then numerically tested on three configurations for which another model, the Turbulent Flame Closure model, is known to give very good agreement. Some effects of numerics are commented and results for both models are compared. While based on very different approaches the two models lead to substantially similar results. In this way, we have shown that the local brush width can effectively be used, giving an additional degree of freedom for premixed turbulent combustion modelling.     

Evaluation of the approximated diffusion flamelet concept using fuels with different chemical complexity

The ability of flamelet models to reproduce turbulent combustion in devices such as diesel engines or gas turbines has enhanced the usage of these approaches in Computational Fluid Dynamics (CFD) simulations. The models based on turbulent look-up tables generated from counterflow laminar diffusion flames (DF model) permit drastic reduction of the computational cost of the CFD calculation. Nevertheless, for complex molecular fuels, such as n-heptane, the oxidation process involves hundreds of species and the calculation of the transport equations together with the ODE system that models the chemical kinetics for the DF solution becomes unaffordable for industrial devices where hundreds of flamelets are required. In this context, new hypotheses have to be introduced in order to reduce the computational cost maintaining the coherence of the combustion process. Recently, a new model known as Approximated Diffusion Flamelet (ADF) has been proposed with the aim of solving the turbulent combustion for complex fuels in a reduced time. However, the validity of this model is still an open question and has to be verified in order to justify subsequent CFD calculations. This work assesses the ADF model and its ability to reproduce accurately the combustion process and its main parameters for three fuels with different chemical complexity and boundary conditions by its comparison with the DF model. Results show that although some discrepancies arise, the ADF model has the ability to correctly describe the ignition delay and the combustion structure in the auto-ignition zone that is the most relevant one for industrial processes.     

Large eddy simulation of non-premixed pulverized coal combustion in corner-fired furnace for various excess air ratios

Large-eddy simulations (LES) were carried out to study the effects of burning atmosphere on the coal combustion process in a corner-fired furnace. The LES for the turbulent gas was coupled with the discrete phase model (DPM) for coal particles trajectories and the non-premixed mixture fraction probability density function (MF-PDF) combustion model for pulverized coal combustion. The coal combustion processes, including the flame characteristics, burning coal behaviors and NO_x pollutant emissions, for different burning atmospheres are analyzed qualitatively and quantitatively. The heat and momentum transfer between burning coal and turbulent gas are greatly enhanced by the corner-fired flow. With a given particle size, the char particles present a similar distribution in the whole chamber. For a fuel-rich atmosphere, the concentration is obviously much higher and exhibits much higher spatial variability than the other two conditions. The coal combustion efficiency decreases in oxygen-rich and fuel-rich burning atmospheres, but the flame stability is more affected at the fuel-rich atmosphere by the lack of oxygen. NO pollutant is obviously reduced at the fuel-rich atmospheres, and the NO pollutant emissions are more affected by the reducing atmosphere than the low temperature. These findings may provide insight into strategies to design and monitor tangentially-fired pulverized coal boilers.     

Large eddy simulation of an ethylene-air turbulent premixed V-flame

Large eddy simulation (LES) using a dynamic eddy viscosity subgrid scale stress model and a fast-chemistry combustion model without accounting for the finite-rate chemical kinetics is applied to study the ignition and propagation of a turbulent premixed V-flame. A progress variable c-equation is applied to describe the flame front propagation. The equations are solved two dimensionally by a projection-based fractional step method for low Mach number flows. The flow field with a stabilizing rod without reaction is first obtained as the initial field and ignition happens just upstream of the stabilizing rod. The shape of the flame is affected by the velocity field, and following the flame propagation, the vortices fade and move to locations along the flame front. The LES computed time-averaged velocity agrees well with data obtained from experiments.     

层焰系统奇点结构的研究

在一定的假设条件下,将气体燃烧的物理模型简化为层焰系统．借助于热力学理论和相关的守恒定律建立了层焰系统的数学模型并对其进行了深入地分析．利用常微分方程的定性理论和方法,结合燃烧的实际需要,对应于参数的不同取值,研究了位于瑞利线上奇点的个数和位置,确定了暴燃区和爆炸区内不同奇点的定性结构及其稳定性,给出了在反应速度-滞止焓平面及燃烧速度-滞止焓平面内层焰系统轨线的相图．     

Progress in the Self-Similar Turbulent Flame premixed combustion model

This paper is devoted to premixed combustion modeling in turbulent flow. First, we briefly remind the main features of the Self-Similar Turbulent Flame model that was more extensively developed in a former paper. Then, we carefully describe some improvements of the model. The determination of the turbulent flame velocity is based on the observed self-similarity of the turbulent flame and uses the local flame brush width as a fundamental parameter, which must be retrieved. With respect to the former version, we now derive more rigorously how the density variation has to be taken into account in the width retrieving function. We reformulate the diffusion term as a classical flux divergence term. We enforce the compatibility of the model for the limit of weak turbulence. We include a contracting effect of the source term, thus allowing to give a stationary mono-dimensional asymptotic solution with a finite width. We also include in a preliminary form, a stretch factor, which proves to be useful for controlling the flame behavior close to the flame holder and near the walls. The model implementation in the Star-CD CFD code is then tested on three different flame configurations. Finally, we shortly discuss the model improvements and the simulation results.     

A mixing model providing statistics of the conditional scalar dissipation rate

Joint composition probability density function (PDF) methods are used for the numerical simulation of turbulent reactive flows. Here, other than in classical Reynolds averaged Navier–Stokes (RANS) or large eddy simulation (LES) approaches, the highly non-linear chemical source term appears in closed form. On the other hand, mixing models are required for the closure of the molecular diffusion term. In the present work, the joint statistics of the scalar and the scalar dissipation rate provided by the parameterized scalar profile (PSP) mixing model are validated. The goal is to combine the PDF method with a flamelet approach, where the scalar dissipation rate plays a crucial role in determining the contribution of the chemical source term. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Application of a non-asymptotic approach to prediction of the propagation of a flame through a fuel and/or oxidant droplet cloud

The mathematical analysis of laminar premixed spray propagation has generally been based on exploiting the inverse of the large chemical activation energy as an appropriate parameter for asymptotic analysis. In the current work we apply a modification of a recently suggested non-asymptotic approach for gaseous flames which makes use of a different approximation. In it, only the Arrhenius exponential term in the reaction rate expression is approximated using a step function chosen so that the two functions are in proximity in an integral sense. Application of this approach is more amenable and is shown to yield a simple formula for the burning velocity of a flame propagating through a cloud of fuel and/or oxidant droplets, for the fuel rich off-stoichiometric case in which the only reactant present in the chemical reaction term is the deficient oxidant which appears linearly. Results computed with the new analytical solutions are presented and a comparison is made with the predictions using the usual large activation energy approach. In addition, a double spray is considered for the first time in which both liquid oxidant and liquid fuel feature as sprays of droplets in the unburned pre-mixture. Such a situation arises in rocket engines in which two initially separate spray streams mix in a turbulent shear flow so that locally one dimensionally propagating double spray premixed flames are created. The analysis leads to an analytical expression for the laminar burning velocity dependent on the spray- and gas-related parameters. Typical thermal and velocity maps in parametric space are presented.     

One-dimensional, two-fluid modelling of turbulent premixed flames

The “two-fluid” mathematical model for turbulent combustion is applied to a one-dimensional, premixed, stabilized ducted flame. The flame is assumed to consist of two interspersed fluids (“reactants” and “products”), each characterized by its own properties and interacting through the exchange of mass, heat, and momentum. The distributions of pressure, densities, velocities, and volume fractions across the duct were successfully simulated. From a parametric study on the effects of the empirical constants involved in the interfluid relations, the significant dependence of the system on the parameters that characterize the mass transfer rate and the relative effect of mass transfer to momentum transfer was confirmed. The application of the model to transient states proved its ability to predict system oscillations.     

Parameterized Scalar Profile Mixing Model for Turbulent Combustion Simulations

The composition fields in turbulent reacting flows are affected by turbulent transport (macromixing), molecular diffusion (micromixing), and chemical reactions. In the joint velocity-composition probability density function transport equation the highly non-linear macromixing and chemical reaction terms appear in closed form. This is a considerable advantage over moment closure methods. Micromixing on the other hand requires modeling and especially for turbulent combustion accurate mixing models are crucial. Our approach to model the mixing of scalars, e.g. species mass fractions or temperature, is based on considering one-dimensional parameterized scalar profiles (PSP). Here, an extension of the PSP mixing model to inhomogeneous flows is presented. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Wall effect on determination of laminar burning velocity in a constant volume bomb using a quasi-dimensional model

In experiment, two optical and pressure-based methods are frequently used to evaluate laminar burning velocity of a combustible mixture. In the currently reported work, the pressure-based method was utilized to find the laminar burning velocity using the measurement of pressure variations during the combustion process in a spherical bomb and analyzing them through a multi-zone quasi-dimensional model. To check the results of the method, isooctane–air mixtures were used at equivalence ratios of 0.85 and 1.0 and initial pressures of 95 and 150 kPa with 343 K initial temperature. The time history of the bomb pressure during the combustion event, initial pressure and temperature, fuel type, and equivalence ratio were applied as input to a Fortran program written by the author based on the multi-zone combustion model; and, flame radius-time, flame speed, and laminar burning velocity at different pressures and temperatures were evaluated assuming spherical flame growth. The obtained results were compared with those of some other researchers and a reasonable agreement was observed. The wall effect on the laminar burning velocity at the end of the combustion process was clearly highlighted and a reliable range of burning velocity was distinguished. The results showed that the evaluated laminar burning velocity was not reliable at the late part of the combustion process due to possible local contact of flame front and the bomb wall, the wall effect on the reacting species, flow to small crevices, and the boundary layer effect.     

Mixing Model Based on Parameterized Scalar Profiles

The composition fields in turbulent reacting flows are affected by turbulent transport (macromixing), molecular diffusion (micromixing), and chemical reactions. In the joint velocity-composition probability density function transport equation the highly nonlinear macromixing and chemical reaction terms appear in closed form. This is a considerable advantage over second moment closure methods. Micromixing on the other hand requires modelling and especially for turbulent combustion accurate mixing models are crucial. In this paper we present an approach to model the mixing of scalars, e.g. species mass fractions or temperature, based on considering one-dimensional parameterized scalar profiles (PSP). (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)