Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reactive burn model”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Numerical modeling of plasma assisted deflagration to detonation transition in a microscale channel

Here, this work numerically studies the plasma assisted deflagration to detonation transition (DDT) of H 2 /O 2 mixtures in a microscale channel with detailed chemistry and transport. The results show that the DDT onset time is non-monotonically dependent on the discharge pulse number. The DDT is accelerated with small pulse numbers, whereas retarded with large ones. Two different DDT regimes, respectively at a small and large plasma discharge number, via acoustic choking of the burned gas and plasma-enhanced reactivity gradient without acoustic choking, are observed. Without plasma discharge, pronounced pressure and temperature gradients in front of the flame are generated by acoustic compression after the choking of the burned gas, triggering DDT via autoignition. With small plasma pulse numbers, the plasma-generated species enhance the ignition kinetics and lead to an increased reactivity in the boundary layer. After the choking of the burned gas, the plasma-enhanced reactivity advances the sequence of autoignition near the wall, strengthens ignition-shock wave coupling, and accelerates DDT. However, with a large discharge pulse number, a direct autoignition initiating DDT can occur without the acoustic choking of the burned gas due to the strongly accelerated reactivity and elevated temperature. In this case, DDT onset is retarded because the elevated temperature increases sonic velocity and the increased reactivity accelerates fuel oxidation in front of the flame, decelerating the formation of a leading shock and subsequent pressure buildup ahead of the flame. The present modeling reveals that no matter with or without plasma discharge, DDT is initiated by autoignition in thermal, pressure, and reactivity gradient fields via the Zel'dovich gradient mechanism. The acoustic choking of the burned gas may not be the necessary condition of DDT with strong plasma-enhanced reactivity gradient. This work provides an answer to the experimentally observed non-monotonic DDT onset time by plasma, which provides guidance to control DDT in advanced detonation engines and fire safety of hydrogen-fueled catalytic reactors in microchannels by non-equilibrium plasma discharge.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Type Ia supernovae deflagration-to-detonation transition explosions powered by the Zel’dovich reactivity gradient mechanism

ABSTRACT Our aim in this work is to identify and explain the necessary conditions required for an energetic explosion of a Chandrasekhar-mass white dwarf. We construct and analyse weakly compressible turbulence models with nuclear burning effects for carbon/oxygen plasma at a density expected for the deflagration-to-detonation transition (DDT) to occur. We observe the formation of carbon deflagrations and transient carbon detonations at early times. As turbulence becomes increasingly inhomogeneous, sustained carbon detonations are initiated by the Zel’dovich reactivity gradient mechanism. The fuel is suitably preconditioned by the action of compressive turbulent modes with wavelength comparable to the size of resolved turbulent eddies; no acoustic wave is involved in this process. Oxygen detonations are initiated, aided either by reactivity gradients or by collisions of carbon detonations. The observed evolutionary time-scales are found to be sufficiently short for the above process to occur in the expanding, centrally ignited massive white dwarf. The inhomogeneous conditions produced prior to the DDT might be of consequence for the chemical composition of the outer ejecta regions of Type Ia supernovae from the single degenerate channel, and offer the potential for validation of the proposed model.

Brooker, E.↗
DOE OSTI2021

Modeling Volatility-Based Aerosol Phase State Predictions in the Amazon Rainforest

