Search NASA⌕ Search

SEARCH · Search NASA

Results for “reactivation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Coupled geomechanical investigation of depletion-induced fault reactivation

Fault reactivation during subsurface fluid production pose significant challenges to safe and sustainable resource extraction. Here, this study presents a three-dimensional coupled geomechanical framework to investigate the processes driving fault reactivation, capturing the interactions between reservoir dynamics and geomechanical responses. Verification against theoretical estimations based on linear poroelasticity confirms the model's capacity in representing reservoir background stress responses. However, the study reveals that relying solely on background stress states can underestimate or overestimate fault reactivation potential, emphasizing the importance of including localized stress perturbations such as differential compaction and stress redistribution. Applied to a fault (M1) inspired by the geological characteristics of the Groningen field, the model shows slip initiation at 2965 m depth with 16.0 MPa depletion, aligning with field observations where seismicity occurred at approximately 3 km depth after 15.8 MPa depletion. Parametric studies reveal: (1) inelastic reservoir compaction delays fault reactivation and mitigates fault slip by reducing stress concentration, (2) higher intermediate in-situ stress magnitudes decrease the Coulomb Failure Stress (CFS) increase rate and reduce fault slip, (3) larger fault offsets amplify shear stress near the offset zone, promoting earlier reactivation and longer rupture propagation, and (4) fault permeability significantly influences pressure diffusion, with low-permeability faults leading to sharper stress changes and earlier fault destabilization. These insights highlight the critical role of geological and mechanical parameters in fault reactivation and provide a predictive framework for mitigating induced seismicity risks.

Coupled geomechanical modeling↗
DOE OSTI2020

Different Kinetic Reactivities of Electrons in Distinct TiO 2 Nanoparticle Trap States

Electrons added to TiO 2 and other semiconductors often occupy trap states, whose reactivity can determine the catalytic and stoichiometric chemistry of the material. We previously showed that reduced aqueous colloidal TiO 2 nanoparticles have two distinct classes of thermally equilibrated trapped electrons, termed Red/e – and Blue/e – . Presented here are parallel optical and electron paramagnetic resonance (EPR) kinetic studies of the reactivity of these electrons with solution-based oxidants. Optical stopped-flow measurements monitoring the reactions of TiO 2 /e – with substoichiometric oxidants showed a surprising pattern: an initial fast (seconds) decrease in TiO 2 /e – absorbance followed by a secondary, slow (minutes) increase in the broad TiO 2 /e – optical feature. The analysis revealed that the fast decrease is due to the preferential oxidation of the Red/e – trap states and the slow increase results from the re-equilibration of electrons from Blue/e - to Red/e - states. This kinetic model was confirmed by freeze-quench EPR measurements. Quantitative analysis of the kinetic data demonstrated that Red/e – react ~5 times faster than Blue/e – with the nitroxyl radical oxidant 4-methoxy-2,2,6,6-tetramethyl-1-piperidinyloxyl (4-MeO-TEMPO). Similar reactivity patterns were also observed in oxidations of TiO 2 /e – by O 2 , which like 4-MeO-TEMPO is a proton-coupled electron transfer (PCET) oxidant, and by the pure electron transfer (ET) oxidant potassium triiodide (KI 3 ). Furthermore, this suggests that the faster intrinsic reactivity of one trap state over another on the seconds–minutes time scale is likely a general feature of reduced TiO 2 reactivity. This differential trap-state reactivity is likely to influence the performance of TiO 2 in photochemical/electrochemical devices, and it suggests an opportunity for tuning catalysis.

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

ZND Verification Tests for Reactive Burn Models in FLAG

