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 829 records · Page 46

Hybrid propulsion for launch vehicle boosters: A program status update

Results obtained in studying the origin and suppression of large-amplitude pressure oscillations in a 24 in. diameter hybrid motor using a liquid oxygen/hydroxylterminated polybutadiene/polycyclopentadiene propellant system are discussed. Tests conducted with liquid oxygen flow rates varying from 10 to 40 lbm/sec were designed to gauge the effectiveness of various vaporization chamber flow fields, injector designs, and levels of heat addition in suppressing high-frequency longitudinal mode oscillations. Longitudinal acoustic modes did not arise in any tests. However, initial testing revealed the presence of high-amplitude, sinusoidal, nonacoustic oscillations persisting throughout the burn durations. Analysis showed this to be analogous to chug mode instability in liquid rocket engines brought about by a coupling of motor combustion processes and the liquid oxygen feed system. Analytical models were developed and verified by test data to predict the amplitude and frequency of feed-system-coupled combustion pressure oscillations. Subsequent testing showed that increasing the feed system impedance eliminated the bulk mode instability. This paper documents the work completed to date in performance of the Hybrid Propulsion Technology for Launch Vehicle Boosters Program (NAS8-39942) sponsored by NASA's George C. Marshall Space Flight Center.

Carpenter, R. L.↗

Heat and mass transfer in combustion - Fundamental concepts and analytical techniques

Fundamental combustion phenomena and the associated flame structures in laminar gaseous flows are discussed on physical bases within the framework of the three nondimensional parameters of interest to heat and mass transfer in chemically-reacting flows, namely the Damkoehler number, the Lewis number, and the Arrhenius number which is the ratio of the reaction activation energy to the characteristic thermal energy. The model problems selected for illustration are droplet combustion, boundary layer combustion, and the propagation, flammability, and stability of premixed flames. Fundamental concepts discussed include the flame structures for large activation energy reactions, S-curve interpretation of the ignition and extinctin states, reaction-induced local-similarity and non-similarity in boundary layer flows, the origin and removal of the cold boundary difficulty in modeling flame propagation, and effects of flame stretch and preferential diffusion on flame extinction and stability. Analytical techniques introduced include the Shvab-Zeldovich formulation, the local Shvab-Zeldovich formulation, flame-sheet approximation and the associated jump formulation, and large activation energy matched asymptotic analysis. Potentially promising research areas are suggested.

Law, C. K.↗

CFD propels NASP propulsion progress

The most complex aerothermodynamics encountered in the National Aerospace Plane (NASP) propulsion system are associated with the fuel-mixing and combustion-reaction flows of its combustor section; adequate CFD tools must be developed to model shock-wave systems, turbulent hydrogen/air mixing, flow separation, and combustion. Improvements to existing CFD codes have involved extension from two dimensions to three, as well as the addition of finite-rate hydrogen-air chemistry. A novel CFD code for the treatment of reacting flows throughout the NASP, designated GASP, uses the most advanced upwind-differencing technology.

Povinelli, Louis A.↗

LSPRAY: Lagrangian Spray Solver for Applications With Parallel Computing and Unstructured Gas-Phase Flow Solvers

Sprays occur in a wide variety of industrial and power applications and in the processing of materials. A liquid spray is a phase flow with a gas as the continuous phase and a liquid as the dispersed phase (in the form of droplets or ligaments). Interactions between the two phases, which are coupled through exchanges of mass, momentum, and energy, can occur in different ways at different times and locations involving various thermal, mass, and fluid dynamic factors. An understanding of the flow, combustion, and thermal properties of a rapidly vaporizing spray requires careful modeling of the rate-controlling processes associated with the spray's turbulent transport, mixing, chemical kinetics, evaporation, and spreading rates, as well as other phenomena. In an attempt to advance the state-of-the-art in multidimensional numerical methods, we at the NASA Lewis Research Center extended our previous work on sprays to unstructured grids and parallel computing. LSPRAY, which was developed by M.S. Raju of Nyma, Inc., is designed to be massively parallel and could easily be coupled with any existing gas-phase flow and/or Monte Carlo probability density function (PDF) solver. The LSPRAY solver accommodates the use of an unstructured mesh with mixed triangular, quadrilateral, and/or tetrahedral elements in the gas-phase solvers. It is used specifically for fuel sprays within gas turbine combustors, but it has many other uses. The spray model used in LSPRAY provided favorable results when applied to stratified-charge rotary combustion (Wankel) engines and several other confined and unconfined spray flames. The source code will be available with the National Combustion Code (NCC) as a complete package.

