Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Flamelet model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Flamelet Model Application for Non-Premixed Turbulent Combustion

The current Final Report contains results of the study which was performed in Scientific Research Center 'ECOLEN' (Moscow, Russia). The study concerns the development and verification of non-expensive approach for modeling of supersonic turbulent diffusion flames based on flamelet consideration of the chemistry/turbulence interaction (FL approach). Research work included: development of the approach and CFD tests of the flamelet model for supersonic jet flames; development of the simplified procedure for solution of the flamelet equations based on partial equilibrium chemistry assumption; study of the flame ignition/extinction predictions provided by flamelet model. The performed investigation demonstrated that FL approach allowed to describe satisfactory main features of supersonic H 2/air jet flames. Model demonstrated also high capabilities for reduction of the computational expenses in CFD modeling of the supersonic flames taking into account detailed oxidation chemistry. However, some disadvantages and restrictions of the existing version of approach were found in this study. They were: (1) inaccuracy in predictions of the passive scalar statistics by our turbulence model for one of the considered test cases; and (2) applicability of the available version of the flamelet model to flames without large ignition delay distance only. Based on the results of the performed investigation, we formulated and submitted to the National Aeronautics and Space Administration our Project Proposal for the next step research directed toward further improvement of the FL approach.

Secundov, A.↗

Laminar flamelet modeling of turbulent diffusion flames

In modeling turbulent combustion, decoupling the chemistry from the turbulence is of great practical significance. In cases in which the equilibrium chemistry model breaks down, laminar flamelet modeling (LFM) is a promising approach to decoupling. Here, the validity of this approach is investigated using direct numerical simulation of a simple chemical reaction in two-dimensional turbulence.

Mell, W. E.↗

Study of Turbulent Premixed Flame Propagation using a Laminar Flamelet Model

The laminar flamelet concept in turbulent reacting flows is considered applicable to many practical combustion systems (Linan & Williams 1993). For turbulent premixed combustion, the laminar flamelet regime is valid when turbulent Karlovitz number is less than unity, which is equivalent to stating that the characteristic thickness of the flame is less than that of a Kolmogorov eddy; this is known as the Klimov-Williams criterion (Williams 1985). In such a case, the flame maintains its laminar structure, and the effect of turbulent flow is merely to wrinkle and strain the flame front. The propagating wrinkled premixed flame can then be described as an infinitesimally thin surface dividing the unburnt fresh mixture and the burnt product.

Im, H. G.↗

A Priori Analysis of a Compressible Flamelet Model using RANS Data for a Dual-Mode Scramjet Combustor

In an effort to make large eddy simulation of hydrocarbon-fueled scramjet combustors more computationally accessible using realistic chemical reaction mechanisms, a compressible flamelet/progress variable (FPV) model was proposed that extends current FPV model formulations to high-speed, compressible flows. Development of this model relied on observations garnered from an a priori analysis of the Reynolds-Averaged Navier-Stokes (RANS) data obtained for the Hypersonic International Flight Research and Experimentation (HI-FiRE) dual-mode scramjet combustor. The RANS data were obtained using a reduced chemical mechanism for the combustion of a JP-7 surrogate and were validated using avail- able experimental data. These RANS data were then post-processed to obtain, in an a priori fashion, the scalar fields corresponding to an FPV-based modeling approach. In the current work, in addition to the proposed compressible flamelet model, a standard incompressible FPV model was also considered. Several candidate progress variables were investigated for their ability to recover static temperature and major and minor product species. The effects of pressure and temperature on the tabulated progress variable source term were characterized, and model coupling terms embedded in the Reynolds- averaged Navier-Stokes equations were studied. Finally, results for the novel compressible flamelet/progress variable model were presented to demonstrate the improvement attained by modeling the effects of pressure and flamelet boundary conditions on the combustion.

Quinlan, Jesse R.↗

Statistics for laminar flamelet modeling

