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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 433 records · Page 24

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Probabilistic constrained Bayesian inversion for transpiration cooling

Abstract To enable safe operations in applications such as rocket combustion chambers, the materials require cooling to avoid material damage. Here, transpiration cooling is a promising cooling technique. Numerous studies investigate possibilities to simulate and evaluate the complex cooling mechanism. One naturally arising question is the amount of coolant required to ensure a safe operation. To study this, we introduce an approach that determines the posterior probability distribution of the Reynolds number using an inverse problem and constraining the maximum temperature of the system under parameter uncertainties. Mathematically, this chance inequality constraint is dealt with by a generalized polynomial chaos expansion of the system. The posterior distribution will be evaluated by different Markov chain Monte Carlo based methods. A novel method for the constrained case is proposed and tested among others on two‐dimensional transpiration cooling models.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Multi-dimensional modeling of mixture preparation in a direct injection engine fueled with gaseous hydrogen

With the recent advances of direct injection (DI) technology, introducing hydrogen into the combustion chamber through DI is being considered as a viable approach to circumvent backfire and pre-ignition encountered in early generations of hydrogen engines. As part of a broader vision to develop a robust numerical model to study hydrogen spark ignition (SI) combustion in internal combustion (IC) engines, the present numerical investigation focuses on mixture preparation in a hydrogen DI SI engine. This study is carried out with a single hole injector with gaseous hydrogen injected at 100 bar injection pressure. Simulations are carried out for high and low tumble configurations and validated against optical data acquired from planar laser induced fluorescence (PLIF) measurements. Varying mesh configurations are investigated for the impact on in-cylinder mixture distribution. A particular emphasis is placed on the effect of nozzle geometry and mesh orientation near the wall. Overall, the computational model is found to predict the mixture distribution in the combustion cylinder reasonably well. The results showed that the alignment of mesh with the flow direction is important to achieve good agreement between numerical analysis and optical measurement data.

08 HYDROGEN↗

DOE OSTI2026

Irradiation Driven Restructuring of Nanocrystalline ThO 2 and Th 1–x U x O 2 Thin Films

Irradiation induced structural changes of actinide oxide materials is a key consideration in their development and use as nuclear fuels. This study reported on the synthesis of ThO 2 and Th 1–x U x O 2 (x = 0.15, 0.50) thin films, fabricated using electrospray-assisted solution combustion synthesis, and their responses to ion irradiation. Krypton ion irradiations, up to a fluence of 1 × 10 16 ions/cm 2 , were carried out to simulate radiation damage induced by fission products in a reactor environment. Structural and chemical changes induced by irradiation were analyzed using high-resolution scanning transmission electron microscopy (STEM), energy-dispersive X-ray spectroscopy (EDS), and electron energy-loss spectroscopy (EELS). It was determined that the extent and nature of irradiation-induced damage are strongly correlated with the uranium content. ThO 2 films were most susceptible to radiation-induced damage, with significant cavity formation and delamination from the substrate at high fluence. Of the compositions studied, Th 0.85 U 0.15 O 2 films showed the highest stability, characterized by moderate grain growth and the absence of voids or severe defect structures. In contrast, Th 0.5 U 0.5 O 2 films accumulated extensive damage, including the formation of a nanocrystalline central region. EELS analysis indicated that oxygen displacement is the primary driver of structural degradation in Th 0.5 U 0.5 O 2 films. α-particle spectroscopy confirmed minimal actinide loss across all compositions, underscoring the mechanical robustness of the films. These findings provide insight into the irradiation-induced damage mechanisms in Th O2 and Th 1–x U x O 2 systems, supporting their development as potential materials for nuclear fuels and irradiation-tolerant thin film targets in nuclear physics measurements.

Th1−xUxO2↗

DOE OSTI2021

Application of an automated machine learning-genetic algorithm (AutoML-GA) coupled with computational fluid dynamics simulations for rapid engine design optimization

In recent years, the use of machine learning-based surrogate models for computational fluid dynamics (CFD) simulations has emerged as a promising technique for reducing the computational cost associated with engine design optimization. However, such methods still suffer from drawbacks. One main disadvantage is that the default machine learning (ML) hyperparameters are often severely suboptimal for a given problem. This has often been addressed by manually trying out different hyperparameter settings, but this solution is ineffective in case of a high-dimensional hyperparameter space. Besides this problem, the amount of data needed for training is also not known a priori. In response to these issues that need to be addressed, the present work describes and validates an automated active learning approach, AutoML-GA, for surrogate-based optimization of internal combustion engines. In this approach, a Bayesian optimization technique is used to find the best machine learning hyperparameters based on an initial dataset obtained from a small number of CFD simulations. Subsequently, a genetic algorithm is employed to locate the design optimum on the ML surrogate surface. In the vicinity of the design optimum, the solution is refined by repeatedly running CFD simulations at the projected optima and adding the newly obtained data to the training dataset. It is demonstrated that AutoML-GA leads to a better optimum with a lower number of CFD simulations, compared to the use of default hyperparameters. The proposed framework offers the advantage of being a more hands-off approach that can be readily utilized by researchers and engineers in industry who do not have extensive machine learning expertise.

