Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reactive transport model”

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

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

At least 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Implementation of Perturbation Theory and Sensitivity Capabilities in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor Multiphysics analysis application, jointly developed by Argonne and Idaho National Laboratories under the DOE-NE NEAMS program. This fiscal year, capabilities for reactivity and sensitivity evaluation using perturbation methods were implemented and verified. The First Order Perturbation Method (FOPT) was employed to compute reactivity worth resulting from small perturbations in input parameters, while the Generalized Perturbation Theory (GPT) was used to evaluate sensitivities of a range of response types, including reaction rate ratio, k-eigenvalue, neutron generation time, and effective delayed neutron fraction. These perturbation methods enable users to quantify how response quantities change due to a perturbation in a input parameter without explicitly performing an additional transport simulation for each perturbed state. In particular, the GPT formulation accounts for indirect effects arising from flux changes by solving generalized inhomogeneous equations, for which a Neumann series-based iterative solution method was developed and implemented in Griffin. The implemented reactivity and sensitivity evaluation capabilities were verified using two test problems: an infinite homogeneous system and a two-dimensional hexagonal core. The results showed excellent agreement with reference solutions obtained by a direct method based on finite difference approximation as well as GPT-based results from the PERSENT code, confirming the accuracy of both reactivity and sensitivity evaluations. Additionally, preliminary uncertainty quantification (UQ) results were obtained by combining the sensitivity values computed using GPT and external covariance data, demonstrating that the implemented sensitivity results can be reliably used for uncertainty calculations. To further demonstrate the generality and practical strength of the implementation, the sensitivity evaluation capability was successfully applied to the Empire microreactor with a geometrically complex design that poses significant modeling challenges. The results confirm that Griffin enables sensitivity evaluations even for irregular and highly heterogeneous reactor configurations, thereby establishing a foundation for UQ applications in advanced reactor designs and analyses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2007

Summertime Influence of Asian Pollution in the Free Troposphere over North America

We analyze aircraft observations obtained during INTEX-A (1 July 14 - August 2004) to examine the summertime influence of Asian pollution in the free troposphere over North America. By applying correlation analysis and Principal Component Analysis (PCA) to the observations between 6-12 km, we find dominant influences from recent convection and lightning (13 percent of observations), Asia (7 percent), the lower stratosphere (7 percent), and boreal forest fires (2 percent), with the remaining 71 percent assigned to background. Asian airmasses are marked by high levels of CO, O3, HCN, PAN, acetylene, benzene, methanol, and SO4(2-). The partitioning of reactive nitrogen species in the Asian plumes is dominated by peroxyacetyl nitrate (PAN) (approximately 600 pptv), with varying NO(x)/HNO3 ratios in individual plumes consistent with different plumes ages ranging from 3 to 9 days. Export of Asian pollution in warm conveyor belts of mid-latitude cyclones, deep convection, and lifting in typhoons all contributed to the five major Asian pollution plumes. Compared to past measurement campaigns of Asian outflow during spring, INTEX-A observations display unique characteristics: lower levels of anthropogenic pollutants (CO, propane, ethane, benzene) due to their shorter summer lifetimes; higher levels of biogenic tracers (methanol and acetone) because of a more active biosphere; as well as higher levels of PAN, NO(x), HNO3, and O3 (more active photochemistry possibly enhanced by injection of lightning NO(x)). The high delta O3/delta CO ratio (0.76 mol mol(exp -1)) of Asian plumes during INTEX-A is due to a combination of strong photochemical production and mixing with stratospheric air along isentropic surfaces. The GEOS-Chem global chemical transport model captures the timing and location of the Asian plumes remarkably well. However, it significantly underestimates the magnitude of the enhancements.

Liang, Q.↗
DOE OSTI2025

Constraining the Low-Temperature Oxidation Mechanism of n -Hexanol through the Detection and Identification of C 6 Elusive Intermediates

