Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reactivity Coefficients”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Effect of CO 2 -brine-rock reactions on pore architecture and permeability in dolostone: Implications for CO 2 storage and EOR

Geologic carbon sequestration (GCS) is considered a feasible technology for storing substantive volumes of greenhouse gases in subsurface geological formations. In the reservoir, far from carbon dioxide (CO 2 ) injection wells or in post-injection scenarios, diffusion dominates over advection. This condition conjoins with spatially distributed geochemical reactions to induce heterogeneous changes in pore architecture, i.e. pore body and throat sizes or surface roughness. These changes can affect CO 2 transport properties and storage capacity. In this work, we investigated mineral dissolution and precipitation in dolomite samples saturated with a CO 2 -saturated brine at 93 °C and 34.5 MPa, aged without flow. Two rock types samples, i.e. intergranular- and vuggy-dominant, were selected to investigate changes in pore size, porosity and permeability under reactive conditions. Mineral dissolution and precipitation were characterized using scanning electron microscopy. Changes in pore size were quantified via time-domain nuclear magnetic resonance (TD-NMR) transverse relaxation time (T 2 ) and diffusion coefficient (D) distributions. We show that mineral dissolution likely occurs in highly permeable pathways. These observations are confirmed through analysis of (T 2 ) and diffusion coefficient (D) distributions. In contrast to results during CO 2 -enriched brine continuous injection, mineral precipitation was observed in micropores. The leftward shift of the T 2 peaks, corresponding to micropores, also evidenced mineral precipitation in lowpermeability zones. However, microscale alterations resulted only in a subtle increase in porosity and permeability. Results in this study shed light on effects of geochemical reactions on alteration of rock properties in diffusion-dominated regions during CO 2 storage.

58 GEOSCIENCES↗
DOE OSTI2026

Pulse radiolysis and transient absorption spectra of aqueous solutions of sodium sulfamate

Chemical kinetics for the reactions of sulfamate ions (NH 2 SO 3 − ) with the primary products of water and nitric acid radiolysis were measured in aqueous solutions at ambient temperature. Using time-resolved electron pulse radiolysis techniques with a custom multichannel detection system, we examined the reactivity of NH 2 SO 3 − with the hydroxyl radical ( • OH), hydrogen atom (H • ), and nitrate radical (NO 3 • ). The sulfamate ion was found to react with • OH and H • with second-order rate coefficients of k • OH = (5.60 ± 0.04) × 10 6 M −1 s −1 and k H • = (7.96 ± 0.10) × 10 6 M −1 s −1 , respectively, and with NO 3 • with a rate coefficient of k NO 3 • = (1.67 ± 0.06) × 10 7 M −1 s −1 . The reactions of NH 2 SO 3 − with • OH and H • resulted in the formation of two transient radical species, one with maximum absorbance at 300 nm and a second with maxima at both 300 nm and 600 nm. These spectra are tentatively assigned to • NH 2 SO 3 and • NHSO 3 − , respectively. By measuring the absorbance of these radicals as a function of pH, the radical pK a was determined to be 9.5 ± 0.1. Overall, this work has implications for the longevity and performance of ferrous sulfamate, Fe(NH 2 SO 3 ) 2 , as a plutonium reductant in the reprocessing of used nuclear fuel.