Organic aerosol (OA) is a complex matrix of various constituents—fresh (primary organic aerosols—POA) and aged via oxidation (secondary organic aerosols—SOA), generated from biogenic, anthropogenic, and biomass burning sources. The viscosity of OA can be critical in influencing new particle formation, reactive uptake processes that impact evaporation-growth kinetics, and the lifetime of particles in the atmosphere. We report this work utilizes a well-defined relationship between volatility and viscosity for pure compounds, which we incorporated within the Weather Research and Forecasting Model coupled to chemistry (WRF-Chem) to simulate the phase state and viscosity of bulk OA during the dry-to-wet transition season (September–October) in the Amazon rainforest during 2014. Our simulations indicate spatial and temporal heterogeneity in aerosol phase state often not captured by global-scale models. We show the strong role of water associated with organic aerosol ($w_s$) as the dominant factor that can be used to quantitatively estimate OA viscosity. Analysis of WRF-Chem simulations across the entire atmospheric column indicates a strong inverse log-linear relationship between $w_s$ and OA viscosity with a correlation coefficient approaching 1, in the background and biomass burning-influenced conditions. At high altitudes where relative humidity (RH) and temperatures are low, our simulations indicate that OA exists in a semisolid-/solid-like phase state, consistent with previous studies. OA hygroscopicity is strongly correlated (ca. –0.8) with OA viscosity at RH ca. 30–50%, but this RH range is found mostly at low OA concentrations and the middle troposphere (ca. 6–10 km altitudes) in our simulated domain. OA hygroscopicity is uncorrelated with viscosity at higher-RH (near surface) and lower-RH (upper troposphere) regimes. At the urban site near surface, where day–night differences in RH are significant, RH is found to drive the phase state. At the background forested site near surface, where day–night RH differences are small, biomass burning-influenced OA is semisolid and a significant OA associated with background conditions is liquid-like. Simulations indicate a long tail of OA viscosity frequency distributions extending in the semisolid/solid regimes over background biogenic-influenced conditions due to the role of low-volatility OA components such as monoterpene oxidation products.

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

Control-oriented core-SOL-divertor model to address integrated burn and divertor control challenges in ITER

The real-time regulation of a burning plasma’s temperature and density, or burn control, will be necessary to produce high fusion power in future tokamaks like ITER. This is made more challenging due to the plasma’s nonlinear characteristics and the interdependence between the core-plasma and edge-plasma regions. For example, a raising plasma temperature leads to increasing reactivity and therefore to more alpha-particle heating, which further increases temperature. Furthermore, a raise of the fusion power increases the heat flow through the scrape-off-layer (SOL), which can compromise the integrity of the divertor without proper safeguards. For control design, a model-based approach is attractive because it can directly incorporate the nonlinear, coupled, burning-plasma dynamics into the design. To facilitate this design approach, a control-oriented core-SOL-divertor (CSD) model is presented in this work. In this CSD model, a core-plasma model captures the nonlinear dynamics of the core’s density and temperature, and a SOL-divertor model defines the plasma conditions at the separatrix and divertor including the heat load on the target plates. The core-plasma and SOL-divertor models are coupled through the exchange of various variables. In particular, the SOL-divertor model yields the separatrix temperature and the influx of recycled particles into the core-plasma. Further, these variables influence the power and particle balances captured by the core-plasma model. In return, the core-plasma model determines the intensity of the heat and particles fluxes across the separatrix, and this outflow strongly impacts the SOL-divertor model. Therefore, the power and density of the core-plasma, which can be readily modulated through external heating systems and pellet injection, can be viewed as control knobs for the SOL-divertor region in addition to the gas puffing. In simulations of the CSD model, it is demonstrated how external actuation can be utilized to meet burn control and divertor control objectives simultaneously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

SCALE Modeling of the Sodium Cooled Fast-Spectrum Advanced Burner Test Reactor

This report documents the modeling and simulation of a sodium-cooled fast reactor (SFR) as part of a U.S. Nuclear Regulatory Commission–sponsored project to assess the modeling and simulation capabilities for accident progression, source term, and consequence analysis for advanced reactor technologies with the Oak Ridge National Laboratory code SCALE and the Sandia National Laboratories (SNL) code MELCOR. Based on publicly available benchmark specifications, a fully heterogeneous 3D SCALE model of the 250 MWth Advanced Burner Test Reactor (ABTR) was developed to demonstrate SCALE’s capabilities for full-core reactivity analysis, fuel inventory prediction, and decay heat analysis of an SFR. The benchmark specifications contain modeling details for the ABTR core at the beginning of equilibrium cycle (BOEC) at operating conditions; they were derived from a 2006 preconceptual design report produced by Argonne National Laboratory. The ABTR was designed to demonstrate reactor-based transmutation of transuranics, that is, to “burn” transuranics recovered from light-water reactor (LWR) spent fuel. The ABTR’s fuel is designed to operate in 4 month cycles using uranium/transuranic (U/TRU) metallic fuel, with a TRU content of approximately 20%, at a conversion ratio of approximately 0.6. Various reactivity calculations were performed with SCALE for the ABTR and, where possible, compared with results available in the open literature. Additionally, SCALE was used to perform a full-core depletion calculation over the 4 month cycle to obtain the nuclide inventory at the end of equilibrium cycle (EOEC). These nuclide inventories, decay heat, power profiles, and reactivity feedback coefficients at EOEC represent the initial conditions for analyzing severe accident scenarios with MELCOR.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