Alcohol-based fuels are currently considered to be viable energy carriers for the transportation sector. Consequently, a comprehensive mechanistic understanding of the low-temperature oxidation of alcohols is essential for application in advanced low-temperature compression engines. Here, in this work, a multidimensional approach involving experimental investigations, kinetic modeling, and theoretical calculations was used to provide new insights into the low-temperature oxidation mechanism of a C 6 alcohol, n -hexanol (CH 3 (CH 2 ) 5 OH), through the detection and identification of elusive C 6 intermediates. The oxidation of n-hexanol was investigated in a jet-stirred reactor under stoichiometric conditions (ϕ = 1.0), an initial fuel concentration of 2%, a residence time of 2 s, a temperature range between 500 and 660 K, and a pressure of 700 Torr. The reactants, intermediates, and final products were detected and identified by means of molecular-beam mass spectrometry coupled with single-photon ionization employing tunable synchrotron-generated vacuum ultraviolet radiation. Chemical kinetic simulations were performed using a previously published kinetic model (Togbé et al., Energy Fuels 2010, 11, 5859−5875) to predict the reactivity of n-hexanol and elucidate the predominant formation pathways of the observed low-temperature species. Experimental photoionization efficiency curves in conjunction with ab initio calculations, enabled the identification of important low-temperature species, such as C 6 unsaturated alcohols, C 6 olefinic hydroperoxides, C 6 cyclic ethers, C 6 diones, and C 6 ketohydroperoxides. The results of this study provide valuable insight into the mechanism of the low-temperature oxidation chemistry of n -hexanol, contributing to the development of kinetic models for the low-temperature oxidation of n-hexanol and other long-chain linear alcohols.

alcohols↗
DOE OSTI2022

Advanced Distributed Wind Turbine Controls Series: Part 2-Wind Energy in Isolated Grids – Microgrids, Infrastructure Resilience, and Advanced Controls Launchpad (MIRACL)

In an isolated grid, wind turbines are typically deployed to provide energy to maximize energy production, reduce diesel-fuel consumption, reduce carbon emissions, and reduce costs for energy and fuel transportation. However, in addition to solely providing energy to the power system, wind turbines contain rotating masses and inverter-based controls that can enable various reliability and resilience services through advance controls. As part of the Microgrids, Infrastructure Resilience, and Advanced Controls Launchpad (MIRACL) this paper demonstrates, through desktop simulations, advanced wind turbine controls that can be employed to support higher contributions of wind in isolated grids, and to demonstrate ways that wind can play a role in supporting stability of an isolated grid. This isolated grid used in these desktop simulations is comprised of a wind turbine (600 kW), solar PV (430 kW), battery energy storage system (1 MW/1MWh), a simulated diesel generator (2 MW) and various types of loads (critical, dynamic). We developed a model of the subsystems in MATLAB/Simulink and validated them with available data from real-world components on NREL's Flatirons Campus. These validated models are then configured for various case studies. We compare the output of the desktop simulation with a baseline case with the diesel generator. Active and reactive power control of the wind turbine can help improve frequency and voltage responses in the isolated grid, respectively. By utilizing a small integrated battery energy storage system in the DC-link of the wind turbine, we also demonstrate that wind turbines can help blackstart a critical load comparable to its rated power and support other renewables (e.g. solar PV) come online and pick up an additional load. This report illustrates some of these reliability and resilience services a wind turbine can provide in an isolated grid.

17 WIND ENERGY↗
DOE OSTI2020

Direct kinetic measurements and theoretical predictions of an isoprene-derived Criegee intermediate

Isoprene has the highest emission into Earth’s atmosphere of any nonmethane hydrocarbon. Atmospheric processing of alkenes, including isoprene, via ozonolysis leads to the formation of zwitterionic reactive intermediates, known as Criegee intermediates (CIs). Direct studies have revealed that reactions involving simple CIs can significantly impact the tropospheric oxidizing capacity, enhance particulate formation, and degrade local air quality. Methyl vinyl ketone oxide (MVK-oxide) is a four-carbon, asymmetric, resonance-stabilized CI, produced with 21 to 23% yield from isoprene ozonolysis, yet its reactivity has not been directly studied. We present direct kinetic measurements of MVK-oxide reactions with key atmospheric species using absorption spectroscopy. Direct UV-Vis absorption spectra from two independent flow cell experiments overlap with the molecular beam UV-Vis-depletion spectra reported recently [M. F. Vansco, B. Marchetti, M. I. Lester, J. Chem. Phys. 149, 44309 (2018)] but suggest different conformer distributions under jet-cooled and thermal conditions. Comparison of the experimental lifetime herein with theory indicates only the syn -conformers are observed; anti -conformers are calculated to be removed much more rapidly via unimolecular decay. We observe experimentally and predict theoretically fast reaction of syn -MVK-oxide with SO 2 and formic acid, similar to smaller alkyl-substituted CIs, and by contrast, slow removal in the presence of water. We determine products through complementary multiplexed photoionization mass spectrometry, observing SO 3 and identifying organic hydroperoxide formation from reaction with SO 2 and formic acid, respectively. The tropospheric implications of these reactions are evaluated using a global chemistry and transport model.

