Search NASA⌕ Search

SEARCH · Search NASA

Results for “Subsurface”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Atomic Dynamics of Multi‐Interfacial Migration and Transformations

Redox-induced interconversions of metal oxidation states typically result in multiple phase boundaries that separate chemically and structurally distinct oxides and suboxides. Directly probing such multi-interfacial reactions is challenging because of the difficulty in simultaneously resolving the multiple reaction fronts at the atomic scale. Using the example of CuO reduction in H 2 gas, a reaction pathway of CuO → monoclinic m-Cu 4 O 3 → Cu 2 O is demonstrated and identifies interfacial reaction fronts at the atomic scale, where the Cu 2 O/m-Cu 4 O 3 interface shows a diffuse-type interfacial transformation; while the lateral flow of interfacial ledges appears to control the m-Cu 4 O 3 /CuO transformation. Together with atomistic modeling, it is shown that such a multi-interface transformation results from the surface-reaction-induced formation of oxygen vacancies that diffuse into deeper atomic layers, thereby resulting in the formation of the lower oxides of Cu 2 O and m-Cu 4 O 3 , and activate the interfacial transformations. In conclusion, these results demonstrate the lively dynamics at the reaction fronts of the multiple interfaces and have substantial implications for controlling the microstructure and interphase boundaries by coupling the interplay between the surface reaction dynamics and the resulting mass transport and phase evolution in the subsurface and bulk.

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

Identifying preferential flow from soil moisture time series: Review of methodologies

Abstract Identifying and quantifying preferential flow (PF) through soil—the rapid movement of water through spatially distinct pathways in the subsurface—is vital to understanding how the hydrologic cycle responds to climate, land cover, and anthropogenic changes. In recent decades, methods have been developed that use measured soil moisture time series to identify PF. Because they allow for continuous monitoring and are relatively easy to implement, these methods have become an important tool for recognizing when, where, and under what conditions PF occurs. The methods seek to identify a pattern or quantification that indicates the occurrence of PF. Most commonly, the chosen signature is either (1) a nonsequential response to infiltrated water, in which soil moisture responses do not occur in order of shallowest to deepest, or (2) a velocity criterion, in which newly infiltrated water is detected at depth earlier than is possible by nonpreferential flow processes. Alternative signatures have also been developed that have certain advantages but are less commonly utilized. Choosing among these possible signatures requires attention to their pertinent characteristics, including susceptibility to errors, possible bias toward false negatives or false positives, reliance on subjective judgments, and possible requirements for additional types of data. We review 77 studies that have applied such methods to highlight important information for readers who want to identify PF from soil moisture data and to inform those who aim to develop new methods or improve existing ones. Core Ideas Soil moisture data can be used to identify the occurrence of preferential flow (PF) and its initiating conditions. Various data‐analysis methods to identify PF differ in susceptibility to error, bias, and subjectivity. These methods can utilize vast amounts of data from soil moisture monitoring networks to develop understanding of when, where, and under what conditions PF occurs. Newly developed methods may lead to better accuracy and reliability, and reduce the need for subjective judgments. Plain Language Summary Preferential flow through soil occurs when a large amount of water is suddenly available, as during an intense storm. This type of flow moves rapidly through the soil in distinct narrow pathways rather than moving evenly throughout the body of soil, with major consequences for groundwater resources, ecosystems, spreading of contaminants, and other vital concerns. Methods of detecting preferential flow have been developed that utilize measurements of soil water content made by sensors installed at various depths. This measurement technology has been widely implemented, many locations now having datasets years in length, and various methods have been developed for using these to identify preferential flow. The various methods are based on different features in the soil moisture records and vary in their advantages and shortcomings. In this review, we explain and evaluate these methods, highlighting important information for their implementation to identify preferential flow from soil moisture data and for efforts to develop new methods or improve existing ones.

Nimmo, John R↗
DOE OSTI2024

Simulated plant-mediated oxygen input has strong impacts on fine-scale porewater biogeochemistry and weak impacts on integrated methane fluxes in coastal wetlands