The ZND theory, named after Zeldovich, von Neumann, and Doering, provides a simple model for one-dimensional ideal steady-state detonation. It assumes that the detonation wave front starts with a shock that is a discontinuous jump and is followed by a finite-length reaction zone. Reactive burn (also called reactive flow) models are based on ZND theory, as they model the shock initiation and detonation process with a finite reaction rate. The ZND wave propagation test is essentially the only available test case where an analytic solution exists for verification of reactive burn models in numerical codes. However, there are extensions and variants of the ZND test that have been devised for verification of multidimensional flows. The objective of this work is to provide verification of the reactive burn models currently implemented in the Lagrangian hydrocode FLAG and investigate the influence of mesh resolution, artificial viscosity models, and the Arbitrary Lagrangian-Eulerian (ALE) Euler relaxer on the simulation results. The burn models of interest are the Wescott-Stewart-Davis (WSD) model, the Scaled Uniform Reactive Flow (SURF) model (specifically with the SURFplus model extension), and the Arrhenius shock temperature state dependent WSD (AWSD) model. Previously, Ralph Menikoff has used ZND tests for verification of the SURF and SURFplus models in the Eulerian hydrocode xRAGE. The ZND tests here are somewhat different than the approach by Menikoff. In particular, we use a piston-driven ZND detonation wave (via a prescribed constant velocity boundary condition) in a Lagrangian framework whereas Menikoff had a ZND wave followed by a invariant rarefaction wave. The xRAGE simulations were carried out on uniform grids and adaptive mesh refinement (AMR) grids. Although AMR was recently implemented in FLAG and now fully functional for 2D simulations, it will not be evaluated with ZND tests at this time. Some work has been done previously for validating the reactive burn models in FLAG. For example, the validation studies of SURF with shock-to-detonation (SDT) tests, cylinder tests, and gap-stick tests. Further validation of the AWSD, WSD, and SURF models is described in for SDT, multi-shock, cylinder, and corner-turning tests. Recently, a large parameter study with approximately twenty SDT tests was performed to validate the AWSD, WSD, SURF, and SURFplus models while also investigating mesh resolution and artificial viscosity settings. To the best of our knowledge, the current work represents the first documented verification of these burn models in FLAG. The reader should note that many of the tables and figures in this report use units of cm/g/µs, which are the default units for FLAG. However, some lengths and velocities are expressed in µm and mm/µs, respectively, which are typical units for detonation analyses.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2022

Analysis of Reactive Power Load Modeling Techniques for PV Impact Studies [Slides]

The increasing availability of advanced metering infrastructure (AMI) data has led to significant improvements in load modeling accuracy. However, since many AMI devices were installed to facilitate billing practices, few utilities record or store reactive power demand measurements from their AMI. When reactive power measurements are unavailable, simplifying assumptions are often applied for load modeling purposes, such as applying constant power factors to the loads. The objective of this work is to quantify the impact that reactive power load modeling practices can have on distribution system analysis, with a particular focus on evaluating the behaviors of distributed photovoltaic (PV) systems with advanced inverter capabilities. Quasi-static time-series simulations were conducted after applying a variety of reactive power load modeling approaches, and the results were compared to a baseline scenario in which real and reactive power measurements were available at all customer locations on the circuit. Overall, it was observed that applying constant power factors to loads can lead to significant errors when evaluating customer voltage profiles, but that performing per-phase time-series reactive power allocation can be utilized to reduce these errors by about 6x, on average, resulting in more accurate evaluations of advanced inverter functions.

14 SOLAR ENERGY↗
DOE OSTI2021

Characterizing Reactive Transport Behavior in a Three-Dimensional Discrete Fracture Network

While several studies have linked network and in-fracture scale properties to conservative transport behavior in subsurface fractured media, studies on reactive transport cases remain relatively underdeveloped. In this study, we explore the behavior of an irreversible kinetic reaction during the interaction of two solute plumes, one consisting of species A and the other species B. When the plumes converge, these species react kinetically to form a new species C via A+B→kC. This reactive system is studied using a three-dimensional discrete fracture network (DFN) model coupled with reactive Lagrangian particle tracking. We find that the interplay of network topology and chemical properties of the reactive solutes controls reactive transport processes. The network topology drives species A and B together, and the chemical properties dictate whether and how quickly a reaction occurs. Results demonstrate that reactions are most likely to occur in high-velocity fractures that make up the network backbone. The interplay between species’ chemical properties and transport is characterized by a non-dimensional Damköhler (Da) number. We show that the spatial distribution of reactions is sensitive to Da, which subsequently influences late-time tailing behavior in outlet breakthrough time distributions. The results of this study provide initial insights into how an irreversible reaction occurs during transport in a fracture network, using a methodology that can be applied to study reactive transport in a wide range of fractured media environments and contexts.

