Search NASA⌕ Search

SEARCH · Search NASA

Results for “direct sampling”

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

Zinc isotope constraints on the cycling of carbon in the Bermuda mantle source

Volatile recycling and storage in the mantle transition zone (MTZ) is important for the refertilization of the upper mantle and is associated with the generation of high-µ (HIMU, where µ is 238 U/ 204 Pb) mantle. One way to probe the MTZ and the processes associated with mantle convection is to sample lavas that originate from the shallow mantle and were contaminated by upwelling from the MTZ, such as at the previously proposed shallow plume of Bermuda. Here we present the first δ 66 Zn isotopic compositions of Bermuda silica-undersaturated and silica-saturated lavas to explore the origin of the carbon-rich lithologies and the genesis of the large seamount found in the western North Atlantic Ocean. Contrasting with global δ 66 Zn data sets, our results (δ 66 Zn between 0.24 ± 0.04 and 0.41 ± 0.04) do not support direct sampling of recycled marine carbonates in the Bermuda HIMU mantle. Instead, we show that δ 66 Zn fractionation toward higher values is associated with magmatic processes and incorporation of carbon sourced from deep fluids associated with the formation of carbonatites. These carbon-rich fluids are likely sourced from the metasomatic reactions between the subducted cold slab of the Iapetus oceanic lithosphere ca. 500 Ma and the thickened continental lithospheric mantle of Pangea. Melting of this metasomatized mantle was triggered by the arrival of the Farallon slab to the eastern North American margin in the late Cenozoic via shallow convection.