Raju, Manthena S.↗

Burning and Suppression of Solids–II (BASS-II) Summary Report

The Burning and Suppression of Solids (BASS) experiment hardware is a small flow duct that provides containment for small scale burning of solid samples within the Microgravity Science Glovebox (MSG) aboard the International Space Station (ISS). A video camera with a data overlay and a 35-mm still camera record the combustion events. The controls (ignition, fan speed, etc.) are operated by an astronaut while the principal investigator team monitors the experiment from the ground and communicates directly with the astronaut. For the first time on ISS, BASS–II utilized MSG working volume dilution with gaseous N2. We developed a perfectly stirred reactor model to determine the N2 flow time and flow rate to obtain the desired reduced O2 concentration in the working volume for each test. We calibrated the model with the Compound Specific Analyzer-Combustion Products (CSA-CP) O2 readings offset using the Major Constituents Analyzer reading of the ISS ambient atmosphere data for that day. This worked out extremely well for operations, and added a new vital variable, ambient O2 level, to our test matrices. The main variables tested in BASS–II were ambient O2 concentration, ventilation flow velocity, and fuel type, thickness, and geometry. BASS–II also utilized the onboard CSA-CP for O2 and CO readings, and the Carbon Dioxide Monitor for CO2 readings before and after each test. Readings from these sensors allow us to evaluate the completeness of the combustion. The O2 and CO2 readings before and after each test were analyzed and compared very well to stoichiometric ratios for a one-step gas-phase reaction. The CO versus CO2 followed a linear trend for some datasets, but not for all the different geometries of fuel and flow tested. We calculated the heat release rates during each test from the O2 consumption and burn times, using the constant 13.1 kJ of heat released per gram of O2 consumed. The results showed that most of the tests had heat release rates well below 100 W. Lastly, the global equivalence ratio for the tests is estimated to be fuel rich, 1.3 on average using mass loss and O2 consumption data.

Microgravity↗

Vapour phase details in the oscillatory combustion of propellants A porous plate analogue

A perforated porous plate burner is designed to scale up the vapor phase details in composite propellant combustion. A fuel vapor is passed through the pores and an oxidizer vapor is passed through the discrete perforations drilled in the plate. Discussion of the scale modeling of the vapor phase details in oscillatory combustion of composite propellants leads to the following conclusions: (1) the concept of a perforated porous plate analog for the composite propellant vapor phase combustion appears to be valid during oscillatory combustion as well; (2) the flame standoff distance can be conveniently determined with motion picture photography; (3) the flame standoff distance varies with the flow velocity during oscillatory combustion much the same way as during time-independent combustion; and (4) the overall structure of the combustion zone (flame) does not appear to vary too much from its time-independent structure during oscillatory combustion

Ramohalli, K.↗

Composition and Climate Impacts of Increasing Launches to Low Earth Orbit