Owoyele, Opeoluwa↗

DOE OSTI2020

Benchmark of the KGMf with a coupled Boltzmann equation solver

The Kinetic Global Model framework (KGMf) is an open-source general-purpose global model (spatially averaged) simulation code developed to explore the reaction kinetics and pathways in plasma discharge systems. It contains species continuity and electron energy balance equations, with a time-dependent evaluated electron energy distribution function (EEDF) for electron impact reactions. The EEDF is utilized to determine the rate coefficients for electron impact reactions, which can have profound impact on the temporal evolutions of plasma parameters. Previously, the EEDF was commonly assumed as Maxwellian or an analytical function of the effective electron temperature. In this work, the KGMf is coupled with a Boltzmann equation (BE) solver to self-consistently compute the EEDF. The EEDF evolution frequency is determined based on relative changes of the reduced electric field. The KGMf is benchmarked with the ZDPlasKin code based on high-pressure low-temperature argon plasma discharge cases. Additionally, the temporal evolutions of reduced electric field, electron temperature, EEDF, reaction rates, and species densities, are obtained and compared under different discharge conditions, showing good agreement between the KGMf and the ZDPlasKin simulations. The application of the KGMf for predicting breakdown times in high power microwave discharges is also presented, which shows qualitative agreement with particle-in-cell simulations. The KGMf can be further applied for more complicated plasma discharge systems, where the reaction kinetics are more intricate (e.g., plasma-assisted combustion systems).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Atmospheric Radiative and Oceanic Biological Productivity Responses to Increasing Anthropogenic Combustion‐Iron Emission in the 1850–2010 Period

Abstract Anthropogenic emission is an important component of the present‐day iron cycle yet its long‐term impacts on climate are poorly understood. Iron mineralogy strongly affects its radiative and oceanic interactions and was unrepresented in previous studies. We perform simulations using a mineralogy‐based inventory and an atmospheric transport model and estimate the 1850–2010 global mean direct radiative forcing (DRF) to be +0.02 to +0.10 W/m 2 . We estimate that the CO 2 sequestration of 0.2–13 ppmv over the last 150 years due to enhanced phytoplankton productivity by anthropogenic iron deposition causes an avoided CO 2 forcing of −0.002 to −0.16 W/m 2 . While globally small, these impacts can be higher in specific regions; the anthropogenic DRF is +0.5 W/m 2 over areas with more coal combustion and metal smelting, and anthropogenic soluble iron sustains >10% of marine net primary productivity in the high‐latitude North Pacific Ocean, a region vulnerable to stratification due to climate change.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Effects of non-thermal termolecular reactions on detonation development in hydrogen (H 2 )/methane (CH 4 ) - air mixtures

The binary fuel blend of H 2 /CH 4 is one of the most promising hydrogen-enriched hydrocarbon fuels in spark-ignition (SI) engines. Yet, the undesirable phenomenon of super-knock, which can severely and in-stantaneously damage an SI engine, limits its widespread adoption. Moreover, there is still a lack of con-sensus on the precise mechanism by which this phenomenon occurs i.e. via flame acceleration or spon-taneous ignition, despite numerous previous investigations. At the same time, recent studies [M. P. Burke, S. J. Klippenstein, Nat. Chem. 9 (2017) 1078 -1082, Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515-522] have demonstrated a high probability of occurrence of non-thermal reactions in premixed flames of such H 2 /CH 4 fuel blends with air due to the presence of non-trivial amounts of highly reactive radicals including H, O and OH apart from O 2 . The present study focuses on the evolution of an initial deflagration front to a detonation wave in H 2 /CH 4 - air mixtures under SI engine relevant conditions through fully resolved, constant volume 1D simulations with and without non-thermal reactivity. Non-thermal reactions were included in the macroscopic kinetics model as chemically termolecular reactions facilitated by the H + CH 3 and H + OH radical-radical recombination and the H + O 2 radical-molecule association reactions. Further, the nonthermal reactions result in a corresponding decrease in the reaction fluxes of the incipient recombination/association reactions. Therefore, an additional set of simulations were performed by applying corrections to the respective incipient recombination/association rate constants using the methodology demonstrated by Tao et al. [Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515-522]. Compared to the baseline case, the onset of spontaneous ignition in the end-gas region was observed to be delayed in the presence of non-thermal termolecular reactions. Concurrently, the developing detonation was observed to be significantly stronger. In contrast, applying corrections to the recombi-nation/association rate constants resulted in a completely different behavior. Specifically, detonation was observed to occur due to self acceleration of the primary flame in the absence of spontaneous ignition in the end-gas region. Sensitivity analysis was performed to quantify the effects of non-thermal reactions on the duration of heat release rate and thereby the mechanism of detonation formation. In addition, chemical explosive mode analysis (CEMA) was performed to identify the dominant species/reactions re-sponsible for the observed results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Large eddy simulation of a supersonic lifted hydrogen flame: Impacts of Lewis, turbulent Schmidt and Prandtl numbers

