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 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

SSME thrust chamber modeling with Navier Stokes equations

The capability of predicting two-dimensional, compressible and reacting flow in the combustion chamber and nozzle of the Space Shuttle Main Engine (SSME) is demonstrated. A nonorthogonal body fitted coordinate system has been used to represent the combustor and nozzle geometry. The Navier-Stokes equations are solved for the entire thrust chamber with the k-epsilon turbulence model accounting for compressibility and large pressure gradients effects. Results of the computational test cases reveal all expected features of the transonic nozzle flows including location of sonic line, internal shock and boundary layer build-up. Calculated performance parameters such as thrust, flow rate, and specific impulse are also in reasonble agreement with available data. The results show promising potential of solving full Navier-Stokes equations with heat transfer and two-phase combustion in truly comprehensive modeling of rocket engines.

Przekwas, A. J.↗
NASA NTRS1993

CFD simulation of coaxial injectors

The development of improved performance models for the Space Shuttle Main Engine (SSME) is an important, ongoing program at NASA MSFC. These models allow prediction of overall system performance, as well as analysis of run-time anomalies which might adversely affect engine performance or safety. Due to the complexity of the flow fields associated with the SSME, NASA has increasingly turned to Computational Fluid Dynamics (CFD) techniques as modeling tools. An important component of the SSME system is the fuel preburner, which consists of a cylindrical chamber with a plate containing 264 coaxial injector elements at one end. A fuel rich mixture of gaseous hydrogen and liquid oxygen is injected and combusted in the chamber. This process preheats the hydrogen fuel before it enters the main combustion chamber, powers the hydrogen turbo-pump, and provides a heat dump for nozzle cooling. Issues of interest include the temperature and pressure fields at the turbine inlet and the thermal compatibility between the preburner chamber and injector plate. Performance anomalies can occur due to incomplete combustion, blocked injector ports, etc. The performance model should include the capability to simulate the effects of these anomalies. The current approach to the numerical simulation of the SSME fuel preburner flow field is to use a global model based on the MSFC sponsored FNDS code. This code does not have the capabilities of modeling several aspects of the problem such as detailed modeling of the coaxial injectors. Therefore, an effort has been initiated to develop a detailed simulation of the preburner coaxial injectors and provide gas phase boundary conditions just downstream of the injector face as input to the FDNS code. This simulation should include three-dimensional geometric effects such as proximity of injectors to baffles and chamber walls and interaction between injectors. This report describes an investigation into the numerical simulation of GH2/LOX coaxial injectors. The following sections will discuss the physical aspects of injectors, the CFD code employed, and preliminary results of a simulation of a single coaxial injector for which experimental data is available. It is hoped that this work will lay the foundation for the development of a unique and useful tool to support the SSME program.

Landrum, D. Brian↗
DOE OSTI2020

Field-emission plasma enhancement of H₂–O₂ micro-combustion

In this paper, a detailed three-dimensional computational model is developed to assess the potential of micron-scale field emission dielectric barrier discharge (FE-DBD) plasma actuators in non-premixed microburners with channel heights in the range of L z = 0.25–0.75 mm. Results for H 2 –O 2 combustion with Reynolds number Re = 100 are studied for different configurations of plasma actuator arrays. Without plasma actuation, results agree with corresponding experiments, although in this case the burner typically suffers from incomplete combustion due to thermal quenching, radical quenching, and slow diffusive mixing. For the range of conditions simulated, plasma actuation increases combustion completeness from 25%–70% to 40%–80% by enhancing mixing, generating radicals, and Joule heating. Ionic wind resulting from the FE-DBD plasma disrupts the otherwise simple laminar flow in the burner, accelerating the growth of a mixing layer and enhancing heat transfer to channel walls. The larger, lower temperature flames reduce thermal gradients (and thereby thermal stresses) in the walls. Radical generation due to electron-impact dissociation extends the reaction zone in the diffusion flame and increases combustion completeness. In all cases, the additional heat release with plasma actuation (100–190 W) is significantly larger than the power needed to sustain the plasma (≲0.1 W), suggesting that plasma actuation may provide a significant advantage in energy generation using realistic portable micro-combustion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Experimental and Simulation Analysis of Binary Mixtures of Biomass and Inert Material