We present simulations of an Earth system model in which launch vehicles inject gaseous and particulate combustion products into the stratosphere. We considered two plausible scenarios representative for the year 2050, one with 1000 launches per year and another with 10 times more, all of which assuming heavy lift (80 tons to low Earth orbit; LEO) launches using methane as a proxy for liquefied natural gas (LNG) fuel. An industry-standard plume flowfield model was used to predict emission of gaseous species including carbon monoxide, nitrogen oxides, and water vapor. Black carbon was introduced assuming vacuum emission equal 25% of equivalent kerosene engine. The gas emissions showed marginal or insignificant changes in atmospheric composition, while the black carbon emission was found to be the most influential component. Feedbacks involve tropopause and stratospheric warming, followed by a moistening of the stratosphere via tropospheric water vapor intrusions, together with a reduction in global albedo. The fact that LNG-fueled rockets could be a concern from aerosol emissions in our simulations, despite burning more efficiently than conventional fuels, underscores the need to better understand actual black carbon emissions from rocket engines in order to accurately predict the global impacts of launch vehicles on future climate.

Earth system model↗

Turbulent flame propagation in partially premixed flames

Turbulent premixed flame propagation is essential in many practical devices. In the past, fundamental and modeling studies of propagating flames have generally focused on turbulent flame propagation in mixtures of homogeneous composition, i.e. a mixture where the fuel-oxidizer mass ratio, or equivalence ratio, is uniform. This situation corresponds to the ideal case of perfect premixing between fuel and oxidizer. In practical situations, however, deviations from this ideal case occur frequently. In stratified reciprocating engines, fuel injection and large-scale flow motions are fine-tuned to create a mean gradient of equivalence ratio in the combustion chamber which provides additional control on combustion performance. In aircraft engines, combustion occurs with fuel and secondary air injected at various locations resulting in a nonuniform equivalence ratio. In both examples, mean values of the equivalence ratio can exhibit strong spatial and temporal variations. These variations in mixture composition are particularly significant in engines that use direct fuel injection into the combustion chamber. In this case, the liquid fuel does not always completely vaporize and mix before combustion occurs, resulting in persistent rich and lean pockets into which the turbulent flame propagates. From a practical point of view, there are several basic and important issues regarding partially premixed combustion that need to be resolved. Two such issues are how reactant composition inhomogeneities affect the laminar and turbulent flame speeds, and how the burnt gas temperature varies as a function of these inhomogeneities. Knowledge of the flame speed is critical in optimizing combustion performance, and the minimization of pollutant emissions relies heavily on the temperature in the burnt gases. Another application of partially premixed combustion is found in the field of active control of turbulent combustion. One possible technique of active control consists of pulsating the fuel flow rate and thereby modulating the equivalence ratio (Bloxsidge et al. 1987). Models of partially premixed combustion would be extremely useful in addressing all these questions related to practical systems. Unfortunately, the lack of a fundamental understanding regarding partially premixed combustion has resulted in an absence of models which accurately capture the complex nature of these flames. Previous work on partially premixed combustion has focused primarily on laminar triple flames. Triple flames correspond to an extreme case where fuel and oxidizer are initially totally separated (Veynante et al. 1994 and Ruetsch et al. 1995). These flames have a nontrivial propagation speed and are believed to be a key element in the stabilization process of jet diffusion flames. Different theories have also been proposed in the literature to describe a turbulent flame propagating in a mixture with variable equivalence ratio (Muller et al. 1994), but few validations are available. The objective of the present study is to provide basic information on the effects of partial premixing in turbulent combustion. In the following, we use direct numerical simulations to study laminar and turbulent flame propagation with variable equivalence ratio.

Poinsot, T.↗

Parabolized Navier-Stokes analysis of ducted turbulent mixing problems with finite-rate chemistry

A hybrid implicit/explicit approach for analyzing 2D ducted turbulent mixing problems with finite-rate chemistry is presented. The approach combines a fully-implicit parabolic mixing algorithm, an explicit viscous-characteristic based shock-capturing algorithm, and linearized implicit chemical kinetic algorithm. The resultant model provides spatial marching solutions of the parabolized Navier-Stokes (PNS) equations for supersonic combustion and viscous nozzle flowfields. Specialized procedures are incorporated to deal with the near wall sublayer and with small central pockets of subsonic flow. A parabolic option is available which can be utilized in the direct or inverse mode for design applications. The ability of the model to treat shock waves and wave/mixing layer interactions is assessed via comparisons of predictions with those of well established explicit shock-capturing Euler (SCIPPY) and PNS (SCIPVIS) models. Applications to a variety of supersonic combustion flowfield problems involving tangential or moderately inclined fuel injection, and, to a viscous nozzle flow problem, are presented. This paper serves to exhibit overall capabilities of the model developed and no comparisons with supersonic combustion data are presented. Such comparisons serve mainly to validate the modeling of the turbulence and turbulence/chemical interactions and will be the subject of a future paper.