Parametric large eddy simulations (LES) of a supersonic lifted hydrogen flame are reported. The emphases are on two aspects: impacts of (1) Lewis number (Lei of the ith species) and (2) turbulent Schmidt and Prandtl numbers (⁠Sc t and Pr t ⁠) on supersonic turbulent flame and flow structures. Five cases are considered: species-specific Le i , Sc t = Sr t = 1.0 (C 0 ); unity Le i , Sc t = Pr t = 1.0 (C 1 ); species-specific Sc t = 0.5, Pr t = 1.0 (C 2 ); species-specific Le i , Sc t = 1.0, Pr t = 0.5 (C 3 ); and species-specific Le i , Sc t = Pr t = 0.5 (C 4 ). Numerical results of instantaneous and/or time-averaged species mole fractions, mixture fraction, heat release rate, flame base location, and mixed modes of premixed and diffusion combustion are compared between cases C 0 and C 1 . Differences in auto-ignition locations and strengths and flame structures and stabilization specify the impacts of Lewis number. They are triggered by different predictions of species mass and thermal diffusions at fuel-coflow and/or coflow-ambient air mixing layers. These differences are rationalized by a scale analysis of mass/thermal diffusion and convection for case C 0 , which suggests the relatively low but non-negligible former against the latter. Cases C 0 and C 2 –C 4 barely see differences in terms of instantaneous and/or time-averaged temperature, velocity, and mixed combustion modes except for further downstream areas where combustion occurs. Both Sc t and Pr t impose less significant influences than Lewis number, as sub-grid scale mass/thermal diffusion is subordinate to its resolved counterpart according to their scale analysis for case C 4 .

42 ENGINEERING↗

DOE OSTI2021

Fundamental insights on ignition and combustion of natural gas in an active fueled pre-chamber spark-ignition system

Pre-chamber spark-ignition (PCSI), either fueled or non-fueled, is a leading concept with the potential to enable diesel-like efficiency in medium-duty (MD) and heavy-duty (HD) natural gas (NG) engines. However, the inadequate scientific base and simulation tools to describe/predict the underlying processes governing PCSI systems is one of the key barriers to market penetration of PCSI for MD/HD NG engines. To this end, experiments were performed in a heavy-duty, optical, single-cylinder engine fitted with an active fueled PCSI module. Furthermore, the spatial and temporal progress of ignition and subsequent combustion of lean-burn natural gas using PCSI system were studied using optical diagnostic imaging and heat release analysis based on main-chamber and pre-chamber pressure measurements.

42 ENGINEERING↗

DOE OSTI2025

A High-Order Discontinuous Galerkin Spectral Element Method for Compressible Reacting Flows

High-order methods have recently been shown to be an effective tool for high-fidelity flow computations like direct numerical simulations and large eddy simulations due to their strong balance between accuracy and computational cost. In this work, a high-order discontinuous Galerkin spectral element method (DGSEM) is developed to solve the chemically reactive Euler equations encountered in high-speed combustion. To handle the disparate length and time scales associated with these equations, we develop a novel method which combines the spectral accuracy of the SEM with the flexibility of DG approach. Thus, the framework is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. The numerical method is implemented within the spectral element solver Nek5000. Validation cases are conducted for both non-reactive and reactive discontinuous flows to demonstrate the solver capability. In particular, canonical one-dimensional and two-dimensional detonation simulations are performed and the high-order numerical results are validated against available literature data.

computational fluid dynamics (CFD)↗

DOE OSTI2020

Range cost-effectiveness of plug-in electric vehicle for heterogeneous consumers: An expanded total ownership cost approach