58 GEOSCIENCES↗
DOE OSTI2020

Impacts of advanced diesel combustion operation and fuel formulation on soot nanostructure and reactivity

Advanced diesel combustion, accomplished via a single pulse fuel injection and high levels of exhaust gas recirculation (referred to as “PCCI”, partially-premixed charge compression ignition), is shown to be a path to reduce oxides of nitrogen and particulate matter simultaneously. This is well established in the literature. Less established is the extent to which such dilute combustion processes influence soot formation and affect soot that is emitted from diesel engines under such combustion modes. This work focuses on characterization of the nanostructure and oxidative reactivity of soot generated by a light-duty turbodiesel engine operating under a PCCI combustion mode, a dilute, low-temperature combustion process. Previous work on a type of PCCI combustion, referred to as high-efficiency clean combustion (HECC), showed soot samples having a fullerenic nanostructure, characterized by high levels of tortuosity of the fringe layers as seen in transmission electron micrograph images and as quantified using an image processing algorithm. Thermogravimetric analysis of the HECC mode soot samples showed that they displayed higher rates of oxidation than soot samples from a conventional diesel combustion mode. The present work returns to PCCI combustion, considering the timing of the main fuel injection, and the effects of operating on fuels rich in n-alkanes, particularly a fuel produced from a low temperature Fischer-Tropsch process (LTFT) and a renewable diesel fuel (RD) produced via hydrodeoxygenation of a plant oil. PCCI combustion conditions yield soot that shows higher reactivity compared to soot from the conventional combustion mode regardless of fuel type. LTFT and RD fuels produce soots with lower reactivity compared to ULSD. Soots produced from PCCI combustion have higher surface oxygen concentration and higher proportion of amorphous carbon. In addition, TEM images show that PCCI soots from all three fuels have smaller primary particle and particle aggregate sizes, and smaller graphene layers. These properties explain the higher reactivity of soot from PCCI combustion. The less reactive soots, which are produced from LTFT and RD fuel under conventional combustion, show internal burning during oxidation. However, soots with higher reactivity which are produced from late injection PCCI combustion and ULSD show shrinking core oxidation, most likely because of their overall amorphous structure.

42 ENGINEERING↗
DOE OSTI2020

Calculation of the First Moment of Energy Using D-T Reactivity Formalisms Under the Maxwell-Boltzmann Distribution (Pt. 1)

Nuclear fusion science is an example of a scientific field with a rich history of expert involvement and scientific publications, which together, form an expert-knowledge base. One example of this history is the utilization of published reaction rates from a variety of authors. Investigators for Deuterium-Tritium (D-T) ion fusion can choose from using frequently cited methods: the Bosch and Hale reactivity, thermonuclear reaction rates from Caughlan and Fowler, and the reactivity evaluation from Miley, Towner & Ivich which forms the basis of the Naval Research Lab (NRL) formulary. There are other choices available. Each of the reactivity formulations considered here, are based upon the Maxwell- Boltzmann velocity distribution for D-T fusion ion reactants. Numerical methods for computer codes simulating hot, energetic plasmas, include tabulations of the reactivity, and the first moment of energy. The purpose of this study is to answer the question: what is the first moment of energy, and how has it been formulated? This report is part of a series of reports by the authors on D-T fusion-reaction formalism. The present focus is on defining the mathematical relationship of the first moment of D-T fusion ion kinetic energy, $\langle$E$\rangle$, with the fusion crosssection, fusion reactivity and its derivative with respect to ion-temperature. Three variants of the first moment $\langle$E$\rangle$ are analytically developed and explored: 1) constant cross-section, 2) a normalized first moment, and 3) a particular function of the first moment. Much of the mathematical development is relegated to the Appendices, for a concise presentation of topics and results. A research component of this discussion is our comparison of some of the reactivity and its derivative results between evaluations (1) versus (2) and (3), which the authors have not found in any other publication to date. This investigation is useful for validation, verification (V&V) and uncertainty quantification (UQ) as it compares the work of several authors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Calculation Of The First Moment Of Energy Using D-T Reactivity Formalisms Under The Maxwell-Boltzmann Distribution--Part II