Statistical information required to support modeling of turbulent premixed combustion by laminar flamelet methods is extracted from a database of the results of Direct Numerical Simulation of turbulent flames. The simulations were carried out previously by Rutland (1989) using a pseudo-spectral code on a three dimensional mesh of 128 points in each direction. One-step Arrhenius chemistry was employed together with small heat release. A framework for the interpretation of the data is provided by the Bray-Moss-Libby model for the mean turbulent reaction rate. Probability density functions are obtained over surfaces of the constant reaction progress variable for the tangential strain rate and the principal curvature. New insights are gained which will greatly aid the development of modeling approaches.

Cant, R. S.↗

Unsteady Three-Dimensional Simulation of a Shear Coaxial GO2/GH2 Rocket Injector with RANS and Hybrid-RAN-LES/DES Using Flamelet Models

Historically, the analysis and design of liquid rocket engines (LREs) has relied on full‐scale testing and one‐dimensional empirical tools. The testing is extremely expensive and the one‐dimensional tools are not designed to capture the highly complex, and multi‐dimensional features that are inherent to LREs. Recent advances in computational fluid dynamics (CFD) tools have made it possible to predict liquid rocket engine performance, stability, to assess the effect of complex flow features, and to evaluate injector‐driven thermal environments, to mitigate the cost of testing. Extensive efforts to verify and validate these CFD tools have been conducted, to provide confidence for using them during the design cycle. Previous validation efforts have documented comparisons of predicted heat flux thermal environments with test data for a single element gaseous oxygen (GO2) and gaseous hydrogen (GH2) injector. The most notable validation effort was a comprehensive validation effort conducted by Tucker et al. [1], in which a number of different groups modeled a GO2/GH2 single element configuration by Pal et al [2]. The tools used for this validation comparison employed a range of algorithms, from both steady and unsteady Reynolds Averaged Navier‐Stokes (U/RANS) calculations, large‐eddy simulations (LES), detached eddy simulations (DES), and various combinations. A more recent effort by Thakur et al. [3] focused on using a state‐of‐the‐art CFD simulation tool, Loci/STREAM, on a two‐dimensional grid. Loci/STREAM was chosen because it has a unique, very efficient flamelet parameterization of combustion reactions that are too computationally expensive to simulate with conventional finite‐rate chemistry calculations. The current effort focuses on further advancement of validation efforts, again using the Loci/STREAM tool with the flamelet parameterization, but this time with a three‐dimensional grid. Comparisons to the Pal et al. heat flux data will be made for both RANS and Hybrid RANSLES/ Detached Eddy simulations (DES). Computation costs will be reported, along with comparison of accuracy and cost to much less expensive two‐dimensional RANS simulations of the same geometry.

Westra, Doug G.↗

Modeling of turbulent chemical reaction

Viewgraphs are presented on modeling turbulent reacting flows, regimes of turbulent combustion, regimes of premixed and regimes of non-premixed turbulent combustion, chemical closure models, flamelet model, conditional moment closure (CMC), NO(x) emissions from turbulent H2 jet flames, probability density function (PDF), departures from chemical equilibrium, mixing models for PDF methods, comparison of predicted and measured H2O mass fractions in turbulent nonpremixed jet flames, experimental evidence of preferential diffusion in turbulent jet flames, and computation of turbulent reacting flows.

Chen, J.-Y.↗

A Priori Analysis of Flamelet-Based Modeling for a Dual-Mode Scramjet Combustor

An a priori investigation of the applicability of flamelet-based combustion models to dual-mode scramjet combustion was performed utilizing Reynolds-averaged simulations (RAS). For this purpose, the HIFiRE Direct Connect Rig (HDCR) flowpath, fueled with a JP-7 fuel surrogate and operating in dual- and scram-mode was considered. The chemistry of the JP-7 fuel surrogate was modeled using a 22 species, 18-step chemical reaction mechanism. Simulation results were compared to experimentally-obtained, time-averaged, wall pressure measurements to validate the RAS solutions. The analysis of the dual-mode operation of this flowpath showed regions of predominately non-premixed, high-Damkohler number, combustion. Regions of premixed combustion were also present but associated with only a small fraction of the total heat-release in the flow. This is in contrast to the scram-mode operation, where a comparable amount of heat is released from non-premixed and premixed combustion modes. Representative flamelet boundary conditions were estimated by analyzing probability density functions for temperature and pressure for pure fuel and oxidizer conditions. The results of the present study reveal the potential for a flamelet model to accurately model the combustion processes in the HDCR and likely other high-speed flowpaths of engineering interest.