Plug-in electric vehicles (PEV) appears to have sales momentum in major personal vehicle markets but are still at the early market stage. Opportunities to accelerate PEV adoption can be discovered through comprehensive total cost of ownership (TCO) analysis. Understanding the cost-effective electric ranges of PEVs for consumers, manufacturers, and the society is critical for any discussion of PEV mass markets. This study expanded the traditional TCO approach by (1) fully considering heterogeneous consumer driving patterns, (2) quantifying the charging inconvenience and range anxiety cost of battery-electric vehicles (BEVs), and (3) monetizing both tangible and intangible PEV policies. Uncertainties were handled through Monte Carlo simulation. The results suggest that BEVs with an electric range of 250–350 km have the lowest TCO in cities with government-enacted purchase limitations, and internal combustion engine vehicles (ICEVs) have the lowest TCO in cities without purchase limitations, even when considering PEV subsidies. The lowest TCO for some consumer groups is obtained by BEVs with an electric range of 400–450 km, especially in northern China, where the weather is colder. The cost-effective all-electric range for BEVs in each city in 2025 will decrease due to improved battery performance in cold environments and an expanded charging infrastructure. Finally, based on TCO, plug-in hybrid electric vehicles (PHEVs) are currently more suitable for drivers with a high average daily mileage or a large mileage variance. However, by 2025, BEVs with a long driving range may become a more cost-effective choice for these drivers.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2021

Diastereomers and Low-Temperature Oxidation

Diastereomers have historically been ignored when building kinetic mechanisms for combustion. Low-temperature oxidation kinetics, which continues to gain interest in both combustion and atmospheric communities, may be affected by the inclusion of diastereomers in radical chain-branching pathways. In this work, key intermediates and transition states lacking stereochemical specification in an existing diethyl ether low-temperature oxidation mechanism were replaced with their diastereomeric counterparts. Rate coefficients for reactions involving diastereomers were computed with ab initio transition state theory master equation calculations. The presence of diastereomers increased rate coefficients by factors of 1.2-1.6 across various temperatures and pressures. Ignition delay simulations incorporating these revised rate coefficients indicate that the diastereomers enhanced the overall reactivity of the mechanism by almost 15% and increased the peak ketohydroperoxide concentration by 30% in the negative temperature coefficient region at combustion-relevant pressures. Furthermore, these results provide an illustrative indication of the important role of stereomeric effects in oxidation kinetics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

The Pele Simulation Suite for Reacting Flows at Exascale

In this work, we present the Pele suite of software tools for compressible and incompressible reacting flows. The Pele suite leverages several different libraries, notably AMReX and SUNDIALS, to achieve performance portability on heterogeneous computing architectures across the supercomputing landscape. The Pele suite is comprised of PeleC, a compressible reacting flow block-structured adaptive mesh refinement solver, PeleLMeX, a low-Mach number reacting flow block-structured adaptive mesh refinement solver, Pele-Physics, a library for transport, thermodynamics, finite rate chemistry, soot, spray and radiation physics. The objective of this paper is (i) to present the code development efforts necessary to achieve highly effective and scalable applications for exascale machines and (ii) to detail the performance results of the Combustion-Pele project applications on Oak Ridge National Laboratory's Frontier. We show good weak and strong scaling results for both PeleC and PeleLMeX up to more than 50 billion cells on more than 4096 Frontier graphics processing unit nodes. We also present a capability demonstration simulation of a dual-fuel pulse compression ignition engine (six adaptive mesh refinement levels, and 60 billion cells or 2.1 trillion degrees of freedom) on Frontier, to date one of the largest simulations performed on the first exascale-class supercomputer.

adaptive mesh refinement↗

DOE OSTI2020

Development of a Physics-Based Combustion Model for Engine Knock Prediction

The objective of this project is to improve the prediction of engine knock by developing a new combustion modeling framework. Engine knock is a limiting factor to constrain the increase of fuel efficiency for spark ignition (SI) engines in most passenger cars. Efforts to increase fuel efficiency, increasing the compression ratio or downsizing, lead to the increase in the tendency of the knock occurrence. The knock is an undesired ignition of the end-gas, unburned fuel/air mixture ahead of the spark-ignited premixed flame, resulting in rapid in-cylinder pressure rises and engine damages. The combustion modeling framework developed in this project can consider turbulence-chemistry interactions during end-gas ignition, while using a reasonably detailed chemical mechanism developed for ignition and combustion reactions under engine relevant conditions, and the subtle characteristics of spark-ignited flame propagation. It is developed in the context of large eddy simulation (LES), which can capture stochastic in-cylinder processes. The developed model is incorporated into a commercial software for engine simulation, CONVERGE CFD, as a user defined function, and validated. Engine knock and knock-free experiments as well as direct numerical simulation (DNS) of end-gas ignition in homogeneous turbulence are performed to help model development and provide data sets for model validation. With further validation, the developed model is expected to advance the predictive capability for engine knock simulations and thus contribute to improving the fuel efficiency.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2020

