Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reactivity”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

An Evaluation of Actinide Reactivity with CO 2 , O 2 , and O 2 /He Gases using Inductively Coupled Plasma Tandem Mass Spectrometry: Application to Simultaneous Measurement of 241 Am/ 241 Pu Ratios in Unseparated Complex Matrices

Accurate actinide measurements are critical within the field of nuclear science. Traditional methods for actinide quantification require time-consuming sample processing prior to analysis. There is a need for rapid analytical techniques that still maintain a high degree of accuracy. In this work, actinide reactivity was assessed for multiple oxygen-containing reaction gases using quadrupole inductively coupled plasma tandem mass spectrometry (Q-ICP-MS/MS) to evaluate actinide analysis in complex sample matrices without analyte-matrix separation. A novel method was developed to measure 241 Am/ 241 Pu in complex sample matrices using O 2 /He reaction gas with no matrix removal or analyte pre-concentration. This inline method reduces matrix-derived polyatomic interferences that complicate traditional ICP-MS analyses by mass-shifting to 241 Am 16 O + and 241 Pu 16 O 2 + , allowing Am and Pu to be mass separated for simultaneous analysis. While mass shifting is efficient, a small portion of Am + (<1.3%) and Pu + (<1.4%) react to from AmO 2 + and PuO + , respectively. Therefore, a mass balance approach was used, in combination with reactivity determined from 242 Pu and 243 Am standard solutions, to correct for residual 241 PuO + and 241 AmO 2 + . The method was validated by measuring 241 Am/ 241 Pu in Pu isotope standards CRM-136 and CRM-137 (separated in March/April 1970 and February 2022, respectively) in both neat solutions and complex matrices containing diluted soil (NIST SRM 2711a, >1000 µg·g -1 ). Method detection limits of 15.9 and 9.6 fg·g -1 were determined for 241 Am and 241 Pu, respectively, and 241 Am/ 241 Pu ratios were measured with accuracies within <3.5%. In conclusion, this work presents the first direct analysis of 241 Am/ 241 Pu in unseparated complex matrices, advancing capabilities for rapid actinide measurements.

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

Vibrational Quantum-State-Controlled Reactivity in the O 2 + + C 3 H 4 Reaction

Quantum-state-controlled reactivity is a long-standing goal in the field of physical chemistry. In this work, we explore the vibrational-state-dependent behavior of the ion–molecule reaction between O 2 + in distinct vibrational states and two isomers of C 3 H 4 , allene (H 2 C 3 H 2 ) and propyne (H 3 C 3 H). While most products are formed regardless of the vibrational state of O 2 + , the branching ratios are influenced by vibrational excitation, and a new product, C 2 O + , appears exclusively in the excited-state reactions. This selective formation of C 2 O + demonstrates that vibrational excitation can effectively activate a reaction pathway, providing direct evidence of quantum-state control in the reactivity. These results represent an important step toward the goal of quantum-state-controlled chemistry in molecular systems.

Chemical reactions↗
DOE OSTI2020

Controlling P–C/C–H Bond Cleavage in Nickel Bis(diphosphine) Complexes: Reactivity Scope, Mechanism, and Computations