Mazza, Sarah E. [Smith College, Northampton, MA (U↗
DOE OSTI2023

An attribution of the low single-scattering albedo of biomass burning aerosol over the southeastern Atlantic

Abstract. Aerosol over the remote southeastern Atlantic is some of the most sunlight-absorbing aerosol on the planet: the in situ free-tropospheric single-scattering albedo at the 530 nm wavelength (SSA530 nm) ranges from 0.83 to 0.89 within ORACLES (ObseRvations of Aerosols above CLouds and their intEractionS) aircraft flights from late August–September. Here we seek to explain the low SSA. The SSA depends strongly on the black carbon (BC) number fraction, which ranges from 0.15 to 0.4. Low organic aerosol (OA)-to-BC mass ratios of 8–14 and modified combustion efficiency values >0.975 point indirectly to the dry, flame-efficient combustion of primarily grass fuels, with back trajectories ending in the miombo woodlands of Angola. The youngest aerosol, aged 4–5 d since emission, occupied the top half of a 5 km thick plume sampled directly west of Angola with a vertically consistent BC:ΔCO (carbon monoxide) ratio, indicating a homogenization of the source emissions. The younger aerosol, transported more quickly off of the continent by stronger winds, overlaid older, slower-moving aerosol with a larger mean particle size and fraction of BC-containing particles. This is consistent with ongoing gas condensation and the coagulation of smaller non-BC particles upon the BC-containing particles. The particle volumes and OA:BC mass ratios of the older aerosol were smaller, attributed primarily to evaporation following fragmentation, instead of dilution or thermodynamics. The CLARIFY (CLoud–Aerosol–Radiation Interaction and Forcing: Year 2017) aircraft campaign sampled aerosols that had traveled further to reach the more remote Ascension Island. CLARIFY reported higher BC number fractions, lower OA:BC mass ratios, and lower SSA yet larger mass absorption coefficients compared to this study's. Values from one ORACLES 2017 flight, held midway to Ascension Island, are intermediate, confirming the long-range changes. Overall the data are most consistent with continuing oxidation through fragmentation releasing aerosols that subsequently enter the gas phase, reducing the OA mass, rather than evaporation through dilution or thermodynamics. The data support the following best fit: SSA530nm=0.801+0055⋅(OA:BC) (r=0.84). The fires of southern Africa emit approximately one-third of the world's carbon; the emitted aerosols are distinct from other regional smoke emissions, and their composition needs to be represented appropriately to realistically depict regional aerosol radiative effects.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Angle-Resolved Polarized Raman Study of Layered Cr 2 Se 3

The polarization-resolved Raman spectra of two-dimensional Cr 2 Se 3 synthesized via chemical vapor deposition (CVD) and chemical vapor transport (CVT) techniques were investigated in detail. The samples were characterized using X-ray diffraction (XRD), transmission electron microscopy (TEM), and energy-dispersive X-ray spectroscopy (EDS). A distinct polarization dependence was observed in the Raman intensity of all the Cr-Cr, Cr-Se, and Se-Se modes in both samples. The observed angle-dependent Raman intensities of each peak could be related to the crystal structure-specific Raman tensor. XRD results of the bulk Cr 2 Se 3 sample synthesized via CVT confirm its trigonal crystal structure, and the Raman peaks can be fitted using the Raman tensors for the A g and E g modes for both the parallel and crossed polarizations. However, for the Cr 2 Se 3 samples directly grown on Si/SiO 2 substrates by CVD, it was necessary to assume the triclinic crystal structure in order to explain the polarized Raman dependence of all the peaks in both parallel and crossed polarization directions. Furthermore, this is the first experimental result suggesting the existence of triclinic Cr 2 Se 3 crystal structure, which has been theoretically predicted in the Materials Project database.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Formation of PAHs, phenol, benzofuran, and dibenzofuran in a flow reactor from the oxidation of ethylene, toluene, and n-decane

The necessity to reveal the formation mechanism of not only polycyclic aromatic hydrocarbons (PAHs) but also oxygenated PAHs (OPAHs) during combustion is increasing. Although many studies on PAHs have been conducted, fundamental studies investigating OPAH formation are still limited. Phenol, benzofuran, and dibenzofuran were selected as OPAHs in this study. We experimentally investigated fuel-rich oxidation in a flow reactor at atmospheric pressure, mean gas temperatures from 1050 to 1350 K, residence times from 0.2 to 1.5 s, and equivalence ratios from 3.0 to 12.0. Ethylene, toluene, and n-decane were used as fuels. Three kinds of OPAHs as well as 23 kinds of PAHs including monocyclic structures were quantitatively measured through a direct sampling using gas chromatography mass spectrometry. The results showed that concentrations of PAHs and OPAHs were strongly affected by the fuel type, with increasing concentrations in the following order: ethylene < n-decane < toluene. A chemical kinetic model for OPAH formation was developed based on a recent PAH growth model. The model showed that the predicted concentrations of OPAHs and PAHs were in reasonable agreement with the measured data in this study. Some modifications were made to the previous model based on recent literature studies and on the comparison of the simulated and measured results. The effect of the fuel type on the formation of PAHs and OPAHs was investigated through kinetic analysis using the model to discuss their reaction pathways.

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

A theoretical kinetic study of ĊH 3 + ṄH 2 : From electronic structure to NH 3 /CH 4 combustion modelling implications

Carbon–nitrogen interaction reactions play an important role in governing the reactivity of ammonia blended fuels. However, there remains uncertainties regarding their detailed reaction pathways and rate constants, hampering the development of high-fidelity chemical kinetic models. In this study, the kinetics of ĊH 3 + ṄH 2 , a key C–N interaction reaction in ammonia/methane blend combustion have been investigated. The potential energy surface has been explored using the high-level ANL0F method, yielding highly accurate stationary point energies that agree with ATcT values within 0.1 kcal mol –1 . Variable reaction coordinate transition state theory is used to treat the barrierless association and decomposition reaction channels, based on directly sampled radical-radical interaction energies at the CASPT2-F12(2e,2o)/cc-pVTZ-F12 level of theory. The minimum transitional mode numbers of states obtained are then coupled with the RRKM/master equation to calculate temperature- and pressure-dependent rate constants. Our a priori calculations capture available experimental measurements from the literature very well. The calculated rate constants have been incorporated into an NH 3 /CH 4 chemical kinetic model currently under development at the University of Galway. The effect of the updated kinetic data for ĊH 3 + ṄH 2 on model predicted NH 3 /CH 4 fuel reactivity is elucidated.

ab initio↗
DOE OSTI2025

Enhancing approximate modular Bayesian inference by emulating the conditional posterior

In modular Bayesian analyses, complex models are composed of distinct modules, each representing different aspects of the data or prior information. In this context, fully Bayesian approaches can sometimes lead to undesirable feedback between modules, compromising the integrity of the inference. The “cut-distribution” prevents unwanted influence between modules by “cutting” feedback. The direct sampling (DS) algorithm is standard practice for approximating the cut-distribution, but it can be computationally intensive, especially when the number of imputations required is large. An enhanced method is proposed, the Emulating the Conditional Posterior (ECP) algorithm, which leverages emulation to increase the number of imputations. Through numerical experiment it is demonstrated that the ECP algorithm outperforms the traditional DS approach in terms of accuracy and computational efficiency, particularly when resources are constrained. Here, it is also shown how the DS algorithm can be improved using ideas from design of experiments. Some practical recommendations are given for algorithm choice in modular Bayesian analyses.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Stratigraphic Identification with Airborne Electromagnetic Methods at the Hanford Site, Washington

Stratigraphic units can influence the fate and transport of subsurface contaminants within groundwater. Units having coarse-grained sediments act as preferential flow pathways, and therefore can accelerate the transport of contaminants to reach human and ecological receptors. At legacy waste sites, detailed knowledge of subsurface stratigraphy can be used for effective monitoring and remediation planning to help minimize risk to human health and the environment. Airborne electromagnetic (AEM) methods can non-invasively provide information on kilometer-scale or larger subsurface stratigraphic features and fill informational gaps in directly sampled data from sparsely located boreholes. In this paper, we present inversion results of a 412 line-km frequency-domain AEM survey to delineate subsurface stratigraphic features at the Hanford Site, located in southeastern Washington State. The inversion was performed using a massively parallel 3D electromagnetic modeling and inversion code, where the modeling is based on solving frequency-domain Maxwell’s equations using an unstructured-mesh finite-element method and the inversion employs a Gauss-Newton optimization scheme. The results are compared to an underlying geologic framework model (GFM), built by interpolating contact depths of stratigraphic units interpreted from site borehole datasets. In areas with good borehole coverage, the inversion results show a good match with the GFM to a depth of about 60 m. Outside of these areas, the inversion results exhibit inconsistencies from the assumptions made to create the GFM, demonstrating that the AEM survey results can be used to improve the understanding of the geological conceptual model.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Single-Cell Metabolic Profiling: Metabolite Formulas from Isotopic Fine Structures in Heterogeneous Plant Cell Populations

Characterization of the metabolic heterogeneity in cell populations requires the analysis of single cells. Most current methods in single-cell analysis rely on cell manipulation, potentially altering the abundance of metabolites in individual cells. Small sample volume and the chemical diversity of metabolites are additional challenges in single-cell metabolomics. Here, we describe the combination of fiber-based laser ablation electrospray ionization (f-LAESI) with 21 tesla Fourier transform ion cyclotron resonance mass spectrometry (21TFTICR-MS) for in situ single-cell metabolic profiling in plant tissue. Single plant cells infected by bacteria were selected and sampled directly from the tissue without cell manipulation through mid-infrared ablation with a fine optical fiber tip for ionization by f-LAESI. Finally, ultrahigh performance 21T-FTICR-MS enabled the simultaneous capture of isotopic fine structures (IFSs) for 47 known and 11 unknown compounds, thus elucidating their elemental compositions from single cells, and providing information on metabolic heterogeneity in the cell population.

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

Monte Carlo MP2-F12 for Noncovalent Interactions: The C 60 Dimer

A scalable stochastic algorithm is presented that can evaluate explicitly correlated (F12) second-order many-body perturbation (MP2) energies of weak, noncovalent, intermolecular interactions. It first transforms the formulas of the MP2 and F12 energy differences into a short sum of high-dimensional integrals of Green’s functions in real space and imaginary time. Furthermore, these integrals are then evaluated by the Monte Carlo method augmented by parallel execution, redundant-walker convergence acceleration, direct-sampling autocorrelation elimination, and control-variate error reduction. By sharing electron-pair walkers across the supermolecule and its subsystems spanned by the joint basis set, the statistical uncertainty is reduced by one to 2 orders of magnitude in the MP2 binding energy corrected for the basis-set incompleteness and superposition errors. The method predicts the MP2-F12/aug-cc-pVDZ binding energy of 19.1 ± 4.0 kcal mol –1 for the C 60 dimer at the center distance of 9.748 Å.

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

Selective One-Dimensional 13 C– 13 C Spin-Diffusion Solid-State Nuclear Magnetic Resonance Methods to Probe Spatial Arrangements in Biopolymers Including Plant Cell Walls, Peptides, and Spider Silk

Two-dimensional (2D) and 3D through-space 13 C– 13 C homonuclear spin-diffusion techniques are powerful solid-state nuclear magnetic resonance (NMR) tools for extracting structural information from 13 C-enriched biomolecules, but necessarily long acquisition times restrict their applications. In this work, we explore the broad utility and underutilized power of a chemical shift-selective one-dimensional (1D) version of a 2D 13 C– 13 C spin-diffusion solid-state NMR technique. The method, which is called 1D dipolar-assisted rotational resonance (DARR) difference, is applied to a variety of biomaterials including lignocellulosic plant cell walls, microcrystalline peptide fMLF, and black widow dragline spider silk. 1D 13 C– 13 C spin-diffusion methods described here apply in select cases in which the 1D 13 C solid-state NMR spectrum displays chemical shift-resolved moieties. Furthermore, this is analogous to the selective 1D nuclear Overhauser effect spectroscopy (NOESY) experiment utilized in liquid-state NMR as a faster (1D instead of 2D) and often less ambiguous (direct sampling of the time domain data, coupled with increased signal averaging) alternative to 2D NOESY. Selective 1D 13 C– 13 C spin-diffusion methods are more time-efficient than their 2D counterparts such as proton-driven spin diffusion (PDSD) and dipolar-assisted rotational resonance. The additional time gained enables measurements of 13 C– 13 C spin-diffusion buildup curves and extraction of spin-diffusion time constants TSD, yielding detailed structural information. Specifically, selective 1D DARR difference buildup curves applied to 13 C-enriched hybrid poplar woody stems confirm strong spatial interaction between lignin and acetylated xylan polymers within poplar plant secondary cell walls, and an interpolymer distance of ~0.45–0.5 nm was estimated. Additionally, Tyr/Gly long-range correlations were observed on isotopically enriched black widow spider dragline silks.

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

Application of sensitivity analysis in DYMOND/Dakota to fuel cycle transition scenarios

The ability to perform sensitivity analysis has been enabled for the nuclear fuel cycle simulator DYMOND through its coupling with the design and analysis toolkit Dakota. To test and demonstrate these new capabilities, a transition scenario and multi-parameter study were devised. The transition scenario represents a partial transition from the US nuclear fleet to a closed fuel cycle with small modular LWRs and fast reactors fueled by reprocessed used nuclear fuel. Four uncertain parameters in this transition were studied – start date of reprocessing, total reprocessing capacity, the nuclear energy demand growth, and the rate at which the fast reactors are deployed – with respect to their impact on four response metrics. The responses – total natural uranium consumed, maximum annual enrichment capacity required, total disposed mass, and total cost of the nuclear fuel cycle – were chosen based on measures known to be of interest in transition scenarios and to be significantly impacted by the varying parameters. Furthermore, analysis of this study was performed both from the direct sampling and through surrogate models developed in Dakota to calculate the global sensitivity measures Sobol’ indices. This example application of this new capability showed that the most consequential parameter to most metrics was the share of new build capacity that is fast reactors. However, for the cost metric, the scaling factor of the energy demand growth was significant and had synergistic behavior with the fast reactor new build share.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Resonance Scattering Treatment with the Windowed Multipole Formalism

A new method for directly sampling the resonance upscattering effect is presented. Alternatives have relied on inefficient rejection sampling techniques or large tabular storage of relative velocities. None of these approaches, which require pointwise energy data, are particularly well suited to the windowed multipole cross-section representation. The new method, called multipole analytic resonance scattering, overcomes these limitations by inverse transform sampling from the target relative velocity distribution where the cross section is expressed in the multipole formalism. The closed-form relative speed distribution contains a novel special function we deem the incomplete Faddeeva function, and we present the first results on its efficient numerical evaluation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Magnetic tunnel junction random number generators applied to dynamically tuned probability trees driven by spin orbit torque

Abstract Perpendicular magnetic tunnel junction (pMTJ)-based true-random number generators (RNGs) can consume orders of magnitude less energy per bit than CMOS pseudo-RNGs. Here, we numerically investigate with a macrospin Landau–Lifshitz-Gilbert equation solver the use of pMTJs driven by spin–orbit torque to directly sample numbers from arbitrary probability distributions with the help of a tunable probability tree. The tree operates by dynamically biasing sequences of pMTJ relaxation events, called ‘coinflips’, via an additional applied spin-transfer-torque current. Specifically, using a single, ideal pMTJ device we successfully draw integer samples on the interval [0, 255] from an exponential distribution based on p -value distribution analysis. In order to investigate device-to-device variations, the thermal stability of the pMTJs are varied based on manufactured device data. It is found that while repeatedly using a varied device inhibits ability to recover the probability distribution, the device variations average out when considering the entire set of devices as a ‘bucket’ to agnostically draw random numbers from. Further, it is noted that the device variations most significantly impact the highest level of the probability tree, with diminishing errors at lower levels. The devices are then used to draw both uniformly and exponentially distributed numbers for the Monte Carlo computation of a problem from particle transport, showing excellent data fit with the analytical solution. Finally, the devices are benchmarked against CMOS and memristor RNGs, showing faster bit generation and significantly lower energy use.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Determination of proton PDF uncertainties with Markov chain Monte Carlo

We present an analysis of parton distribution functions (PDFs) of the proton using Markov chain Monte Carlo (MCMC) methods. The MCMC approach naturally implements Bayes’ theorem and, thus, provides a means to directly sample the underlying probability distribution—in this case, the probability distribution of the PDF parameters. This allows for a straightforward propagation of the resulting uncertainties into any PDF-dependent observable, preserving their simple probabilistic interpretation. In our analysis we include a broad set of deep inelastic scattering data from HERA, BCDMS and NMC experiments along with the Drell-Yan, 𝑊 and 𝑍 boson data from LHC and Tevatron experiments, which combined with theoretical calculations at next-to-next-to-leading order in QCD allow for realistic determination of PDFs. The main focus of this analysis is to explore alternative methods for PDF uncertainty estimation that are more firmly grounded in statistical principles. We show that the flexibility of the Bayes framework, allowing one, e.g., to account for non-Gaussianity or inconsistencies of datasets, is crucial to extract realistic uncertainties when such assumptions are not fulfilled. We also demonstrate that MCMC allows one to determine the Δ⁢𝜒 2 value corresponding to a given confidence level in the sample, which can, in turn, be used as a statistically well-founded tolerance criterion used in the Hessian method, thus addressing one of its main long-standing drawbacks.

Risse, Peter Clemens [Universität Münster (Germany↗
DOE OSTI2026

Spatially Accelerated Winding Numbers for Curved Geometry

The generalized winding number (GWN) is a scalar field that supports robust containment queries on curved geometry, including non-watertight, overlapping, and nested boundary representations. While queries can be easily parallelized over samples, direct evaluation on parametric curves and surfaces remains costly for large and complex models. Fast, state-of-the-art GWN approaches leverage a spatial index to approximate the GWN, typically coupled with a Taylor expansion which approximates the GWN contribution for far clusters of geometric primitives. However, such methods operate only on discrete inputs such as triangle meshes and point clouds, and would introduce containment errors near boundaries if applied to curved input. We extend support for fast GWN evaluation over arbitrary collections of NURBS curves in 2D and trimmed NURBS patches in 3D via a Bounding Volume Hierarchy that stores efficiently precomputed moment data in the hierarchy nodes. When querying the hierarchy, approximations for far clusters are used alongside direct evaluation for nearby NURBS primitives, achieving sub-linear complexity while preserving the geometric features in the vicinity of the query point. Central to our performance improvements is an adaptive subdivision strategy for NURBS primitives during a preprocessing phase, creating better spatial partitions while retaining the same accuracy for containment decisions as a direct evaluation. We demonstrate the performance and accuracy of our approach across a large collection of 2D and 3D datasets.

Computer science↗
DOE OSTI2024

Environmental associations of Ophidiomyces ophidiicola , the causative agent of ophidiomycosis in snakes

Emerging pathogenic fungi have become a topic of conservation concern due to declines observed in several host taxa. One emerging fungal pathogen, Ophidiomyces ophidiicola, is well documented as the causative agent of ophidiomycosis, otherwise known as snake fungal disease (SFD). O. ophidiicola has been found to cause disease in a variety of snake species across the United States, including the eastern massasauga (Sistrurus catenatus), a federally threatened rattlesnake species. Most work to date has involved detecting O. ophidiicola for diagnosis of infection through direct sampling of snakes, and attempts to detect O. ophidiicola in the abiotic environment to better understand its distribution, seasonality, and habitat associations are lacking. We collected topsoil and groundwater samples from four macrohabitat types across multiple seasons in northern Michigan at a site where Ophidiomyces infection has been confirmed in eastern massasauga. Using a quantitative PCR (qPCR) assay developed for O. ophidiicola, we detected Ophidiomyces DNA in topsoil but observed minimal to no detection in groundwater samples. Detection frequency did not differ between habitats, but samples grouped seasonally showed higher detection during mid-summer. We found no relationships of detection with hypothesized environmental correlates such as soil pH, temperature, or moisture content. Furthermore, the distribution of Ophidiomyces positive samples across the site was not linked to estimated space use of massasaugas. Our data suggests that season has some effect on the presence of Ophidiomyces. Differences in presence between habitats may exist but are likely more dependent on the time of sampling and currently uninvestigated soil or biotic parameters. These findings build on our understanding of Ophidiomyces ecology and epidemiology to help inform where and when snakes may be exposed to the fungus in the environment.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Development of Mockups and Instrumentation for Spent Fuel Drying Tests

The purpose of this report is to provide updates on the experimental components, methodology, and instrumentation under development for use in advanced studies of realistic drying operations conducted on surrogate spent nuclear fuel. Validation of the extent of water removal in a dry spent nuclear fuel storage system based on drying procedures used at nuclear power plants is needed to close existing technical gaps. Operational conditions leading to incomplete drying may have potential impacts on the fuel, cladding, and other components in the system. Water remaining in canisters upon completion of drying procedures can lead to cladding corrosion, embrittlement, and breaching, as well as fuel degradation. Additional information is needed on the drying process efficacy to help evaluate the potential impacts of water retention on extended longterm dry storage. A general lack of data suitable for model validation of commercial nuclear canister drying processes necessitates additional, well-designed investigations. Smaller-scale tests that incorporate relevant physics and well-controlled boundary conditions are essential to provide insight and guidance to the simulation of prototypic systems undergoing drying processes. This report describes the implementation of moisture monitoring equipment on a pressurized, submersible system employing a single waterproof, electrically heated spent fuel rod simulator as a demonstration of analytical capabilities during a drying process. A mass spectrometer with specially designed inlets was used to monitor moisture and other gases at 150 kPa to 800 kPa for a test simulating a forced helium dehydration procedure and below 1 torr for tests mimicking a vacuum drying process. The dew point data from the mass spectrometer was found to be in good agreement with a solid-state moisture probe. A distinct advantage of the mass spectrometer system was the capability to directly sample from the hightemperature (>200 °C) head space expected in a prototypic scale experiment where a solid-state moisture probe would suffer considerable loss of accuracy or fail altogether. The operational and analytical experiences gained from this test series are poised to support an expansion to assembly-scale tests at prototypic length. These assemblies are designed to feature prototypic assembly hardware, advanced diagnostics for in situ internal rod pressure monitoring, and failed fuel rod simulators with engineered cladding defects to challenge the drying system with waterlogged fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

A radioisotope - enabled reactive transport model for deep vadose zone carbon

In mountainous regions, which constitute the principle source of recharge to major rivers and regional aquifers, infiltration occurs through fractured, partially saturated, weathered bedrock that acts as a boundary layer between saturated aquifers and surface soil. Commonly this deep vadose zone (DVZ) is many meters thick, and yet its role in regulating the generation, retention and mobility of reactive solutes, including nutrients, contaminants and weathering products, is largely unknown. In particular, many of the key reactions that drive the formation of the weathered DVZ and the quality of water moving through it are redox processes, regulated by the availability of organic carbon and oxygen below the soil layer. The role of the DVZ is thus also poorly constrained in the context of carbon stocks and mobility, particularly in lithologies that are naturally high in organic carbon, such as shales. The overarching hypothesis of this study is that upland regions developed in geologic settings with abundant petrogenic carbon store and actively cycle carbon in the weathered DVZ below the soil and above the water table at rates that are significant and currently unconstrained. In order to quantify this cycling, the current study combines novel instrumentation techniques allowing new direct sampling of DVZ systems with advanced numerical reactive transport simulations of carbon transport and transformation. Critically, these simulations will explicitly treat the three isotopes of carbon (the abundant 12C, the stable rare 13C and the radioactive 14C) in a unified framework, thus clearly parsing between the contributions of modern surface derived carbon and lithologic carbon sources in integrated measurements of fluid and gas phase fluxes. This novel model capability will be applied to test the role of DVZ carbon cycling as a regulator of water quality and geological weathering in two complementary field sites both located in organic carbon rich shale lithologies. The first is the Eel River Critical Zone Observatory (ERCZO) in Mendocino County, California, and the second is the Lawrence Berkeley National Laboratory Watershed Function Scientific Focus Area (SFA) in the East River watershed, near Crested Butte, Colorado. At the ERCZO, a novel vadose zone monitoring system has been installed in a 20 m thick, partially saturated, weathered shale hillslope, and preliminary data already indicate substantial CO2 flux generated many meters below the soil surface. At the SFA field site, an instrumented hillslope transect indicates a more complex multi-dimensional fluid and solute transport regime, which will serve as a key test of the calibrated models. Collectively, this project will advance understanding of the cycling of carbon belowground and in relation to transport pathways across the poorly constrained DVZ characteristic of primary water recharge areas. The key product of this work will be enhanced isotope simulation capabilities that are robust and publicly available for application across a broad diversity of systems.

58 GEOSCIENCES↗