Nuclear fusion science is an example of a scientific field with a rich history of expert involvement and scientific publications, which together, form an expert-knowledge base. One example of a nuclear fusion formalism is the utilization of published reaction rates from a variety of authors. Investigators for Deuterium- Tritium (D-T) ion fusion can choose from using frequently cited methods: the Bosch and Hal reactivity, thermonuclear reaction rates from Caughlan and Fowler, and the reactivity evaluation from Miley, Towner & Ivich which forms the basis of the Naval Research Lab (NRL) formulary. There are other choices available. Each of the reactivity formulations considered here, are based upon the Maxwell-Boltzmann velocity distribution for D-T fusion ion reactants. Numerical methods for computer codes simulating hot, energetic plasmas, include tabulations of the reactivity, and the first moment of energy. This report continues the step toward building understanding of nuclear fusion reactivity formalisms. It is part of a series of reports with the same goal, [5-10] and is the continuation of the Part I paper for defining the mathematical relationship of the first moment of D-T fusion ion kinetic energy, <$E$>, with the fusion cross-section, fusion reactivity and its derivative with ion-temperature. In Part I, three variants of the first moment <$E$> were analytically developed and explored: 1) constant cross-section, 2) a normalized first moment, and 3) a particular function of the first moment from Brysk. In Part II, attention is given to the definition of <$K$>, originally described as a ratio of moments from Brysk, and its relationship to the first moment definitions from Part I. One measure of the progress made in these documents is the identification that Brysk’s ratio of the second moment to the first moment ratio, <$K$>, does not correspond to his provided solution of the first moment of energy. Another measure of (our) progress from this work is the comparison of first moment variants. That comparison includes confirming the importance of cross sections defined in terms of energy. The analytical relationships we have developed among important physics quantities are useful tools in validation and verification (V&V). For example, we can calculate the kinetic energy as a mean or as a first moment, or as a function of the first moment. These analytically-determined values can be compared directly with numerically-determined values, supplied to the authors, representing <$E$>.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Numerical Studies of the On-Line Reactivity Monitoring Techniques in Subcritical Systems: Revisit of the Area-Ratio Method, the Neutron Source-Modulation Method, and the Flux-to-Current Ratio Method

In the subcritical systems, the subcriticality level must be monitored continuously to ensure the safe operation of the systems. Techniques such as the neutron source-modulation method and the flux-to-current ratio method have been proposed for on-line reactivity measurements of the subcritical systems. However, the accuracy of these methods is unknown in the actual applications. In this report, the fundamental concepts of these techniques are visited, with numerical simulation techniques developed to validate these on-line reactivity monitoring methodologies. Specifically, the Monte Carlo numerical simulation results show that the neutron source-modulation method is not accurate in determining the reactivity of the subcritical systems. The flux-to-current ratio method can provide on-line reactivity measurements, and its spatial effects have been studied and well understood through the numerical simulations of the fuel loading of the neutron source facility of Ukraine. In addition, a modified version of the area-ratio method is proposed for the on-line measurements of the absolute reactivity of the subcritical systems. Monte Carlo numerical simulations showed that this method can provide the absolute reactivity of the subcritical systems within ~5% error, if the neutron detectors are placed close to the fuel assembly interface with the reflector zone or in the reflector zones away from the neutron source.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Exploration of Fuel Property Impacts on the Combustion of Late Post Injections Using Binary Blends and High-Reactivity Ether Bioblendstocks