Sinha, N.↗

Chemical kinetic modeling of propane oxidation behind shock waves

The stoichiometric combustion of propane behind incident shock waves was studied experimentally and analytically over a temperature range from 1700 K to 2600 K and a pressure range from 1.2 to 1.9 atm. Measurements of the concentrations of carbon monoxide (CO) and carbon dioxide (CO2) and the product of the oxygen atom and carbon dioxide concentrations (O)(CO) were made after passage of the incident shock wave. A kinetic mechanism was developed which, when used in a computer program for a flowing, reacting gas behind an incident shock wave predicted experimentally measured results quite well. Ignition delay times from the literature were also predicted quite well. The kinetic mechanism consisted of 59 individual kinetic steps.

Mclain, A. G.↗

Kinetics and mechanism of soot formation in hydrocarbon combustion

The focus of this work was on kinetic modeling. The specific objectives were: detailed modeling of soot formation in premixed flames, elucidation of the effects of fuel structure on the pathway to soot, and the development of a numerical technique for accurate modeling of soot particle coagulation and surface growth. Those tasks were successfully completed and are briefly summarized.

Frenklach, Michael↗

Flammability of Heterogeneously Combusting Metals

Most engineering materials, including some metals, most notably aluminum, burn in homogeneous combustion. 'Homogeneous' refers to both the fuel and the oxidizer being in the same phase, which is usually gaseous. The fuel and oxidizer are well mixed in the combustion reaction zone, and heat is released according to some relation like q(sub c) = delta H(sub c)c[((rho/rho(sub 0))]exp a)(exp -E(sub c)/RT), Eq. (1) where the pressure exponent a is usually close to unity. As long as there is enough heat released, combustion is sustained. It is useful to conceive of a threshold pressure beyond which there is sufficient heat to keep the temperature high enough to sustain combustion, and beneath which the heat is so low that temperature drains away and the combustion is extinguished. Some materials burn in heterogeneous combustion, in which the fuel and oxidizer are in different phases. These include iron and nickel based alloys, which burn in the liquid phase with gaseous oxygen. Heterogeneous combustion takes place on the surface of the material (fuel). Products of combustion may appear as a solid slag (oxide) which progressively covers the fuel. Propagation of the combustion melts and exposes fresh fuel. Heterogeneous combustion heat release also follows the general form of Eq.(1), except that the pressure exponent a tends to be much less than 1. Therefore, the increase in heat release with increasing pressure is not as dramatic as it is in homogeneous combustion. Although the concept of a threshold pressure still holds in heterogeneous combustion, the threshold is more difficult to identify experimentally, and pressure itself becomes less important relative to the heat transfer paths extant in any specific application. However, the constants C, a, and E(sub c) may still be identified by suitable data reduction from heterogeneous combustion experiments, and may be applied in a heat transfer model to judge the flammability of a material in any particular actual-use situation. In order to support the above assertions, two investigations are undertaken: 1) PCT data are examined in detail to discover the pressure dependence of heterogeneous combustion experiment results; and 2) heterogeneous combustion in a PCT situation is described by a heat transfer model, which is solved first in simplified form for a simple actual-use situation, and then extended to apply to PCT data reduction (combustion constant identification).

Jones, Peter D.↗

CFD modeling of turbulent flows around the SSME main injector assembly using porosity formulation