09 BIOMASS FUELS↗
DOE OSTI2023

Tuned Reactivity at the Lithium Metal–Argyrodite Solid State Electrolyte Interphase

Thin intermetallic Li 2 Te–LiTe 3 bilayer (0.75 µm) derived from 2D tellurene stabilizes the solid electrolyte interphase (SEI) of lithium metal and argyrodite (LPSCl, Li 6 PS 5 Cl) solid-state electrolyte (SSE). Tellurene is loaded onto a standard battery separator and reacted with lithium through single-pass mechanical rolling or transferred directly to SSE surface by pressing. State-of-the-art electrochemical performance is achieved, e.g., symmetric cell stable for 300 cycles (1800 h) at 1 mA cm -2 and 3 mAh cm -2 (25% DOD, 60 µm foil). Cryo-stage focused ion beam (Cryo-FIB) sectioning and Raman mapping demonstrate that the Li 2 Te–LiTe 3 bilayer impedes SSE decomposition. The unmodified Li–LPSCl interphase is electrochemically unstable with a geometrically heterogeneous reduction decomposition reaction front that extends deep into the SSE. Decomposition drives voiding in Li metal due to its high flux to the reaction front, as well as voiding in the SSE due to the associated volume changes. Analysis of cycled SSE found no evidence for pristine (unreacted) lithium metal filaments/dendrites, implying failure driven by decomposition phases with sufficient electrical conductivity that span electrolyte thickness. In conclusion, DFT calculations clarify thermodynamic stability, interfacial adhesion, and electronic transport properties of interphases, while mesoscale modeling examines interrelations between reaction front heterogeneity (SEI heterogeneity), current distribution, and localized chemo-mechanical stresses.

25 ENERGY STORAGE↗
DOE OSTI2024

Criticality Analysis of FSV Spent Nuclear Fuel in the DOE Standard Canister

The U.S. Department of Energy (DOE) is responsible for managing over 300 types of spent nuclear fuel (SNF). To manage this large variety of fuel types, DOE plans to employ standardized canisters for the transportation, long-term storage, and eventual disposal of SNF. Idaho National Laboratory is currently supporting DOE’s SNF Packaging Demonstration Project, in which Fort Saint Vrain (FSV) fuel assemblies will be loaded into a DOE Standard Canister. This paper presents criticality calculations demonstrating that all four or five FSV fuel assemblies loaded into the DOE Standard Canister will remain subcritical in any expected normal or credible abnormal conditions. Previous criticality analyses were performed for one FSV fuel assembly and 12 Peach Bottom Core 2 fuel elements loaded into a DOE Standard Canister. This paper covers the criticality analysis performed for loading both four and five FSV fuel assemblies into a DOE Standard Canister. Various intact and degraded mode configurations were modeled in conducting the criticality calculations. This analysis encompassed three different configurations: (1) a single DOE Standard Canister loaded into a concrete storage overpack, (2) seven DOE Standard Canisters loaded into a concrete storage overpack, and (3) nine DOE Standard Canisters loaded into a concrete storage overpack. The overpack dimensions were varied for each of the three configurations, and transport, storage, and disposal scenarios were analyzed for each configuration. For transport scenarios, a pair of degradation cases was analyzed. In the first case, the fuel compacts became degraded and were removed from the fuel block, then deposited at the bottom of a horizontally placed canister, thereby simulating a drop event. The canister was considered to remain intact. In the second case, the spacing between horizontally placed canisters in a nine-canister overpack was reduced such that the canisters were piled on top of each other, simulating a drop event. For this case, no degradation of the canister internals or fuel was considered. For storage scenarios, the water moderator location in the system was varied to enable identification of the most reactive configurations. Dry and wet conditions were analyzed for the fuel materials, canister, and overpack. For disposal scenarios, two degradation cases were analyzed. In the first, the stainless-steel internals of the canister degraded to either hematite or goethite under both dry and wet conditions. In the second case, degraded FSV fuel formed a uranium-water slurry that filled the coolant/void holes. None of the cases presented exceeded the application specific upper subcritical limit.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2007