Methane (CH 4 ) emissions from wetland ecosystems are controlled by redox conditions in the soil, which are currently underrepresented in Earth system models. Plant-mediated radial oxygen loss (ROL) can increase soil O 2 availability, affect local redox conditions, and cause heterogeneous distribution of redox-sensitive chemical species at the root scale, which would affect CH 4 emissions integrated over larger scales. In this study, we used a subsurface geochemical simulator (PFLOTRAN) to quantify the effects of incorporating either spatially homogeneous ROL or more complex heterogeneous ROL on model predictions of porewater solute concentration depth profiles (dissolved organic carbon, methane, sulfate, sulfide) and column integrated CH 4 fluxes for a tidal coastal wetland. From the heterogeneous ROL simulation, we obtained 18% higher column averaged CH 4 concentration at the rooting zone but 5% lower total CH 4 flux compared to simulations of the homogeneous ROL or without ROL. This difference is because lower CH 4 concentrations occurred in the same rhizosphere volume that was directly connected with plant-mediated transport of CH 4 from the rooting zone to the atmosphere. Sensitivity analysis indicated that the impacts of heterogeneous ROL on model predictions of porewater oxygen and sulfide concentrations will be more important under conditions of higher ROL fluxes or more heterogeneous root distribution (lower root densities). Despite the small impact on predicted CH 4 emissions, the simulated ROL drastically reduced porewater concentrations of sulfide, an effective phytotoxin, indicating that incorporating ROL combined with sulfur cycling into ecosystem models could potentially improve predictions of plant productivity in coastal wetland ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

An integrated approach to derive relative permeability from capillary pressure

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability curve can be determined. Relative permeability is a numerical description of the interaction between two or more fluids and the porous media. It is a critical parameter for various tools that characterize subsurface multiphase flow systems, such as numerical simulation for carbon sequestration, oil and gas development, and groundwater contamination remediation. Therefore, it is critical to get a good statistical distribution of relative permeability in the porous media under study. Empirical formula for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility that is required to match laboratory-derived relative permeability curves. By expanding the existing methods for calculating relative permeability from capillary pressure data, it is possible to create both two and three-phase relative permeability curves. Mercury intrusion capillary pressure (MICP) data from the Morrow 'B' Sandstone coupled with interfacial tension and contact angle measurements were used to create a suite of relative permeability curves. Furthermore, these curves were then calibrated to a small sample of existing laboratory curves to elucidate common fitting parameters for the formation that were then used to create relative permeability curves from MICP data that does not have an associated laboratory-measured relative permeability curve.

58 GEOSCIENCES↗
DOE OSTI2026

In-situ sensor monitoring of multi-class gas porosity formation in laser powder bed fusion using convolutional neural network

In-situ monitoring of defect formation remains a significant challenge in the laser powder bed fusion (LPBF) process. Recent advances have enabled real-time defect detection with machine learning and in-situ sensing technologies; however, most studies focus on binary classification of keyhole pores, limiting nuanced multi-class pore differentiation and formation mechanisms. This work introduces a multi-class pore detection framework (no pore, small pores < 15 µm, and large pores > 15 µm) by leveraging photodiode sensor data alongside high-fidelity synchrotron X-ray imaging. The 15 µm threshold is selected to distinguish between two fundamentally different defect mechanisms, following the physical size-mechanism boundary established by prior high-resolution synchrotron X-ray characterization of Al6061 LPBF. Distinguishing these classes is critical because large keyhole pores are structurally detrimental, whereas small gas pores are often benign, requiring different process control strategies. Thermal emission monitoring data collected simultaneously with high-speed X-ray imaging at the Stanford Synchrotron Radiation Lightsource (SSRL), are correlated with subsurface melt pool dynamics to establish ground truth. Continuous Wavelet Transform (CWT) with optimized parameters converts the photodiode time-series signals into time–frequency images, facilitating feature extraction. Convolutional Neural Networks (CNN) are then applied for real-time multi-class pore classification in an average inference time of 1 ms per signal window. It achieves 79% accuracy and an Area Under the Receiver Operating Characteristic curve (AUC ROC) score of 0.89 with five-fold cross-validation. The results demonstrate that coupling CWT-based feature engineering with CNN architecture enables reliable multi-class pore detection in Al6061 builds using affordable in-situ sensors. This approach advances scalable and affordable quality assurance in additive manufacturing by moving beyond binary defect detection toward more nuanced classification of porosity mechanisms with in-situ sensors and machine learning.

Laser powder bed fusion, Multi-class pores, In-sit↗
DOE OSTI2026