The synthesis of heteroleptic [Ni(P 2 N 2 )(diphosphine)][BF 4 ] 2 complexes and the cleavage of P–C and C–H bonds of the P 2 N 2 ligand in those complexes are reported here. The products are five-coordinate complexes in which Ni–C and P–H bonds have formed to give a cyclic moiety containing Ni–CH$=$NR 2 . The reactivity of [Ni(P 2 N 2 )(diphosphine)][BF 4 ] 2 complexes is influenced by the rigidity of the diphosphine, the steric effect of the substituents, and length of the carbon linker of the diphosphine ligands. Diphosphine ligands bearing a rigid backbone (e.g., dmpbz, 1,2-bis(dimethylphosphino)benzene) or aromatic substituents (e.g., dppe, 1,2-bis(diphenylphosphino)ethane) react with [Ni(P t Bu 2 N Bn 2 )(CH 3 CN) 2 ][BF 4 ] 2 to give P–C/C–H bond cleavage products. Both [Ni(P t Bu 2 N Bn 2 )(dmpe)(MeCN)][BF 4 ] 2 and [Ni(P t Bu 2 N Bn 2 )(dmpm)(MeCN)][BF 4 ] 2 (dmpm = 1,2-bis(dimethylylphosphino)methane) were prepared by the reaction of [Ni(P t Bu 2 N Bn 2 )(CH 3 CN) 2 ][BF 4 ] 2 with the corresponding diphosphine ligands. [Ni(P t Bu 2 N Bn 2 )(dmpe)(MeCN)][BF 4 ] 2 readily undergoes P–C/C–H bond cleavage in nitromethane. In sharp contrast, [Ni(P t Bu 2 N Bn 2 )(dmpm)][BF 4 ] 2 is stabilized by dmpm, a diphosphine with small bite angle, and does not show P–C/C–H bond cleavage reactivity. Computational results show that for complexes bearing less bulky diphosphine ligands, such as dmpm, the barriers for the rate-determining transition states are in some examples higher than 30 kcal/mol with the M06 functional, higher than those for complexes bearing more rigid or more bulky ligands, consistent with experimental studies. The calculated barriers for the first transition state correlated with increased values of the dihedral angle formed by the two NiP 2 planes.

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

Ultrafast Reactive Laser Sintering of Highly Conductive Garnet-Type LLZTO Solid Electrolytes

Rapid and scalable fabrication of garnet-type solid electrolytes remains a major challenge for the practical deployment of lithium metal batteries. Here, we report reactive laser sintering (RLS) as an ultrafast and potentially scalable strategy for fabricating garnet-type Li 6.4 La 3 Zr 1.4 Ta 0.6 O 12 (LLZTO) solid electrolytes. RLS of LLZTO enables simultaneous reaction and densification, achieving ∼95% relative density while minimizing lithium loss and suppressing secondary phase formation. Compared to conventional furnace sintering, RLS promotes enhanced grain growth and improved densification, leading to improved ionic conductivity (0.36 ± 0.08 mS cm −1 ) while maintaining comparable activation energies for Li + transport. Structural characterization by X-ray diffraction (XRD), Raman spectroscopy, and solid-state 6 Li/ 7 Li NMR confirms the formation of cubic garnet LLZTO with homogeneous microscale elemental distribution. In addition, nanoindentation measurements demonstrate that RLS preserves the mechanical properties of the garnet framework despite ultrafast localized thermal processing. By integrating simultaneous reaction and densification with tunable microstructural control, reactive laser sintering provides a promising manufacturing pathway for high-performance garnet solid electrolytes toward next-generation solid-state batteries.

CO2 laser↗
DOE OSTI2025

Chemical Recycling of Polycaprolactones via Reactive Melt Processing

Chemical recycling is a promising technology for the deconstruction of waste plastics into monomers or other chemical intermediates, which can be converted into other value-added products or repolymerized into polymers. In this work, polycaprolactone and poly(4-propylcaprolactone) were chemically recycled by exploiting ring-closing depolymerization (RCD). To demonstrate the utility of this approach, reactive distillation experiments were initially conducted in a heated round-bottom flask with a distillation head to explore the effect of temperature and catalyst loading on RCD. Inspired by industrial devolatilization equipment, analogous experiments were conducted using a twin-screw melt compounder that was modified to simultaneously heat and mix polymer in the presence of a catalyst while continuously removing the evolved monomer vapors with a vacuum system and collecting them in a cold trap. Yield and selectivity of the recovered monomer were characterized by 1 H nuclear magnetic resonance spectroscopy and gas chromatography-mass spectrometry. Importantly, the more scalable twin-screw melt compounder approach produced values of monomer selectivity (∼90%) and yield (∼80%) similar to those of the round-bottom flask experiments. The recovered monomers were also successfully repolymerized into the same neat polymer with comparable molar masses without purification steps or adding initiator. Here, the results of this study promote the use of reactive devolatilization extrusion as a scalable approach to polyester depolymerization.

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