The Relationship of Loss, Mean Age of Air and the Distribution of CFCs to Stratospheric Circulation and Implications for Atmospheric Lifetimes

Man-made molecules called chlorofluorcarbons (CFCs) are broken apart in the stratosphere by high energy light, and the reactive chlorine gases that come from them cause the ozone hole. Since the ozone layer stops high energy light from reaching low altitudes, CFCs must be transported to high altitudes to be broken apart. The number of molecules per volume (the density) is much smaller at high altitudes than near the surface, and CFC molecules have a very small chance of reaching that altitude in any particular year. Many tons of CFCs were put into the atmosphere during the end of the last century, and it will take many years for all of them to be destroyed. Each CFC has an atmospheric lifetime that depends on the amount of energy required to break them apart. Two of the gases that were made the most are CFC13 and CF2C12. It takes more energy to break apart CF2C12 than CFC13, and its lifetime is about 100 years, nearly twice as long as the lifetime for CFC13. It is hard to figure out the lifetimes from surface measurements because we don't know exactly how much was released into the air each year. Atmospheric models are used to predict what will happen to ozone and other gases as the CFCs decrease and other gases like C02 continue to increase during the next century. CFC lifetimes are used to predict future concentrations and all assessment models use the predicted future concentrations. The models have different circulations and the amount of CFC lost according to the model may not match the loss that is expected according to the lifetime. In models the amount destroyed per year depends on how fast the model pushes air into the stratosphere and how much goes to high altitudes each year. This paper looks at the way the model circulation changes the lifetimes, and looks at measurements that tell us which model is more realistic. Some models do a good job reproducing the age-of-air, which tells us that these models are circulating the stratospheric air at the right speed. These same models also do a good job reproducing the amount of CFCs in the lower atmosphere where they were measured by instruments on NASA's ER-2, a research plane that flies in the lower stratosphere. The lifetime for CFC13 that is calculated using the models that do the best job matching the data is about 25% longer than most people thought. This paper shows that using these measurements to decide which models are more realistic helps us understand why their predictions are different from each other and also to decide which predictions are more likely.

Douglass, A. R.↗
DOE OSTI2021

Dissolution of periodic arrays of grains: Upscaling of pore-scale simulations with fast reactions

Upscaling methods are frequently used to derive transport equations at the macroscopic scale from more fundamental equations formulated at the pore scale. These methods typically give a suitable structure for the macroscopic equations and can also provide explicit expressions for the constitutive parameters, such as permeability and dispersion coefficients. Introducing chemical reactions complicates upscaling in at least two important ways. First, the interplay between chemical reactions and transport processes introduces a new length scale, which can be much smaller than the convective or dispersive length scales. A small reactive length scale breaks one of the key assumptions in upscaling; that there is a significant separation in length between the pore-scale and macro-scale processes. The second complication is that if reactions take place at mineral surfaces (dissolution or precipitation) then the pore space itself is evolving in time. In this paper we suggest ways in which these difficulties can be approached, based on analysis of pore-scale simulation data. First, we noticed that the concentration field in successive unit cells has an almost identical spatial variation, with a single scaling factor for each unit cell that is proportional to the incoming reactant flux. Using pore-scale simulations to determine the mass transfer coefficient in a few unit cells, we can calculate the concentration field in the whole domain, even when dissolution is entirely transport limited. Second, we have noticed a time-dependent mapping of the grain shapes from different unit cells. From these observations, we can deduce constitutive relations where the only time-varying parameter is the porosity. Here, we show that a model based on these ideas can quantitatively account for the pore-scale simulation data.

58 GEOSCIENCES↗
DOE OSTI2020

Colloid-facilitated transport of 238 Pu, 233 U and 137 Cs through fractured chalk: Laboratory experiments, modelling, and implications for nuclear waste disposal