In this study, the impacts of fuel volatility and reactivity on combustion stability and emissions were studied in a light-duty single-cylinder research engine for a three-injection catalyst heating operation strategy with late post-injections. Here, N-heptane and blends of farnesane/2,2,4,4,6,8,8-heptamethylnonane were used to study the impacts of volatility and reactivity. The effect of increased chemical reactivity was also analysed by comparing the baseline #2 diesel operation with a pure blend of mono-ether components (CN > 100) representative of potential high cetane oxygenated bioblendstocks and a 25 vol.% blend of the mono-ether blend and #2 diesel with a cetane number (CN) of 55. At constant reactivity, little to no variation in combustion performance was observed due to differences in volatility, whereas increased reactivity improved combustion stability and efficiency at late injection timings. Fuels with higher reactivity were found to reduce engine-out hydrocarbon and carbon monoxide emissions while also achieving stable combustion at post-injection timings later than those achievable with #2 diesel fuel. The pure ether blend had the latest achievable post-injection timing of +30.5 CAD while still maintaining stable combustion (coefficient of variation of gross-indicated mean effective pressure < 5%). With post-injection timing adjusted to achieve a matched exhaust temperature of 300 °C, the ether-diesel CN 55 blend was observed to have slightly higher thermal efficiency in comparison to the baseline #2 diesel fuel. The results also indicate that cetane number may serve as a good indicator of combustion characteristics at late injection timings used for aftertreatment thermal management operation.

09 BIOMASS FUELS↗
DOE OSTI2023

Molecular determinants of cross-reactivity and potency by VH3-33 antibodies against the Plasmodium falciparum circumsporozoite protein

IGHV3-33 -encoded antibodies are prevalent in the human humoral response against the Plasmodium falciparum circumsporozoite protein (PfCSP). Among VH3-33 antibodies, cross-reactivity between PfCSP major repeat (NANP), minor (NVDP), and junctional (NPDP) motifs is associated with high affinity and potent parasite inhibition. However, the molecular basis of antibody cross-reactivity and the relationship with efficacy remain unresolved. Here, we perform an extensive structure-function characterization of 12 VH3-33 anti-PfCSP monoclonal antibodies (mAbs) with varying degrees of cross-reactivity induced by immunization of mice expressing a human immunoglobulin gene repertoire. We identify residues in the antibody paratope that mediate cross-reactive binding and delineate four distinct epitope conformations induced by antibody binding, with one consistently associated with high protective efficacy and another that confers comparably potent inhibition of parasite liver invasion. Our data show a link between molecular features of cross-reactive VH3-33 mAb binding to PfCSP and mAb potency, relevant for the development of antibody-based interventions against malaria.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Influence of physical and chemical hydrology on bioremediation of a U-contaminated aquifer informed by reactive transport modeling incorporating 238 U/ 235 U ratios

Microbially-catalyzed reductive immobilization of aqueous uranium (U) as a solid phase has been proposed as a U remediation technique. Both laboratory and field experiments have demonstrated that this reduction reaction alters the 238 U/ 235 U ratio, producing a 238 U-enriched U(IV) solid. In contrast, other major U reactive transport processes fractionate these isotopes much less. This suggests the potential to quantify the extent of bioreduction occurring in groundwater containing U using the 238 U/ 235 U ratio, which would substantially improve upon current practices largely relying on U concentration measurements alone. Current reactive transport models for uranium dynamics include only concentration measurements, which are strongly influenced by highly coupled reactive transformation and aqueous transport processes. The complex physical and chemical behavior of U potentially compromises such quantitative analysis of U storage and release. Here we report the first numerical reactive transport model which explicitly incorporates variations in the 238 U/ 235 U ratio of U and demonstrates improved interpretation of the principal chemical reactions and groundwater transport processes affecting the subsurface mobility and distribution of this widespread contaminant. Recent U bioreduction studies performed in a contaminated aquifer in Rifle, Colorado, USA applied Rayleigh distillation models to interpret U stable isotope fractionation observed as a result of acetate amendment. These simplified models were unable to resolve the spatiotemporal pattern of U isotope fractionation recorded in the aqueous solutes. Here, we employ the multi-component, isotope-enabled CrunchTope reactive transport software to interpret these measured U isotope ratios, and demonstrate accurate reproduction of observed trends in both geochemistry and 238 U/ 235 U ratios for two consecutive years of field experiments. Overall, our results indicate that accurately modeling both the U concentration and isotope ratio distributions greatly constrains the parameter space of the model. We find that the transport properties of U in the Rifle aquifer are governed by the presence of low-permeability regions, which the isotopes are uniquely sensitive to. When U reduction is spatially constrained by these low-permeability regions, the shift in the 238 U/ 235 U ratio becomes more muted. Accurate modeling of observed U isotope ratios thus provides a powerful means to better understand bioremediation, and the current study serves to advance the application of this novel method.