Hot gas turbulent flow distribution around the main injector assembly of the Space Shuttle Main Engine (SSME) and LOX flow distribution through the LOX posts have a great effect on the combustion phenomena inside the main combustion chamber. In order to design a CFD model to be an effective engineering analysis tool with good computational turn-around time (especially for 3-D flow problems) and still maintain good accuracy in describing the flow features, the concept of porosity was employed to describe the effects of blockage and drag force due to the presence of the LOX posts in the turbulent flow field around the main injector assembly of the SSME. Two-dimensional numerical studies were conducted to identify the drag coefficients of the flows, both through tube banks and round the shielded posts, over a wide range of Reynolds numbers. Empirical, analytical expressions of the drag coefficients as a function of local flow Reynolds number were then deduced. The porosity model was applied to the turbulent flow around the main injector assembly of the SSME, and analyses were performed. The 3-D CFD analysis was divided into three parts: LOX dome, hot gas injector assembly, and hydrogen cavity. The numerical results indicate that the mixture ratio at the downstream of injector face was close to stoichiometric around baffle elements.

Cheng, Gary C.↗

Aerodynamic properties of turbulent combustion fields

Flow fields involving turbulent flames in premixed gases under a variety of conditions are modeled by the use of a numerical technique based on the random vortex method to solve the Navier-Stokes equations and a flame propagation algorithm to trace the motion of the front and implement the Huygens principle, both due to Chorin. A successive over-relaxation hybrid method is applied to solve the Euler equation for flows in an arbitrarily shaped domain. The method of images, conformal transformation, and the integral-equation technique are also used to treat flows in special cases, according to their particular requirements. Salient features of turbulent flame propagation in premixed gases are interpreted by relating them to the aerodynamic properties of the flow field. Included among them is the well-known cellular structure of flames stabilized by bluff bodies, as well as the formation of the characteristic tulip shape of flames propagating in ducts. In its rudimentary form, the mechanism of propagation of a turbulent flame is shown to consist of: (1) rotary motion of eddies at the flame front, (2) self-advancement of the front at an appropriate normal burning speed, and (3) dynamic effects of expansion due to exothermicity of the combustion reaction. An idealized model is used to illustrate these fundamental mechanisms and to investigate basic aerodynamic features of flames in premixed gases. The case of a confined flame stabilized behind a rearward-facing step is given particular care and attention. Solutions are shown to be in satisfactory agreement with experimental results, especially with respect to global properties such as the average velocity profiles and reattachment length.

Hsiao, C. C.↗

Simulation of Liquid Rocket Engine Failure Propagation Using Self-Evolving Scenarios

Traditional probabilistic risk assessment approaches often require failure scenarios to be explicitly defined through event sequences that are then quantified as part of the integrated analysis. This approach becomes difficult when failure propagation paths change as a function of the system operation. Additionally, if the propagation paths represent interactions among even a modest number of components, the scenario count becomes combinatorially intractable. This paper presents an alternate approach for quantifying the probability of failure propagation in such a case. Rather than explicitly defining scenario sequences, simple physical models are created for each of the components. In this way, only the physical states and rules of component interactions must be defined, rather than event sequences for each individual scenario. Initiating failures are introduced into the system, either randomly or as defined by relative likelihood, and the failures cascade through the system via the interaction rules. This process is repeated using Monte Carlo methods and, as a result, the most probable scenarios “self-evolve” in terms of both sequence path and frequency. This approach was applied to failures occurring in the engine compartment of a space launch vehicle with four liquid rocket engines and four high-pressure helium tanks. Each engine was modeled with key components, such as turbomachinery, combustion chamber, propellant lines, and additional support systems. Three test cases were conducted with different high-energy engine failures. End results of interest included an additional engine-out failure and tank burst, which represent the loss-of-mission (LOM) and loss-of-crew (LOC) failure environments, respectively. Observations show that almost every scenario outcome is unique and that many scenarios involve complex chain reactions that are difficult to predict. This validates the usefulness of the modeling approach in assessing the overall risks to the crew during a launch vehicle abort.

Rocket Engine↗

Ignition and combustion characteristics of metallized propellants, phase 2