A comparison of long-term 37 Ar and 39 Ar transport and implications for 39 Ar as a persistent indicator of underground nuclear explosions

Measurements were first made in 2016 of legacy 39 Ar gas from the decades-old Barnwell underground nuclear explosion site within the Nevada National Security Site. This has raised the question of the ubiquity of the radionuclide 39 Ar following underground nuclear explosions in general locations. This article presents a high-level evaluation of the scientific case for using 39 Ar as a long-term indicator of underground nuclear explosions. Assessment is made based on a broad set of assumptions to reasonably estimate the persistence of 39 Ar in the subsurface environment following an explosion to establish potential detectability years to decades after the originating event.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Methods to Observe Tribological Failures in Self-Mated Steel Contacts

Scuffing, a type of wear found in highly stressed or poorly lubricated contacts, is characterized by a rapid increase in friction and severe plastic deformation of the near-surface material. Scuffing has proven difficult to study because it initiates unpredictably, progresses rapidly, and typically develops within an inaccessible contact interface. Although there have been successful in-situ studies of scuffing in real-time, the transparent counter body needed for these studies changes the interactions between the surfaces and the lubricant, which affects the scuffing process in unknown ways. This paper describes the development of X-ray-compatible tribometry to study the scuffing of self-mated steels in-situ and in real-time. The method uses a crossed cylinders configuration with a thin (500 μm thick) stationary component and a small (≈200 μm) contact width to maximize X-ray interactions with atoms within the stress field generated by the contact. The resulting instrument and method are used to benchmark the scuffing response of self-mated 52,100 steel under tribologically challenging ‘oil-off’ lubrication conditions. The results demonstrate reliable scuffing in this configuration despite the relatively small contact areas and loads used. Following scuffing, gross plastic deformation was observed on both surfaces along with significant subsurface grain refinement and flow only on the stationary surface, which experienced constant contact. Interestingly, high friction initiated at specific locations of the migratory surface, which experienced intermittent contact, and then propagated across the track over time, suggesting that local conditions of the migratory surface dominated friction leading into the failure event.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Validation of Hole-Drilling Residual Stress Measurements in Workpieces of Various Thickness

A recent revision to the ASTM E837 standard for near-surface residual stress measurement by the hole-drilling method describes a new thickness-dependent stress calculation procedure applicable to “thin” and “intermediate” workpieces for which strain versus depth response depends on workpiece thickness. This new calculation procedure differs from that of the prior standard, which applies only to thick workpieces with strain versus depth response independent of thickness. Herein we assess the new calculation procedures by performing hole-drilling residual stress measurements in samples with a range of thickness. Near-surface residual stress is measured in a thick aluminum plate containing near-surface residual stress from a uniform shot peening treatment, and in samples of different thickness removed from the plate at the peened surface. A finite element (FE) model is used to assess consistency between measured residual stress across the range of sample thickness. Measured residual stress varies with sample thickness, with thinner samples exhibiting smaller near-surface compressive stress and a larger gradient of subsurface stress. These trends are consistent with both observed bending (curvature) of the removed samples and the trend in FE-calculated expected residual stress. The measured and expected residual stresses are in good agreement for samples of intermediate thickness, but the agreement decreases with sample thickness. Measured residual stress is invariant with gage circle diameter. The new thickness-dependent stress calculation procedure for hole-drilling provides meaningful improvement compared to thick-workpiece calculations.

42 ENGINEERING↗
DOE OSTI2024

Splitting tensile strength of shale cores: intact versus fractured and sealed with ureolysis-induced calcium carbonate precipitation (UICP)

Ureolysis-induced calcium carbonate precipitation (UICP) is a biomineral solution where the urease enzyme converts urea and calcium into calcium carbonate. The resulting biomineral can bridge gaps in fractured shale, reduce undesired fluid flow, limit fracture propagation, better store carbon dioxide, and potentially enhance well efficiency. The mechanical properties of shale cores were investigated using a modified Brazilian indirect tensile strength test. An investigation of intact shale using Eagle Ford and Wolfcamp cores was conducted at varying temperatures. Results show no significant difference between shale types (average tensile strength = 6.19 MPa). Eagle Ford displayed higher strength at elevated temperature, but temperature did not influence Wolfcamp. Comparatively, cores with a single, lengthwise heterogeneous fracture were sealed with UICP and further tested for tensile strength. UICP was delivered via a flow-through method which injected 20–30 sequential patterns of ureolytic microorganisms and UICP-promoting fluids into the fracture until permeability reduced by three orders of magnitude or with an immersion method which placed cores treated with guar gum and UICP-promoting fluids into a batch reactor, demonstrating that guar gum is a suitable inclusion and may reduce the number of flow-through injections required. Tensile results for both delivery methods were variable (0.15–8 MPa), and in some cores the biomineralized fracture split apart, possibly due to insufficient sealing and/or heterogeneity in the composite UICP-shale cores. Notably in other cores the biomineralized fracture remained intact, demonstrating more cohesion than the surrounding shale, indicating that UICP may produce a strong seal for subsurface application.