The influence of montmorillonite colloids on the mobility of 238 Pu, 233 U and 137 Cs through a chalk fracture was investigated to assess the transport potential for radioactive waste. Radioisotopes of each element, along with the conservative tracer tritium, were injected in the presence and absence of montmorillonite colloids into a naturally fractured chalk core. In parallel, batch experiments were conducted to obtain experimental sorption coefficients (K d , mL/g) for both montmorillonite colloids and the chalk fracture material. Breakthrough curves were modelled to determine diffusivity and sorption of each radionuclide to the chalk and the colloids under advective conditions. Uranium sorbed sparingly to chalk (log K d = 0.7 ± 0.2) in batch sorption experiments. 233 U(VI) breakthrough was controlled primarily by the matrix diffusion and sorption to chalk (15 and 25% recovery with and without colloids, respectively). Cesium, in contrast, sorbed strongly to both the montmorillonite colloids and chalk (batch log K d = 3.2 ± 0.01 and 3.9 ± 0.01, respectively). The high affinity to chalk and low colloid concentrations overwhelmed any colloidal Cs transport, resulting in very low 137 Cs breakthrough (1.1–5.5% mass recovery). Batch and fracture transport results, and the associated modelling revealed that Pu migrates both as Pu (IV) sorbed to montmorillonite colloids and as dissolved Pu(V) (7% recovery). Transport experiments revealed differences in Pu(IV) and Pu(V) transport behavior that could not be quantified in simple batch experiments but are critical to effectively predict transport behavior of redox-sensitive radionuclides. Finally, a brackish groundwater solution was injected after completion of the fracture flow experiments and resulted in remobilization and recovery of 2.2% of the total sorbed radionuclides which remained in the core from previous experiments. In general, our study demonstrates consistency in sorption behavior between batch and advective fracture transport. Finally, the results suggest that colloid-facilitated radionuclide transport will enhance radionuclide migration in fractured chalk for those radionuclides with exceedingly high affinity for colloids.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2008

Modeling of the Ablation of Fibrous Materials in the Knudsen Regime

During atmospheric entry of planetary probes, the thermal protection system (TPS) of the probe is exposed to high temperatures under low pressures. In these conditions, carbonous TPS materials undergo gasification in the Knudsen regime leading to mass loss and wall recession called ablation. This work aims to improve the understanding of materiaVenvironment interactions through a study of the coupling between carbon dioxide transport in the Knudsen regime, heterogeneous oxidation of carbon, and sutface recession. A 3D Monte-Carlo simulation tool is used for this study. The fibrous architecture of the materiils, consisting of high porosity random array of carbon fibers, is numerically reproduced on a 3D Cartesian grid. Mass transport in the Knudsen regime from the boundary layer to the surface, and inside this porous material is simulated by random walk. A reaction probability is used to simulate the heterogeneous oxidation reaction. The surface recession is followed by front tracking using a simplified marching cube approach. The output data of the simulations are ablation velocity and dynamic evolution of the material porosity. A parametric study is carried out to analyze the material behavior as a function of Knudsen number for the porous media (length of the mean free path compared to the mean pore diameter) and the intrinsic reactivity of the carbon fibers. The results enable extrapolation of laboratory experimental data to actual entry conditions.

Lachaud, J.↗
DOE OSTI2021

Repetitive diacetyl vapor exposure promotes ubiquitin proteasome stress and precedes bronchiolitis obliterans pathology

Bronchiolitis obliterans (BO) is a devastating lung disease seen commonly after lung transplant but also seen following severe respiratory tract infection or chemical inhalation exposure. Diacetyl (DA; 2,3-butanedione), highly reactive alpha-diketone is known to cause BO when inhaled, however, the mechanisms of how inhalation exposure leads to BO development remains poorly understood. In the current work, we combined two clinically relevant models for studying the pathogenesis of DA-induced BO: (1) an in vivo rat model of repetitive DA vapor exposures with recovery and (2) an in vitro model of primary human airway epithelial cells exposed to pure DA vapors. Rats exposed to 5 consecutive days 200 parts-per-million DA 6 hrs per day had worsening survival, persistent hypoxemia, poor weight gain, and histologic evidence of BO 14 days after DA exposure cessation. At the end of exposure, increased expression of the ubiquitin stress protein, ubiquitin-C, accumulated within DA-exposed rat lung homogenates, localized primarily to the airway epithelium, the primary site of BO development, and did not occur in air-exposed controls. Lung proteasome activity increased concurrently after DA exposure. In primary human airway cultures, global proteomics analysis identified 519 significantly modified proteins in DA-exposed samples relative to controls with the ubiquitin proteasome system, endosomal reticulum transport, and the response to unfolded protein pathways being upregulated while cell-cell adhesion and oxidation-reduction pathways being downregulated in DA-exposed samples. Altogether, repetitive DA vapor exposure resulted in abundant protein damage, accumulation of polyubiquinated proteins, and ubiquitin proteasome stress prior to the development of chemical-induced BO pathology.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2012