Conrad, Jacy K. [Idaho National Laboratory (INL), ↗
DOE OSTI2022

On the effect of scalar flux weighing of linearly anisotropic scattering matrices in few-group transport calculations

The majority of the codes available for homogenized group constant generation for deterministic transport calculations apply the approximation of scalar flux weighting during energy group condensation of higher-order anisotropic scattering matrices. In this paper, we discuss the effect of scalar flux weighting of the linearly anisotropic scattering matrices in the frame of S P{sub 3} and S{sub 12} calculations performed for a two-dimensional VVER-440 reactor benchmark. To compare group constants generated for two neutron energy groups, an infinite pin cell was homogenized with Serpent 2 and ERANOS ECCO. Serpent 2 applies scalar flux, while ERANOS ECCO performs current weighting of the linearly anisotropic scattering matrices during energy group condensation. For analyzing the effect of the various weighting options, three simple reactor models were built assuming different core sizes using standard rectangular assemblies with 15*15 fuel pins. Diffusion, SP{sub 3} and S{sub 12} calculations were performed for the 3 models using group constants generated with Serpent 2 and ERANOS ECCO. The effect of scalar flux weighting of linearly anisotropic scattering matrices is shown by comparing the decrease in reactivity due to the decreased reactor size, as well as the assembly power distribution to reference results obtained with Serpent 2 Monte Carlo calculations. Neglecting higher than linearly anisotropic scattering and indirect application of diffusion coefficients in higher-order transport calculations is advised if angular flux-moment spectra weighted higher-order scattering matrices cannot be generated. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Multiphase Processing of the Water-Soluble and Insoluble Phases of Biomass Burning Organic Aerosol

Biomass burning is one of the most significant sources of organic aerosol in the atmosphere. Biomass burning organic aerosol (BBOA) has been observed to undergo liquid– liquid phase separation (LLPS) to give core–shell morphology with the hydrophobic phase encapsulating the hydrophilic phase, potentially impacting the evolution of light-absorbing components, i.e., brown carbon (BrC), through multiphase processes. Here, we demonstrate how multiphase processing differs between the watersoluble (i.e., hydrophilic) and insoluble (i.e., hydrophobic) phases of BBOA in terms of reactive uptake of ozone in a coated-wall flow tube. Effects of relative humidity (RH) and ultraviolet (UV) irradiation were investigated. Experimental timeseries were used to inform simulations using multilayer kinetic modeling. Among non-irradiated thin films, the uptake coefficient was greatest for the water-soluble phase at 75% RH (3 × 10 –5 , corresponding to a diffusion coefficient of BrC, D BrC , of 3 × 10 –9 cm 2 s –1 ) and least for the same phase at 0% RH (1 × 10 –5 , corresponding to D BrC of 1 × 10 –10 cm 2 s –1 ). The uptake coefficient for the water-insoluble phase fell between these two (about 1.5 × 10 –5 ), regardless of RH, and the corresponding D BrC increased only slightly (8 × 10 –10 cm 2 s –1 at 0% RH to 9 × 10 –10 cm 2 s –1 at 75% RH). The uptake coefficients of both phases at 0% RH decreased significantly after UV irradiation, consistent with a transition from viscous liquid to solid and supported by qualitative microscopy observations. Modeling multiphase ozone oxidation of primary BrC components in the atmosphere demonstrated, first, that LLPS may extend the lifetime of water-soluble BBOA encapsulated by water-insoluble species by a factor of 1.5 at moderate to high RH and, also, that UV irradiation may extend the lifetime of both phases by more than a factor of 2.5.

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

Data-Driven Chance-Constrained Design of Voltage Droop Control for Distribution Networks: Preprint

This paper addresses the design of local control methods for voltage control in distribution networks with high level of distributed energy resources (DERs). The designed control methods adapt the active and reactive power output of distributed energy resources proportional to the deviation of the local measured voltage magnitudes from a reference voltage, which is referred to as droop control. Thus, the design focuses on determining the droop characteristics which satisfy network-wide voltage magnitude constraints. The uncertainty and variability of DERs renders the design of optimal droop controls very challenging. Hence, this paper proposes chance constraints to limit the risk from intermittent DERs, by designing droop control coefficients that guarantee the satisfaction of network operational constraints with a specific probability. In addition, the proposed approach relies entirely on historical data rather than assuming knowledge of the probability distributions that characterize the uncertainty of DERs. The efficacy of the proposed method is demonstrated on a 37-bus distribution feeder.

chance-constrained optimization↗
DOE OSTI2022

Investigating the kinetic effect of prenol on iso -octane auto-ignition by means of an experimental and modeling study

A detailed experimental and kinetic modeling study was dedicated to understand the reported octane hyperboosting effect of prenol, by means of the measurement of the ignition delay times of its blends with iso -octane, and measurement of the mole fraction profiles of the fuels and intermediates inside the ULille rapid compression machine. Here, these results show that prenol addition leads to a reduction of the first-stage ignition phenomena and negative temperature coefficient behavior, which is only qualitatively captured by the model and is consistent with knock resistance improvement. It is suggested that this behavior is caused by two different factors. The first originates from gas-phase reactivity of prenol, and spans from the formation of unreactive unsaturated species through resonance-stabilized radicals, thereby constituting a competitive pathway for the radical pool generated by iso -octane. The second is of catalytic nature and cannot be captured by means of gas-phase kinetic modeling, but could also play an important role in the behavior of prenol in internal combustion engines.

42 ENGINEERING↗
DOE OSTI2024

Applying a Phase-Separation Parameterization in Modeling Secondary Organic Aerosol Formation from Acid-Driven Reactive Uptake of Isoprene Epoxydiols under Humid Conditions

Secondary organic aerosol (SOA) from acid-driven reactive uptake of isoprene epoxydiols (IEPOX) contributes up to 40% of organic aerosol (OA) mass in fine particulate matter. Our previous work showed substantial conversions of particulate inorganic sulfates to surface-active organosulfates (OSs) by IEPOX decreases aerosol acidity and creates a viscous organic-rich shell that poses as a diffusion barrier, inhibiting additional reactive uptake of IEPOX. In this study, to account for this “self-limiting” effect, a phase-separation box model was developed to evaluate parameterizations of IEPOX reactive uptake against time-resolved chamber measurements of IEPOX-SOA tracers, including 2-methyltetrols (2-MT) and methyltetrol sulfates (MTS), at ~ 50% relative humidity. The phase-separation model was most sensitive to mass accommodation coefficient, IEPOX diffusivity in the organic shell, and ratio of the third-order reaction rate constants forming 2-MT and MTS (k MT /k MTS ). In particular, k MT /k MTS had to be lower than 0.1 to bring model predictions of 2-MT and MTS in closer agreement with chamber measurements, while prior studies reported values larger than 0.71. The model-derived rate constants favor more particulate MTS formation due to 2-MT likely off-gassing at ambient-relevant OA loadings. Incorporating this parametrization into chemical transport models is expected to predict lower IEPOX-SOA mass and volatility due to the predominance of OSs.

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

Tantalum pentoxide: a new material platform for high-performance dielectric metasurface optics in the ultraviolet and visible region

Dielectric metasurfaces, composed of planar arrays of subwavelength dielectric structures that collectively mimic the operation of conventional bulk optical elements, have revolutionized the field of optics by their potential in constructing high-efficiency and multi-functional optoelectronic systems on chip. The performance of a dielectric metasurface is largely determined by its constituent material, which is highly desired to have a high refractive index, low optical loss and wide bandgap, and at the same time, be fabrication friendly. Here, we present a new material platform based on tantalum pentoxide (Ta 2 O 5 ) for implementing high-performance dielectric metasurface optics over the ultraviolet and visible spectral region. This wide-bandgap dielectric, exhibiting a high refractive index exceeding 2.1 and negligible extinction coefficient across a broad spectrum, can be easily deposited over large areas with good quality using straightforward physical vapor deposition, and patterned into high-aspect-ratio subwavelength nanostructures through commonly-available fluorine-gas-based reactive ion etching. We implement a series of high-efficiency ultraviolet and visible metasurfaces with representative light-field modulation functionalities including polarization-independent high-numerical-aperture lensing, spin-selective hologram projection, and vivid structural color generation, and the devices exhibit operational efficiencies up to 80%. Our work overcomes limitations faced by scalability of commonly-employed metasurface dielectrics and their operation into the visible and ultraviolet spectral range, and provides a novel route towards realization of high-performance, robust and foundry-manufacturable metasurface optics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Energy-Storage Fed Smart Inverters for Mitigation of Voltage Fluctuations in Islanded Microgrids

The continuous integration of intermittent low-carbon energy resources makes islanded microgrids vulnerable to voltage fluctuations. Besides, different dynamic response of synchronous-based and inverter-based distributed generation (DG) units can result in an instantaneous power imbalance between supply and demand during transients. As a result, the ac-bus voltage of microgrid starts oscillating which might have severe consequences such as blackouts. This paper modifies the conventional control scheme of battery energy storage systems (BESSs) to participate in improving the dynamic behavior of islanded microgrids by mitigating the voltage fluctuations. A piecewise linear-elliptic (PLE) droop is proposed and employed in BESS to achieve an enhanced voltage profile by injecting/absorbing reactive power during transients. In this way, the conventional inverter implemented in BESS turns into a smart inverter to cope with fast transients. Using the proposed approach in this paper, any linear droop curve with a specified coefficient can be replaced by a PLE droop curve. Compared with linear droop, an enhanced dynamic response is achieved by utilizing the proposed PLE droop. Case study results are presented using PSCAD/EMTDC to demonstrate the superiority of the proposed approach in improving the dynamic behavior of islanded microgrids.

Pilehvar, Mohsen S.↗
DOE OSTI2020

Assessment of Existing Transportation Packages for Use with HALEU

Commercial light water reactor operators and fuel vendors in the United States are pursuing changes to fuel that include increased 235 U enrichment. Economic studies generally anticipate maximum near-term fuel assembly designs with up to 8 wt.% 235 U. Many next-generation nuclear reactor designs require high-assay low-enriched uranium (HALEU) (19.75 wt.% > 235 U > 5 wt.%) fuel. One necessary element for the commercial-scale use of HALEU is the ability to safely transport large quantities of enriched fuel material in multiple forms. However, there is uncertainty as to whether subcriticality requirements can be satisfied with existing package designs and whether existing critical benchmark experiment data are sufficient to support criticality safety code validation for HALEU transportation applications. This study assesses the potential to use currently licensed transportation packages for the transportation of increased enrichment unirradiated U fuel forms. The assessment uses selected package designs that represent the five categories of fuel form-boiling water reactor pins and assemblies, pressurized water reactor pins and assemblies, UF 6 , U-metal and tristructural isotropic (TRISO) particles, and UO 2 pellets or powder-and focuses on demonstrating subcriticality and identifying benchmark critical experiments appropriate for use in criticality computer code validation. Key quantities of interest that relate to subcriticality are limiting conditions (e.g., optimum moderation), package or package array k eff , package capacity, and package transportation array size. The SCALE TSUNAMI-3D and TSUNAMI Indices and Parameters codes are used for sensitivity and uncertainty calculations and for identification of candidate critical benchmark experiments for code validation. The primary metric for identifying candidate benchmarks is the similarity coefficient, c k . For each fuel form category, a representative package is evaluated. Results provided for each package evaluation include enrichment and packaging limits (e.g., maximum transportation array size as a function of enrichment) and benchmark critical experiment similarity coefficients. Results indicate that there are viable means for increasing enrichments into the HALEU range across the spectrum of fuel forms with differing increase amounts available for different packages. Sources of subcriticality margin to offset increased enrichment reactivity include reduced transportation array size, reduced fissile mass, burnable absorber credit, and safety analysis margin harvesting. For all packages except the DN-30, numerous critical benchmark experiment candidates for validation were identified.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Electrodeposition of Magnonic V(tetracyanoethylene) 2 Thin Films

Molecule-based magnetic materials have been identified as promising candidates for application in magnonic technologies, owing not only to their solution processability but also because they can exhibit narrow ferromagnetic resonance (FMR) linewidths and low Gilbert damping coefficients—crucial prerequisites for the transmission of coherent magnons over macroscopic distances. In particular, V(TCNE) 2 , a compound with a threedimensional network structure composed of vanadium(II) centers linked by tetracyanoethylene (TCNE •− ) radical anions, displays magnonic properties comparable to yttrium iron garnet, the quintessential magnonic material in the field. However, existing solution and chemical vapor deposition methods for synthesizing V(TCNE) 2 require the use of highly reactive zero-valent molecular vanadium precursors, stymying research on this important material. Herein, we report a facile electrochemical method for the deposition of thin films of V(TCNE) 2 using readily obtainable and stable divalent vanadium precursors and TCNE •− anions generated by electrochemical reduction. Magnetization measurements reveal that the films exhibit ferrimagnetic ordering above room temperature, consistent with V(TCNE) 2 films synthesized via other methods. Moreover, the electrodeposited films exhibit narrow FMR linewidths as low as 17.5 G and a low Gilbert damping coefficient of 1.1 × 10 −3 , values that are on par with some currently integrated metallic magnonic materials. More generally, these results demonstrate that electrodeposition can provide a straightforward means of generating highperformance magnonic materials using readily available molecular precursors.

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

Atomistic thermodynamics and kinetics of dicalcium silicate dissolution

Low-temperature cement manufacturing has garnered academic and industrial attention for its low environmental footprints. However, the sluggish hydration kinetics of the resultant cement affects their early-age strength development. This motivates fundamental studies to unravel the mechanistic picture of the dissolution process and discover science-informed pathways to accelerate hydration. Standard atomistic simulations seldomly exceed a microsecond making them impractical to study slow dissolution processes. We are using rare event sampling techniques, we provide the mechanistic picture of Ca 2+ ion dissolution from a kink site on the dicalcium silicate surface. The Ca 2+ ion dissolution is comprised of two sequential stages: breaking restraints from the kink sites to form a ledge adatom and detaching from the ledge/terrace adatom sites into the solution. The first and second stages feature free energy barriers of ~63 kJ/mol and ~ 29 kJ/mol respectively, making the first stage the rate-limiting step of the entire Ca 2+ dissolution kinetics. Using the reactive flux method, the rate and equilibrium constants for each reaction step are calculated, which yield the Ca 2+ ion activity of ~1.03 × 10 -5 . The diffusion calculations indicate that the surface effects lower the self-diffusion coefficient of Ca 2+ ions at the solid-water interface.

36 MATERIALS SCIENCE↗
DOE OSTI2024

COVID-19–Related School Closures, United States, July 27, 2020–June 30, 2022

As part of a multiyear project that monitored illness-related school closures, we conducted systematic daily online searches during July 27, 2020–June 30, 2022, to identify public announcements of COVID-19–related school closures (COVID-SCs) in the United States lasting ≥1 day. We explored the temporospatial patterns of COVID-SCs and analyzed associations between COVID-SCs and national COVID-19 surveillance data. COVID-SCs reflected national surveillance data: correlation was highest between COVID-SCs and both new PCR test positivity (correlation coefficient [r] = 0.73, 95% CI 0.56–0.84) and new cases (r = 0.72, 95% CI 0.54–0.83) during 2020–21 and with hospitalization rates among all ages (r = 0.81, 95% CI 0.67–0.89) during 2021–22. The numbers of reactive COVID-SCs during 2020–21 and 2021–22 greatly exceeded previously observed numbers of illness-related reactive school closures in the United States, notably being nearly 5-fold greater than reactive closures observed during the 2009 influenza (H1N1) pandemic.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Solvent selection for biphasic extraction of 5-hydroxymethylfurfural via multiscale modeling and experiments

We introduce a comprehensive conceptual framework for selecting solvents for reactive extraction in biphasic organic-water systems and demonstrate it for the separation of HMF (5-hydroxylmethylfurfural), a platform chemical produced in the acid-catalyzed dehydration of hexoses. We first perform in silico screening of ~2500 solvents, from the ADFCRS-2018 database using the ADF COSMO-RS implementation, and classification, based on the solvent partition coefficient. We then determine experimentally the partition coefficients for HMF, fructose, and products of HMF rehydration (levulinic acid (LA), and formic acid (FA)), the mutual water-organic solvent solubilities, and the separation factors in >50 select solvents spanning multiple homologous series at room temperature and a typical reaction temperature with in situ sampling. We find that COSMO-RS is excellent for screening purposes (typical error in most cases within a factor of ~2). Increased temperatures lead to significant reduction in partitioning, and room temperature measurements are clearly inadequate for solvent selection. Upon down selecting classes of solvents based on separation performance, we perform experimental thermal stability and reaction compatibility studies of a small set of solvents at relevant reactive-extraction temperatures. We discover that many substituted phenols exhibit an order-of-magnitude increase in partitioning compared to conventional solvents due chiefly to hydrogen bond interactions and show the necessary stability but retain a significant fraction of water and LA, factors that need to be considered in technoeconomic analysis. In contrast, anilines, aldehydes, and acids are good to excellent regarding separation but incompatible with this specific reaction medium. Furthermore, this multifaceted framework can be extended to other biomass-derived products and processes.

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

Direct Demonstration of Unified Brønsted-Evans-Polanyi Relationships for Proton-Coupled Electron Transfer Reactions on Transition Metal Surfaces

Brønsted-Evans-Polanyi (BEP) relationships, which relate elementary reaction barriers to reaction thermodynamics, have long been discussed in electrochemical science, but demonstration of their existence in heterogeneous electrocatalysis, across many different catalyst surfaces and voltages, is lacking. In this contribution, the BEP model is demonstrated to describe both multiple catalyst surfaces and variable voltages in electrochemical environments. For proton-coupled electron transfer (PCET) to surface nitrogen (N*) and nitric oxide (NO*), reaction energies and activation barriers are calculated using Density Functional Theory (DFT) on a parallel plate capacitor model at three different potentials and for nine different transition metal surfaces. Linear BEP relations that describe all potentials and catalyst surfaces are obtained for these elementary reactions, and the slopes of the correlations are shown to be directly related to the fractional coordinate of the transition states (FCTS) of the reactions. The results, which are explained in terms of Marcus Theory, prove a direct equivalence between unified BEP coefficients, describing both variable catalyst surfaces and voltages, and electrochemical symmetry factors and provide a straightforward means of estimating this quantify for PCET reactions on transition metal surfaces. The resulting relationships, in turn, could lead to predictions of electrocatalytic reactivity trends of enhanced accuracy and efficiency.

Electrochemistry↗
DOE OSTI2022

Replicating HCCI-like autoignition behavior: What gasoline surrogate fidelity is needed?

This work seeks to characterize the fidelity needed in a gasoline surrogate with the intent to replicate the complex autoignition behavior exhibited within advanced combustion engines, and specifically Homogeneous Charge Compression Ignition (HCCI). A low-temperature gasoline combustion (LGTC) engine operating in HCCI mode and a rapid compression machine (RCM) are utilized to experimentally quantify fuel reactivity, through autoignition and preliminary heat release characteristics. Fuels considered include a research grade E10 U.S. gasoline (RD5-87), three multi-component surrogates (PACE-1, PACE-8, PACE-20), and a binary surrogate (PRF88.4). Each fuel was studied at lean/HCCI-like conditions covering a wide range of temperatures and pressures that are representative of naturally aspirated to high boost engine operation. Detailed chemical kinetic modeling is also undertaken using a recently updated gasoline surrogate kinetic model to simulate the RCM experiments and to provide chemical insight into surrogate-to-surrogate differences. The LGTC engine experiments demonstrate nearly identical reactivity between PACE-20 and RD5-87 across studied conditions, while faster phasing is seen for both PACE-1 and PACE-8 due to their stronger intermediate- and low-temperature heat release (ITHR/LTHR) at naturally aspirated and boosted conditions, respectively. The RCM experiments reveal typical low-temperature, negative temperature coefficient (NTC) and intermediate-temperature autoignition behaviors at all pressure conditions for RD5-87, which are qualitatively reproduced by all surrogates. Quantitative discrepancies in both autoignition and preliminary heat release are observed for all surrogates, while their ability to replicate RD5-87 autoignition behavior follows the order of PACE-20 > PACE-1 > PACE-8 > PRF88.4. Excellent mapping is obtained between the LGTC engine and the RCM, where the engine pressure-time trajectories can be characterized by the regimes represented by the RCM autoignition isopleths. The kinetic model performs commendably when simulating both autoignition and preliminary heat release of PACE-20, while typically overpredicting ignition delay times for PACE-1, PACE-8 and PRF88.4 at high-pressure and low-temperature/NTC conditions. Sensitivity and rate of production (ROP) analyses highlight surrogate-to-surrogate differences in the governing chemical kinetics where n-pentane initiates rapid OH branching at a faster rate and an earlier timing for PACE-20 than iso-pentane does for PACE-1 and PACE-8, making it computationally more reactive than the other surrogates. The current study highlights the need to include non-standardized properties, such as the lean/HCCI-like autoignition characteristics, in addition to ASTM properties (e.g., RON, MON) as metrics of fuel reactivity and targets to be matched when formulating high-fidelity surrogates that fully capture gasoline advanced combustion behavior such as HCCI-like autoignition.

42 ENGINEERING↗
DOE OSTI2024

Integration and validation of some modules for modelling of high-speed chemically reactive flows in two-phase gas-droplet mixtures

Three modules are integrated into the built-in OpenFOAM rhoCentralFoam solver towards accurate and efficient modelling of high-speed chemically reactive flows in two-phase gas-droplet mixtures within the OpenFOAM 10.0 framework. The first module is the mixture-averaged diffusion model. The second module is the built-in OpenFOAM Lagrangian solver coupled with optimised droplet drag coefficient and convective heat transfer coefficient sub-models. The last module is a sparse stiff chemistry solver based on dynamic adaptive hybrid integration (AHI-S). The optimised droplet sub-models are first verified in correct implementation for subsequent simulations in this work. Further, they show good accuracy against experimental and analytical data in the modelling of ammonia droplet acceleration and cooling in the flowing and/or low-temperature air. The accuracy and efficiency gains related to the mixture-averaged diffusion model and the AHI-S chemistry solver are examined by simulating 1-D detonation propagation in ammonia droplet-free/laden ammoniaoxygen mixtures. Numerical results of detonation propagation speed, gaseous temperature, density, and species distributions around the induction zone show good agreement with experimental data and analytical solutions. Compared to the built-in OpenFOAM diffusion model, the mixture-averaged diffusion model provides different numerical predictions of pulsating instabilities in detonation propagation. It shows better accuracy in depicting the detonation structure within the droplet-free section attributed to improved multi-component diffusion modelling. Compared to the built-in OpenFOAM solver EulerImplicit (backward Euler), the AHI-S chemistry solver reduces the computational cost by around 50%. It achieves satisfactory accuracy in calculating detonation propagation speed within the droplet-free section with the optimal efficiency when the safety factor, β, equals 0.5.

42 ENGINEERING↗
DOE OSTI2026

Ab Initio Modeling of Aqueous Methanol Mixtures at DFT-SCAN Level Using Machine Learning Interatomic Potentials

Abstract Methanol–water mixtures find use in many applications, particularly catalytic energy conversion processes. Their importance has motivated numerous computational studies, most of which employed molecular dynamics based on classical force fields. These enable simulations of large systems on long time scales but do not reliably describe reactive dynamics involving bond breaking and bond formation. In contrast, ab initio molecular dynamics (AIMD) based on density functional theory (DFT) is generally more reliable for such applications but has a high computational cost, which discourages systematic studies of alcohol-water mixtures. To remedy this, we trained a machine learning interatomic potential capable of probing the properties of aqueous methanol mixtures at the DFT level using the SCAN functional. Our results show that SCAN qualitatively reproduces multiple key experimental features arising from the amphiphilic nature of methanol, including density, diffusion coefficients, X-ray structure factors, and Kirkwood–Buff integrals. We also find that structural correlations between water molecules are somewhat overestimated, leading to a stronger preferential association than that predicted by experiments. However, increasing the temperature by 30 K mitigates this effect and also recovers the correct mobilities of both methanol and water. These results indicate that SCAN provides an accurate description of methanol–water mixtures, making it a reliable choice for investigating the reactive dynamics in such systems.

Park, Sanghyun J. [Princeton University , , , ,]↗