58 GEOSCIENCES↗
DOE OSTI2026

From bulk to surface: Structure and dynamics of amorphous alumina from deep potential molecular dynamics

Understanding the atomic-scale structure and dynamics of amorphous oxide surfaces is essential for interpreting their chemical reactivity, mechanical stability, and interfacial behavior, yet direct experimental characterization remains challenging. We employ Deep Potential (DP) molecular dynamics to generate large-scale, ab initio -quality models of amorphous Al 2 O 3 bulk glasses and melt-quenched free surfaces, enabling a quantitative analysis of both structure and relaxation dynamics with statistical confidence inaccessible to direct ab initio simulation. The trained DP model reproduces experimental liquid and glass structure, captures the cooling-rate dependence of the bulk glass transition, and corrects systematic biases in the polyhedral populations predicted by widely used classical force fields. At the free surface, mass density recovers to bulk values over ~10 Å, while local coordination requires a slightly wider subsurface region to fully converge. The outermost layer is oxygen-enriched, exhibits altered polyhedral connectivity with contracted Al–O bonds, and hosts a broad population of under-coordinated motifs (notably AlO 3 and OAl 2 ) whose abundances are governed by glass stability. These under-coordinated surface motifs exhibit distinct vibrational signatures and occur as locally paired Lewis acid and Brønsted base sites consistent with bond-valence compensation, yet remain spatially dispersed rather than aggregating into extended clusters. Despite this pronounced structural heterogeneity, surface relaxation and the glass-transition temperature remain comparable to their bulk counterparts, suggesting that the disordered surface is kinetically stable once formed. Together, these results establish a molecular-level picture of amorphous alumina surfaces and demonstrate the capability of machine-learned potentials to resolve structure–property relationships in disordered oxide interfaces.

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

Gas pore correlations in laser powder bed fusion of Al6061

Additive manufacturing (AM) of metal materials based on powder bed fusion technology is widely used now in many industries. A known limitation of this type of manufacturing is the formation of gas pores in the bulk material. Here we present a combined X-ray imaging and mid-infrared pyrometry study of pore formation in side-by-side tracks of Al6061 for different processing conditions using both in-situ and post-processing analysis. By carefully quantifying the distributions and correlations of pore positions, we show that an existing pore in one track often catalyzes the formation of another pore in an adjacent track. In a raster scan strategy commonly used to construct bulk material, this phenomenon has the result of forming subsurface perforations, or lines of pores transverse to the scanning direction in a rastered patch. If controlled, this effect can be eliminated to improve the yield strength of the build, or exploited to create programmable failures for specific purposes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Real-time tracking and analysis of gas bubble dynamics in laser powder bed fusion using in-situ X-ray characterization and machine learning

Porosity defects remain a significant challenge in the laser powder bed fusion (LPBF) process, adversely affecting the mechanical properties and reliability of additively manufactured components. Here, this study investigates the real-time formation and trajectory of gas bubbles during LPBF of Al6061 alloy using advanced in-situ X-ray characterization and machine learning. The unsupervised Gaussian mixture model and particle tracking algorithm developed are able to precisely track and quantify the properties of gas bubbles and keyhole pores. Our analysis identified five distinct types of gas bubble formation and movement patterns, emphasizing the diverse origins and behaviors of these defects. It enables precise quantification of trajectories, velocities, and morphological changes of gas bubbles, offering a granular view of the subsurface dynamics within the melt pool. Additionally, we explored keyhole-induced pore dynamics, revealing the critical role of keyhole oscillation and collapse for the formation of both large and small gas pores. It defines four different regions of gas bubble movement within the melt pool, providing a clearer understanding of how local fluid dynamics affect pore behavior. The results underscore the importance of integrating in-situ experimental observation and automated machine learning to develop a more robust predictive model for defect formation in LPBF.

