Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Combustion modeling”

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 271 records · Page 15

Application of Flux Conserved Modeling to an Unsteady Combustion Driven Pseudo-Shock

Quasi-steady correlation-based and flux conserved models of the pseudo-shock pressure profile are experimentally assessed in the Air Force Research Lab RC-18 Sub-Scale DirectConnect Isolator Combustion Facility. A bifurcated normal pseudo-shock is established in the 2.67 aspect ratio isolator test section by a downstream cavity-stabilized combustor. The Waltrup-Billig correlation, Smart flux conserved model, and Fiévet data driven flux conserved models are compared against measured pressure profiles. The literature models are found to perform poorly in a high-amplitude forced environment of RC-18, where shock speeds upwards of 20% the isolator Mach number occur. Two modifications to the quasi-steady Fiévet data driven flux conserved model are proposed including incorporating an estimate of pseudo-shock static pressure recovery and exercising the model in a shock stationary reference frame. It is shown that these modifications allow the quasi-steady flux conserved modeling approach to accurately predict pseudo-shock pressure profiles over the discrete phases of high-amplitude forced dynamics.

Shock Train↗

Partial equilibrium model for predicting concentration of CO in combustion

A simple kinetic model, based upon the concept of partial equilibrium, is developed for predicting carbon monoxide concentrations in such steady-flow hydrocarbon-air combustion systems as gas-turbine combustors. The only two kinetic constraints used in the model are on the rate of change of the total number of gaseous particles in the system, and on the rate of change of the CO concentration. The accuracy of the model is verified by comparison with experimental results obtained by burning kerosene with air in a 130,000 Btu/hr atmospheric-pressure steady-flow burner. Burnt-gas cooling rates of 1,000,000 K/sec, about the magnitude encountered in the dilution zone of gas-turbine combustors, were obtained with a compact water-cooled heat exchanger. Good agreement between measured CO concentrations and the values predicted by the partial-equilibrium model was obtained.

Morr, A. R.↗

Forced cocurrent smoldering combustion

An analytical model of cocurrent smoldering combustion through a very porous solid fuel is developed. Smoldering is initiated at the top of a long radially insulated uniform fuel cylinder, so that the smolder wave propagates downward, opposing an upward-forced flow of oxidizer, with the solid fuel and the gaseous oxidizer entering the reaction zone from the same direction (hence, cocurrent). Radiative heat transfer was incorporated using a diffusion approximation, and smoldering was modeled using a one-step reaction mechanism. The results indicate that, for a given fuel, the final temperature depends only on the initial oxygen mass flux, increasing logarithmically with the mass flux. The smolder velocity is linearly dependent on the initial oxygen mass flux, and, at a fixed value of the flux, increases with initial oxygen mass fraction. The mathematical relationship determining the conditions for steady smolder propagation is presented.

Dosanjh, Sudip S.↗

Multigrid Method for Modeling Multi-Dimensional Combustion with Detailed Chemistry

A highly accurate and efficient numerical method is developed for modeling 3-D reacting flows with detailed chemistry. A contravariant velocity-based governing system is developed for general curvilinear coordinates to maintain simplicity of the continuity equation and compactness of the discretization stencil. A fully-implicit backward Euler technique and a third-order monotone upwind-biased scheme on a staggered grid are used for the respective temporal and spatial terms. An efficient semi-coarsening multigrid method based on line-distributive relaxation is used as the flow solver. The species equations are solved in a fully coupled way and the chemical reaction source terms are treated implicitly. Example results are shown for a 3-D gas turbine combustor with strong swirling inflows.

Zheng, Xiaoqing↗

JANNAF 25th Airbreathing Propulsion Subcommittee, 37th Combustion Subcommittee and 1st Modeling and Simulation Subcommittee Joint Meeting