Overview of the 2010 Carbonaceous Aerosols and Radiative Effects Study (CARES)

Substantial uncertainties still exist in the scientific understanding of the possible interactions between urban and natural (biogenic) emissions in the production and transformation of atmospheric aerosol and the resulting impact on climate change. The US Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) program's Carbonaceous Aerosol and Radiative Effects Study (CARES) carried out in June 2010 in Central Valley, California, was a comprehensive effort designed to improve this understanding. The primary objective of the field study was to investigate the evolution of secondary organic and black carbon aerosols and their climate-related properties in the Sacramento urban plume as it was routinely transported into the forested Sierra Nevada foothills area. Urban aerosols and trace gases experienced significant physical and chemical transformations as they mixed with the reactive biogenic hydrocarbons emitted from the forest. Two heavily-instrumented ground sites - one within the Sacramento urban area and another about 40 km to the northeast in the foothills area - were set up to characterize the evolution of meteorological variables, trace gases, aerosol precursors, aerosol size, composition, and climaterelated properties in freshly polluted and "aged" urban air. On selected days, the DOE G-1 aircraft was deployed to make similar measurements upwind and across the evolving Sacramento plume in the morning and again in the afternoon. The NASA B-200 aircraft, carrying remote sensing instruments, was also deployed to characterize the vertical and horizontal distribution of aerosols and aerosol optical properties within and around the plume. This overview provides: (a) the scientific background and motivation for the study, (b) the operational and logistical information pertinent to the execution of the study, (c) an overview of key observations and initial findings from the aircraft and ground-based sampling platforms, and (d) a roadmap of planned data analyses and focused modeling efforts that will facilitate the integration of new knowledge into improved representations of key aerosol processes and properties in climate models.

aerosols↗
DOE OSTI2022

Quantifying Subsurface Flow and Solute Transport in a Snowmelt‐Recharged Hillslope With Multiyear Water Balance

Abstract Quantifying flow and transport from hillslopes is vital for understanding water quantity and quality in rivers, but remains obscure because of limited subsurface measurements. Using measured hydraulic conductivity K profiles and water balance over a single year to calibrate a transmissivity feedback model for a hillslope in the East River watershed (Colorado) proved unsatisfactory for predicting flow over the subsequent years. Well‐constrained field‐scale K were obtained by optimizing subsurface flux predictions over years having large differences in recharge, and by including estimates of interannual transfer of excess snowmelt recharge. Water and solute exports during high snowmelt recharge occur predominantly via shallow groundwater flow through weathered rock and soil because of their enlarged transmissivities under saturated conditions. Conversely, these shallow pathways are less active in snow drought years when the water table remains deeper within the weathering zone. Hillslope soil water monitoring showed that rainfall does not infiltrate deeply during summer and fall months, and revealed water losses consistent with model ET predictions. By combining water table‐dependent fluxes with pore water chemistry in different zones, time‐dependent rates of solute exports become predictable. As an example, calibrated K were combined with dissolved nitrogen concentrations in pore waters to show the snowmelt‐dependence of reactive nitrogen exported from the hillslope, further supporting the recent finding that the weathering zone is the dominant source of reactive nitrogen at this site. Subsurface export predictions can now be obtained for wide ranges of recharge based on measurements of water table elevation and profiles of pore water chemistry.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Materials genome innovation for computational software (magics) center