Experimental and analytical investigations focusing on aluminum/hydrocarbon gel droplet secondary atomization and its effects on gel-fueled rocket engine performance are being conducted. A single laser sheet sizing/velocimetry diagnostic technique, which should eliminate sizing bias in the data collection process, has been designed and constructed to overcome limitations of the two-color forward-scatter technique used in previous work. Calibration of this system is in progress and the data acquisition/validation code is being written. Narrow-band measurements of radiant emission, discussed in previous reports, will be used to determine if aluminum ignition has occurred in a gel droplet. A one-dimensional model of a gel-fueled rocket combustion chamber, described in earlier reports, has been exercised in conjunction with a two-dimensional, two-phase nozzle code to predict the performance of an aluminum/hydrocarbon fueled engine. Estimated secondary atomization effects on propellant burnout distance, condensed particle radiation losses to the chamber walls, and nozzle two phase flow losses are also investigated. Calculations indicate that only modest secondary atomization is required to significantly reduce propellant burnout distances, aluminum oxide residual size, and radiation heat losses. Radiation losses equal to approximately 2-13 percent of the energy released during combustion were estimated, depending on secondary atomization intensity. A two-dimensional, two-phase nozzle code was employed to estimate radiation and nozzle two phase flow effects on overall engine performance. Radiation losses yielded a one percent decrease in engine Isp. Results also indicate that secondary atomization may have less effect on two-phase losses than it does on propellant burnout distance and no effect if oxide particle coagulation and shear induced droplet breakup govern oxide particle size. Engine Isp was found to decrease from 337.4 to 293.7 seconds as gel aluminum mass loading was varied from 0-70 wt percent. Engine Isp efficiencies, accounting for radiation and two phase flow effects, on the order of 0.946 were calculated for a 60 wt percent gel, assuming a fragmentation ratio of five.

Mueller, D. C.↗

Theoretical studies of the ignition and combustion of silane-hydrogen-air mixtures

A chemical kinetic mechanism is proposed for the combustion of silane-hydrogen-oxygen-nitrogen mixtures in the initial temperature range from 800K to 1250K and pressure range from 0.5 to 1.35 atm. The mechanism yields results which are in agreement with published ignition delay times obtained from shock tube experiments. Comparisons between the results obtained using the proposed mechanism and that of an alternative mechanism reveal that the former predicts appreciably shorter ignition delay times, but a flame blowout envelope which is shifted so as to decrease the stable flame region. Over much of the thermodynamic range examined, the mechanism predicts long reaction times. A three step global mechanism is proposed which closely models the ignition phase of SiH4 - H2 - air combustion; however, the reaction phase is less well reproduced by the global model. The necessity for additional experimental data to further assess the proposed models is stressed.

Chinitz, W.↗

Rocket Engine Oscillation Diagnostics

Rocket engine oscillating data can reveal many physical phenomena ranging from unsteady flow and acoustics to rotordynamics and structural dynamics. Because of this, engine diagnostics based on oscillation data should employ both signal analysis and physical modeling. This paper describes an approach to rocket engine oscillation diagnostics, types of problems encountered, and example problems solved. Determination of design guidelines and environments (or loads) from oscillating phenomena is required during initial stages of rocket engine design, while the additional tasks of health monitoring, incipient failure detection, and anomaly diagnostics occur during engine development and operation. Oscillations in rocket engines are typically related to flow driven acoustics, flow excited structures, or rotational forces. Additional sources of oscillatory energy are combustion and cavitation. Included in the example problems is a sampling of signal analysis tools employed in diagnostics. The rocket engine hardware includes combustion devices, valves, turbopumps, and ducts. Simple models of an oscillating fluid system or structure can be constructed to estimate pertinent dynamic parameters governing the unsteady behavior of engine systems or components. In the example problems it is shown that simple physical modeling when combined with signal analysis can be successfully employed to diagnose complex rocket engine oscillatory phenomena.

Nesman, Tom↗