Volume I, the first of three volumes, is a compilation of 24 unclassified/unlimited-distribution technical papers presented at the Joint Army-Navy-NASA-Air Force (JANNAF) 25th Airbreathing Propulsion Subcommittee, 37th Combustion Subcommittee and 1st Modeling and Simulation Subcommittee (MSS) meeting held jointly with the 19th Propulsion Systems Hazards Subcommittee. The meeting was held 13-17 November 2000 at the Naval Postgraduate School and Hyatt Regency Hotel, Monterey, California. Topics covered include: a Keynote Address on Future Combat Systems, a review of the new JANNAF Modeling and Simulation Subcommittee, and technical papers on Hyper-X propulsion development and verification; GTX airbreathing launch vehicles; Hypersonic technology development, including program overviews, fuels for advanced propulsion, ramjet and scramjet research, hypersonic test medium effects; and RBCC engine design and performance, and PDE and UCAV advanced and combined cycle engine technologies.

Fry, Ronald S.↗

Hydrodynamic Instability in an Extended Landau/Levich Model of Liquid-Propellant Combustion at Normal and Reduced Gravity

The classical Landau/Levich models of liquid-propellant combustion, despite their relative simplicity, serve as seminal examples that correctly describe the onset of hydrodynamic instability in reactive systems. Recently, these two separate models have been combined and extended to account for a dynamic dependence, absent in the original formulations, of the local burning rate on the local pressure and temperature fields. The resulting model admits an extremely rich variety of both hydrodynamic and reactive/diffusive instabilities that can be analyzed either numerically or analytically in various limiting parameter regimes. In the present work, a formal asymptotic analysis, based on the realistic smallness of the gas-to-liquid density ratio, is developed to investigate the combined effects of gravity and other parameters on the hydrodynamic instability of the propagating liquid/gas interface. In particular, an analytical expression is derived for the neutral stability boundary A(p)(k), where A(p) is the pressure sensitivity of the burning rate and k is the wavenumber of the disturbance. The results demonstrate explicitly the stabilizing effect of gravity on long-wave disturbances, the stabilizing effect of viscosity (both liquid and gas) and surface tension on short-wave perturbations, and the instability associated with intermediate wavenumbers for negative values of A(p). In the limiting case of weak gravity, it is shown that hydrodynamic instability in liquid-propellant combustion is a long-wave instability phenomenon, whereas at normal gravity, this instability is first manifested through O(1) wavenumber disturbances. it is also demonstrated that, in general, surface tension and the viscosity of both the liquid and gas phases each produce comparable stabilizing effects in the long-wavenumber regime, thereby providing important modifications to previous analyses in which one or more of these effects were neglected.

Margolis, S. B.↗

Hydrodynamic Instability in an Extended Landau/Levich Model of Liquid-Propellant Combustion at Normal and Reduced Gravity

The classical Landau/Levich models of liquid-propellant combustion, despite their relative simplicity, serve as seminal examples that correctly describe the onset of hydrodynamic instability in reactive systems. Recently, these two separate models have been combined and extended to account for a dynamic dependence, absent in the original formulations, of the local burning rate on the local pressure and temperature fields. The resulting model admits an extremely rich variety of both hydrodynamic and reactive/diffusive instabilities that can be analyzed either numerically or analytically in various limiting parameter regimes. In the present work, a formal asymptotic analysis, based on the realistic smallness of the gas-to-liquid density ratio, is developed to investigate the combined effects of gravity and other parameters on the hydrodynamic instability of the propagating liquid/gas interface. In particular, an analytical expression is derived for the neutral stability boundary A(sub p)(k), where A(sub p) is the pressure sensitivity of the burning rate and k is the wavenumber of the disturbance. The results demonstrate explicitly the stabilizing effect of gravity on long-wave disturbances, the stabilizing effect of viscosity (both liquid and gas) and surface tension on short-wave perturbations, and the instability associated with intermediate wavenumbers for critical negative values of A(sub p). In the limiting case of weak gravity, it is shown that hydrodynamic instability in liquid-propellant combustion is a long-wave instability phenomenon, whereas at normal gravity, this instability is first manifested through O(1) wavenumber disturbances. It is also demonstrated that, in general, surface tension and the viscosity of both the liquid and gas phases each produce comparable stabilizing effects in the large-wavenumber regime, thereby providing important modifications to previous analyses in which one or more of these effects were neglected.

