Turbulence models for gas turbine combustors
This viewgraph presentation covers gas turbine combustor flow physics, turbulence model investigations, turbulent combustion modeling, and present status and future needs.
Brankovic, Andreja↗
SEARCH · Search NASA
Results for “turbulent combustion”
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.
Recent progress in the joint velocity-scalar PDF method
This viewgraph presentation discusses joint velocity-scalar PDF method; turbulent combustion modeling issues for gas turbine combustors; PDF calculations for a recirculating flow; stochastic dissipation model; joint PDF calculations for swirling flows; spray calculations; reduced kinetics/manifold methods; parallel processing; and joint PDF focus areas.
Anand, M. S.↗
Heat release effects on the instability of parallel shear layers
The influence of time-dependent heat addition on the linear instablity of shear layers is of considerable interest in understanding the dynamic behavior of reacting flows and combustion-turbulence interactions. The approach is based upon the Bernoulli enthalpy aeroacoustics theory, which utilizes the specific enthalpy and specific entropy as the primary thermodynamic variables. In addition, velocity oscillations are split into Helmoholtz decomposition theorem.
Hegde, U.↗
Characteristics of transitional and turbulent jet diffusion flames in microgravity
This paper presents the ground-based results obtained to date in preparation of a proposed space experiment to study the role of large-scale structures in microgravity transitional and turbulent gas-jet diffusion flames by investigating the dynamics of vortex/flame interactions and their influence on flame characteristics. The overall objective is to gain an understanding of the fundamental characteristics of transitional and turbulent gas-jet diffusion flames. Understanding of the role of large-scale structures on the characteristics of microgravity transitional and turbulent flames will ultimately lead to improved understanding of normal-gravity turbulent combustion.
Bahadori, Yousef M.↗
Computational fluid dynamic analysis of hybrid rocket combustor flowfields
Computational fluid dynamic analyses of the Navier-Stokes equations coupled with solid-phase pyrolysis, gas-phase combustion, turbulence and radiation are performed to study hybrid rocket combustor flowfields. The computational study is closely co-ordinated with a companion experimental program using a planar slab burner configuration with HTPB as fuel and gaseous oxygen. Computational predictions agree reasonably well with measurement data of fuel regression rates and surface temperatures. Additionally, most of the parametric trends predicted by the model are in general agreement with experimental trends. The computational model is applied to extend the results from the lab-scale to a full-scale axisymmetric configuration. The numerical predictions indicate that the full-scale configuration burns at a slower rate than the lab-scale combustor under identical specific flow rate conditions. The results demonstrate that detailed CFD analyses can play a useful role in the design of hybrid combustors.
Venkateswaran, S.↗
LES, DNS, and RANS for the Analysis of High-Speed Turbulent Reacting Flows
A filtered density function (FDF) method suitable for chemically reactive flows is developed in the context of large eddy simulation. The advantage of the FDF methodology is its inherent ability to resolve subgrid scales (SGS) scalar correlations that otherwise have to be modeled. Because of the lack of robust models to accurately predict these correlations in turbulent reactive flows, simulations involving turbulent combustion are often met with a degree of skepticism. The FDF methodology avoids the closure problem associated with these terms and treats the reaction in an exact manner. The scalar FDF approach is particularly attractive since it can be coupled with existing hydrodynamic computational fluid dynamics (CFD) codes.
Colucci, P. J.↗
Dynamics and Structure of Weakly-Strained Flames In Normal- and Micro-Gravity
Strained laminar flames have been systematically studied, as the understanding of their structure and dynamic behavior is of relevance to turbulent combustion. Most of these studies have been conducted in opposed-jet, stagnation-type flow configurations. Flame studies in stagnation flows also allow for the determination of fundamental flame properties such as laminar flame speeds and extinction states that can be also conveniently modeled in detail. Studies at high strain rates are important in quantifying and understanding the response of vigorously-burning flames under conditions in which the transport time scales become comparable to the chemical time scales. Studies of weakly-strained flames can be of particular interest for all stoichiometries. For example, the laminar flame speeds for any equivalence ratio, phi, can be accurately determined by using the counterflow technique only if measurements are obtained at very low strain rates. Furthermore, near-limit flames can be only stabilized by weak strain rates. Previous studies have shown that weakly-burning flames are particularly sensitive to chain mechanisms, thermal radiation, and unsteadiness. The stabilization and study of weakly-strained flames is complicated by the presence of buoyancy that can render the flames unstable to the point of extinction. Such instabilities are caused either by the induced natural convection or the fact that higher density fluid can find itself on top of a lower density fluid. Thus, the use of microgravity (mu-g) becomes essential in order to provide meaningful insight into this important combustion regime. In view of the foregoing considerations, the main objectives of the program are to: (1) Experimentally determine the laminar flame speed at near-zero strain rates; (2) Experimentally determine the extinction limits of near-limit flames; (3) Experimentally determine the response of near-limit flames to unsteadiness and heat loss; (4) Introduce Digital Particle Image Velocimetry (DPIV) in microgravity; (5) Conduct detailed numerical simulations of the experiments; and (6) Provide physical insight into the underlying physico-chemical mechanisms.
Zhang, Hai↗
Vortex/Flame Interactions in Microgravity Pulsed Jet Diffusion Flames
The problem of vortex/flame interaction is of fundamental importance to turbulent combustion. These interactions have been studied in normal gravity. It was found that due to the interactions between the imposed disturbances and buoyancy induced instabilities, several overall length scales dominated the flame. The problem of multiple scales does not exist in microgravity for a pulsed laminar flame, since there are no buoyancy induced instabilities. The absence of buoyant convection therefore provides an environment to study the role of vortices interacting with flames in a controlled manner. There are strong similarities between imposed and naturally occurring perturbations, since both can be described by the same spatial instability theory. Hence, imposing a harmonic disturbance on a microgravity laminar flame creates effects similar to those occurring naturally in transitional/turbulent diffusion flames observed in microgravity. In this study, controlled, large-scale, axisymmetric vortices are imposed on a microgravity laminar diffusion flame. The experimental results and predictions from a numerical model of transient jet diffusion flames are presented and the characteristics of pulsed flame are described.
Bahadori, M. Y.↗
Quantitative Studies on the Propagation and Extinction of Near-Limit Premixed Flames under Normal- and Micro-gravity
Strained laminar flames have been systematically studied, as the understanding of their structure and dynamic behavior is of relevance to turbulent combustion.. Most of these studies have been conducted in opposed-jet, stagnation-type flow configurations. Studies at high strain rates are important in quantifying and understanding the response of vigorously burning flames and determine extinction states. Studies of weakly strained flames can be of particular interest for all stoichiometries. For example, the laminar flame speeds, S(sup o)(sub u), can be accurately determined by using the counterflow technique only if measurements are obtained at very low strain rates. Furthermore, near-limit flames are stabilized by weak strain rates. Previous studies have shown that near-limit flames are particularly sensitive to chain mechanisms, thermal radiation, and unsteadiness. The stabilization and study of weakly strained flames is complicated by the presence of buoyancy that can render the flames unstable to the point of extinction. Thus, the use of microgravity (mu-g) becomes essential in order to provide meaningful insight into this important combustion regime. In our past studies the laminar flame speeds and extinction strain rates were directly measured at ultra-low strain rates. The laminar flame speeds were measured by having a positively strained planar flame undergoing a transition to a negatively strained Bunsen flame and by measuring the propagation speed during that transition. The extinction strain rates of near-limit flames were measured in mu-g. Results obtained for CH4/air and C3H8/air mixtures are in agreement with those obtained by Maruta et al.
Egolfopoulos, F. N.↗
Flame-Vortex Interactions Imaged in Microgravity - To Assess the Theory Flame Stretch
The goals of this research are to: 1) Assess the Theory of Flame Stretch by operating a unique flame-vortex experiment under microgravity conditions in the NASA Glenn 2.2 Second Drop Tower (drops to identify operating conditions have been completed); 2) Obtain high speed shadowgraph images (500-1000 frames/s) using the drop rig (images were obtained at one-g, and the NASA Kodak RO camera is being mounted on the drop rig); 3) Obtain shadowgraph and PIV images at 1-g while varying the effects of buoyancy by controlling the Froude number (completed); 4) Numerically model the inwardly-propagating spherical flame that is observed in the experiment using full chemistry and the RUN 1DL code (completed); 5) Send images of the flame shape to Dr. G. Patniak at NRL who is numerically simulating the entire flame-vortex interaction of the present experiment (data transfer completed); and 6) Assess the feasibility of obtaining PIV velocity field images in the drop rig, which would be useful (but not required) for our assessment of the Theory of Flame Stretch (PIV images were obtained at one-g using same low laser power that is available from fiber optic cable in drop tower). The motivation for the work is to obtain novel measurement needed to develop a physically accurate model of turbulent combustion that can help in the control of engine pollutants. The unique experiment allows, for the first time, the detailed study of a negatively-curved (negatively stretched) flame, which is one of the five fundamental types of premixed flames. While there have been studies of flat flames, positively-curved (outwardly-propagating) cases and positively-strained (counterflow) cases, this is the first detailed study of a negatively-curved (inwardly-propagating) flame. The first set of drops in the 2.2 Second Drop Tower showed that microgravity provides more favorable conditions for achieving inwardly-propagating flames (IPFs) than 1-g. A vortex interacts with a flame and creates a spherical pocket, which burns inwardly. Shadowgraphs at 1000 frames/sec quantify the Markstein number and flame speed. A Low-Laser Power PIV System was developed and is being added to the drop package. Numerical computations were required to explain why the Markstein numbers measured for the inwardly-propagating flames differ from those of outward propagating flames; this is an important research issue in the assessment of the Theory of Flame Stretch. The RUN-1DL code (developed by Prof. B. Rogg) was run for IPF and OPFs with complex methane and propane chemistry. Results confirmed that Ma for the IPFs are larger than for OPFs as was observed experimentally. Physical reasons for these new findings about the Theory of Flame Stretch are being determined from the experiments and the computations. Several journal papers have been published; the drop package is described in the AIAA Journal, while the one-g results appear in three other journal papers.
Driscoll, James F.↗
Carbon Nanostructure: Its Evolution During its Impact Upon Soot Growth and Oxidation
The proposed work is a ground-based study to define and quantify soot nanostructural changes in response to growth conditions, thermal and oxidative treatments and to quantify their impact upon further oxidation and growth of highly ordered carbon materials. Experimental data relating soot oxidation rates to multiple oxidizing species concentrations will directly test for additive or synergistic soot oxidation rates. Such validation is central for assessing the applicability of individual soot oxidation rates and designing oxidative strategies for controlling soot loadings in and emissions from turbulent combustion processes. Through these experiments, new insights into soot nanostructure evolution during and its impact upon oxidation by O2 and OH will be realized. It is expected that the results of this effort will spawn new research directions in future microgravity and 1g environments. Issues raised by positive or even negative demonstration of the hypotheses of this proposal have direct bearing on modelling and controlling soot formation and its destruction in nearly every combustion process producing soot.
Source record↗
Efficient Low Dissipative High Order Schemes for Multiscale MHD Flows
Accurate numerical simulations of complex multiscale compressible viscous flows, especially high speed turbulence combustion and acoustics, demand high order schemes with adaptive numerical dissipation controls. Standard high resolution shock-capturing methods are too dissipative to capture the small scales and/or long-time wave propagations without extreme grid refinements and small time steps. An integrated approach for the control of numerical dissipation in high order schemes for the compressible Euler and Navier-Stokes equations has been developed and verified by the authors and collaborators. These schemes are suitable for the problems in question. Basically, the scheme consists of sixth-order or higher non-dissipative spatial difference operators as the base scheme. To control the amount of numerical dissipation, multiresolution wavelets are used as sensors to adaptively limit the amount and to aid the selection and/or blending of the appropriate types of numerical dissipation to be used. Magnetohydrodynamics (MHD) waves play a key role in drag reduction in highly maneuverable high speed combat aircraft, in space weather forecasting, and in the understanding of the dynamics of the evolution of our solar system and the main sequence stars. Although there exist a few well-studied second and third-order high-resolution shock-capturing schemes for the MHD in the literature, these schemes are too diffusive and not practical for turbulence/combustion MHD flows. On the other hand, extension of higher than third-order high-resolution schemes to the MHD system of equations is not straightforward. Unlike the hydrodynamic equations, the inviscid MHD system is non-strictly hyperbolic with non-convex fluxes. The wave structures and shock types are different from their hydrodynamic counterparts. Many of the non-traditional hydrodynamic shocks are not fully understood. Consequently, reliable and highly accurate numerical schemes for multiscale MHD equations pose a great challenge to algorithm development. In addition, controlling the numerical error of the divergence free condition of the magnetic fields for high order methods has been a stumbling block. Lower order methods are not practical for the astrophysical problems in question. We propose to extend our hydrodynamics schemes to the MHD equations with several desired properties over commonly used MHD schemes.
Sjoegreen, Bjoern↗
Characteristics of supercritical transitional temporal mixing layers
Several direct numerical simulation realizations of supercritical, three-dimensional, temporal mixing layers are used to investigate aspects of interest to turbulent combustion.
supercritical temporal mixing layers↗
Error Representation in Time For Compressible Flow Calculations
Time plays an essential role in most real world fluid mechanics problems, e.g. turbulence, combustion, acoustic noise, moving geometries, blast waves, etc. Time dependent calculations now dominate the computational landscape at the various NASA Research Centers but the accuracy of these computations is often not well understood. In this presentation, we investigate error representation (and error control) for time-periodic problems as a prelude to the investigation of feasibility of error control for stationary statistics and space-time averages. o These statistics and averages (e.g. time-averaged lift and drag forces) are often the output quantities sought by engineers. o For systems such as the Navier-Stokes equations, pointwise error estimates deteriorate rapidly which increasing Reynolds number while statistics and averages may remain well behaved.
Barth, Timothy J.↗
A Model of Reduced Kinetics for Alkane Oxidation Using Constituents and Species for N-Heptane
The reduction of elementary or skeletal oxidation kinetics to a subgroup of tractable reactions for inclusion in turbulent combustion codes has been the subject of numerous studies. The skeletal mechanism is obtained from the elementary mechanism by removing from it reactions that are considered negligible for the intent of the specific study considered. As of now, there are many chemical reduction methodologies. A methodology for deriving a reduced kinetic mechanism for alkane oxidation is described and applied to n-heptane. The model is based on partitioning the species of the skeletal kinetic mechanism into lights, defined as those having a carbon number smaller than 3, and heavies, which are the complement of the species ensemble. For modeling purposes, the heavy species are mathematically decomposed into constituents, which are similar but not identical to groups in the group additivity theory. From analysis of the LLNL (Lawrence Livermore National Laboratory) skeletal mechanism in conjunction with CHEMKIN II, it is shown that a similarity variable can be formed such that the appropriately non-dimensionalized global constituent molar density exhibits a self-similar behavior over a very wide range of equivalence ratios, initial pressures and initial temperatures that is of interest for predicting n-heptane oxidation. Furthermore, the oxygen and water molar densities are shown to display a quasi-linear behavior with respect to the similarity variable. The light species ensemble is partitioned into quasi-steady and unsteady species. The reduced model is based on concepts consistent with those of Large Eddy Simulation (LES) in which functional forms are used to replace the small scales eliminated through filtering of the governing equations; in LES, these small scales are unimportant as far as the overwhelming part of dynamic energy is concerned. Here, the scales thought unimportant for recovering the thermodynamic energy are removed. The concept is tested by using tabular information from the LLNL skeletal mechanism in conjunction with CHEMKIN II utilized as surrogate ideal functions replacing the necessary functional forms. The test reveals that the similarity concept is indeed justified and that the combustion temperature is well predicted, but that the ignition time is over-predicted, a fact traced to neglecting a detailed description of the processes leading to the heavies chemical decomposition. To palliate this deficiency, functional modeling is incorporated into this conceptual reduction in addition to the modeling the evolution of the global constituent molar density, the enthalpy evolution of the heavies, the contribution to the reaction rate of the unsteady lights from other light species and from the heavies, the molar density evolution of oxygen and water, and the mole fractions of the quasisteady light species. The model is compact in that there are only nine species-related progress variables. Results are presented showing the performance of the model for predicting the temperature and species evolution. The model reproduces the ignition time over a wide range of equivalence ratios, initial pressure, and initial temperature.
Harstad, Kenneth G.↗
Time-Filtered Navier-Stokes Approach and Emulation of Turbulence-Chemistry Interaction
This paper describes the time-filtered Navier-Stokes approach capable of capturing unsteady flow structures important for turbulent mixing and an accompanying subgrid model directly accounting for the major processes in turbulence-chemistry interaction. They have been applied to the computation of two-phase turbulent combustion occurring in a single-element lean-direct-injection combustor. Some of the preliminary results from this computational effort are presented in this paper.
Liu, Nan-Suey↗
Development of a computerized analysis for solid propellant combustion instability with turbulence
A multi-dimensional numerical model has been developed for the unsteady state oscillatory combustion of solid propellants subject to acoustic pressure disturbances. Including the gas phase unsteady effects, the assumption of uniform pressure across the flame zone, which has been conventionally used, is relaxed so that a higher frequency response in the long flame of a double-base propellant can be calculated. The formulation is based on a premixed, laminar flame with a one-step overall chemical reaction and the Arrhenius law of decomposition with no condensed phase reaction. In a given geometry, the Galerkin finite element solution shows the strong resonance and damping effect at the lower frequencies, similar to the result of Denison and Baum. Extended studies deal with the higher frequency region where the pressure varies in the flame thickness. The nonlinear system behavior is investigated by carrying out the second order expansion in wave amplitude when the acoustic pressure oscillations are finite in amplitude. Offset in the burning rate shows a negative sign in the whole frequency region considered, and it verifies the experimental results of Price. Finally, the velocity coupling in the two-dimensional model is discussed.
Chung, T. J.↗
Experimental investigation of the combustion zone in a turbulent flame /eksperimental'noe issledovanie zony goreniya turbulentnogo plameni/
Structure of combustion zone in turbulent flame - measurement of ionization phenomena
COMBUSTION↗