CH 4 Activation over Perovskite Catalysts: True Density and Reactivity of Active Sites

The high thermal stability of perovskites has drawn attention toward their applications for catalytic CH 4 activation and conversion, typically occurring at high temperatures. The reaction rates of perovskite catalysts for CH 4 combustion, however, trail behind those of noble metal catalysts. Ways to optimize the performance of perovskite catalysts are destined to trial-and-error approaches unless their complex reconstructed surfaces are correlated with fundamental kinetic parameters. Discerning the intrinsic activity of surface catalytic sites and the density of those sites is crucial to rationally envision complex metal oxides with enhanced catalytic performance. Here, the present work presents a detailed kinetic analysis of catalytic CH 4 combustion over a set of seven perovskites (SrTiO 3 , SrZrO 3 , SrFeO 3 , LaFeO 3 , LaInO 3 , LaCoO 3 , LaMnO 3 ) with various surface terminations. Steady-state isotopic transient kinetic analysis was employed to measure turnover frequency (TOF) and density of surface intermediates (N) under operando conditions. Top surface characterization elucidated performance-structure relationships between near-monolayer surface composition and intrinsic reactivity of the catalysts. By using a chemical etching procedure to expose Fe-sites at the top surface of LaFeO 3 (LaFeO 3 , HNO 3 ), its TOF was increased 4-fold, compared with the unmodified sample, although N on the surface of LaFeO 3 , HNO 3 decreased. Density functional theory simulations corroborated that surface Fe-termination and La-Fe termination offer lower energetic barriers for CH 4 activation when compared with La-termination. In general, surface reconstruction is shown as a tool to tune TOF and N to improve reaction rates. This work fills a gap in current kinetic studies of perovskites through a careful assessment and discussion of the density and intrinsic reactivity of active sites for methane combustion over well-characterized reconstructed perovskite surfaces.

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

A Thioether-Ligated Cupric Superoxide Model with Hydrogen Atom Abstraction Reactivity

The central role of cupric superoxide intermediates proposed in hormone and neurotransmitter biosynthesis by noncoupled binuclear copper monooxygenases like dopamine-β-monooxygenase has drawn significant attention to the unusual methionine ligation of the Cu M (“Cu B ”) active site characteristic of this class of enzymes. The copper–sulfur interaction has proven critical for turnover, raising still-unresolved questions concerning Nature’s selection of an oxidizable Met residue to facilitate C–H oxygenation. We describe herein a model for Cu M , [( TMG N 3 S)Cu I ] + ([1] + ), and its O 2 -bound analog [( TMG N 3 S)Cu II (O 2 •– )] + ([1·O 2 ] + ). The latter is the first reported cupric superoxide with an experimentally proven Cu–S bond which also possesses demonstrated hydrogen atom abstraction (HAA) reactivity. Introduction of O2 to a precooled solution of the cuprous precursor [1]B(C 6 F 5 ) 4 (-135 °C, 2-methyltetrahydrofuran (2-MeTHF)) reversibly forms [1·O 2 ]B(C 6 F 5 ) 4 (UV/vis spectroscopy: λ max 442, 642, 742 nm). Resonance Raman studies (413 nm) using 16 O 2 [ 18 O 2 ] corroborated the identity of [1·O 2 ] + by revealing Cu–O (446 [425] cm –1 ) and O–O (1105 [1042] cm –1 ) stretches, and extended X-ray absorption fine structure (EXAFS) spectroscopy showed a Cu–S interatomic distance of 2.55 Å. HAA reactivity between [1·O 2 ] + and TEMPO–H proceeds rapidly (1.28 × 10 –1 M –1 s –1 , -135 °C, 2-MeTHF) with a primary kinetic isotope effect of kH/kD = 5.4. Comparisons of the O 2 -binding behavior and redox activity of [1] + vs [2] + , the latter a close analog of [1]+ but with all N atom ligation (i.e., N 3 S vs N 4 ), are presented.

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

Enhanced CO 2 Reactive Capture and Conversion Using Aminothiolate Ligand–Metal Interface