Quinlan, Jesse R.↗

The evolution equation for the flame surface density in turbulent premixed combustion

The mean reaction rate in flamelet models for turbulent premixed combustion depends on two basic quantities: a mean chemical rate, called the flamelet speed, and the flame surface density. Our previous work had been primarily focused on the problem of the structure and topology of turbulent premixed flames, and it was then determined that the flamelet speed, when space-averaged, is only weakly sensitive to the turbulent flow field. Consequently, the flame surface density is the key quantity that conveys most of the effects of the turbulence on the rate of energy release. In flamelet models, this quantity is obtained via a modeled transport equation called the Sigma-equation. Past theoretical work has produced a rigorous approach that leads to an exact but unclosed formulation for the turbulent Sigma-equation. In the exact Sigma-equation, it appears that the dynamical properties of the flame surface density are determined by a single parameter, namely the turbulent flame stretch. Unfortunately, the turbulent flame stretch as well as the flame surface density is not available from experiments, and, in the absence of experimental data, little is known on the validity of the closure assumptions used in current flamelet models. Direct Numerical Simulation (DNS) is the alternative approach to get basic information on these fundamental quantities. In the present work, three-dimensional DNS of premixed flames in isotropic turbulent flow is used to estimate the different terms appearing in the Sigma-equation. A new methodology is proposed to provide the source and sink terms for the flame surface density, resolved both temporally and spatially throughout the turbulent flame brush. Using this methodology, our objective is to extract the turbulent flame stretch from the DNS data base and then perform extensive comparisons with flamelet models. Thanks to the detailed information produced by the DNS-based analysis, it is expected that this type of comparison will not only underscore the shortcomings of current models, but also suggest ways to improve them.

Trouve, Arnaud↗

The evolution equation for the flame surface density in turbulent premixed combustion

One central ingredient in flamelet models for turbulent premixed combustion is the flame surface density. This quantity conveys most of the effects of the turbulence on the rate of energy release and is obtained via a modeled transport equation, called the Sigma-equation. Past theoretical work has produced a rigorous approach that leads to an exact, but unclosed, formulation for the turbulent Sigma-equation. In this exact Sigma-equation, it appears that the dynamical properties of the flame surface density are determined by a single parameter, namely the turbulent flame stretch. Unfortunately, the flame surface density and the turbulent flame stretch are not available from experiments and, in the absence of experimental data, little is known on the validity of the closure assumptions used in current flamelet models. Direct Numerical Simulation (DNS) is the obvious, complementary approach to get basic information on these fundamental quantities. Three-dimensional DNS of premixed flames in isotropic turbulent flow is used to estimate the different terms appearing in the Sigma-equation. A new methodology is proposed to provide the source and sink terms for the flame surface density, resolved both temporally and spatially throughout the turbulent flame brush. Using this methodology, the effects of the Lewis number on the rate of production of flame surface area are described in great detail and meaningful comparisons with flamelet models can be performed. The analysis reveals in particular the tendency of the models to overpredict flame surface dissipation as well as their inability to reproduce variations due to thermo-diffusive phenomena. Thanks to the detailed information produced by a DNS-based analysis, this type of comparison not only underscores the shortcomings of current models but also suggests ways to improve them.

Trouve, A.↗

Progress Toward Affordable High Fidelity Combustion Simulations Using Filtered Density Functions for Hypersonic Flows in Complex Geometries