Functional layered material (LM) architectures will dominate nanomaterials science in this century. We have developed theory, modeling, simulation, and software and data tools that enhance understanding and AI guide synthesis, enable characterization of complex structures, and improve capabilities in the predictive design and growth of LMs. Research at the Center has focused on: Computational synthesis and characterization: AI guided synthesis and experimental synthesis of stacked LMs with tailored properties via optimized chemical vapor deposition (CVD) growth and liquid-phase exfoliation; study defects, edges, grain boundaries, wrinkling of atomic layers and their effects on chemical, mechanical, electrical, and optical properties. Far-from-equilibrium processes: Joint experimental and simulation based probe of electronic processes with NAQMD and ultrafast X-ray free-electron laser (XFEL) and ultrafast electron diffraction (UED) facilities at Stanford. Experimentally validate NAQMD by ultrafast electron diffraction and X-ray spectroscopy studies of structural and excited state dynamics, shape fluctuations, and phonon dynamics. Scalable software: Simulation engines for desktop-to-exascale platforms using low-overhead, linear-scaling QMD algorithms; divide-conquer-recombine NAQMD with electronic excitations; extended-Lagrangian reactive molecular dynamics (RMD), machine learning (ML) based neural-network quantum molecular dynamics (NNQMD), and super-state accelerated molecular dynamics (AMD) and kinetic Monte Carlo codes; thermal and electrical transport software; and design 3D architectures of LMs with desired functionality using scalable software. Distribution of software and data, and training: Software and simulation-experimental data generated within the Center are distributed to the materials science community via Berkeley Materials Project (MP) framework. We have also organized three workshops for software distribution and training at USC (Nov. 2017, Mar. 2018) and Gaithersburg, MD (Nov. 2018) to train researchers, with the last one in focused on underrepresented groups, in collaboration with Howard University which is one of the largest HBCUs. The Center supported a total of 46 personnel and 6 undergraduate students. These include 14 faculty, 11 postdoctoral research associates, 20 graduate research assistants, and mentored 6 undergraduate students. This resulted in the publications of 63 research papers that include 46 publications on Reactive and Quantum Dynamics Simulations, 13 publications on Machine Learning for Quantum Materials, and 4 publications on Quantum Computing.

2D Materials↗
NASA NTRS1997

Reactive Nitrogen and its Correlation with O3 and CO Over the Pacific in Winter and Early Spring

Measurements of NO, NO(y), O3, and CO were made during NASA's Global Tropospheric Experiment/Pacific Exploratory Mission-West B (GTE/PEM-West B) carried out over the western Pacific in February and March 1994. NO(x) was calculated from NO using a photostationary state model ((NO(x)(sub mc)). Correlations between these species are presented, and some insights into the sources of NO(x) and NO(y) are described. The boundaries between the lower, middle, and upper troposphere have been defined at potential temperatures of 311 K and 328 K, which correspond to the geometric altitudes of about 5 and 9 km at 30degN. Enhancements in the mixing ratios of NO(y) and CO were observed in the lower and middle troposphere. A positive correlation was found between these two species suggesting that the high NO(y) values were due to anthropogenic emissions over the continental surface. On the other hand, O3 increased little with increase in CO. As a result, NO(y)/O3 ratios were higher in air more influenced by pollution. NO(y), values in 55 and 28% of the air masses sampled in the lower and middle troposphere, respectively, were higher than the clean free tropospheric NO(y)-O3 range when O3 values simultaneously observed were used. High (NOx)mc/NOy ratios between 0.15 and 0.3 were found in the boundary layer with relatively low mixing ratios of CO and NOy during the three flights. These air masses were transported from a higher altitude (approximately 5 km) and a higher latitude (approximately 50degN) within a few days. The peroxyacetyl nitrate (PAN)/NO(y) ratios were generally high (approximately 0.4) in these air masses, and the thermal decomposition of PAN was a probable source of NO(x). In the middle troposphere the (NO(x))mc mixing ratio did not generally increase with NO(y) or CO, suggesting that the transport of air masses affected by anthropogenic emissions did not increase the NO(x) level significantly. In the upper troposphere, very minor effects from the continental surface sources were seen in the CO mixing ratio. By contrast, NO(y) values in 33% of the air masses were higher than those expected when stratospheric air intrusion is assumed to be a single source of NO(y) based on NO(y)-O3 correlation analyses. This result suggests significant free tropospheric NO(y) sources, namely exhaust from the aircraft and NO production by lightning activity. In fact, spikes in the (NO(x))(sub m)c mixing ratios were observed near the aircraft corridor south of Tokyo at an altitude of 10 km. These two free tropospheric NO(x) sources were considered to be important in determining the levels of the upper tropospheric NO(x) and NO(y) during PEM-West B.