Metallic catalyst modification by organic ligands is an emerging catalyst design in enhancing the activity and selectivity of electrocatalytic carbon dioxide (CO 2 ) reactive capture and reduction to value-added fuels. However, a lack of fundamental science on how these ligand–metal interfaces interact with CO 2 and key intermediates under working conditions has resulted in a trial-and-error approach for experimental designs. With the aid of density functional theory calculations, we provided a comprehensive mechanism study of CO 2 reduction to multicarbon products over aminothiolate-coated copper (Cu) catalysts. Our results indicate that the CO 2 reduction performance was closely related to the alkyl chain length, ligand coverage, ligand configuration, and Cu facet. The aminothiolate ligand–Cu interface significantly promoted initial CO 2 activation and lowered the activation barrier of carbon–carbon coupling through the organic (nitrogen (N)) and inorganic (Cu) interfacial active sites. Experimentally, the selectivity and partial current density of the multicarbon products over aminothiolate-coated Cu increased by 1.5-fold and 2-fold, respectively, as compared to the pristine Cu at –1.16 V RHE , consistent with our theoretical findings. Furthermore, this work highlights the promising strategy of designing the ligand–metal interface for CO 2 reactive capture and conversion to multicarbon products.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Unraveling Reactivity Differences: Room-Temperature Ring-Opening Metathesis Polymerization (ROMP) versus Frontal ROMP

Here, in this study, we explore the distinct reactivity patterns between frontal ring-opening metathesis polymerization (FROMP) and room-temperature solventless ring-opening metathesis polymerization (ROMP). Despite their shared mechanism, we find that FROMP is less sensitive to inhibitor concentration than room-temperature ROMP. By increasing the initiator-to-monomer ratio for a fixed inhibitor/initiator quantity, we find reduction in the ROMP background reactivity at room temperature (i.e., increased resin pot life). At elevated temperatures where inhibitor dissociation prevails, accelerated frontal polymerization rates are observed because of the concentrated presence of the initiator. Surprisingly, the strategy of employing higher initiator loading enhances both pot life and front speeds, which leads to FROMP rates exceeding prior reported values by over 5 times. This counterintuitive behavior is attributed to an increase in the proximity of the inhibitor to the initiator within the bulk resin and to whether the temperature favors coordination or dissociation of the inhibitor. A rapid method was developed for assessing resin pot life, and a straightforward model for active initiator behavior was established. Modified resin systems enabled direct ink writing of robust thermoset structures at rates much faster than previously possible.

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

Reactive Transport Modeling of Mineral Precipitation and Carbon Trapping in Discrete Fracture Networks

In this study we use numerical experiments to analyze reactive flow and transport behavior in discrete intersecting fracture networks, focusing on (a) how reaction-induced changes in physical and chemical properties affect flow connectivity and (b) how fracture networks developed in the Earth's critical zone contribute to carbon sequestration via mineral weathering reactions. In the first part of the study, we used two-dimensional reactive flow and transport simulations to analyze the impacts of mixing in a natural discrete fracture network. We concluded that reaction-induced changes can substantially alter the flow connectivity, especially at fracture intersections. The second set of simulations considered the problem of natural weathering of fractured mafic and ultramafic rocks in the partially saturated Earth's critical zone as a function of infiltration rates, fracture permeability, and partially saturated flow parameters. As a model system, we considered an incongruent reaction network with dissolution of forsterite and precipitation of magnesite. The behavior is complex in terms of the rate-controlling processes because of the multicomponent nature of the system as shown by the grid Peclet number: the CO 2 behavior is gas diffusion-controlled in the partially saturated zone, while the rate of water flow via the Damkӧhler number controls Mg 2+ transport through the fracture network. The amounts of carbon that can be trapped are modest, but the naturally fractured domain considered here provides a useful “base case” against which various engineered solutions can be compared.

58 GEOSCIENCES↗
DOE OSTI2026

Identifying mechanistic differences between co-fed CO 2 hydrogenation and reactive CO 2 capture using Ru and Pd dual function materials