Margolis, Stephen B.↗

Stress and deformation modeling of multiple rotary combustion engine trochoid housings

This paper documents the development of the capability to produce finite element models of alternate trochoid housing configurations. The effort needed to produce these models is greatly reduced by the use of a newly developed specialized finite element preprocessor which is described. The results of static stress comparisons conducted on a Mazda trochoid housing are presented. Planned future development of this modeling capability to operational situations is also presented.

Lychuk, W. M.↗

Frame-Transfer Gating Raman Spectroscopy for Time-Resolved Multiscalar Combustion Diagnostics

Accurate experimental measurement of spatially and temporally resolved variations in chemical composition (species concentrations) and temperature in turbulent flames is vital for characterizing the complex phenomena occurring in most practical combustion systems. These diagnostic measurements are called multiscalar because they are capable of acquiring multiple scalar quantities simultaneously. Multiscalar diagnostics also play a critical role in the area of computational code validation. In order to improve the design of combustion devices, computational codes for modeling turbulent combustion are often used to speed up and optimize the development process. The experimental validation of these codes is a critical step in accepting their predictions for engine performance in the absence of cost-prohibitive testing. One of the most critical aspects of setting up a time-resolved stimulated Raman scattering (SRS) diagnostic system is the temporal optical gating scheme. A short optical gate is necessary in order for weak SRS signals to be detected with a good signal- to-noise ratio (SNR) in the presence of strong background optical emissions. This time-synchronized optical gating is a classical problem even to other spectroscopic techniques such as laser-induced fluorescence (LIF) or laser-induced breakdown spectroscopy (LIBS). Traditionally, experimenters have had basically two options for gating: (1) an electronic means of gating using an image intensifier before the charge-coupled-device (CCD), or (2) a mechanical optical shutter (a rotary chopper/mechanical shutter combination). A new diagnostic technology has been developed at the NASA Glenn Research Center that utilizes a frame-transfer CCD sensor, in conjunction with a pulsed laser and multiplex optical fiber collection, to realize time-resolved Raman spectroscopy of turbulent flames that is free from optical background noise (interference). The technology permits not only shorter temporal optical gating (down to <1 s, in principle), but also higher optical throughput, thus resulting in a substantial increase in measurement SNR.

Nguyen, Quang-Viet↗

Effect of model selection on combustor performance and stability using ROCCID

The ROCket Combustor Interactive Design (ROCCID) methodology is an interactive computer program that combines previously developed combustion analysis models to calculate the combustion performance and stability of liquid rocket engines. Test data from a 213 kN (48,000 lbf) Liquid Oxygen (LOX)/RP-1 combustor with a O-F-O (oxidizer-fuel-oxidizer) triplet injector were used to characterize the predictive capabilities of the ROCCID analysis models for this injector/propellant configuration. Thirteen combustion performance and stability models have been incorporated into ROCCID, and ten of them, which have options for triplet injectors, were examined in this study. Calculations using different combinations of analysis models, with little or no anchoring, were carried out on a test matrix of operating conditions matching those of the test program. Results of the computer analyses were compared to test data, and the ability of the model combinations to correctly predict combustion stability or instability was determined. For the best model combination(s), sensitivity of the calculations to fuel drop size and mixing efficiency was examined. Error in the stability calculations due to uncertainty in the pressure interaction index (N) was examined. The recommended model combinations for this O-F-O triplet LOX/RP-1 configuration are proposed.

Giuliani, James E.↗

Non-isotropic turbulence effects on spray combustion

A numerical model for the prediction of local properties of statistically stationary spray-combusting flows is evaluated by comparison with experimental data. To appraise the relative performance of turbulence models, computations were carried out by the k-epsilon model and the algebraic stress model. The present numerical results show the qualitative agreement with experimental data. In terms of overall local flow properties, the algebraic stress model improves a degree of conformity to the experimental data due to its ability to introduce the nonisotropic turbulence effects. Two swirl numbers are considered to investigate the influence of swirl on the droplet evaporation and trajectories, and the effects of droplet/turbulence interactions in flow properties. It is found that the large swirl produces a higher evaporationn rate, and more intensive turbulent mixing and burning. The discrepancies observed in the results are attributed mainly to uncertainties in the initial spray size and velocity distributions, the droplet/wall impingement interaction, the combustion model with the fast chemistry and the turbulence models dealing with the strong streamline curvature and complex interactions between the dispersed droplets and the continuous gas-phase flows.