238U/235U↗
DOE OSTI2020

Modeling of fluid injection-induced fault reactivation using coupled fluid flow and mechanical interface model

The present study is aimed at developing a numerical model to reproduce coupled hydro-mechanical processes associated with fault reactivation by fluid injection in low permeability rock, as part of the DECOVALEX-2019 project Task B. We proposed a modeling approach for simulating the processes using the TOUGH-FLAC simulator, and modeled a fault reactivation experiment conducted at Mont Terri Rock Laboratory in Switzerland. The first step of the study involved benchmark calculations considering a simplified fault plane and geometry. Fluid flow along a fault was modeled using elements of aperture-sized thickness on the basis of Darcy's law and the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by zero-thickness interface elements in FLAC3D upon which a slip and/or separation is allowed. A methodology to connect a TOUGH2 volume element to a FLAC3D interface element was developed for handling the hydro-mechanical interactions on the fault during fluid injection. Two different fault models for describing the evolutions of hydraulic aperture by elastic fracture opening and failure-induced aperture increase were considered in the benchmark calculations. In the coupling process, the changes in geometrical features and hydrological properties induced by mechanical deformation were continuously updated. The transient responses of the fault and host rock to stepwise pressurization were examined during the simulation. The hydro-mechanical behavior, including the injection flow rate, pressure distribution around the borehole, stress conditions, and displacements in normal and shear directions were monitored in the surrounding rock and along the fault. The results of benchmark calculations suggest that the developed model reasonably represents the hydro-mechanical behavior of a fault and the surrounding rock. This modeling approach was applied to the fault reactivation experiment of the Mont Terri Rock Laboratory. In this interpretive modeling, a parametric study was conducted to examine the effects of input parameters regarding in situ stress and fault properties on the hydro-mechanical responses of the fault to water injection. Then, an optimal parameter set to reproduce the field experiment results was chosen by trial-and-error. The injection flow rate and pressure response during fault reactivation closely matched those obtained at the site, which indicates the capability of the model to appropriately capture the progressive pathway evolution during fault reactivation tests at the site. Finally, the anchor displacements were overestimated by the model, but a fair agreement was obtained in terms of the order of magnitude and the variation tendency.

58 GEOSCIENCES↗
DOE OSTI2022

Impacts of engine lubrication oil-derived ash on soot oxidative reactivity on a catalytic gasoline particulate filter

Characterizing soot oxidation kinetics is crucial for understanding how gasoline particulate filters (GPFs) perform both in terms of filtration efficiency and pressure drop. The most common method for measuring soot oxidative reactivity is thermogravimetric analysis (TGA). Because TGA is an offline method, it may inaccurately predict how soot oxidizes on a GPF, especially if a catalytic washcoat is present that may enhance oxidation rates through surface interactions as the soot loads on the filter. In this work, a novel in-situ soot oxidative reactivity measurement method was developed. The method involved loading a catalytic GPF under high temperature conditions conducive to soot oxidation and measuring the filtration efficiency of 100 nm particles. A correlation was made between filtration efficiency and loaded soot mass on the filter to allow calculation of the soot oxidation rate. The method was evaluated in experiments using a 2.0 L gasoline direct injection (GDI) engine with three oils of varying additive packages, including an oil with a high zinc dialkyldithiophosphate (ZDDP) concentration, a non-additive, pure poly-alpha olefin (PAO) oil, and an oil with a high concentration of calcium sulfonate. Calculated specific soot oxidation rates ranged from ~0.06–0.6 min –1 and showed strong dependence on GPF inlet gas temperature and lubrication oil additive type. The results clearly demonstrated a catalytic effect of calcium-containing ash particles. Reactivity of soot produced by the engine running with the high calcium containing lubrication oil was increased. Similarly, the results indicated that the catalytic washcoat on the filter increased soot oxidation rates on-filter, especially at low soot loading. Conversely, the ZDDP oil additive exhibited a reactivity inhibiting effect, resulting in lower soot oxidation rates. Finally, this work represents the first known in-situ soot oxidative reactivity measurements on a GPF and elucidates the effect of lubrication oil additives and catalytic washcoat on oxidation rate.