Development of a Physics-Based Combustion Model for Engine Knock Prediction

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2023

Aqueous Bromide Discharges from U.S. Coal-Fired Power Plants: Points of Origin, Concentration Ranges, and Effluent Treatment Costs

Bromide discharges from coal-fired power plants have received increased attention from regulatory bodies due to their contribution to the formation of disinfection byproducts (DBPs) in downstream drinking water treatment plants. Here, this paper characterizes the relative contributions of bromide from coal feedstocks and bromine-based mercury control processes, estimates the distribution of bromide concentrations at 85 active coal-fired power plants across the United States (U.S.) with wet flue gas desulfurization units, and estimates the cost of bromide removal from wastewater discharge using year 2020 data. Bromide discharges are estimated at the plant level using a combination of the reported coal rank and composition combusted, estimates of bromide addition in mercury control techniques under multiple halogen addition scenarios, and the estimated flue gas desulfurization (FGD) wastewater flow rate. The median, simulated plant-level estimation of total FGD wastewater flow is 18.3 gallons/min at a bromide concentration of 319 mg/L, equivalent to ~11.6 tonnes/year of bromide discharges to the environment. Next, we evaluated the expected cost of employing the best available technology (BAT) to control bromide discharges in FGD wastewater to prevent contributions to DBP formation. Treatment would need to remove more than 99.8% of bromide to reach the 0.2 mg/L voluntary incentive program (VIP) limit. The total cost of treatment depends on whether disposal is on- or off-site; the average costs for all plants combined come to an average of $110 million ($95.2/kgal) in 2021 U.S. dollars for on-site disposal, or $134 million ($115/kgal) for off-site disposal.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2021

Examination of probability distribution of mixture fraction in LES/FDF modelling of a turbulent partially premixed jet flame

An accurate prediction of the probability density function (PDF) of the mixture fraction is crucial to the prediction of combustion since mixing plays an important role in turbulent non-premixed and partially premixed flames. This work provides an assessment of the large-eddy simulation (LES)/filtered density function (FDF) method for the prediction of the PDF of the mixture fraction. The advantage of the LES/FDF method is that it provides the full predictions of the statistical distribution of scalars including the mixture fraction. The predictive accuracy of the method for the PDF is yet to be fully validated. Assessing the prediction of the PDF of the mixture fraction, a conserved scalar, is an important starting point. The Sydney/Sandia inhomogeneous inlet jet flame is used as a test case. A quick comparison shows that the LES/FDF predicted PDF shapes of the mixture fraction deviate significantly from the commonly presumed Beta-PDF as well as from the experimental data in the flame. Here, to examine the source of the discrepancy, we clarify the different PDF definitions used in the comparison among the predictions, measurements, and the presumed shape PDFs. The discrepancy observed from the comparison is largely reconciled by clarifying the difference between the PDFs that are examined. The PDF of the resolved mixture fraction is shown to be close to the Beta-PDF in both the measurements and predictions, while the PDF directly deduced from the LES/FDF particles deviates significantly from the Beta-PDF. A multimodal PDF analysis and a pseudo convergence analysis are conducted to provide plausible evidence to support the predicted multimodal PDF shapes. The sub-filter scale FDF is shown to be close to the Beta-PDF too through the construction of a synthesized PDF, which supports the common presumed Beta-PDF assumption used in the presumed PDF methods when combined with LES.

42 ENGINEERING↗

DOE OSTI2025

Evolution at the Edge: Real-Time Evolution for Neuromorphic Engine Control

Neuromorphic computing systems are attractive for real-time control at the edge because of their low power operation, real-time processing capabilities and their potential ability to do online learning. In this work, we describe an approach for performing real-time evolution of spiking neural networks for neuromorphic systems at the edge called Neuromorphic Optimization using Dynamic Evolutionary Systems or NODES. We apply this approach to real-time combustion engine control and develop an engine-specific hardware platform for NODES called FireBox. We demonstrate how the real-time evolution approach works in simulation and the performance of networks trained in simulation on the physical engine.

Maldonado Puente, Bryan [ORNL] (ORCID:000000033880↗