A second order yield-temperature relation for accurate inference of burn-averaged quantities in multi-species plasmas

Measured yields and ion temperatures inferred from the fusion product energy spectra can be used as metrics for the performance of an ICF implosion. This can be to infer species separation, thermal decoupling, flows or other effects that can cause the inferred ion temperatures to deviate from the true underlying thermal temperature and the yield ratio to deviate from the expected value. Direct inference of the impact of these effects on observed temperatures and yields can be difficult to uncover due to underlying dependence on the shape and time evolution of the temperature and density proles of the fusing plasma. Due to differences in the temperature dependence of the reactivities, different fusion products are emitted from different regions and times within the implosion. In order to properly account for this, a second order analytic expression relating the apparent temperatures and yield ratios is developed. This expression can be coupled to models of yield and/or temperature altering effects to infer their burn-averaged impact on an implosion. The second order expression shows significant improvement over lower order expressions in synthetic data studies. Demonstrations of its applications to synthetic data coupled with models of ion thermal decoupling and radial flows are presented. In the case of thermal decoupling both first and second order expressions show reasonable levels of accuracy. To consistently infer the amplitude of radial flow with <10% error the second order equation is required.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A second order yield-temperature relation for accurate inference of burn-averaged quantities in multi-species plasmas

Measured yields and ion temperatures inferred from the fusion product energy spectra can be used as metrics for the performance of an ICF implosion. This can be used to infer species separation, thermal decoupling, flows, or other effects that can cause the inferred ion temperatures to deviate from the true underlying thermal temperature and the yield ratio to deviate from the expected value. Direct inference of the impact of these effects on observed temperatures and yields can be difficult to uncover due to the underlying dependence on the shape and time evolution of the temperature and density profiles of the fusing plasma. Due to differences in the temperature dependence of the reactivities, different fusion products are emitted from different regions and times within the implosion. In order to properly account for this, a second-order analytical expression relating the apparent temperatures and yield ratios is developed. This expression can be coupled to models of yield and/or temperature altering effects to infer their burn-averaged impact on an implosion. The second-order expression shows significant improvement over lower-order expressions in synthetic data studies. Demonstrations of its applications to synthetic data coupled with models of ion thermal decoupling and radial flows are presented. In the case of thermal decoupling, both first and second-order expressions show reasonable levels of accuracy. To consistently infer the amplitude of radial flow with a <10% error, the second-order equation is required.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Experimental and Kinetic Modeling Study of 3-Methyl-2-butenol (Prenol) Oxidation

Longer chain alcohols with 4–5 carbon atoms are attractive alternative fuels as they can be derived from biological sources and since their combustion leads to lower exhaust gas levels of NO x and soot compared to commercial fossil fuels. The auto-ignition behavior of fuels that contain both a hydroxyl group and a C=C double bond in their molecular structure is not well established in the literature. Understanding the influence of these functional groups on the ignition behavior of fuels is critical in the development of tailor-made fuels for advanced combustion engines. In this study, ignition delay times of an unsaturated alcohol, 3-methyl-2-butenol (prenol), are measured using a high-pressure shock tube and a rapid compression machine at pressures of 15 and 30 bar at equivalence ratios of 0.5, 1.0, and 2.0 in “air” in the temperature range 600–1400 K. Furthermore, a detailed kinetic model is developed and validated using the new experimental data in this study. In addition, speciation data in a jet-stirred reactor, ignition delay times, and laminar burning velocities available in the literature were also used to validate the new kinetic model. Fuel flux and sensitivity analyses are performed using this new model to determine the important fuel consumption pathways and critical reactions that affect the reactivity of prenol.

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