A literature review was conducted to identify experiments on various scales of biomass combustion specific to circulating fluidized bed (CFB) combustors. Typical feedstock and inert bed material characteristics for biomass combustion in a CFB were gathered. This information was used to design laboratory-scale, non-reacting experiments on biomass-inert bed material mixtures in a 1-in. cylindrical test rig in the Multiphase Flow Analysis Laboratory (MFAL) in Morgantown, WV. The biomass chosen was a hardwood. Pellets from the hardwood were milled and sieved to obtain three size ranges—small, medium, and large. Three mass loadings of each of the biomass samples were fluidized with glass beads as the inert bed material to study the mixing characteristics under varying fluidizing gas velocities. Hydrodynamic data including pressure drop measurements at select vertical locations and video recordings of the mixture behavior were generated in this effort.

09 BIOMASS FUELS↗
DOE OSTI2021

Fuel Property Effects on Knock Propensity and Thermal Efficiency in a Direct-Injection Spark-Ignition Engine

Engine knock remains one of the major barriers to further improvement in thermal efficiency of Direct-Injection Spark-Ignition (DISI) engines. While Research Octane Number and Motor Octane Number are often used as standard rating methods for knock resistance of fuels, the impacts of other fuel properties on knock propensity in modern engines such as heat of vaporization (HoV) and laminar flame speed (LFS) require better understanding in order to co-optimize fuels and engine designs to achieve higher thermal efficiency and lower CO2 emission. In the present study, computational fluid dynamics (CFD) is used to model a boosted DISI engine with a focus on knock prediction and fuel property effects. A level-set G-equation model is employed to capture turbulent premixed combustion, and is coupled with a transported Livengood-Wu (L-W) integral approach to predict autoignition in the unburnt region. A criterion associated with the L-W integral is developed to accurately predict knock onset and knock-limited spark-advance. This model is then applied to a sensitivity analysis of HoV and LFS on knock tendency and thermal efficiency. The pressure-temperature trajectory framework is applied and extended to study the fuel effects on auto-ignition process in the engine. An existing efficiency-based merit function, which is derived from experiments for boosted SI engines, is evaluated and improved based on the current CFD results.

DISI↗
DOE OSTI2021

An Experimental and Numerical Investigation to Characterize the Low-Temperature Heat Release in Stoichiometric and Lean Combustion

This work reports on an experimental and modeling study on the low-temperature heat release (LTHR) characteristics for three RON 90 binary blends (n-heptane blended with isooctane, toluene and ethanol) in a Cooperative Fuel Research (CFR) engine at lean and stoichiometric conditions that are representative of homogeneous charge compression ignition (HCCI) and spark-ignition (SI) end-gas combustion conditions, respectively. An analysis of the end-gas temperature-pressure (T-P) trajectories was performed to identify the intake conditions leading to similar T-P trajectories between the two lambdas for each fuel blend. A heat release analysis was then conducted for the identified cases, where fuel-to-fuel differences in LTHR were identified and found to be sensitive to the operating condition. Simulations were conducted for these cases using a recently updated chemical kinetic model and a 0-D engine model, where good qualitative and reasonable quantitative agreements in LTHR were obtained. Sensitivity analysis was also performed directly on the rates of LTHR, to understand the controlling chemical reactions of LTHR, providing further insights into the fuel-to-fuel differences. The results demonstrate the significant promoting effect of n-heptane on LTHR rates, while inhibiting effects were seen for ethanol and toluene. Also highlighted was the importance of H atom abstraction reactions from the chemistry of each fuel component, which could lead to contradictory fuel behavior depending on the locations of the H site of the abstraction reaction due to the different ensuing pathways for the primary fuel radicals.

LTHR, RON 90 binary blends, sensitivity analysis↗
NASA NTRS1984

Turbulent transport models for scramjet flowfields

Turbulence modeling approaches were examined from the standpoint of their capability to predict the complex flowfield features observed in scramjet combustions. Thus, for example, the accuracy of each turbulence model, with respect to the prediction of recirculating flows, was examined. It was observed that for large diameter ratio axisymmetric sudden expansion flows, a choice of turbulence model was not critical because of the domination of their flowfields by pressure forces. For low diameter ratio axisymmetric sudden expansions and planar backward-facing steps flows, where turbulent shear stresses are of greater significance, the algebraic Reynolds stress approach, modified to increase its sensitivity to streamline curvature, was found to provide the best results. Results of the study also showed that strongly swirling flows provide a stringent test of turbulence model assumptions. Thus, although flows with very high swirl are not of great practical interest, they are useful for turbulence model development. Finally, it was also noted that numerical flowfields solution techniques have a strong interrelation with turbulence models, particularly with the turbulent transport models which involve source-dominated transport equations.