Dual function materials (DFMs) enable reactive carbon capture (RCkeC), an intensified approach to carbon dioxide capture and utilization for cost and energy input reductions. Yet, there is a fundamental lack of understanding of mechanisms around CO 2 adsorption and subsequent conversion on these materials, hindering further development. Herein, we investigated several supported alkaline metal oxides for their CO 2 adsorption characteristics to find that Na/Al 2 O 3 had the highest CO 2 adsorption capacity, accompanied by a variety of CO 2 adsorption geometries as identified by in situ DRIFTS and computational modeling. The addition of catalytic metals (Ru, Pd) increased the adsorption capacity of Na/Al 2 O 3 without altering binding modes. In the subsequent reactive desorption step, acetate and formate intermediates were observed. Notably, this mechanistic investigation identified that the formation of acetate species was unique to RCC on a DFM, as these species were not observed in co-fed hydrogenation over the DFM or RCC over a Na-free catalyst.

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

Active learning for robust, high-complexity reactive atomistic simulations

Machine learned reactive force fields based on polynomial expansions have been shown to be highly effective for describing simulations involving reactive materials. Nevertheless, the highly flexible nature of these models can give rise to a large number of candidate parameters for complicated systems. In these cases, reliable parameterization requires a well-formed training set, which can be difficult to achieve through standard iterative fitting methods. In this paper, we present an active learning approach based on cluster analysis and inspired by Shannon information theory to enable semi-automated generation of informative training sets and robust machine learned force fields. The use of this tool is demonstrated for development of a model based on linear combinations of Chebyshev polynomials explicitly describing up to four-body interactions, for a chemically and structurally diverse system of C/O under extreme conditions. We show that this flexible training database management approach enables development of models exhibiting excellent agreement with Kohn–Sham density functional theory in terms of structure, dynamics, and speciation.

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

Quantum reactive scattering calculations for the cold and ultracold Li + LiNa → Li 2 + Na reaction

A first-principles based quantum dynamics study of the Li + LiNa(ν = 0, j = 0) → Li 2 (ν', j') + Na reaction is reported for collision energies spanning the ultracold (1 nK) to cold (1 K) regimes. A full-dimensional ab initio potential energy surface for the ground electronic state of Li 2 Na is utilized that includes an accurate treatment of the long-range interactions. The Li + LiNa reaction is barrierless and exoergic and exhibits a deep attractive potential well that supports complex formation. Thus, significant reactivity occurs even for collision temperatures approaching absolute zero. The reactive scattering calculations are based on a numerically exact time-independent quantum dynamics methodology in hyperspherical coordinates. Total and rotationally resolved rate coefficients are reported at 56 collision energies and include all contributing partial waves. Several shape resonances are observed in many of the rotationally resolved rate coefficients and a small resonance feature is also reported in the total rate coefficient near 50 mK. Furthermore, of particular interest, the angular distributions or differential cross sections are reported as a function of both the collision energy and scattering angle. Unique quantum fingerprints (bumps, channels, and ripples) are observed in the angular distributions for each product rotational state due to quantum interference and shape resonance contributions. The Li + LiNa reaction is under active experimental investigation so that these intriguing features could be verified experimentally when sufficient product state resolution becomes feasible for collision energies below 1 K.

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

A physics-informed deep learning description of Knudsen layer reactivity reduction

A physics-informed neural network (PINN) is used to evaluate the fast ion distribution in the hot spot of an inertial confinement fusion target. The use of tailored input and output layers to the neural network is shown to enable a PINN to learn the parametric solution to the Vlasov–Fokker–Planck equation in the absence of any synthetic or experimental data. As an explicit demonstration of the approach, the specific problem of Knudsen layer fusion yield reduction is treated. Here, the predictions from the Vlasov–Fokker–Planck PINN are used to provide a non-perturbative solution of the fast ion tail in the vicinity of the hot spot, thus allowing the spatial profile of the fusion reactivity to be evaluated for a range of collisionalities and hot spot conditions. Excellent agreement is found between the predictions of the Vlasov–Fokker–Planck PINN and the results from traditional numerical solvers with respect to both the energy and spatial distribution of fast ions and the fusion reactivity profile, demonstrating that the Vlasov–Fokker–Planck PINN provides an accurate and efficient means of determining the impact of Knudsen layer yield reduction across a broad range of plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Electric fields imbue enzyme reactivity by aligning active site fragment orbitals