Koike, M.↗
DOE OSTI2023

Geomechanical in situ testing of fault reactivation in argillite repositories

Abstract. Pressurization of natural faults as a result of repository-induced effects can lead to their reactivation and permeability generation in case such features are present near disposal tunnels. Potential driving forces for such pressurization are the temperature increase caused by heat-producing high-level radioactive waste and the generation of hydrogen and other gases due to corrosion of engineered materials. We are primarily concerned about pressurization and fault activation in host rocks and faults that have very low natural permeability for pore pressure increases to dissipate, such as the argillite rocks currently investigated in Switzerland, France and other countries. This presentation discusses a series of in situ experiments of fault activation by fluid injection conducted in the argillite rock (Opalinus clay) at the Mont Terri underground research laboratory in Switzerland. A multi-model monitoring test bed was installed at Mont Terri that includes distributed fiber optics for strain, temperature, and acoustics; local fault pore pressure and three-dimensional displacement sensors; active seismic imaging; and passive seismic monitoring. The fault experiments (and their subsequent analysis via hydromechanical modeling) provide a new fundamental understanding of the coupling between pore pressure, fault deformation and permeability generation as a function of time and allow exploring how seismic and aseismic events may impact the integrity of faulted argillite host rock. Our experimental observations demonstrate that significant flow (and transport) can occur along the initially impermeable argillite fault when rupture is activated. However, the fault permeability decreases to almost its pre-activation value when fluid injection ceases and fluid pressure drops. Dilatant slip on the fault plane alone does not explain the observed increase in fault permeability; pressure-induced fault opening also plays a role, favored by the softness of the shale along with the fact that the structure of the fault zone prevents fluids from diffusing into the adjacent damage zone. Rupture initiation and permeability generation is initially aseismic, associated only with an increase in noise level and emerging tremors. Micro-earthquakes are initiated later in the experiments and typically occur away from the fluid-pressurized area. Hydromechanical models show that stress transferred from the initial aseismic deformation can build up to stress criticality and later induce seismic rupture. After presenting the experimental results, we will close the presentation with an outlook to future experimental campaigns using the fault test bed at Mont Terri. We are currently planning a controlled thermal stimulation of the fault, where instead of fluid injection as a trigger mechanism we will heat up the nearby rock volume and measure potential effects on fault stability. Lessons learned from our past and future experiments will help inform the safety assessment of geologic disposal in argillite host rock.

Birkholzer, Jens T.↗
DOE OSTI2020

Feasibility of using reactive silicate particles with temperature-responsive coatings to enhance the security of geologic carbon storage

The large-scale deployment of geologic carbon storage will require minimizing the risk of CO 2 leakage. Here we report on a novel strategy to seal leakage pathways by injecting functionalized reactive calcium silicate (CaSiO 3 ) particles. The particles consist of a mineral silicate core within a temperature-sensitive Poly(N,N-dimethylaminoethyl methacrylate) coating that limits contact between the core and the surrounding CO 2 (aq) above a switching temperature making them relatively unreactive deep in the subsurface. If the CO 2 leaks to shallower depths and cooler temperatures, the coatings would extend and enable the silicate core to react with CO 2(aq) to form solid carbonate precipitates that reduce the permeability of the fracture or flow path. A combined laboratory- and simulation-based investigation was conducted to explore key early-stage research questions related to the feasibility of this leakage management strategy. Sand columns were injected with either coated or uncoated CaSiO 3 particles and the permeability was measured before and after exposure to CO 2 . Experiments were conducted at 35 °C and 70 °C, a range that spans the switching temperature for the coatings. Results show that at the cooler temperature the coating effectively enabled reaction and permeability reduction. At higher temperatures, the coated particles limited mineral carbonation and the permeability decrease was minimal compared to columns with uncoated particles. Here, the morphology and location of the precipitates provides insight about permeability evolution. Model simulations of 2D radial particle transport from an injection point provided insight into the effects of particle size, coating properties, and concentrations of particles in several representative rock formations.

54 ENVIRONMENTAL SCIENCES↗