Automated power-following control for nuclear thermal propulsion startup and shutdown using MOOSE-based applications

This paper describes an investigation of various automated control strategies applied to a full-core multiphysics Griffin/Bison/RELAP-7 model of a prototypical nuclear thermal propulsion system. In all cases, control is achieved by actuating control drums based on the demanded power and predicted quantities from the numerical model. One key finding is that hybrid proportional integral derivative controllers – a novel type of controller that uses both power and reactivity predicted signals – can demonstrate a level of performance rivaling that of period-generated control. The former requires parameter tuning, while the latter mostly necessitates providing reactivity coefficients and temperature rates of change, which could be very challenging to accurately measure in real-time. In addition, decay heat plays an important role in determining cooling requirements during cooldown phases. A decay heat model, accounting for burn time and throttling, was derived, and is incorporated within the model to simulate the steady-state and shutdown phases and satisfyingly follow the power demand. However, temperature overshoots and non-constant specific impulse during throttling will necessitate further improvements.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

An experimental and kinetic modeling study of cyclopentane and dimethyl ether blends

Cyclopentane is a suitable naphthene, or cycloalkane, in a palette for multi-component gasoline surrogate fuels due to its presence in market fuels and its relevance to alkyl substituted cyclopentanes also present. However, the previous oxidation studies of cyclopentane have primarily focused on neat mixtures. Blending cyclopentane with dimethyl ether in this work therefore serves to inform our understanding of, and improve predictive models for, multi-component mixtures. Here, the auto-ignition of cyclopentane/dimethyl ether blends was studied in a high-pressure shock tube and in a rapid compression machine. A wide range of temperatures (650 – 1350 K) and elevated pressures of 20 and 40 bar were studied at equivalence ratios of 0.5, 1.0 and 2.0 in air for two blending ratios (30/70 and 70/30 mole% cyclopentane/di-methyl ether mixtures). A detailed kinetic model for cyclopentane was revised to capture the measured ignition delay times and apparent heat release rates in this study. Literature ignition delay time, jet-stirred reactor, and laminar burning velocity measurements of neat cyclopentane were used as additional validation. Improvements to the kinetic model were based on recent literature studies related to sub-models including cyclopentene and cyclopentadiene which allowed the removal of previous local rate-constant optimizations. Low temperature reactivity of cyclopentane was found to be controlled by the branching ratio between concerted elimination of HO͘ 2 and the strained formation of Q&#775OOH radicals in agreement with previous studies. In this study, the low branching ratio of Q&#775OOH formation increases the influence of a competing consumption pathway for cyclopentyl-peroxy (CPTO͘ 2 J) radicals. The sensitivity of the simulated ignition delay times to the formation of cyclopentyl hydroperoxide (CPTO 2 H), from CPTO͘ 2 J and HO͘ 2 , is discussed. The current model is used to analyze the influence of dimethyl ether on the reactivity of cyclopentane in the context of previous literature studies of dimethyl ether binary blends with ethanol and toluene.

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

Modeling Novel Aqueous Particle and Cloud Chemistry Processes of Biomass Burning Phenols and Their Potential to Form Secondary Organic Aerosols