In-situ X-ray imaging↗
DOE OSTI2024

Numerical investigation of multiphase flow through self-affine rough fractures

Multiphase flow through fractures has great significance in subsurface energy recovery and gas storage applications. Different fracture and flow properties affect flow through a fracture which is difficult to control in laboratory experiments. Here, we perform lattice Boltzmann simulations in an ensemble of synthetically generated fractures. Drainage simulations are performed at different capillary numbers, wettability, and viscosity ratios. We track the invading front and quantify breakthrough saturations and show that roughness and wettability have a strong effect on fluid invasion through a complex fracture. Invading a more viscous fluid results in more stable displacement regardless of the capillary number while at very low capillary numbers, fluid migration is dependent on the inherent structure of the fracture. We develop a fluid displacement phase diagram in a single rough fracture and compare our results from that in the literature. Finally, we extend the phase diagrams across multiple fractures and demonstrate the importance of natural fracture features of roughness and wettability in identifying stable versus unstable displacement regimes during multiphase flow through rough fractures. Our work presents an end-to-end numerical pathway for testing on experimental data and expanding numerical data sets for testing combinations of different physical phenomenon and make valuable predictions on fluid flow through rough fractures.

02 PETROLEUM↗
DOE OSTI2025

The rigorous upscaling of advection-dominated transport in heterogeneous porous media via the Method of Finite Averages

Systems involving advection-dominated transport through heterogeneous porous and fractured media are ubiquitous in subsurface engineering applications. However, upscaling such systems continues to challenge rigorous modeling efforts, particularly when advection is stronger than diffusion at fine spatial scales (i.e., when the Péclet number is greater than one at length scales that characterize a system’s unit-cells, representative elementary volumes, or averaging regions). Here, in this work, we propose and validate a strategy for extending the Method of Finite Averages (MoFA), a rigorous upscaling methodology for heterogeneous porous media, to upscale transport systems experiencing stronger advection than diffusion at fine scales (i.e., fine-scale Péclet numbers greater than one). We detail the strategy, the physical conditions under which it can be applied while retaining a priori modeling error guarantees, and implement the strategy to obtain a MoFA model for advective-diffusive transport that accommodates advective physics at fine spatial scales. We then perform two numerical experiments considering systems with system-scale Péclet numbers of 300 and 1000 — which correspond to fine-scale Péclet numbers of 30 and 100, respectively — to verify that the error guarantees are met under the strategy. After, we conduct a numerical study to demonstrate the strategy’s advantages over the original MoFA methodology. The results suggest that rigorously-upscaled transport models for heterogeneous porous media experiencing advective physics at finer spatial scales can be derived through MoFA and resolved orders of magnitude faster than their pore-scale counterparts. The results also suggest that the presented strategy is limited to modeling shallow concentration gradients when there are large differences between the time scales related to advection and a system’s temporally-varying boundary conditions. This limitation hinders the strategy’s practicality in modeling more advective systems, and as such, opportunity exists for developing additional strategies that accommodate rapidly-varying boundary conditions — and consequentially, steeper concentration gradients — while modeling advective systems with MoFA.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Underground hydrogen storage resource assessment for the Cook Inlet, Alaska

Underground hydrogen storage will be essential to a hydrogen economy. This work focuses on the challenge of identifying and screening candidate storage systems, given the unique behavior of hydrogen in the subsurface. Here, we describe a resource assessment methodology and apply it to Alaska’s Cook Inlet region. Alaska provides an interesting case study because of its abundant renewable energy resources, relatively low energy demand, and isolated electrical grid. The assessment framework considers each site’s ability to (1) store a specific volume, (2) physically-contain the stored gas, and (3) limit biogeochemical activity. We estimate that reservoirs in the Cook Inlet area could theoretically store a total of 286 TWh (or 8.6 million tonnes [Mt]) in hydrogen working gas in 92 pools. This is likely sufficient to meet both local hydrogen demand and support an array of exportable products. We further identify seven pools that may be especially well-suited for hydrogen storage. Broadly, this work demonstrates a framework for regional resource assessments. On a finer scale, this work supports an early demonstration of porous–media hydrogen storage in the United States.