Sindir, M. M.↗
NASA NTRS1995

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.↗
NASA NTRS1984

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.↗
DOE OSTI2022

Closures for multicomponent reacting flows based on dispersion analysis

This work presents algebraic closure models associated with advective transport and nonlinear reactions in a Reynolds-averaged Navier-Stokes context for a system of species subject to binary reactions and transport by advection and diffusion. Expanding upon analysis originally developed for non-reactive transport in the context of Taylor dispersion of scalars, this work extends the modified gradient diffusion model explicated by Peters [N. Peters, Turbulent Combustion, Cambridge Monographs on Mechanics (Cambridge University Press, Cambridge, 2000)] and based on work by Corrsin [S. Corrsin, The reactant concentration spectrum in turbulent mixing with a first-order reaction, J. Fluid Mech. 11, 407 (1961)] beyond single-component transport phenomena and involving nonlinear reactions. The presented model forms, from this weakly nonlinear extension of the original dispersion theory, lead to an analytic expression for the eddy diffusivity matrix that explicitly captures the influence of the reaction kinetics on the closure operators. Furthermore, we demonstrate that the derived model form directly translates between flow topologies through a priori and a posteriori testing of a binary species system subject to homogeneous isotropic turbulence. Using two- and three-dimensional direct numerical simulations involving laminar and turbulent flows, it is shown that this framework improves prediction of mean quantities compared to previous results. Lastly, the presented model form, collapses to the earlier gradient diffusion and its modified version derived by Corrsin in the limits of nonreactive species and linear reactions, respectively.

42 ENGINEERING↗
DOE OSTI2021

Modeling Nanosecond-Pulsed Spark Discharge and Flame Kernel Evolution

Dilute combustion, either using exhaust gas recirculation or with excess air, is considered a promising strategy to improve the thermal efficiency of internal combustion engines. However, the dilute air-fuel mixture, especially under intensified turbulence and high-pressure conditions, poses significant challenges for ignitability and combustion stability, which may limit the attainable efficiency benefits. In-depth knowledge of the flame kernel evolution to stabilize ignition and combustion in a challenging environment is crucial for effective engine development and optimization. To date, a comprehensive understanding of ignition processes that result in the development of fully predictive ignition models usable by the automotive industry does not yet exist. Spark-ignition consists of a wide range of physics that includes electrical discharge, plasma evolution, joule-heating of gas, and flame kernel initiation and growth into a self-sustainable flame. In this study, an advanced approach is proposed to model spark-ignition energy deposition and flame kernel growth. To decouple the flame kernel growth from the electrical discharge, a nanosecond-pulsed high-voltage discharge is used to trigger spark-ignition in an optically accessible small ignition test vessel with a quiescent mixture of air and methane. Initial conditions for the flame kernel, including its thermodynamic state and species composition, are derived from a plasma-chemical equilibrium calculation. The geometric shape and dimension of the kernel are characterized using a multi-dimensional thermal plasma solver. Here, the proposed modeling approach is evaluated using a high-fidelity computational fluid dynamics procedure to compare the simulated flame kernel evolution against flame boundaries from companion Schlieren images.

42 ENGINEERING↗
NASA NTRS1990

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.↗
NASA NTRS1998

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.↗
NASA NTRS2021

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↗
NASA NTRS1979

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.↗
DOE OSTI2022

Exploring the potential of hydrogen in decarbonizing China's light-duty vehicle market