Kim, Y. M.↗

Effect of model selection on combustor performance and stability predictions using ROCCID

The ROCket Combustor Interactive Design (ROCCID) methodology is an interactive computer program that combines previously developed combustion analysis models to calculate the combustion performance and stability of liquid rocket engines. Test data from 213 kN (48,000 lbf) Liquid Oxygen (LOX)/RP-1 combustor with an O-F-O (oxidizer-fuel-oxidizer) triplet injector were used to characterize the predictive capabilities of the ROCCID analysis models for this injector/propellant configuration. Thirteen combustion performance and stability models were incorporated into ROCCID, and ten of them, which have options for triplet injectors, were examined. Calculations using different combinations of analysis models, with little or no anchoring, were carried out on a test matrix of operating combinations matching those of the test program. Results of the computer analyses were compared to test data, and the ability of the model combinations to correctly predict combustion stability or instability was determined. For the best model combination(s), sensitivity of the calculations to fuel drop size and mixing efficiency was examined. Error in the stability calculations due to uncertainty in the pressure interaction index (N) was examined. The recommended model combinations for this O-F-O triplet LOX/RP-1 configuration are proposed.

Giuliani, James E.↗

Thermophysics Characterization of Kerosene Combustion

A one-formula surrogate fuel formulation and its quasi-global combustion kinetics model are developed to support the design of injectors and thrust chambers of kerosene-fueled rocket engines. This surrogate fuel model depicts a fuel blend that properly represents the general physical and chemical properties of kerosene. The accompanying gaseous-phase thermodynamics of the surrogate fuel is anchored with the heat of formation of kerosene and verified by comparing a series of one-dimensional rocket thrust chamber calculations. The quasi-global combustion kinetics model consists of several global steps for parent fuel decomposition, soot formation, and soot oxidation and a detailed wet-CO mechanism to complete the combustion process. The final thermophysics formulations are incorporated with a computational fluid dynamics model for prediction of the combustion efficiency of an unielement, tripropellant combustor and the radiation of a kerosene-fueled thruster plume. The model predictions agreed reasonably well with those of the tests.

Wang, Ten-See↗

3-D CFD Simulation and Validation of Oxygen-Rich Hydrocarbon Combustion in a Gas-Centered Swirl Coaxial Injector using a Flamelet-Based Approach

Injector design is a critical part of the development of a rocket Thrust Chamber Assembly (TCA). Proper detailed injector design can maximize propulsion efficiency while minimizing the potential for failures in the combustion chamber. Traditional design and analysis methods for hydrocarbon-fuel injector elements are based heavily on empirical data and models developed from heritage hardware tests. Using this limited set of data produces challenges when trying to design a new propulsion system where the operating conditions may greatly differ from heritage applications. Time-accurate, Three-Dimensional (3-D) Computational Fluid Dynamics (CFD) modeling of combusting flows inside of injectors has long been a goal of the fluid analysis group at Marshall Space Flight Center (MSFC) and the larger CFD modeling community. CFD simulation can provide insight into the design and function of an injector that cannot be obtained easily through testing or empirical comparisons to existing hardware. However, the traditional finite-rate chemistry modeling approach utilized to simulate combusting flows for complex fuels, such as Rocket Propellant-2 (RP-2), is prohibitively expensive and time consuming even with a large amount of computational resources. MSFC has been working, in partnership with Streamline Numerics, Inc., to develop a computationally efficient, flamelet-based approach for modeling complex combusting flow applications. In this work, a flamelet modeling approach is used to simulate time-accurate, 3-D, combusting flow inside a single Gas Centered Swirl Coaxial (GCSC) injector using the flow solver, Loci-STREAM. CFD simulations were performed for several different injector geometries. Results of the CFD analysis helped guide the design of the injector from an initial concept to a tested prototype. The results of the CFD analysis are compared to data gathered from several hot-fire, single element injector tests performed in the Air Force Research Lab EC-1 test facility located at Edwards Air Force Base.