It is broadly recognized that intramolecular electric fields, produced by the protein scaffold and acting on the active site, facilitate enzymatic catalysis. This field effect can be described by several theoretical models, each of which is intuitive to varying degrees. In this contribution, we show that a fundamental effect of electric fields is to generate electrostatic potentials that facilitate the energetic alignment of reactant frontier orbitals. We apply this model to demystify the impact of electric fields on high-valent iron–oxo heme proteins: catalases, peroxidases, and peroxygenases/monooxygenases. Specifically, we show that this model easily accounts for the observed field-induced changes to the spin distribution within peroxidase active sites and explains the transition between epoxidation and hydroxylation pathways seen in Cytochrome P450 active site models. Thus, for the intuitive interpretation of the chemical effect of the field, the strategy involves analyzing the response of the orbitals of active site fragments, and their energetic alignment. We note that the energy difference between fragment orbitals involved in charge redistribution acts as a measure for the chemical hardness/softness of the reactive complex. This measure, and its sensitivity to electric fields, offers a single parameter model from which to quantitatively assess the effects of electric fields on reactivity and selectivity. Thus, the model provides an additional perspective to describe electrostatic preorganization and offers ways for its manipulation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Global Sensitivity Analysis of a Reactive Transport Model for Mineral Scale Formation During Hydraulic Fracturing

Injection of water-based hydraulic fracturing fluid (HFF) into tight shale gas/oil formations can increase formation permeability and enhance production rates, but this process frequently causes mineral scale formation that can occlude pore space and hinder flow. To identify the most important factors that control the formation of mineral scales, we applied a novel global sensitivity analysis method—distance-based generalized sensitivity analysis (DGSA)—to a reactive transport model (RTM) that was previously built and calibrated to simulate precipitation of barite [BaSO4] and iron (hydr)oxide [Fe(OH) 3 ] in shale matrices and on fracture surfaces. Reactive transport simulations were run with model parameters randomly sampled based on assigned uncertainties. Modeling results for barite and Fe(OH)3 formation were clustered using machine-learning algorithms. A list of ranked critical input parameters was obtained after statistical quantification of cumulative distribution functions of input parameters. We found that barite formation is most sensitive to the rate of sulfate ion generation, which is determined by the pyrite dissolution rate coefficient and oxidant availability. In addition, barite formation is sensitive to the initial amounts of barite in HFF and shale, followed by barite thermodynamics/kinetics. For Fe(OH) 3 formation, the ranked factors are Fe(OH)3 precipitation rate coefficients, initial HFF pH, initial Fe(OH) 3 amount in HFF, and oxidant availability. Overall, our results provide insights into managing mineral scale formation during hydraulic fracturing to enhance production. Meanwhile, this study serves as an example of global sensitivity analysis of RTMs using the efficient, straightforward, and open-source DGSA method.

58 GEOSCIENCES↗
DOE OSTI2021

Reactive Flows in Porous Media: Challenges in Theoretical and Numerical Methods

We review theoretical and computational research, primarily from the past 10 years, addressing the flow of reactive fluids in porous media. The focus is on systems where chemical reactions at the solid–fluid interface cause dissolution of the surrounding porous matrix, creating nonlinear feedback mechanisms that can often lead to greatly enhanced permeability. Here, we discuss insights into the evolution of geological forms that can be inferred from these feedback mechanisms, as well as some geotechnical applications such as enhanced oil recovery, hydraulic fracturing, and carbon sequestration. Until recently, most practical applications of reactive transport have been based on Darcy-scale modeling, where averaged equations for the flow and reactant transport are solved. We summarize the successes and limitations of volume averaging, which leads to Darcy-scale equations, as an introduction to pore-scale modeling. Pore-scale modeling is computationally intensive but offers new insights as well as tests of averaging theories and pore-network models. We include recent research devoted to validation of pore-scale simulations, particularly the use of visual observations from microfluidic experiments.