Significant progress has been made in the development of subgrid scale (SGS) closures based on a filtered density function (FDF) for large eddy simulations (LES) of turbulent reacting flows. The FDF is the counterpart of the probability density function (PDF) method, which has proven effective in Reynolds averaged simulations (RAS). However, while systematic progress is being made advancing the FDF models for relatively simple flows and lab-scale flames, the application of these methods in complex geometries and high speed, wall-bounded flows with shocks remains a challenge. The key difficulties are the significant computational cost associated with solving the FDF transport equation and numerically stiff finite rate chemistry. For LES/FDF methods to make a more significant impact in practical applications a pragmatic approach must be taken that significantly reduces the computational cost while maintaining high modeling fidelity. An example of one such ongoing effort is at the NASA Langley Research Center, where the first generation FDF models, namely the scalar filtered mass density function (SFMDF) are being implemented into VULCAN, a production-quality RAS and LES solver widely used for design of high speed propulsion flowpaths. This effort leverages internal and external collaborations to reduce the overall computational cost of high fidelity simulations in VULCAN by: implementing high order methods that allow reduction in the total number of computational cells without loss in accuracy; implementing first generation of high fidelity scalar PDF/FDF models applicable to high-speed compressible flows; coupling RAS/PDF and LES/FDF into a hybrid framework to efficiently and accurately model the effects of combustion in the vicinity of the walls; developing efficient Lagrangian particle tracking algorithms to support robust solutions of the FDF equations for high speed flows; and utilizing finite rate chemistry parametrization, such as flamelet models, to reduce the number of transported reactive species and remove numerical stiffness. This paper briefly introduces the SFMDF model (highlighting key benefits and challenges), and discusses particle tracking for flows with shocks, the hybrid coupled RAS/PDF and LES/FDF model, flamelet generated manifolds (FGM) model, and the Irregularly Portioned Lagrangian Monte Carlo Finite Difference (IPLMCFD) methodology for scalable simulation of high-speed reacting compressible flows.

Drozda, Tomasz G.↗

Full numerical simulation of coflowing, axisymmetric jet diffusion flames

The near field of a non-premixed flame in a low speed, coflowing axisymmetric jet is investigated numerically using full simulation. The time-dependent governing equations are solved by a second-order, explicit finite difference scheme and a single-step, finite rate model is used to represent the chemistry. Steady laminar flame results show the correct dependence of flame height on Peclet number and reaction zone thickness on Damkoehler number. Forced simulations reveal a large difference in the instantaneous structure of scalar dissipation fields between nonbuoyant and buoyant cases. In the former, the scalar dissipation marks intense reaction zones, supporting the flamelet concept; however, results suggest that flamelet modeling assumptions need to be reexamined. In the latter, this correspondence breaks down, suggesting that modifications to the flamelet modeling approach are needed in buoyant turbulent diffusion flames.

Mahalingam, S.↗

The compositional structure and the effects of exothermicity in an unpremixed planar jet flame

Results are presented of direct numerical simulations (DNS) of a randomly perturbed compressible, spatially developing planar jet under the influence of a finite rate Arrhenius chemical reaction of the type F + O yields Product + Heat. The objectives of the simulations are to assess the compositional structure of the flame and to determine the influence of reaction exothermicity by means of statistical sampling of the data generated by DNS. It is shown that even with this idealized kinetics model, the simulated results exhibit features in accordance with experimental data. These results indicate that the Damkoehler number is an important parameter in determining the statistical composition of the reacting field, and that the results are insensitive to the mechanism by which this parameter is varied. It is demonstrated that as the intensity of mixing is increased and the effect of finite rate chemistry is more pronounced, the magnitudes of the ensemble mean and rms of the product mass fraction decrease, and those of the reactants mass fraction increase. Also, at higher mixing rates the joint probability density functions of the reactants' mass fractions shift towards higher values within the composition domain indicating a lower reactedness. These trends are consistent with those observed experimentally and are very useful in portraying the statistical structure of non-equilibrium diffusion flames. The DNS generated data are also utilized to examine the applicability of the 'laminar diffusion flamelet model' in predicting the rate of the reactant conversion with finite rate chemistry. This examination indicates that the bounds of the product formation scale reasonably well with those obtained by the flamelet model. However, in the range of the relatively low values of the Damkoehler numbers considered, the scatter of the results is substantially more than that to be completely determined by the closure.

Steinberger, C. J.↗

Scalar Dissipation Modeling for Passive and Active Scalars: a priori Study Using Direct Numerical Simulation