Richardson, Brian↗

Finite area combustor theoretical rocket performance

Previous to this report, the computer program of NASA SP-273 and NASA TM-86885 was capable of calculating theoretical rocket performance based only on the assumption of an infinite area combustion chamber (IAC). An option was added to this program which now also permits the calculation of rocket performance based on the assumption of a finite area combustion chamber (FAC). In the FAC model, the combustion process in the cylindrical chamber is assumed to be adiabatic, but nonisentropic. This results in a stagnation pressure drop from the injector face to the end of the chamber and a lower calculated performance for the FAC model than the IAC model.

Gordon, Sanford↗

Numerical Simulation for Droplet Combustion Using Lagrangian Hydrodynamics

A predictive model of spray combustion must incorporate models for the wide variety of physical environments in a practical combustor. In regions where droplets are closely spaced, combustion resembles a diffusion flame; where they are well separated, an envelope or wake flame results. The relative velocity field between the fuel droplets and oxidizer in influences boundary layer development about the droplet, recirculating flow patterns, and droplet shape and stability. A model must encompass these interacting temporal and spatial effects as well as complicated combustor boundaries. The objective of the current work is to develop the triangular gridding method for describing the individual and collective properties of vaporizing and burning fuel droplets.

Fritts, M. J.↗

Improved Modeling of Finite-Rate Turbulent Combustion Processes in Research Combustors

The objective of this thesis is to further develop and test a stochastic model of turbulent combustion in recirculating flows. There is a requirement to increase the accuracy of multi-dimensional combustion predictions. As turbulence affects reaction rates, this interaction must be more accurately evaluated. In this work a more physically correct way of handling the interaction of turbulence on combustion is further developed and tested. As turbulence involves randomness, stochastic modeling is used. Averaged values such as temperature and species concentration are found by integrating the probability density function (pdf) over the range of the scalar. The model in this work does not assume the pdf type, but solves for the evolution of the pdf using the Monte Carlo solution technique. The model is further developed by including a more robust reaction solver, by using accurate thermodynamics and by more accurate transport elements. The stochastic method is used with Semi-Implicit Method for Pressure-Linked Equations. The SIMPLE method is used to solve for velocity, pressure, turbulent kinetic energy and dissipation. The pdf solver solves for temperature and species concentration. Thus, the method is partially familiar to combustor engineers. The method is compared to benchmark experimental data and baseline calculations. The baseline method was tested on isothermal flows, evaporating sprays and combusting sprays. Pdf and baseline predictions were performed for three diffusion flames and one premixed flame. The pdf method predicted lower combustion rates than the baseline method in agreement with the data, except for the premixed flame. The baseline and stochastic predictions bounded the experimental data for the premixed flame. The use of a continuous mixing model or relax to mean mixing model had little effect on the prediction of average temperature. Two grids were used in a hydrogen diffusion flame simulation. Grid density did not effect the predictions except for peak temperature and tangential velocity. The hybrid pdf method did take longer and required more memory, but has a theoretical basis to extend to many reaction steps which cannot be said of current turbulent combustion models.

VanOverbeke, Thomas J.↗

Vapor deposition and condensate flow on combustion turbine blades - Theoretical model to predict/understand some corrosion rate consequences of molten alkali sulfate deposition in the field or laboratory

An analysis is undertaken of aerodynamically- and centrifugally-driven liquid condensate layers on nonisothermal combustion turbines' stator vanes and rotor blades. Attention is given to the quantitative consequences of one possible mechanism for the initiation of 'hot corrosion' in the underlying blade material through a 'fluxing' of the protective oxide coating by the molten salt of the Newtonian condensate film. Illustrative calculations are presented for the condensate streamline pattern and the distributions of the steady-state condensate layer thickness, together with the corresponding oxide dissolution rate, for a test turbine blade.

Rosner, Daniel E.↗