wormhole formation↗
DOE OSTI2025

2D reactive transport model of shale chemical weathering and biogeochemical fluxes along a mountainous hillslope, East River Watershed, Colorado: Input files and simulation results

This data package contains input files and simulation results for a two-dimensional (2D) reactive transport model used to quantitatively analyze the coupled hydrological and biogeochemical processes governing shale weathering and associated biogeochemical fluxes under realistic environmental conditions in the high-elevation East River Watershed. These data support the conclusions presented in Stolze et al. (Water Resources Research, under review), "Model-based interpretation of solute exports and carbon partitioning during shale weathering in a mountainous hillslope". The model simulates atmospheric-subsurface gas exchange, subsurface water flow, and shale weathering processes under dynamic, year-scale conditions along a shale-underlain hillslope located in the East River watershed. The simulations were performed using the PFLOTRAN flow and reactive transport code and executed on the Perlmutter supercomputer to leverage its large-scale parallel computing capabilities. The data package contains two zipped folders, "model_input_files" and "simulation_results", and one readme.txt file. "model_input_files" contains the necessary input files to run the calibrated base-base model presented in Stolze et al. (Water Resources Research, under review). "simulation_results" contains a single hdf5 file ("Output_2D_hillslope_model.h5") which includes the results of simulation performed using the base-case model. This file can be opened with HDFView 3.1.4, Python, or MATLAB. "readme.txt" contains relevant information about the base-case model and provides guidelines on how to run the associated input files provided in the folder "model_input_files". Furthermore, readme.txt provides information regarding the model results provided in "Output_2D_hillslope_model.h5" such as matrix dimensionality and output units. Field datasets used to evaluate model performance were collected at three monitoring wells located along a hillslope transect (PLM1, PLM2, and PLM3). Dissolved ion concentration data were collected from November 2016 to October 2021 for Ca, Mg, DIC, Na, K, SO4 (Dong et al., 2025 - dic_npoc_data_2014_2024.zip - DOI:10.15485/1660459; Williams et al., 2025 - anion_data_2014_2024.zip - DOI:10.15485/1668054; Dong et al., 2025 - cation_data_2014_2024.zip - DOI:10.15485/1668055). Note that we used the files named er_PLM1_xx_yy, er_PLM2_xx_yy, and er_PLM3_xx_yy where xx stands for the name of the aqueous species and yy stands for the depth where the measurements were performed. Soil water content ([0 - 1] m) and water table depth were collected from November 2016 to October 2021 (Wan et al., 2024 - Dynamic_water_table__depthsFig2b.csv and Soil_water_content_Fig4e.csv - DOI:10.15485/2322567). Gaseous CO2 concentration were collected from October 2020 to December 2021(Wan et al., 2024 - Soil_CO2_concentrations_Fig4h.csv - DOI:10.15485/2322567) Gaseous CO2 flux from the subsurface to the atmosphere were collected in the vicinity of PLM2 from October 2019 to May 2022 (Wu et al., 2025). Soil microbial biomass concentration was measured from August 2016 to June 2017 (Sorensen et al., 2019 - 2017_East_River_Pumphouse_Microbial_Biomass__1_.csv - DOI:10.15485/1577267) All field data are published as CSV files compatible with Microsoft Excel, MATLAB, and Python, or as text files. The coordinates of the monitoring wells and the CO2(g) flux sensor in the coordinate system WGS84 are: -PLM1: [38.9197710 ; -106.9492750] -PLM2: [38.9201580 ; -106.9487170] -PLM3: [38.9207843 ; -106.9483668] -PLM4: 38.9210060 ; -106.9479528] -CO2(g) flux sensor: [38.9199180 ; -106.9489906] ------------------------------------------------------------------------------------------- This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231. This research used resources of the National Energy Research Scientific Computing Center (NERSC), a Department of Energy User Facility using NERSC award BER-ERCAP 23980, BER-ERCAP 28550, and BER-ERCAP 33789.

54 ENVIRONMENTAL SCIENCES↗