Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Combustion simulation”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Modeling–Experiment–Theory Analysis of Reactions Initiated from Cl + Methyl Formate

Methyl formate (MF; CH 3 OCHO) is the smallest representative of esters, which are common components of biodiesel. The present study characterizes the thermal dissociation kinetics of the radicals formed by H atom abstraction from MF—CH 3 OCO and CH 2 OCHO—through a combination of modeling, experiment, and theory. For the experimental effort, excimer laser photolysis of Cl 2 was used as a source of Cl atoms to initiate reactions with MF in the gas phase. Time-resolved species profiles of MF, Cl 2 , HCl, CO 2 , CH 3 , CH 3 Cl, CH 2 O, and CH 2 ClOCHO were measured and quantified using photoionization mass spectrometry at temperatures of 400–750 K and 10 Torr. The experimental data were simulated using a kinetic model, which was informed by ab initio-based theoretical kinetics calculations and included chlorine chemistry and secondary reactions of radical decomposition products. Here, we calculated the rate coefficients for the H-abstraction reactions Cl + MF → HCl + CH 3 OCO (R1a) and Cl + MF → HCl + CH 2 OCHO (R1b): k 1a,theory = 6.71 × 10 –15 ·T 1.14 ·exp(—606/T) cm 3 /molecule·s; k 1b,theory = 4.67 × 10 –18 ·T 2.21 ·exp(—245/T) cm 3 /molecule·s over T = 200–2000 K. Electronic structure calculations indicate that the barriers to CH 3 OCO and CH 2 OCHO dissociation are 13.7 and 31.6 kcal/mol and lead to CH 3 + CO 2 (R3) and CH 2 O + HCO (R5), respectively. The master equation-based theoretical rate coefficients are k 3,theory (P = ∞) = 2.94 × 10 9 ·T 1.21 ·exp(—6209/T) s –1 and k 5,theory (P = ∞) = 8.45 × 10 8 ·T 1.39 ·exp(—15132/T) s –1 over T = 300–1500 K. The calculated branching fractions into R1a and R1b and the rate coefficient for R5 were validated by modeling of the experimental species time profiles and found to be in excellent agreement with theory. Additionally, we found that the bimolecular reactions CH 2 OCHO + Cl, CH 2 OCHO + Cl 2 , and CH 3 + Cl 2 were critical to accurately model the experimental data and constrain the kinetics of MF-radicals. Inclusion of the kinetic parameters determined in this study showed a significant impact on combustion simulations of larger methyl esters, which are considered as biodiesel surrogates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

The role of energy transfer and competing bimolecular reactions in characterizing the unimolecular dissociations of allylic radicals