42 ENGINEERING↗
DOE OSTI2025

A quantitative risk assessment framework for fault reactivation in underground hydrogen storage: Coupled simulation and deep learning approach

Underground hydrogen storage (UHS) is emerging as a critical solution for large-scale energy storage. However, like all subsurface fluid injection activities, UHS poses the risk of injection-induced fault reactivation. Accurate risk assessment is essential to ensuring the safety and efficiency of UHS operations. This study presents the development of deep-learning surrogate models for fault reactivation prediction in UHS, trained on a comprehensive database of fully coupled fluid flow-geomechanics simulations. Our findings reveal that analytical models often yield unreliable estimates, with errors up to 54% in the allowable injection pressure, potentially leading to a 40% reduction in UHS operational capacity. The developed surrogate models were incorporated into a quantitative risk assessment (QRA) framework, enabling probabilistic evaluation of fault reactivation risk while accounting for uncertainties in the input variables. Site-specific features, such as horizontal stress gradients, fault’s dip and strike angles, and operational parameters like bottom-hole injection pressure and well-fault distance, were identified as the primary drivers of fault reactivation across various stress regimes. Whereas other hydraulic, geological, and poroelastic reservoir properties were found to have a secondary impact. Notably, we observed that the risk of fault reactivation for a critically oriented fault with a static friction coefficient greater than 0.55 remains below 10% in a normal faulting stress regime. However, the risk significantly increases as the stress regime transitions from normal to strike-slip and ultimately to reverse faulting conditions. These findings underscore the importance of rigorous site characterization and comprehensive QRA evaluations to optimize UHS performance and minimize geomechanical risks.

25 ENERGY STORAGE↗
DOE OSTI2021

Physics-based, neural network force fields for reactive molecular dynamics: Investigation of carbene formation from [EMIM + ][OAc - ]

Reactive molecular dynamics simulations enable detailed understanding of solvent effects on chemical reaction mechanisms and reaction rates. While classical molecular dynamics using reactive force fields allows significantly longer simulation time scales and larger system sizes compared with ab initio molecular dynamics, constructing reactive force fields is a difficult and complex task. In this work, we describe a general approach following the Empirical Valence Bond (EVB) framework for constructing ab initio reactive force fields for condensed phase simulations by combining physics-based methods with neural networks (PB/NN). The physics-based terms ensure correct asymptotic behavior of electrostatic, polarization, and dispersion interactions, and are compatible with existing solvent force fields. Neural networks are utilized for versatile description of short-range orbital interactions within the transition state region, and accurate rendering of vibrational motion of the reacting complex. Herein, we demonstrate our methodology for a simple deprotonation reaction of the 1-ethyl-3-methylimidazolium (EMIM+) cation with acetate to form 1-ethyl- 3-methylimidazol-2-ylidene and acetic acid. Our PB/NN force field exhibits ~ 1 kJ/mol MAE accuracy within the transition state region for the gas-phase complex. To characterize solvent modulation of the reaction profile, we compute potentials of mean force (PMFs) for the gas-phase reaction as well as the reaction within a four ion cluster, and benchmark against ab initio molecular dynamics simulations. We find that the surrounding ionic environment significantly destabilizes formation of the carbene product, and we show that this effect is accurately captured by the reactive force field. By construction, the PB/NN potential may be directly employed for simulations of other solvents/chemical environments without additional parameterization.

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