Phenols emitted from biomass burning contribute significantly to secondary organic aerosol (SOA) formation through the partitioning of semivolatile products formed from gas-phase chemistry and multiphase chemistry in aerosol liquid water and clouds. The aqueous-phase SOA (aqSOA) formed via hydroxyl radical (˙OH), singlet molecular oxygen ( 1 O 2 *), and triplet excited states of organic compounds ( 3 C*), which oxidize dissolved phenols in the aqueous phase, might play a significant role in the evolution of organic aerosol (OA). However, a quantitative and predictive understanding of aqSOA has been challenging. Here, we develop a stand-alone box model to investigate the formation of SOA from gas-phase ˙OH chemistry and aqSOA formed by the dissolution of phenols followed by their aqueous-phase reactions with ˙OH, 1 O 2 *, and 3 C* in cloud droplets and aerosol liquid water. Here, we investigate four phenolic compounds, i.e., phenol, guaiacol, syringol, and guaiacyl acetone (GA), which represent some of the key potential sources of aqSOA from biomass burning in clouds. For the same initial precursor organic gas that dissolves in aerosol/cloud liquid water and subsequently reacts with aqueous phase oxidants, we predict that the aqSOA formation potential (defined as aqSOA formed per unit dissolved organic gas concentration) of these phenols is higher than that of isoprene-epoxydiol (IEPOX), a well-known aqSOA precursor. Cloud droplets can dissolve a broader range of soluble phenols compared to aqueous aerosols, since the liquid water contents of aerosols are orders of magnitude smaller than cloud droplets. Further, our simulations suggest that highly soluble and reactive multifunctional phenols like GA would predominantly undergo cloud chemistry within cloud layers, while gas-phase chemistry is likely to be more important for less soluble phenols. But in the absence of clouds, the condensation of low-volatility products from gas-phase oxidation followed by their reversible partitioning to organic aerosols dominates SOA formation, while the SOA formed through aqueous aerosol chemistry increases with relative humidity (RH), approaching 40% of the sum of gas and aqueous aerosol chemistry at 95% RH for GA. Our model developments of biomass-burning phenols and their aqueous chemistry can be readily implemented in regional and global atmospheric chemistry models to investigate the aqueous aerosol and cloud chemistry of biomass-burning organic gases in the atmosphere.

09 BIOMASS FUELS↗
DOE OSTI2023

A Fully Explicit Integrator for Modeling Astrophysical Reactive Flows

Simulating complex astrophysical reacting flows is computationally expensive—reactions are stiff and typically require implicit integration methods. The reaction update is often the most expensive part of a simulation, which motivates the exploration of more economical methods. In this research note, we investigate how the explicit Runge–Kutta–Chebyshev (RKC) method performs compared to an implicit method when applied to astrophysical reactive flows. These integrators are applied to simulations of X-ray bursts arising from unstable thermonuclear burning of accreted fuel on the surface of neutron stars. We show that the RKC method performs with similar accuracy to our traditional implicit integrator, but is more computationally efficient when run on CPUs.

GPU computing↗
DOE OSTI2025

Scoping Analysis of Pebble-Bed Reactors for the Destruction of the Transuranic Inventory of LWR Spent Nuclear Fuel

With the forecasted increase in the construction and operation of nuclear reactors, there will be a corresponding increase in the quantity of spent nuclear fuel (SNF) that requires long-term storage. In SNF, transuranic isotopes contribute the most to the long-term radiotoxicity of the fuel and pose a proliferation risk. One option that has been explored to address these issues is the removal of the transuranic isotopes from SNF and the conversion of these isotopes into transuranic fuel (TRU fuel). Here, this work sought to determine how effective a micro-modular Pebble-Bed High-Temperature Gas-Cooled Reactor (PB-HTGR); the 10-MW High Temperature Gas-cooled Test Reactor (HTR-10); and a salt-cooled small-modular pebble-bed reactor (PBR), i.e. the generic Fluoride-cooled High-temperature Reactor (gFHR), are at reducing the inventory of transuranic isotopes while still maintaining the intrinsic safety features of the PBR designs, such as negative temperature coefficients of reactivity. Optimized pebble designs utilizing TRU fuel were found for both reactors through the adjustment for the packing fraction of fuel in each pebble. The Axial Zone Equilibrium Modeling (A-ZEM) method was used in this work to help select the optimized pebble design. Once an optimized pebble design was selected and an equilibrium model was produced, the results from the deep burn (DB) HTR-10 and gFHR designs were compared to the results of two models from the literature. While both the DB gFHR and the DB HTR-10 were able to reduce the weapons-usable transuranic inventory, the performance of these reactors did not match that of the small-modular PB-HTGRs in the literature. Therefore, a need was identified for further refinement of the gFHR design using TRU fuel, as the results for this model were more promising than those of the DB HTR-10, which was strongly limited by the high leakage intrinsic to micro-modular PB-HTGRs.

