A robust framework for the simulation of scalar transport in turbulent spray combustion.
Abstract not provided.
Wenzel, Everett↗
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.
Simulation of Engine Expansion for Transparent Nozzle Combustion Research.
Abstract not provided.
Abers, Paul Michael↗
Lagrangian simulations of spray interaction with a surface using a stochastic multi-regime impingement model
Modeling of spray impingement upon the piston or cylinder surface is important for hydrocarbon emission and soot predictions in high-fidelity simulation of internal combustion engines with direct injections. Due to the complexity of the physical process of spray impingement, current computational models heavily rely on the empirical studies which involve large uncertainties. A computational model is presented for spray impingement on solid walls in general computational fluid dynamics (CFD) simulations. The model adopts the Marmanis-Thoroddsen theory for determining the number of secondary droplets and a stochastic method for calculating the velocities of the secondary droplets. The numerical methods and parallel implementation of the impingement model are presented in detail. The model is implemented in Grit, which is a recently developed CPU/GPU performance-portable Lagrangian particle library, for accelerating spray simulations. The model and its implementation are tested for examples of spray impingements on flat and smooth surfaces with different impingement angles.
Ge, Wenjun↗
High-efficiency catalytic reduction of residual oxygen for purification of carbon dioxide streams from high-pressure oxy-combustion systems
Pressurized oxy-combustion is a promising technology for carbon capture, utilization, and storage. For the captured CO 2 to be used for enhanced oil recovery or stored in geological formations, flue gas impurities, including residual O 2 in the CO 2 stream, must be purified to meet the purity specifications. A catalytic approach to reducing residual O 2 with CH 4 was investigated in this study. Five CoMn- and Cu-based catalysts were synthesized or acquired, and a reverse-flow fixed-bed reactor was used to assess their performance for O 2 removal from a simulated oxy-combustion flue gas at 15 bar. The impacts of the operating parameters on O 2 removal, such as temperature, gas hourly space velocity, O 2 /CH 4 ratio, and gas pressure, were investigated. Among the tested catalysts, the two CoMn catalysts were superior in both activity and selectivity, with the reaction lighting off at about 350 °C and achieving 99% O 2 removal at about 500 °C. Finally, the kinetics of the catalytic reaction is discussed, and the Mars–van Krevelen redox mechanism is deemed valid for describing the reaction pathway for the top-performing CoMn catalysts. The catalytic reaction was determined to be first order in CH 4 and zero order in O 2 under the test conditions.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Data-driven selection of stiff chemistry ODE solver in operator-splitting schemes
Most computational fluid dynamics simulations of practical combustion applications employ operator-splitting schemes, where chemistry and transport are separated and integrated with distinct numerical methods. The changes in composition due to chemistry are evaluated by solving ordinary differential equations (ODE) in each cell of the computational domain, which typically dominates the computational cost when detailed chemistry is considered. In this work, a data-driven approach for the selection of chemistry ODE solvers in operator-splitting schemes is presented. Neural networks are used to predict the ODE solvers CPU times and errors for a given thermochemical state. This allows the selection of an optimal ODE solver on a cell-by-cell, timestep-by-timestep basis. The models are trained using a wide set of thermochemical states generated through partially-stirred reactors and flames simulations. The methodology is validated by quantifying the prediction errors, the classification accuracy, and the computational speedup. The model predicts the optimal ODE solver for 70 to 95% of the validation cases and decreases the computional cost by a factor of 3 or more. The generalizability of the methodology to different chemical mechanisms and different fuels is assessed and it is shown that the model’s performance is only slightly degraded and its applicability is significantly enhanced if the inputs to the neural networks are restricted to a small set of thermochemical state variables present in most chemical mechanisms. In conclusion, the models are used in an homogeneous reactor case and a multi-dimensional CFD simulation of a diesel spray at high pressure where a speedup of more than 3 is achieved.
42 ENGINEERING↗
Characterizing Tradeoffs in Memory, Accuracy, and Speed for Chemistry Tabulation Techniques
Chemistry tabulation is a common approach in practical simulations of turbulent combustion at engineering scales. Linear interpolants have traditionally been used for accessing precomputed multidimensional tables but suffer from large memory requirements and discontinuous derivatives. Higher-degree interpolants address some of these restrictions but are similarly limited to relatively low-dimensional tabulation. Artificial neural networks (ANNs) can be used to overcome these limitations but cannot guarantee the same accuracy as interpolants and introduce challenges in reproducibility and reliable training. These challenges are enhanced as the physics complexity to be represented within the tabulation increases. Here, we assess the efficiency, accuracy, and memory requirements of Lagrange polynomials, tensor product B-splines, and ANNs as tabulation strategies. We analyze results in the context of nonadiabatic flamelet modeling where higher dimension counts are necessary. While ANNs do not require structuring of data, providing benefits for complex physics representation, interpolation approaches often rely on some structuring of the table. Interpolation using structured table inputs that are not directly related to the variables transported in a simulation can incur additional query costs. This is demonstrated in the present implementation of heat losses. We show that ANNs, despite being difficult to train and reproduce, can be advantageous for high-dimensional, unstructured datasets relevant to nonadiabatic flamelet models. Furthermore we demonstrate that Lagrange polynomials show significant speedup for similar accuracy compared to B-splines.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
A direct numerical simulation study of the dilution tolerance of propane combustion under spark-ignition engine conditions
Modern spark ignition internal combustion (IC) engines rely on highly diluted fuel-air mixtures to achieve high brake thermal efficiencies. To support this, new engine designs have introduced high stroke-to-bore ratios and cylinder head designs that promote high tumble flow and turbulence intensities. However, mixture dilution through exhaust gas recirculation (EGR) is limited by combustion instabilities manifested in the form of cycle-to-cycle variability. Propane has been observed to have superior EGR dilution tolerance than gasoline, which makes it a very competitive low-carbon fuel for the new IC engines without sacrificing efficiency. Two-dimensional direct numerical simulations (DNS) are performed with detailed chemistry to study and contrast the effect of turbulence intensity and dilution on propane and iso-octane premixed flames at high pressure conditions similar to those in-cylinder. A new reduced mechanism for propane consisting of 53 transported species and 17 quasi-steady state species is developed based on a previously published mechanism and used in these simulations. Three levels of turbulence intensity and two levels of exhaust gas dilution are chosen based on conditions relevant to IC engine operation. The DNS results are analyzed based on the evolution of the flame surface area and the statistics of its driving terms, which are found to be similar for both fuels when there is no dilution but considerably different under high dilution. The analysis of the DNS data provides fundamental insights into the underlying mechanisms for improved stability under dilution.
42 ENGINEERING↗
A unified reaction network on the formation of five-membered ringed polycyclic aromatic hydrocarbons (PAHs) and their role in ring expansion processes through radical–radical reactions
Exploiting a chemical microreactor in combination with an isomer-selective product identification through fragment-free photoionization utilizing tunable vacuum ultraviolet (VUV) light in tandem with the detection of the ionized molecules by a high resolution reflection time-of-flight mass spectrometer (Re-TOF-MS), the present investigation reveals molecular mass growth processes to four distinct polycyclic aromatic hydrocarbons carrying two six- and one five-membered ring (C 13 H 10 ): 3H-cyclopenta[a]naphthalene, 1H-cyclopenta[b]naphthalene, 1H-cyclopenta[a]naphthalene, and fluorene in the gas phase. Temperatures of 973 and 1023 K simulating conditions in combustion settings along with circumstellar envelopes of carbon-rich stars and planetary nebulae. These reactions highlight the importance of methyl-substituted aromatic reactants (biphenyl, naphthalene) which can be converted to the methylene (–CH 2 ˙) motive by hydrogen abstraction or photolysis. Upon reaction with acetylene, methylene-substituted aromatics carrying a hydrogen atom at the ortho position of the ring can be then converted to cyclopentadiene-annulated aromatics thus providing a versatile pathway to five-membered ring aromatics at elevated temperatures.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Measuring Critical Fuel Properties for Simulations to Accelerate Fuel Qualification
NREL's sustainable aviation fuel (SAF) research in the area of fuel properties and end use is focused on accelerating SAF development and deployment. This will be achieved by close collaboration with BETO program researchers on fuel property measurement and development of structure-property relationships. Additionally, fuel properties measured at often extreme engine operating conditions will be measured, and will be used as inputs to exascale ready simulations of aviation combustion turbines to reveal impacts of properties on turbine operation. This research will potentially lead to a widening of the allowable fuel property range for SAF, expanding SAF production process options and production volumes.
ADVANCED PROPULSION SYSTEMS↗
PeleMP: The Multiphysics Solver for the Combustion Pele Adaptive Mesh Refinement Code Suite
Combustion encompasses multiscale, multiphase reacting flow physics spanning a wide range of scales from the molecular scales, where chemical reactions occur, to the device scales, where the turbulent flow is affected by the geometry of the combustor. This scale disparity and the limited measurement capabilities from experiments make modeling combustion a significant challenge. Recent advancements in high-performance computing (HPC), particularly with the Department of Energy's Exascale Computing Project (ECP), have enabled high-fidelity simulations of practical applications to be performed. The major physics submodels, including chemical reactions, turbulence, sprays, soot, and thermal radiation, exhibit distinctive computational characteristics that need to be examined separately to ensure efficient utilization of computational resources. This paper presents the multiphysics solver for the Pele code suite, called PeleMP, which consists of models for spray, soot, and thermal radiation. Here, the mathematical and algorithmic aspects of the model implementations are described in detail as well as the verification process. The computational performance of these models is benchmarked on multiple supercomputers, including Frontier, an exascale machine. Results are presented from production simulations of a turbulent sooting ethylene flame and a bluff-body swirl stabilized spray flame with sustainable aviation fuels to demonstrate the capability of the Pele codes for modeling practical combustion problems with multiphysics. This work is an important step toward the exascale computing era for high-fidelity combustion simulations providing physical insights and data for predictive modeling of real-world devices.
42 ENGINEERING↗
The role of cool-flame fluctuations in high-pressure spray flames, studied using high-speed optical diagnostics and Large-Eddy Simulations
Here this work investigates the low- and high-temperature ignition and combustion processes, applied to the Engine Combustion Network Spray A flame, combining advanced optical diagnostics and large-eddy simulations (LES). Simultaneous high-speed (50 kHz) formaldehyde (CH 2 0) planar laser-induced fluorescence (PLIF) and line-of-sight OH* chemiluminescence imaging were used to measure the low- and high-temperature flame, during ignition as well as during quasi-steady combustion. While tracking the cool flame at the laser sheet plane, the present experimental setup allows detection of distinct ignition spots and dynamic fluctuations of the lift-off length over time, which overcomes limitations for flame tracking when using schlieren imaging. After significant development to improve LES prediction of the low-and high-temperature flame position, both during the ignition processes and quasi-steady combustion, the simulations were analyzed to gain understanding of the mixture variance and how this variance affects formation/consumption of CH 2 0. Analysis of the high-temperature ignition period shows that a key improvement in the LES is the ability to predict heterogeneous ignition sites, not only in the head of the jet, but in shear layers at the jet edge close to the position where flame lift-off eventually stabilizes. The LES analysis also shows concentrated pockets of CH 2 0, in the center of jet and at 20 mm downstream of the injector (in regions where the equivalence ratio is greater than 6), that are of similar length scale and frequency as the experiment (approximately 5–6 kHz). The periodic oscillation of CH 2 0 match the frequency of pressure waves generated during auto-ignition and reflected within the constant-volume vessel throughout injection. The ability of LES to capture the periodic appearance and destruction of CH 2 0 is particularly important because these structures travel downstream and become rich premixed flames that affect soot production.
42 ENGINEERING↗
CFD/FEA Co-Simulation Framework for Analysis of the Thermal Barrier Coating Design and Its Impact on the HD Diesel Engine Performance
Thermal barrier coatings (TBCs) have been investigated both experimentally and through simulation for mixing controlled combustion (MCC) concepts as a method for reducing heat transfer losses and increasing cycle efficiency, but it is still a very active research area. Early studies were inconclusive, with different groups discovering obstacles to realizing the theoretical potential. Nuanced papers have shown that coating material properties, thickness, microstructure, and surface morphology/roughness all can impact the efficacy of the thermal barrier coating and must be accounted for. Adding to the complexities, a strong spatial and temporal heat flux inhomogeneity exists for mixing controlled combustion (diesel) imposed onto the surfaces from the impinging flame jets. In support of the United States Department of Energy SuperTruck II program goal to achieve 55% brake thermal efficiency on a heavy-duty diesel engines, this study sought to develop a deeper insight into the inhomogeneous heat flux from mixing controlled combustion on thermal barrier coatings and to infer concrete guidance for designing coatings. To that end, a co-simulation approach was developed that couples high-fidelity computational fluid dynamics (CFD) modeling of in-cylinder processes and combustion, and finite element analysis (FEA) modeling of the thermal barrier-coated and metal engine components to resolve spatial and temporal thermal boundary conditions. The models interface at the surface of the combustion chamber; FEA modeling predicts the spatially resolved surface temperature profile, while CFD develops insights into the effect of the thermal barrier coating on the combustion process and the boundary conditions on the gas side. The paper demonstrates the capability of the framework to estimate cycle impacts of the temperature swing at the surface, as well as identify critical locations on the piston/thermal barrier coating that exhibit the highest charge temperature and highest heat fluxes. In addition, the FEA results include predictions of thermal stresses, thus enabling insight into factors affecting coating durability. An example of the capability of the framework is provided to illustrate its use for investigating novel coatings and provide deeper insights to guide future coating design.
42 ENGINEERING↗
Large-eddy simulation study on cycle-to-cycle variation of knocking combustion in a spark-ignition engine
The cycle-to-cycle variation in the knock intensity is commonly encountered under abnormal combustion conditions. The severity of these abnormal combustion events can vary significantly, and the efficiency of engines at high loads is limited in practice by heavy knocking phenomena. Since, a thorough analysis of such recurrent but non-cyclic phenomena via experiments alone becomes highly cumbersome, in the present work, a multi-cycle large-eddy simulation study was performed to quantitatively predict cyclic variability in the combustion process and cyclic knock intensity variability in a direct injection spark-ignition engine. To account for the turbulence-chemistry interaction effects on flame propagation, the G-equation combustion model was used. Detailed chemistry was solved outside the flame front with a toluene primary reference fuel skeletal kinetic mechanism. For both the mild knock and heavy knock conditions, the numerical results were validated against experimental measurements. Based on the simulation results, a correlation analysis was performed considering combustion phasing, peak cylinder pressure and maximum amplitude of pressure oscillation. Furthermore, a detailed three-dimensional spatial analysis illustrated the evolution of auto-ignition kernel development and propagation of pressure waves during knocking combustion for three typical cycles with different knock intensities. In this process, it was found that an early occurrence of auto-ignition in the end gas was prone to high knock intensity. Although multiple auto-ignition kernels were observed in different cycles, the degree of coupling between chemical heat release and pressure waves varied, thereby leading to different maximum amplitude of pressure oscillation values.
42 ENGINEERING↗
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↗
The effects of mixture preburning on detonation wave propagation
Pressure gain combustion in the form of continuous detonations can provide a significant increase in the efficiency of a variety of propulsion and energy conversion devices. In this regard, rotating detonation engines (RDEs) that utilize an azimuthally-moving detonation wave in annular systems are increasingly seen as a viable approach to realizing pressure gain combustion. However, practical RDEs that employ non- premixed fuel and oxidizer injection need to minimize losses through a number of mechanisms, including turbulence-induced shock-front variations, incomplete fuel-air mixing, and premature deflagration. In this study, a canonical stratified detonation configuration is used to understand the impact of preburning on detonation efficiency. It was found that heat release ahead of the detonation wave leads to weaker shock fronts, delayed combustion of partially-oxidized fuel-air mixture, and non-compact heat release. Furthermore, large variations in wave speeds were observed, which is consistent with wave behavior in full-scale RDEs. Peak pressures in the compression region or near triple points were considerably lower than the theoretically-predicted values for ideal detonations. Analysis of the detonation structure indicates that this deflagration process is parasitic in nature, reducing the detonation efficiency but also leading to heat release far behind the wave that cannot directly strengthen the shock wave. As a result, this parasitic combustion leads to commensal combustion (heat release far downstream of the wave), indicating that it is the root cause of combustion efficiency losses.
42 ENGINEERING↗
Analysis of combustion acoustic phenomena in compression–ignition engines using large eddy simulation
As computational capabilities continue to grow, exploring the limits of computational fluid dynamics to capture complex and elusive phenomena, which are otherwise difficult to study by experimental techniques, is one of the main targets for the research community. This paper presents a detailed analysis of the physical processes that lead to combustion noise emissions in internal combustion engines. In particular, diesel combustion in a compression-ignition (CI) engine is studied in order to understand the singular behavior of the in-cylinder flow field responsible for the acoustic emissions. The main objective is, therefore, to improve the understanding of the phenomena involved in CI engine noise using large eddy simulations. Several visualization methods are employed to investigate the connection between combustion behavior and its effects on the pressure field. In addition, proper orthogonal decomposition is used to analyze the modal energy distribution among all the acoustic modes. The results show that the acoustic signature is fundamentally conditioned by the intensity of the premixed combustion rather than by the pressure oscillations generated by turbulent fluctuations in the flame surface established during the diffusion stage.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Reactive Molecular Dynamics Simulations and Quantum Chemistry Calculations To Investigate Soot-Relevant Reaction Pathways for Hexylamine Isomers
Sooting tendencies of a series of nitrogen-containing hydrocarbons (NHCs) have been recently characterized experimentally using the yield sooting index (YSI) methodology. This work aims to identify soot-relevant reaction pathways for three selected C6H15N amines, namely, dipropylamine (DPA), diisopropylamine (DIPA), and 3,3-dimethylbutylamine (DMBA) using ReaxFF molecular dynamics (MD) simulations and quantum mechanical (QM) calculations and to interpret the experimentally observed trends. ReaxFF MD simulations are performed to determine the important intermediate species and radicals involved in the fuel decomposition and soot formation processes. QM calculations are employed to extensively search for chemical reactions involving these species and radicals based on the ReaxFF MD results and also to quantitatively characterize the potential energy surfaces. Specifically, ReaxFF simulations are carried out in the NVT ensemble at 1400, 1600, and 1800 K, where soot has been identified to form in the YSI experiment. These simulations account for the interactions among test fuel molecules and pre-existing radicals and intermediate species generated from rich methane combustion, using a recently proposed simulation framework. ReaxFF simulations predict that the reactivity of the amines decrease in the order DIPA > DPA > DMBA, independent of temperature. Both QM calculations and ReaxFF simulations predict that C2H4, C3H6, and C4H8 are the main nonaromatic soot precursors formed during the decomposition of DPA, DIPA, and DMBA, respectively, and the associated reaction pathways are identified for each amine. Both theoretical methods predict that sooting tendency increases in the order DPA, DIPA, and DMBA, consistent with the experimentally measured trend in YSI. This work demonstrates that sooting tendencies and soot-relevant reaction pathways of fuels with unknown chemical kinetics can be identified efficiently through combined ReaxFF and QM simulations. Overall, predictions from ReaxFF simulations and QM calculations are consistent, in terms of fuel reactivity, major intermediates, and major nonaromatic soot precursors.
fuel decomposition↗
C3MechLite: An integrated component library of compact kinetic mechanisms for low-carbon, carbon neutral and zero-carbon fuels
Based on our latest detailed chemical reaction mechanism, C3MechV4.0, we have developed two reduced reaction mechanisms—C3MechLite and C3MechCore—targeting C 0 –C 3 chemical species including NH 3 . C3MechLite (61 species), contains a number of species comparable to GRI-Mech (53 species), that can accurately predict the combustion characteristics of hydrogen, carbon monoxide, ammonia, methane, natural gas, nitrogen oxides, and their mixtures for a wide range of conditions. C3MechCore (118 species) targets a more comprehensive range of C 0 –C 3 fuels, including ammonia, methanol, ethanol, and dimethyl ether. Both mechanisms demonstrate predictive accuracy comparable to C3MechV4.0 for the combustion characteristics of the target fuels. C3MechLite is designed with a component library structure, enabling further reduction in mechanism size depending on the fuel(s) of interest for 2D/3D numerical simulations. Various combinations of component libraries were validated, and the average prediction error remains within 1 % compared to C3MechLite. Furthermore, the mechanism was applied to 3D LES simulations of H 2 lifted flames and was confirmed to reproduce flame characteristics with high accuracy. C3MechLite and its component library structure enable high-fidelity and computationally efficient chemical kinetic mechanisms, paving the way for application in more complex combustion simulations.
Ammonia↗