Structural features of xylan dictate reactivity and functionalization potential for bio-based materials

Plant-based materials have the potential to replace some petroleum-based products, offering compostability and biodegradability as critical advantages. Xylan-rich biomass sources are gaining recognition due to their abundance and underutilization in current industrial applications. Research of potential xylan applications has been complicated by the complex and heterogeneous structure that varies for different xylan feedstocks. Acylation is a broadly used reaction in functionalization of polysaccharides at an industrial scale. However, the efficiency of this reaction varies with the xylan source. To optimize xylan valorization, a systematic understanding of structure–reactivity relationships is essential. This study explores, characterizes, and compares various xylan feedstocks in the acylation process. Xylan feedstocks were analyzed for their chemical composition, degree of polymerization, branching, solubility, and presence of impurities. These features were correlated with xylan glycotypes’ reactivity toward functionalization with succinic anhydride in an optimized DMSO/KOH condition, achieving carboxyl contents of up to 1.46. We used principal component analysis and hierarchical clustering to identify key structural features of xylan that promote its reactivity. Our findings reveal that xylans with higher xylose content and lower degrees of branching exhibit enhanced reactivity, achieving higher carboxyl content and yields. Structural analyses confirmed successful modification, and light scattering analyses showed dramatic changes in the solution properties. Succinylation improves the solubility and film-forming properties of native xylans. This study shows key structure–reactivity relationships in xylan succinylation, establishing that low branching, high xylose content, and reduced lignin impurity enhance chemical functionalization. The results offer a framework for selecting optimal biomass feedstocks and support future efforts in genetic and synthetic biology to design plants with tunable xylan architectures. These findings advance the hemicellulose valorization for applications in coatings and packaging.

Acylation↗
DOE OSTI2022

The Contribution of Proton-Donor pKa on Reactivity Profiles of [FeFe]-hydrogenases

The [FeFe]-hydrogenases are enzymes that catalyze the reversible activation of H 2 coupled to the reduction–oxidation of electron carriers. Members of the different taxonomic groups of [FeFe]-hydrogenases display a wide range of preference, or bias, for H 2 oxidation or H 2 production reactions, despite sharing a common catalytic cofactor, or H-cluster. Identifying the properties that control reactivity remains an active area of investigation, and models have emerged that include diversity in the catalytic site coordination environments and compositions of electron transfer chains. The kinetics of proton-coupled electron transfer at the H-cluster might be expected to be a point of control of reactivity. To test this hypothesis, systematic changes were made to the conserved cysteine residue that functions in proton exchange with the H-cluster in the three model enzymes: CaI, CpII, and CrHydA1. CaI and CpII both employ electron transfer accessory clusters but differ in bias, whereas CrHydA1 lacks accessory clusters having only the H-cluster. Changing from cysteine to either serine (more basic) or aspartate (more acidic) modifies the sidechain pKa and thus the barrier for the proton exchange step. The reaction rates for H 2 oxidation or H 2 evolution were surveyed and measured for model [FeFe]-hydrogenases, and the results show that the initial proton-transfer step in [FeFe]-hydrogenase is tightly coupled to the control of reactivity; a change from cysteine to more basic serine favored H 2 oxidation in all enzymes, whereas a change to more acidic aspartate caused a shift in preference toward H 2 evolution. Overall, the changes in reactivity profiles were profound, spanning 10 5 in ratio of the H 2 oxidation-to-H 2 evolution rates. The fact that the change in reactivity follows a common trend implies that the effect of changing the proton-transfer residue pKa may also be framed as an effect on the scaling relationship between the H-cluster di(thiolmethyl)amine (DTMA) ligand pKa and E m values of the H-cluster. Experimental observations that support this relationship, and how it relates to catalytic function in [FeFe]-hydrogenases, are discussed.

59 BASIC BIOLOGICAL SCIENCES↗