The Chinese government has pledged to achieve overall carbon neutrality by 2060. Currently, the transportation sector contributes to about 10% of total greenhouse gas (GHG) emissions in China. Hence, China has created a well-defined energy vehicle development strategy to reduce GHG emissions from the transportation sector, further expanding into hydrogen vehicle technologies. In this study, the Transportation Energy Analysis Model (TEAM) investigates the potential of hydrogen internal combustion engine vehicles (H2- ICEVs) and fuel cell vehicles (FCEVs) as a reliable pathway towards the government's aspiration of carbon neutrality in the transportation sector. According to TEAM, by adopting FCEVs and H2-ICEVs in the vehicle market, hydrogen demand could reach 25% of the total light-duty transportation energy demand in 2050. Consequently, this will lead to an annual reduction of more than 35 million tons GHG compared to only counting on the electrification pathway in the decarbonization task. Besides, FCEVs would take longer to penetrate the light-duty vehicle market compared to H2-ICEVs, as the current fuel cell technology still requires much improvement to attain a competitive vehicle cost of production.

Life cycle analysis↗
DOE OSTI2023

Simulation of a rapid compression machine for evaluation of ignition chemistry and soot formation using gasoline/ethanol blends

Due to the projected decline of demand for gasoline in light duty engines and the advent of ethanol as a green fuel, the use of gasoline/ethanol blend fuels in heavy duty applications are being investigated as they are projected to have lower cost and lower lifecycle green house gas (GHG) emissions. In heavy duty engines, the primary mode of combustion is mixing controlled combustion where wide range of mixture conditions (equivalence ratio) exist. Soot emissions of these fuels in richer conditions are not well understood. The goal of this research is to evaluate some commercially available soot modeling codes for the particulate matter emissions from gasoline/ethanol fuel blends, especially at fuel rich conditions. A Rapid Compression Machine (RCM) is modeled in a three-dimensional numerical simulation using CONVERGE computational software using a reduced chemical kinetic mechanism with SAGE chemistry solver and a RANS k-ϵturbulence model with a sector model including the creviced piston. The creviced piston is used in the experimental setup to reduce boundary layer effects and to maintain a homogeneous core in the reaction cylinder. Computational fluid dynamics simulations are conducted for different gasoline-ethanol fuel blends from E10 (10% ethanol v/v) to E100. The fuel blend is modeled as a surrogate mixture of toluene, iso-octane, n-heptane for gasoline content, and ethanol. The computational results were validated against experimental results using pressure measurements and laser extinction diagnostics. Different soot models are investigated to evaluate their capability of predicting the sooting tendencies of fuel blends, especially in richer conditions experienced during mixing-controlled combustion. The experimental combustion characteristics such as the ignition delay of different blends of fuel are reasonably well predicted. The Particulate Size Mimic (PSM) model accurately predicts the soot generation characteristics of the different fuels, but the Hiroyasu-NSC model falls short in this regard. For accurate prediction of soot with the PSM model, the thermodynamic conditions during combustion must be accurately modeled. While the current computational modeling tools can produce accurate results for the prediction of particulate matter emissions, there is much work to be done in improving our understanding of the underlying fundamental processes.

Energy & Fuels↗
DOE OSTI2021

The impact of the third O 2 addition reaction network on ignition delay times of neo-pentane

In this work, we studied the oxidation of neo-pentane by combining experiments, theoretical calculations, and mechanistic developments to elucidate the impact of the 3rd O 2 addition reaction network on ignition delay time predictions. The experiments are based on photoionization mass spectrometry in jet-stirred and time-resolved flow reactors allowing for sensitive detection of the keto-hydroperoxide (KHP) and keto-dihydroperoxide (KDHP) intermediates. With neo-pentane exhibiting a unique symmetric molecular structure, which consequently results only in single KHP and KDHP isomers, theoretical calculations of ionization and fragment appearance energies and of absolute photoionization cross sections enabled the unambiguous identification and quantification of the KHP intermediate. Its temperature and time-resolved profiles together with calculated and experimentally observed KHP-to-KDHP signal ratios were compared to simulation results based on a newly developed mechanism that describes the 3rd O 2 addition reaction network. A satisfactory agreement has been observed between the experimental data points and the simulation results, thus adding confidence to the model's overall performance. Finally, this mechanism was used to predict ignition delay times reported previously in shock tube and rapid compression machine experiments (J. Bugler et al., Combust. Flame 163 (2016) 138–156). While the model accurately reproduces the experimental data, simulations with and without the 3rd O 2 addition reaction network included reveal only a negligible effect on the predicted ignition delay times at 10 and 20 atm. According to model calculations, low temperatures and high pressures promote the importance of the 3rd O2 addition reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