Transuranic fuel↗
DOE OSTI2025

Turbulent burning velocity of lean premixed hydrogen/air flames at engine conditions: Effects of turbulence intensity and length scale

For turbulent lean premixed hydrogen flames with strong thermodiffusively instabilities, most previous studies have focused on the influence of turbulence intensity, whereas the role of turbulence length scale is less well understood. Here, this study addresses this gap by conducting direct numerical simulations (DNS) of statistically planar turbulent premixed flames for a lean (ϕ=0.35) hydrogen/air mixture under independently varied turbulence intensity (u') and length scale (l T ) at engine-relevant thermodynamics conditions. Results show that as u' increases, the flame front becomes increasingly wrinkled, forming smaller cellular structures. In contrast, l T variations do not significantly alter the size of these structures. For the turbulent burning velocity (S T ), the normalized S T (i.e., S T /S L , where S L is the laminar flame speed) increases linearly with u', driven by both enhanced flame surface wrinkling (i.e., increased A T /A L ) and enhanced local burning rate (i.e., increased I 0 ). However, increasing l T reduces I 0 , despite a continued increase in A T /A L , resulting in only a marginal increase in S T /S L . To reveal the underlying mechanisms, especially the decreasing trend of I 0 with l T , local flame dynamics analyses are performed. It is found that as l T increases, the interaction between thermodiffusive effects and turbulence weakens due to the reduced tangential strain rate, while the flame curvature remains largely unchanged. This suppresses local reactivity enhancement and thus decreases I 0 , In contrast, an increase in u' enhances the interaction by amplifying both curvature fluctuation and tangential strain rate, leading to increased local reactivity (increased I 0 ). Finally, based on the DNS data, several new scaling models are proposed for the three global properties, S T /S L , A T /A L , and I 0 , and showed improvements compared to existing models. These findings provide new insights into the flame-turbulence interactions in thermodiffusively unstable hydrogen flames. The DNS dataset is also useful for the development of turbulent combustion models applicable to practical engine simulations.

Engine-relevant condition↗
DOE OSTI2022

MODELING AND PARAMETRIC STUDY OF END-GAS AUTOIGNITION TO ALLOW THE REALIZATION OF ULTRA-LOW EMISSIONS, HIGH-EFFICIENCY HEAVY-DUTY SPARK-IGNITED NATURAL GAS ENGINES

Engine knock and misfire are barriers to pathways leading to high-efficiency Spark-Ignited (SI) Natural Gas (NG) engines. The general tendency to knock is highly dependent on engine operating conditions and the fuel reactivity. The problem is further complicated by the low emission limits and the wide range of chemical reactivity in pipeline-quality natural gas. Depending on the region and the source of the natural gas, its reactivity, described by its Methane Number (MN), which is analogous to the Octane Number for liquid SI fuels, can span from 65 to 95. In order to realize diesel-like efficiencies, SI NG engines must be designed to operate at high Brake Mean Effective Pressures (BMEP), near or beyond knock limits, over a wide range of fuel reactivity. This requires a deep understanding of the combustion-engine interactions pertaining to flame propagation and End-Gas Autoignition (EGAI), i.e., the autoignition of the unburned gas (end gas) ahead of the flame front. However, EGAI, if controlled, provides an opportunity to increase SI NG engine efficiency by increasing the combustion rate and the total fraction of burned fuel, mitigating the effects of the slow flame speeds characteristic of natural gas fuels, which generally reduce BMEP and increase unburned hydrocarbon emissions. For this reason, to realize diesel-like efficiencies and ultra-low emissions on SI NG engines, this work proposes the study of the main parameters influencing the modeling and prediction of NG EGAI to allow for its control. In this work, a novel EGAI detection and onset determination method was developed to reliably quantify EGAI for data analysis and engine control. The new method allowed the prediction of EGAI on SI NG engines without the need to use engine- and operating-condition-dependent thresholds and reduced the error in quantifying the fraction of the total energy released by the EGAI event by up to 40%pts. One- and three-dimensional engine models were then developed to study the engine/fuel interactions that lead to NG EGAI and its performance benefits. These models, although having decent agreement with experimental data, showed the need to account for NOx chemistry when predicting NG EGAI due to a consistently later prediction of the EGAI onset (~1.65 crank-angle degrees) and thus, a new reduced chemical mechanism for real NG fuels was developed containing NOx chemistry. The new reduced mechanism improved the EGAI onset prediction agreement to within ±0.5 crank-angle degrees and decreased simulation time during combustion by nearly 50% when using the further reduced AREIS50NOx chemical mechanism. These models were then used to study the role of NG composition on EGAI, evaluate the engine/fuel interactions leading to NG EGAI, and perform engine optimization while leveraging EGAI to increase thermal efficiency. Piston design optimization combined with a Controlled EGAI (C-EGAI) combustion mode allowed a Heavy-Duty (HD) SI NG engine to operate at diesel-like efficiencies, i.e., Brake Thermal Efficiency (BTE) ≥44%. Experimental and modeling data analysis revealed that earlier and faster heat release increases combustion efficiency by an average of 1%pts, increases work transferred to the piston resulting in a decrease in exhaust losses by 50% depending on the engine operating condition while slightly increasing heat losses. Finally, the simulation results revealed an opportunity to further enhance the BTE (up to 50%) by enabling C-EGAI combustion at leaner conditions, λ=1.4-1.6.

Bestel, Diego Bernardi↗
DOE OSTI2023

Evaluation of Spray and Combustion Models for Simulating Dilute Combustion in a Direct-Injection Spark-Ignition Engine

Dilute combustion in spark-ignition engines has the potential to improve thermal efficiency by mitigating knock and by reducing throttling and wall heat losses. However, ignition and combustion processes can become unstable for dilute operation due to a lowered laminar flame speed, resulting in excessive cycle-to-cycle variability (CCV) of the combustion process. To compensate for the slower combustion in less reactive mixtures, a modified intake port geometry can be employed to generate a strong tumble flow in the cylinder and elevate turbulence levels around the spark plug, thereby promoting a faster transition to turbulent deflagration. Consequently, optimizing combustion chamber geometry and operating strategy is crucial to maximizing the benefits of using dilute combustion with enhanced in-cylinder turbulence across a wide range of operating conditions. Computational fluid dynamics (CFD) simulations can be utilized for virtual engine optimization tasks, but this would require the models to be truly predictive regarding the impact of changes to the engine design and operational parameters.In this study, multicycle large-eddy simulations (LES) are performed for a direct-injection spark-ignition engine to investigate the model performance in predicting engine combustion characteristics with respect to changes in the intake configuration. A tumble plate that blocks the lower part of the intake port inlet is used to vary the tumble. A set of CFD models that have been recently developed are employed, which takes into account the drag of nonspherical droplets, flash-boiling behavior of liquid sprays, spray-wall interaction, surrogate formulation of a research-grade E10 gasoline, and fast chemical kinetic solvers. Simulation results are compared to experimental engine data in terms of cylinder pressure, apparent heat release rate, mass fraction burned timing, and flame images. It is found that LES employing the state-of-the-art CFD models are capable of properly predicting the spray processes and reproducing the measured mean cylinder pressure for the case with the tumble plate. On the other hand, the LES over-predicts the combustion rate during the early combustion stage and under-estimates the CCV, and these discrepancies become larger when the tumble plate is removed.

computational fluid dynamics simulation↗
DOE OSTI2021

Efficient connection of reactive shock wave profiles to reaction models

The future of energetic materials development, especially high-performance high explosives, demands new efficient and cost-effective experimental methods to understand the performance of new materials when only limited, development-scale quantities are available. Indeed, perhaps thousands of new molecules that have promising thermodynamic properties have been synthesized, only to be documented in the literature and never evaluated for their potential due to the cost of scaling up to produce quantities needed for traditional performance measurements.

36 MATERIALS SCIENCE↗