Resonantly stabilized radicals (RSR's) such as allylic radicals typically have bond energies exceeding 45 kcal/mol and are therefore reasonably long-lived up to quite high temperatures. Consequently, kinetics of the unimolecular dissociations and bimolecular reactions of these long-lived RSR's assume significance in high temperature chemically reacting systems such as in flames. Our recent analyses (J. Cho et al., Proc. Comb. Inst., 2022) of the uncertainty propagated by the energy transfer parameters for 1-methylallyl (1MA) dissociation to the laminar flame speed, indicates an intricate coupling between the kinetics of 1MA dissociation (chain propagation) and its reaction with H-atoms (chain-termination). Here, this work extends such analyses to other C 3 – C 5 allylic radicals: allyl, 2-methylallyl, 1,1-dimethylallyl, 1,2-dimethylallyl, and 1,3-dimethylallyl. Theoretical kinetics for the C 3 – C 4 allylic radicals were taken from prior literature studies (J.A. Miller et al., J. Phys. Chem., 2008; J. Cho et al., Proc. Comb. Inst., 2022; R.S. Tranter et al., Proc. Comb. Inst., 2017). Potential energy surfaces were newly computed for the dissociations of the larger C5 allylic radicals and the reactions of these radicals with H-atoms. These PES's were then used to calculate rate constants for the relevant unimolecular and bimolecular reactions. Energy transfer has a crucial role to play in these pressure-dependent reactions and therefore the present kinetics calculations relied on a-priori predictions of collisional energy transfer parameters to best characterize the branching between chain propagating unimolecular dissociations of these allylic radicals and chain terminating bimolecular reactions of these RSR's with H-atoms. Prompt dissociations are now established to be a universal feature of free radical dissociations. Here we have also characterized the propensity for prompt dissociations in these more-stable allylic RSR's and their role in combustion simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Understanding and Estimating Error Propagation in Neural Networks for Scientific Data Analysis

Neural networks are increasingly integrated into scientific discovery, where input data reduction and model quantization play a key role in accelerating inference. However, understanding and mitigating the impact of these techniques on output error is critical for ensuring reliable results, particularly in tasks demanding high numerical precision. This paper introduces a comprehensive framework for optimizing neural network inference in scientific computing by combining data reduction and weight quantization while maintaining error-controlled outcomes. We develop theoretical analyses to bound error propagation under these reductions and propose a framework that balances computational performance with error constraints. Evaluation on real-world learning-based combustion simulations and satellite image classification demonstrates that our derived error bounds accurately predict observed errors while enabling significant computational speedup under our framework. This work highlights the potential for further leveraging advancements in modern lossy compression algorithms and hardware accelerators that support lower-precision formats.

He, Weiming [New Jersey Institute of Technology]↗
DOE OSTI2021

Comparison of Radiation Models for a Turbulent Piloted Methane/Air Jet Flame: A Frozen-Field Study

Numerical modeling of radiative transfer in nongray reacting media is a challenging problem in computational science and engineering. The choice of radiation models is important for accurate and efficient high-fidelity combustion simulations. Different applications usually involve different degrees of complexity, so there is yet no consensus in the community. In this paper, the performance of different radiative transfer equation (RTE) solvers and spectral models for a turbulent piloted methane/air jet flame are studied. The flame is scaled from the Sandia Flame D with a Reynolds number of 22,400. Three classes of RTE solvers, namely the discrete ordinates method, spherical harmonics method, and Monte Carlo method, are examined. The spectral models include the Planck-mean model, the full-spectrum k-distribution (FSK) method, and the line-by-line (LBL) calculation. The performances of different radiation models in terms of accuracy and computational cost are benchmarked. The results have shown that both RTE solvers and spectral models are critical in the prediction of radiative heat source terms for this jet flame. The trade-offs between the accuracy, the computational cost, and the implementation difficulty are discussed in detail. The results can be used as a reference for radiation model selection in combustor simulations.

David, Chloe↗
DOE OSTI2021

In Situ Machine Learning for Intelligent Data Capture on Exascale Platforms. Final Report

In many dynamic systems, interesting events occur locally in time and space. Examples of such systems include ignition events in combustion simulations, material fractures in mechanics simulations, and extreme weather events in climate simulations. Due to memory constraints and data I/O costs, current simulation workflows save data at regularly spaced time-steps, at a fixed rate determined before the start of the simulation. Often this mode of operation results in missed events of interest, necessitating a simulation restart from before an event occurred with more frequent data saves. This data saving workflow is grossly inefficient and is already a bottleneck in the computing process. We propose to develop machine learning algorithms that can detect when interesting dynamical events are occurring, triggering data saves. These machine learning algorithms will perform in situ anomaly detection to flag regions with different dynamical properties than those previously recorded. The adaptive data saves would be local in time and space to match the event of interest, thereby enabling a much more efficient workflow that will reduce data I/O costs and data storage memory requirements. The algorithms will be tested on two applications: auto-ignition simulations and climate simulations. A critical component of this project will be developing machine learning algorithms that can be deployed efficiently in situ on HPC platforms with out-of-the-box functionality. The development of in situ machine learning methods to detect anomalous events would enable a more efficient and effective workflow, in which all the relevant data are saved in a single simulation run, without re-starts or scientist intervention.

42 ENGINEERING↗
DOE OSTI2022

ExaLearn – GenTen Tensor Software ECP Milestone

The objective of this milestone was to finish integrating GenTen tensor software with combustion application Pele using the Ascent in situ analysis software, partnering with the ALPINE and Pele teams. Also, to demonstrate the usage of the tensor analysis as part of a combustion simulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Design and Evaluation of Environmental Barrier Coatings for Protection of Ceramic Matrix Composites in Hydrogen-Based Turbines

This presentation at the 2024 FECM Spring Project Review Meeting summarizes ongoing work at NETL aimed at enabling the use of ceramic matrix composite (CMC) materials in next-generation high-efficiency power generation turbines fueled by hydrogen, natural gas, or other fuels. The work broadly involves three aspects (1) phase-field damage modeling of CMCs and the environmental barrier coatings (EBCs) used to protect them (2) computational design of new and improved EBCs and (3) performance testing of EBCs using a novel gas turbine combustion simulation rig developed at NETL. An overview on current status and next steps for each aspect of the project are presented.

ceramic matrix composite (CMC)↗
NASA NTRS1972

The simulation of high pressure hydrogen/oxygen rocket engines.

Initial experiments have established the feasibility of using a backward running detonation wave in a shock tube to produce the correct enthalpy and chemistry of the simulated combustion chamber gases. This technique produces plume Reynolds numbers higher than those achievable by any other known technique. The practical pressure limitations are being explored and detailed surveys of both single and multiple engine plumes and impingement data on a cylindrical body are being obtained. Facility details to allow safe handling of the high pressure hydrogen and oxygen up to a design chamber pressure of 20,000 psi are described.

Hopkins, H.↗
NASA NTRS1992

Operating characteristics at Mach 4 of an inlet having forward-swept, sidewall-compression surfaces

A three-dimensional, forward-swept sidewall-compression inlet has been tested in the Mach 4 Blowdown Facility at the NASA Langley Research Center to examine parameters that affect inlet starting, operability, and performance. The inlet was designed to simultaneously provide good starting and mass capture characteristics through the combination of forward-swept sidewalls and the aft placement of the downstream cowl. Parametrics examined included top surface compression, geometric contraction ratio, sidewall compression angle, boundary-layer thickness on the vehicle undersurface, and cowl position. To simulate combustion effects on inlet operation, the inlet was mechanically back pressured to determine maximum achievable combustor-to-inlet pressure ratio prior to inlet unstart. Inlet operability and performance are shown in terms of inlet starting (pulsed and self-starting) and unstarting characteristics, inlet static pressure distributions, combustor back pressure limits, mass capture, and flow uniformity.

Hudgens, Julie A.↗
NASA NTRS1993

An experimental study of the effects of bodyside compression on forward swept sidewall compression inlets ingesting a turbulent boundary layer

Forward swept sidewall compression inlets have been tested in the Mach 4 Blowdown Facility at the NASA Langley Research Center to study the effects of bodyside compression surfaces on inlet performance in the presence of an incoming turbulent boundary layer. The measurements include mass flow capture and mean surface pressure distributions obtained during simulated combustion pressure increases downstream of the inlet. The kerosene-lampblack surface tracer technique has been used to obtain patterns of the local wall shear stress direction. Inlet performance is evaluated using starting and unstarting characteristics, mass capture, mean surface pressure distributions and permissible back pressure limits. The results indicate that inlet performance can be improved with selected bodyside compression surfaces placed between the inlet sidewalls.

Rodi, Patrick E.↗
NASA NTRS1959

Investigation on the use of a Freely Rotating Rotor at the Cowl Face of a Supersonic Conical Inlet to Reduce Inlet Flow Distortion

An investigation has been made on the use of a freely rotating rotor at the cowl face of a supersonic conical diffuser to determine its effectiveness in reducing inlet flow distortion and the penalty in terms of total-pressure loss imposed by such a device when distortions are negligible. Tests were made with a rotor having an inlet tip diameter of 2.18 inches and a ratio of hub radius to tip radius of 0.52, in conjunction with a conical inlet having a 25 deg semi-vertex cone angle, at a Mach number of 2.1 over an angle-of-attack range of 0 deg to 8 deg. A simplified analysis showing that a supersonic, freely rotating rotor with maximum solidity for noninterference between blades will operate in an undistorted flow with a total-pressure defect of 1 percent or less was experimentally verified. Overall total-pressure distortions of 0.1 to 0.4 and Mach number distortions of 0.4 to 1.4, obtained at 4 deg to 8 deg angle of attack, were reduced about 30 percent and 23 percent, respectively, because of the presence of the rotor, with no measurable total-pressure loss. The rotor increased the peak total-pressure recovery at the simulated combustion chamber 1 1/2 and 3 1/2 percent at 6 deg and 8 deg angles of attack, respectively. This increase is attributed to lower diffusion duct losses as a consequence of a more uniform flow created by the rotor.

Goldberg, Theodore J.↗
NASA NTRS2000

National Combustion Code: Parallel Implementation and Performance

The National Combustion Code (NCC) is being developed by an industry-government team for the design and analysis of combustion systems. CORSAIR-CCD is the current baseline reacting flow solver for NCC. This is a parallel, unstructured grid code which uses a distributed memory, message passing model for its parallel implementation. The focus of the present effort has been to improve the performance of the NCC flow solver to meet combustor designer requirements for model accuracy and analysis turnaround time. Improving the performance of this code contributes significantly to the overall reduction in time and cost of the combustor design cycle. This paper describes the parallel implementation of the NCC flow solver and summarizes its current parallel performance on an SGI Origin 2000. Earlier parallel performance results on an IBM SP-2 are also included. The performance improvements which have enabled a turnaround of less than 15 hours for a 1.3 million element fully reacting combustion simulation are described.

Quealy, A.↗
NASA NTRS2002

Hybrid Propulsion Testing at Marshall Space Flight Center

Hybrid propulsion testing involving eleven and twenty-four inch motors performed at the Marshall Space Flight Center (MSFC) from the early 1990's to the present are discussed. Topics covered include: Solid Propulsion Investigation Program, Joint NASA Industry Research and Development (JIRAD) program, Large Subscale Solid Rocket Combustion Simulator (LSSRCS), Hybrid Propulsion Demonstration Program (HPDP), Hybrid Propulsion for Launch Vehicle Booster (HPTLV), Peroxide Hybrid Upper Stage (PHUS) and Solid Fuel Torch (SFT).

Prince, Andrew S.↗
NASA NTRS2002

National Combustion Code Parallel Performance Enhancements

The National Combustion Code (NCC) is being developed by an industry-government team for the design and analysis of combustion systems. The unstructured grid, reacting flow code uses a distributed memory, message passing model for its parallel implementation. The focus of the present effort has been to improve the performance of the NCC code to meet combustor designer requirements for model accuracy and analysis turnaround time. Improving the performance of this code contributes significantly to the overall reduction in time and cost of the combustor design cycle. This report describes recent parallel processing modifications to NCC that have improved the parallel scalability of the code, enabling a two hour turnaround for a 1.3 million element fully reacting combustion simulation on an SGI Origin 2000.

Quealy, Angela↗
NASA NTRS1995

Affect of Brush Seals on Wave Rotor Performance Assessed

The NASA Lewis Research Center's experimental and theoretical research shows that wave rotor topping can significantly enhance gas turbine engine performance levels. Engine-specific fuel consumption and specific power are potentially enhanced by 15 and 20 percent, respectively, in small (e.g., 400 to 700 hp) and intermediate (e.g., 3000 to 5000 hp) turboshaft engines. Furthermore, there is potential for a 3- to 6-percent specific fuel consumption enhancement in large (e.g., 80,000 to 100,000 lbf) turbofan engines. This wave-rotor-enhanced engine performance is accomplished within current material-limited temperature constraints. The completed first phase of experimental testing involved a three-port wave rotor cycle in which medium total pressure inlet air was divided into two outlet streams, one of higher total pressure and one of lower total pressure. The experiment successfully provided the data needed to characterize viscous, partial admission, and leakage loss mechanisms. Statistical analysis indicated that wave rotor product efficiency decreases linearly with the rotor to end-wall gap, the square of the friction factor, and the square of the passage of nondimensional opening time. Brush seals were installed to further minimize rotor passage-to-cavity leakage. The graph shows the effect of brush seals on wave rotor product efficiency. For the second-phase experiment, which involves a four-port wave rotor cycle in which heat is added to the Brayton cycle in an external burner, a one-dimensional design/analysis code is used in conjunction with a wave rotor performance optimization scheme and a two-dimensional Navier-Stokes code. The purpose of the four-port experiment is to demonstrate and validate the numerically predicted four-port pressure ratio versus temperature ratio at pressures and temperatures lower than those that would be encountered in a future wave rotor/demonstrator engine test. Lewis and the Allison Engine Company are collaborating to investigate wave rotor integration in an existing turboshaft engine. Recent theoretical efforts include simulating wave rotor dynamics (e.g., startup and load-change transient analysis), modifying the one-dimensional wave rotor code to simulate combustion internal to the wave rotor, and developing an analytical wave rotor design/analysis tool based on macroscopic balances for parametric wave rotor/engine analysis.

Source record↗
NASA NTRS2012

Environmental Durability and Stress Rupture of EBC/CMCs

This research focuses on the strength and creep performance of SiC fiber-reinforced SiC ceramic matrix composite (CMC) environmental barrier coating (EBC) systems under complex simulated engine environments. Tensile-strength and stress-rupture testing was conducted to illustrate the material properties under isothermal and thermal gradient conditions. To determine material durability, further testing was conducted under exposure to thermal cycling, thermal gradients and simulated combustion environments. Emphasis is placed on experimental techniques as well as implementation of non-destructive evaluation, including modal acoustic emission and electrical resistivity monitoring, to characterize strength degradation and damage mechanisms. Currently, little is known about the behavior of EBC-CMCs under these conditions; consequently, this work will prove invaluable in the development of structural components for use in high temperature applications.

Appleby, Matthew↗
NASA NTRS2021

High Temperature Solid Particle Erosion in a Melt-Infiltrated SiC/SiC Ceramic Matrix Composite

Ceramic matrix composites (CMCs) are an enabling propulsion material system that offer weight benefits over current Ni-based superalloys, and have higher temperature capabilities that can reduce cooling requirements. Incorporating CMCs into the hot section of gas-turbine engines therefore leads to an increase in engine efficiency. While significant advancements have been made, challenges still remain for current and next-generation gas-turbines; particularly when operating in dust-laden or erosive environments. Solid particles entrained in the gas flow can impact engine hardware resulting in localized damage and material removal due to repeated, cumulative impacts. In this study, the erosion behavior of a melt-infiltrated (MI) silicon carbide fiber-reinforced silicon carbide (SiC/SiC) CMC is investigated at high temperature (1,200 °C) in a simulated combustion environment using 150 µm alumina particles as erodent. Particle impact velocities ranged from 100 to 200 m/s and the angle of impingement varied from 30° to 90°. Erosion testing was also performed on α-SiC to elucidate similarities and differences in the erosion response of the composite compared to that of a monolithic ceramic. Scanning electron microscopy (SEM) was used to study the post-erosion damage morphology and the governing mechanisms of material removal.

Ceramic matrix composites (CMCs), solid particle e↗