Transitional databases from Direct Numerical Simulation (DNS) of three-dimensional mixing layers for single-phase flows and two-phase flows with evaporation are analyzed and used to examine the typical hypothesis that the scalar dissipation Probability Distribution Function (PDF) may be modeled as a Gaussian. The databases encompass a single-component fuel and four multicomponent fuels, two initial Reynolds numbers (Re), two mass loadings for two-phase flows and two free-stream gas temperatures. Using the DNS calculated moments of the scalar-dissipation PDF, it is shown, consistent with existing experimental information on single-phase flows, that the Gaussian is a modest approximation of the DNS-extracted PDF, particularly poor in the range of the high scalar-dissipation values, which are significant for turbulent reaction rate modeling in non-premixed flows using flamelet models. With the same DNS calculated moments of the scalar-dissipation PDF and making a change of variables, a model of this PDF is proposed in the form of the (beta)-PDF which is shown to approximate much better the DNS-extracted PDF, particularly in the regime of the high scalar-dissipation values. Several types of statistical measures are calculated over the ensemble of the fourteen databases. For each statistical measure, the proposed (beta)-PDF model is shown to be much superior to the Gaussian in approximating the DNS-extracted PDF. Additionally, the agreement between the DNS-extracted PDF and the (beta)-PDF even improves when the comparison is performed for higher initial Re layers, whereas the comparison with the Gaussian is independent of the initial Re values. For two-phase flows, the comparison between the DNS-extracted PDF and the (beta)-PDF also improves with increasing free-stream gas temperature and mass loading. The higher fidelity approximation of the DNS-extracted PDF by the (beta)-PDF with increasing Re, gas temperature and mass loading bodes well for turbulent reaction rate modeling.

scalar dissipation modeling↗

Premixed Edge-Flames in Spatially-Varying Straining Flows

Flames subject to temporally and spatially uniform hydrodynamic strain are frequently used to model the local interactions of flame fronts with turbulent flow fields (Williams, 1985; Peters, 1986; Bradley, 1992). The applicability of laminar flamelet models in strongly turbulent flows have been questioned recently (Shay and Ronney, 1998) because in turbulent flows the strain rate (sigma) changes at rates comparable to sigma itself and the scale over which the flame front curvature and sigma changes is comparable to the curvature scale itself. Therefore quasi-static, local models of turbulent strain and curvature effects on laminar flamelets may not be accurate under conditions where the strain and curvature effects are most significant. The purpose of this study is to examine flames in spatially-varying strain and compare their properties to those of uniformly strained flames.

Liu, Jian-Bang↗

DNS and modeling of the interaction between turbulent premixed flames and walls

The interaction between turbulent premixed flames and walls is studied using a two-dimensional full Navier-Stokes solver with simple chemistry. The effects of wall distance on the local and global flame structure are investigated. Quenching distances and maximum wall heat fluxes during quenching are computed in laminar cases and are found to be comparable to experimental and analytical results. For turbulent cases, it is shown that quenching distances and maximum heat fluxes remain of the same order as for laminar flames. Based on simulation results, a 'law-of-the-wall' model is derived to describe the interaction between a turbulent premixed flame and a wall. This model is constructed to provide reasonable behavior of flame surface density near a wall under the assumption that flame-wall interaction takes place at scales smaller than the computational mesh. It can be implemented in conjunction with any of several recent flamelet models based on a modeled surface density equation, with no additional constraints on mesh size or time step.

Poinsot, T. J.↗

Direct simulation of turbulent combustion

Understanding and modeling of turbulent combustion are key-problems in the computation of numerous practical systems. Because of the lack of analytical theories in this field and of the difficulty of performing precise experiments, direct simulation appears to be one of the most attractive tools to use in addressing this problem. The present work can be split into two parts: (1) Development and validation of a direct simulation method for turbulent combustion; (2) Applications of the method to premixed turbulent combustion problems. The goal of part 1 is to define and to test a numerical method for direct simulation of reacting flows. A high level of confidence should be attached to direct simulation results, and this can only be achieved through extensive validation tests. In part 2, direct simulation is used to address some of the many critical problems related to turbulent combustion. At the present time, I have limited this work to premixed combustion and considered only four basic issues: (1) The effect of pressure waves on flame propagation; (2) The interaction between flame fronts and vortices; (3) The influence of curvature on premixed flame fronts; and (4) The validation of flamelet models for premixed turbulent combustion.

Poinsot, T. J.↗