08 HYDROGEN↗
DOE OSTI2024

Wellbore cement alteration and roles of CO 2 and shale during underground hydrogen storage

To mitigate climate change and adopt renewable energy, energy storage is crucial and can be done in the form of hydrogen gas (H 2 ). Subsurface geologic reservoirs are positioned to store H 2 on the largest scales for the longest terms of all potential options. However, H 2 injection may boost reactions that consume hydrogen, generate undesired gases, and alter pore structures of geomedia. To explore the extent of H 2 -associated biotic reactions at a near wellbore location, four experiments were conducted under underground storage conditions with wellbore cement cores and, in most instances, shale samples submerged in synthetic formation brine. Post-reaction gas, aqueous, and solid phase samples were analyzed using olfactory screening and, later, gas chromatography (GC-MS), inductively coupled plasma optical emission spectroscopy (ICP-OES), scanning electron microscopy (SEM), and synchrotron micro-scale x-ray fluorescence (μ-XRF). Within a period of 16 weeks, hydrogen sulfide (H 2 S) was generated in systems containing both H 2 and shale. XRF mapping identified a zone enriched in iron(II) and reduced sulfur along the rim of cement cross sections that was largely associated with CO 2 -induced cement carbonation. Shale did not show noticeable alteration, but there is evidence it contributed to the initial inoculation of the system and provided nutrients for microbes via water-rock interactions. Here, this study considers both rock formations and wellbore cement not previously evaluated concurrently. Findings support understanding and modeling of H 2 -associated biogeochemical reactions during underground hydrogen storage.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Unraveling Hydrogen Induced Geochemical Reaction Mechanisms through Coupled Geochemical Modeling and Machine Learning

Underground hydrogen storage (UHS) provides a promising large-scale, long-term energy storage solution. A reasonable recovery of stored hydrogen is critical for a successful storage scheme. However, in subsurface reservoirs hydrogen is subject to active geochemical reactions that might result in hydrogen loss. In this study, we implemented a geochemical modeling approach coupled with an unsupervised machine learning technique called non-negative matrix factorization (NMF) to unravel the complex brine-rock-H 2 geochemical processes responsible for hydrogen losses, with particular focus on sulfate reduction reactions. NMF is applied to modeled mineral evolution and fluid component profiles to retrieve profiles that can be interpreted to more easily assess competing processes. NMF decouples simulated competing equilibrium reactions. This facilitates separation of overlapping reaction profiles from redox processes, dissolution fronts, and secondary precipitation while considering the effects of simulation parameters such as salinity, temperature, and total H 2 pressure. NMF successfully discriminates these competing effects in nonlinear ways, allowing robust interpretation. In addition, NMF reveals subtle coupled mineral associations and reaction fronts that are invisible to conventional model analysis. This integrated approach strengthens the conceptual understanding of complex nonlinear hydrogen-brine-rock interactions and advances geochemical research on UHS systems to resolve complexities in modeled geochemical systems without the need for direct experiments or prior knowledge. Furthermore, this study highlights the efficacy of combining geochemical modeling with machine learning techniques to enhance the interpretability of the intricate geochemical simulation output through deciphering the overlapping reaction path that cannot be achieved only using conventional analysis of geochemical models alone.

08 HYDROGEN↗
DOE OSTI2024

A new perspective on density and strength loss profiles at the surface of thermally oxidized nuclear graphite

Oxidation of graphite components could influence their designed life in a high-temperature nuclear reactor. The oxidized regions could potentially lower the allowed stress capacity. The American Society of Mechanical Engineers rules for the design and construction of graphite-moderated reactors recommend that subsurface regions that might become excessively damaged by oxidation during reactor operation be identified and excluded from geometry and stress calculations. Identification of oxidation-affected regions is possible, in principle, through complex modeling exercises of reactor behavior during hypothetical accident scenarios coupled with graphite oxidation models, but this procedure may not have the precision needed for informed decisions. Here, this paper proposes an alternate method, based on interpretation of a series of well-designed oxidation experiments, which could augment the designer's tools. The procedure is illustrated by data on oxidation by air of several graphite grades (NBG-18, PCEA, IG-110, R4-650) that are corroborated with independent literature information, when available. The Wichner model for graphite oxidation used for this analysis provides conservative results that could be quickly implemented in the design process.

36 MATERIALS SCIENCE↗