Search NASASearch

SEARCH · Search NASA

Results for “sodium source saving”

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.

364 records · Page 18

SAM Code Enhancements for Fission Product Tracking of Noble Gases and Metals in MSRs

This report documents fiscal year 2026 enhancements to the System Analysis Module (SAM) for modeling fission product transport in liquid-fueled molten salt reactors (MSRs). The work advances three principal areas: noble gas transport, noble metal deposition, and user interface improvements. The noble gas transport capability integrates drift-flux gas transport, Henry’s law two-film interphase mass transfer with pressure-based nucleation suppression, Knudsen-regime pore diffusion into porous graphite with a conjugate salt-graphite interface constraint, built-in material properties, five Sherwood-number mass transfer correlations including three derived from high-fidelity NekRS simulations, and xenon-135 reactivity feedback through SAM’s point-kinetics model. This work also presents a comprehensive verification test suite, including new analytically verified cases for pressure-dependent onset of interphase gas transfer in a stagnant vertical pipe, a postulated FLiBe-graphite Xe extraction permeator, a gravity riser with a fission-product source, and a descending pipe with gas redissolution driven by hydrostatic pressure. A machine learning framework for bubble rise velocity prediction in molten salt systems is developed and benchmarked on molten-salt and diverse aqueous bubble datasets. The best-performing fine-tuned transfer-learning networks achieve an 82% reduction in RMSE relative to the Clift correlation, and is implemented directly in SAM. The noble metal transport capability is developed, including a liquid-wall deposition model and a gas-surface flotation mechanism that transfers insoluble particles entrained by sparging gas to wetted structures. Verification tests and demonstration cases cover the surface deposition, flotation efflux, and flotation shedding. Finally, a new [SpeciesTransport] input structure replaces positional global vectors with selfcontained, order-independent, named species blocks, simplifies the specification of multiphase species and decay chains, and remains fully compatible with existing SAM input files. Together, these developments improve the physical fidelity, verification basis, and usability of SAM for system-level analyses of fissionproduct behavior in MSRs.

Mui, Travis (ORCID:0000000303736470)

PACT Center: Perovskite PV Accelerator for Commercializing Technologies (Final Technical Report)

The Perovskite PV Accelerator for Commercializing Technologies (PACT) center was established in July 2021 as a national resource to accelerate the commercialization of perovskite photovoltaic (PV) technology in the United States. Since its inception, PACT has been led by Sandia National Laboratories (Sandia) in partnership with the National Laboratory of the Rockies (NLR), formerly known as NREL. From FY20-FY23, Los Alamos National Laboratory (LANL), CFV Labs, Black & Veatch (B&V), and the Electric Power Research Institute (EPRI) were part of the project team. LANL brought expertise in perovskite PV device designs and processing, CFV Labs (now GroundWork Renewables) provided initial indoor and outdoor measurement hardware technology, B&V led the initial effort on perovskite PV bankability, and EPRI worked on reviewing testing standards, identifying commercialization gaps, and helping to run PACT’s Industry Advisory Board, a group including representatives from commercial testing labs, independent engineering firms, insurance companies, state regulators, and electric utilities. To source perovskite PV module samples, PACT contracted with the University of North Carolina (UNC), the University of Toledo, the University of Washington, and SLAC/Stanford University to provide a steady stream of research-grade perovskite mini modules, enabling protocol development in advance of commercial module availability. The project period ran from July 1, 2021, through December 31, 2025, including a No Cost Extension. Starting in FY25, the project was continued as a Core Capability in the Lab Call portfolio and continues at a reduced budget with only Sandia and NLR as funded recipients. Notably, starting in FY25 PACT expanded its scope beyond MHP modules to accept all emerging PV mini module technologies for testing, including organic PV (OPV) and all-thin-film tandems, with the aim of supporting commercialization across the broader emerging PV ecosystem. With this change in scope the program was renamed the PV Accelerator for Commercializing Technologies, dropping perovskite from the name.

14 SOLAR ENERGY

Divergent carbon use efficiency-growth rate tradeoff in popular biological growth models

Carbon use efficiency (CUE) is an important trait emerging from processes regulating biological growth. CUE can be computed either based on the growth of structural biomass or total biomass divided by substrate uptake rate. Nonequilibrium thermodynamics and observations suggest that, for an exponentially growing population of cells, structural biomass CUE should first increase, then peak, and finally decrease with specific growth rate; meanwhile, total biomass CUE increases asymptotically with specific growth rate. We compared predictions from six popular models that are often used for plant and microbial growth in existing ecosystem models. We found that, for an exponentially growing population of biological cells, (1) the source-driven Pirt and Compromise models predict that structural biomass CUE increase asymptotically with growth rate; (2) the apparent sink-driven modified Droop model predicts that structural biomass CUE decreases with growth rate; and (3) the sink-driven variable internal storage model and two dynamic energy budget models predict that structural biomass CUE first increases, then peaks, and finally decreases with growth rate. Moreover, the modified Droop model predicts that total biomass CUE is constant with growth rate, while all other five models predict that total biomass CUE increases with growth rate asymptotically. For non-exponential biological growth, we show that there is no static relationship between total biomass CUE or structural biomass CUE with respect to either growth rate or temperature. Therefore, we contend that biological growth models should explicitly represent interactions between substrate acquisition, substate transformation, and maintenance respiration to better capture observed CUE dynamics, and the sink-driven model should be preferred for general ecosystem biogeochemistry modeling.

Tang, Jinyun [Lawrence Berkeley National Laborator

Roadmap and Benchmarking: Privacy in Federated Load Forecasting

Data-driven techniques for energy demand forecasting continue to emerge with promising impacts on distribution grid planning. However, the development of robust and generalizable machine learning models requires that representative high quality training data are available. Distributed energy resources have begun to embed intelligence, gathering large amounts of data on customer demand, behavior, and household devices that are connected to the grid. Though utilities aggregate meter-level demand data for load shaping, demand response, outage management, reliability planning, and billing applications, there lies an inherent privacy concern in sharing consumption data that may identify individual consumer behavioral patterns. Hence, while sharing the data is crucial, the private sensitive customer data must be safeguarded from being exposed or manipulated. In this study, we propose a roadmap for implementing a based privacy preserving framework to support the advancement of data-driven analytics in data-sensitive distributed energy resources environments. The roadmap incorporates federated learning–a distributed training framework, differential privacy–a statistical framework that provides guarantees to safeguard the leakage of sensitive data, secure multiparty computation and homomorphic encryption– techniques for encrypting model gradients and applying secure aggregation on the server. Moreover, we perform baseline experiments on the federated short-term load forecasting (STLF) task using open-source residential load profile datasets, offering insights into the challenges of integrating differential privacy into federated learning.

Abebe, Waqwoya [Oak Ridge National Laboratory (ORN

Chronology of ungrouped achondrites Northwest Africa 11119 and 7325/8486: implications for early evolution of a heterogeneous crust on a differentiated planetesimal

The chronology of ungrouped achondrites provides key insights into the timeline of igneous activity in the early Solar System. Two ungrouped achondrites – Northwest Africa (NWA) 11119 (andesite-dacite) and NWA 8486 and its pair NWA 7325 (olivine gabbro) – are the focus of this study, where their chronologies are investigated using multiple high-resolution techniques. Here we report the lead-lead (207Pb-206Pb) and manganese-chromium (53Mn-53Cr) systematics of NWA 11119, as well as the 207Pb-206Pb systematics of NWA 8486 alongside aluminum-magnesium (26Al-26Mg) systematics for NWA 7325. The U-corrected 207Pb-206Pb ages of NWA 11119 and the combined ages of NWA 7325/8486 are 4566.4 ± 0.8 Ma and 4563.8 ± 1.9 Ma, respectively. Additionally, we report the 53Mn-53Cr age of NWA 11119 to be 4564.4 ± 2.5 Ma and the 26Al-26Mg age of NWA 7325/8486 to be 4563.1 ± 0.3 Ma. The formation of NWA 11119 requires partial melting of a (likely chondritic) source reservoir leading to eruption of Si rich, alkali depleted magmas, while NWA 7325/8486 likely formed from a chemically fractionated reservoir with superchondritic Al/Mg. The clear geochemical and isotopic similarities of these achondrites, combined with the chronology reported here, is suggestive of formation of these two ungrouped achondrites on a common parent body which likely formed in the inner Solar System and experienced early differentiation under reducing conditions. If these achondrites did share a parent body, it would suggest that primary asteroids in the early Solar System commonly produced mineralogically and geochemically heterogeneous crusts. While mineralogical and geochemical heterogeneity is known to exist in the (mostly mafic) crust of asteroid Vesta, our findings show that even more significant crustal heterogeneity (representing felsic and mafic compositions) may exist on other asteroids.

Cartwright, Julia [University of Leicester, UK]

Qualifying LaBr3:Ce+Sr detector performance for the Mu2e experiment at Fermilab Using the ELBE Accelerator

A LaBr3:Ce+Sr detector will be used to measure the stopped muon captures at the Mu2e experiment at Fermilab. It has been benchmarked in a test beam experiment performed at the ELBE electron accelerator located at the Helmholtz-Zentrum Dresden-Rossendorf, Germany. ELBE’s pulsed bremsstrahlung beam line was set to deliver an average γ -ray energy of between 4–5 MeV. The detector response was mapped to match Mu2e beam conditions, including rates up to 1 Mcps, energy flux, and time structure. A radioactive calibration source was used to mimic the characteristic 1808.7 keV γ -ray, emitted during the atomic muon nuclear capture in the Mu2e aluminum stopping target. The detector energy resolution was measured as a function of the average energy flux: up to 1 TeV/s for 0.34 s, the steady-operation beam-on time and up to 4 TeV/s for 5 ms to get a conservative estimate of the effect of high intensity Mu2e beam fluctuations. The PMT gain variation as a function of the beam spill length and average intensity has been parametrized and corrected for. When a PMT gain correction corresponding to the average beam-spill intensity is applied, the residual effect of beam intensity fluctuations around the average degrades the energy resolution, σ E γ /E γ , at 1808.7 keV from 0.66% to 0.83%.

Huang, Shihua [Purdue U., West Lafayette]

Selective Binding and Light-Driven Release of Fluorous PF 6 – and Radioactive 99 TcO 4 – Anions for All-to-Nothing Liquid–Liquid Extraction

The removal of anions from aqueous media using molecular receptors in liquid–liquid extraction is a long-standing strategy to clean up contaminated water sources. Therein, high selectivity is needed to remove specific ions from mixtures of other ions, and high affinity provides extractions at low concentrations. However, the high affinity creates a conundrum by impeding the release of the ions in any stripping steps needed for further processing. To circumvent this problem, light-responsive receptors have been proposed as candidates for turning off the binding, but they are currently untested in liquid–liquid extraction. We tested the feasibility of light-driven release using a cyanostar macrocycle. We demonstrate the selective extraction of PF$^{–}_{6}$ anions over large excesses of competing anions (Cl – , NO$^{–}_{3}$, SO$^{2–}_{4}$) followed by photodriven release for quantitative isolation of the target. Release relies on photoisomerization of the macrocycle’s five stilbenes generating distorted isomers to turn off binding. With modest reversibility, only a single-shot release was demonstrated, akin to photodriven uncaging. These methods were extended to the capture and photodriven release of ReO$^{–}_{4}$ and radioactive 99 TcO$^{–}_{4}$ anions at ∼90% efficiency. Extraction was demonstrated down to the highly dilute 4 ppb levels of the 99 TcO$^{–}_{4}$ anion. Furthermore, this proof-of-concept demonstration verifies the use of a large change in affinity for the all-to-nothing capture and release of target anions between liquid phases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Working Fluid Characterization and Performance Assessment of Subcritical Organic Rankine Cycles Based on the Lee–Kesler Approach for Energy Recovery

Here, a generalized model using the Lee–Kesler approach based on the corresponding states principle is developed to assess the performance of subcritical Organic Rankine Cycles operating with different working fluids. Each fluid is characterized by five parameters: the acentric factor, critical temperature, critical pressure, molar mass, and the ideal-gas ratio of specific heats at the critical temperature. The model was developed using the compressibility factor modified version of the Benedict–Webb–Rubin equation proposed by Lee and Kesler and the enthalpy and entropy functions to calculate thermodynamic state properties. The model was validated by comparing the results calculated with the model and working fluid thermodynamic properties obtained with the CoolProp database. This comparison was conducted for 91 working fluids, obtaining a relative error below 5% for 88 out of the 91 fluids (∼97%). A generalized parametric study was conducted to determine the influence of the pinch point and each fluid parameter on the performance of Organic Rankine Cycle (ORC) systems. It was found that efficiency increases with critical temperature, ideal-gas ratio of specific heats at the critical temperature, and acentric factor, reaching up to 13%. The developed model enables the evaluation of ORC system performance for existing working fluids. It also allows the formulation and evaluation of new fluids to enhance the performance of the ORC while retrieving energy from any kind of source; and likewise, the methodology can be applied to other power generation cycles.

24 POWER TRANSMISSION AND DISTRIBUTION

Development of large-scale, 3D Printed high temperature ceramic material

Through this collaborative effort, a new 3D printing platform for ceramics called laser-induced slip casting (LIS) was explored and developed for improving the processing and manufacture of silicon carbide (SiC), a high temperature ceramic material. This method prints layers of ceramic slips/slurries with subsequent selective-laser heating to dry each layer to build a 3D structure. The completed and dried part is then sintered. This manufacturing process is inherently lower cost than alternative ceramic printing methods such as binder jet technology or stereolithography when it comes to the feedstock material but is more expensive than robocasting or direct ink writing. However, it has the potential to make more controlled parts with less defects compared to robocasting. The largest cost is the heating source for the printer. With this technique, there is potential to make large ceramic parts in a near-net shape. Further, there is no known commercial manufacturing of 3D printed ceramics that creates a large range of different high-density ceramics at large scale. The goals and outcomes for the development of the new printing technology were: 1) assessing the technology with alumina by characterizing coupons, 2) printing and sintering of silicon carbide (SiC), and 3) characterization and properties testing of materials printed and a scale-up of SiC part(s). Through this collaborative project, it was sought to provide the best solution to achieving highly dense and near-net shaped ceramic parts at large scale and reduced cost. The goal was to produce high density sintered materials for applications in defense, energy generating systems, armor, and wear parts using 3D printing methods. Variables including powder particle sizes, dispersant molecular weights (MWs), binders, printing parameters, and post processing were explored as well as the characterization of the physical properties (add specifics here).

36 MATERIALS SCIENCE

The brighter-fatter effect in the JWST MIRI Si:As IBC detectors: I. Observations, impact on science, and modeling

Context.The Mid-Infrared Instrument (MIRI) on board theJames WebbSpace Telescope (JWST) uses three Si:As impurity band conduction (IBC) detector arrays. The output voltage level of each MIRI detector pixel is digitally recorded by sampling up the ramp. For uniform or low-contrast illumination, the pixel ramps become nonlinear in a predictable way, but in areas of high contrast, the nonlinearity curve becomes much more complex. The origin of the effect is poorly understood and currently not calibrated out of the data. Aims.We provide observational evidence of the brighter-fatter effect (BFE) in MIRI conventional and high-contrast coronagraphic imaging, low-resolution spectroscopy, and medium-resolution spectroscopy data, and we investigate the physical mechanism that gives rise to the effect on the MIRI detector pixel raw voltage integration ramps. Methods.We used public data from the JWST/MIRI commissioning and Cycle 1 phase. We also developed a numerical electrostatic model of the MIRI detectors using a modified version of the publicPoisson_CCDcode. Results.We find that the physical mechanism behind the BFE manifesting in MIRI data is fundamentally different to that of charge-coupled devices and photodiode arrays such as the Hawaii-XRG near-infrared detectors used by the NIRISS, NIRCam, and NIRSpec instruments on board JWST. Observationally, the BFE makes the JWST MIRI data yield 10–25% larger point sources and spectral line profiles as a function of the relative level of de-biasing of neighboring detector pixels. This broadening impacts the MIRI absolute flux calibration, time-series observations of faint companions, and point spread function modeling and subtraction. We also find that the intra-pixel 2D profile of the shrinking Si:As IBC detector depletion region directly impacts the accuracy of the pixel ramp nonlinearity calibration model.

Astronomy & Astrophysics

From Pollutant Removal to Renewable Energy: MoS2-Enhanced P25-Graphene Photocatalysts for Malathion Degradation and H2 Evolution

The widespread presence of pesticides—especially malathion—in aquatic environments presents a major obstacle to conventional remediation strategies, while the ongoing global energy crisis underscores the urgency of developing renewable energy sources such as hydrogen. In this context, photocatalytic water splitting emerges as a promising approach, though its practical application remains limited by poor charge carrier dynamics and insufficient visible-light utilization. Herein, we report the design and evaluation of a series of TiO2-based ternary nanocomposites comprising commercial P25 TiO2, reduced graphene oxide (rGO), and molybdenum disulfide (MoS2), with MoS2 loadings ranging from 1% to 10% by weight. The photocatalysts were fabricated via a two-step method: hydrothermal integration of rGO into P25 followed by solution-phase self-assembly of exfoliated MoS2 nanosheets. The composites were systematically characterized using X-ray diffraction (XRD), Raman spectroscopy, transmission electron microscopy (TEM), UV-Vis diffuse reflectance spectroscopy (DRS), and photoluminescence (PL) spectroscopy. Photocatalytic activity was assessed through two key applications: the degradation of malathion (20 mg/L) under simulated solar irradiation and hydrogen evolution from water in the presence of sacrificial agents. Quantification was performed using UV-Vis spectroscopy, gas chromatography–mass spectrometry (GC-MS), and thermal conductivity detection (GC-TCD). Results showed that the integration of rGO significantly enhanced surface area and charge mobility, while MoS2 served as an effective co-catalyst, promoting interfacial charge separation and acting as an active site for hydrogen evolution. Nearly complete malathion degradation (~100%) was achieved within two hours, and hydrogen production reached up to 6000 µmol g−1 h−1 under optimal MoS2 loading. Notably, photocatalytic performance declined with higher MoS2 content due to recombination effects. Overall, this work demonstrates the synergistic enhancement provided by rGO and MoS2 in a stable P25-based system and underscores the viability of such ternary nanocomposites for addressing both environmental remediation and sustainable energy conversion challenges.

Chemistry

RADAI: A Large-Scale Realistic Dataset for Radiation Detection Algorithm Development

Open, realistic datasets are essential for developing and benchmarking radiation detection algorithms, yet they remain scarce. The Radiological Anomaly Detection and Identification (RADAI) project was develop to create datasets that meet the training and testing needs for sophisticated radiation detection algorithms. The RADAI dataset is a large-scale synthetic resource that integrates high-fidelity Monte Carlo simulations with realistic urban scenarios to capture both background variability and source signatures. RADAI models construction-material NORM, people and vehicles, urban clutter, and dynamic environmental effects such as cosmic-ray and rain-induced transients, and they provide list-mode detector data with motion and response modeling suitable for algorithm training and evaluation. The RADAI project resulted in three publicly-released complementary datasets together with an online scoring portal for standardized performance assessment and an open software toolkit that supports data access, augmentation, model development, and evaluation. These resources enable reproducible comparisons across methods and promote rigorous studies at the scale required by contemporary machine learning. By grounding algorithm development in realistic, well-documented conditions, RADAI supports progress toward more robust detection, identification, and localization in complex urban environments.

Ghawaly, James M. [Division of Computer Science an

Calculations of the Moon's Thermal History at Different Concentrations of Radioactive Elements, Taking Into Account Differentiation on Melting

Calculations of the thermal history of the moon were carried out by solving the thermal conductivity equation for the case in which the heat sources are the long-lived radioactive elements Th, U, and K 40 . The concentrations of these elements were adjusted to give four variations of the heat flow: 1.35 x 10 -8 and 0.91 x 10 -6 cal cm -2 s -1 (I and I**), 0.61 x 10 -6 (variant II, the terrestrial mixture of Lyubimova and Starkova), and 0.236 x 10 -6 (variant III, the chondrite model of Urey and MacDonald). In the same calculations, we considered layering of the differentiated material with transport to the surface of the radioactive elements after the temperature of the layer rose to 200 K above the melting temperature, which is considered in five variants that differ in the amount of transported radioactive elements: 100 percent ( n = 1), 80 percent n = 0.8), 40 percent, 20 percent ( n = 0.6, 0.4, 0.2). During fusion the heat capacity and heat conductivity were changed. We considered two variants of an initially cold (273 K) and hot (900 K) Moon. Calculations show that the interior of the Moon was heated to melting during the first 0.7 to 2.3 x 10 9 years. The maximum fusion involved practically the entire Moon to a distance from 15 to 45 km beneath the surface, and started 3.5 to 4.0 x 10 9 years ago (I, I**), or 2.5 to 3.0 x 10 9 years ago (II, III) and continued for 1 to 2 x 10 9 years. Today the Moon is cooling. The current thickness of the solid crust is from 150 to 200 km and the heat flow exceeds the stationary value 1.5 fold. Apparently the most realistic variant is II (terrestrial mixture) for an initially hot Moon, and gives, regardless of the moderate concentration of radioactive elements, a heat flow of 0.9 to 0.95 x 10 -6 cal cm -2 s -1 , which agrees with the radioastronomical measurements of Troitsky and Krotikov and with the Apollo 15 data.

O I Ornatskaya

Impact of Limited Degree of Freedom Drag Coefficients on a Floating Offshore Wind Turbine Simulation

The worldwide effort to design and commission floating offshore wind turbines (FOWT) is motivating the need for reliable numerical models that adequately represent their physical behavior under realistic sea states. However, properly representing the hydrodynamic quadratic damping for FOWT remains uncertain, because of its dependency on the choice of drag coefficients (dimensionless or not). It is hypothesized that the limited degree of freedom (DoF) drag coefficient formulation that uses only translational drag coefficients causes mischaracterization of the rotational DoF drag, leading to underestimation of FOWT global loads, such as tower base fore-aft shear. To address these hydrodynamic modeling uncertainties, different quadratic drag models implemented in the open-source mid-fidelity simulation tool, OpenFAST, were investigated and compared with the experimental data from the Offshore Code Comparison Collaboration, Continued, with Correlation (OC5) project. The tower base fore-aft shear and up-wave mooring line tension were compared under an irregular wave loading condition to demonstrate the effects of the different damping models. Two types of hydrodynamic quadratic drag formulations were considered: (1) member-based dimensionless drag coefficients applied only at the translational DoF (namely limited-DoF drag model) and (2) quadratic drag matrix model (in dimensional form). Based on the results, the former consistently underestimated the 95th percentile peak loads and spectral responses when compared to the OC5 experimental data. In contrast, the drag matrix models reduced errors in estimates of the tower base shear peak load by 7–10% compared to the limited-DoF drag model. The underestimation in the tower base fore-aft shear was thus inferred be related to mischaracterization of the rotational pitch drag and the heave motion/drag by the limited-DoF model.

17 WIND ENERGY

Utilities Perspective on Protection Challenges with High IBR Penetration

Utilities have seen rapid increase in Solar, Wind, and Battery Energy Storage Resources that interconnect to their electric system through Inverters. Inverter-Based Resources (IBRs) have fault current characteristics that are unlike the fault current response of traditional rotating-machine-based generators, which is well known and repeatable. IBR’s non-traditional fault current behavior is due to the IBR control scheme, which is configured to provide a clean AC output but also protect the inverter’s sensitive power electronics devices from damage, one source of which is overcurrent. This results in low fault current magnitude, low or no negative sequence current injection, the variability of sequence component currents, the variability of voltage with respect to current angles, and the lack of inertia. The control scheme also results in a fault current response that can vary between manufactures and between models of the same manufacturer. High penetration of IBRs can adversely affect the protection schemes applied in areas with high penetration of IBRs. With the proliferation of IBRs, utilities are finding out that conventional protection schemes are not adequately equipped to protect the electric systems. This is mainly because the existing protection elements and practices have been designed based on the fault current response of conventional rotating machines. In several cases, the available literature does not provide any clear solution for the issues when the protection scheme does not operate properly near IBRs. This presentation identifies various protection challenges due to IBRs that industry is facing, from the utility perspective. Instead of facing on one issue, we are looking broadly on all the challenges that system protection has experienced with high penetration of IBRs. Based on the IBR response from various utilities during real fault events and gathering perspective from different utility SMEs via questionnaire, the presentation summarizes on gathered data and internal experiences.

24 POWER TRANSMISSION AND DISTRIBUTION

SENTRA: A Modular Computational Graph Framework for Critical Mineral and Materials Supply Chains: Part I: Network Construction Latent-Quantity Estimation, and Temporal Graph Forecasting

Global supply chains for critical minerals and materials are complex, evolving networks of countries, products, production stages, and trade relationships. Existing analytical approaches are limited by fragmented data and static network representations that do not capture the dynamic production dependencies linking raw materials, intermediate products, and final goods across multiple countries. Trade and production statistics provide only a partial view of domestic production, inventories, and material flows, making it difficult to identify indirect sourcing pathways, hidden dependencies, and embedded foreign exposures. This paper introduces the Supply Chain Exposure Network Tracking and Risk Assessment (SENTRA) framework, a modular graph-based computational framework for constructing, analyzing, and forecasting dynamic supply chain networks. As the first paper in a three-part methodological series, it establishes the computational foundation of SENTRA by constructing a temporal attributed multi-relational graph whose nodes represent product–country pairs and whose edges encode observed trade and within-country value-chain relationships. Statistical estimation and constrained optimization recover latent production, final demand, and product input dependency coefficients while enforcing economic accounting constraints. Graph-derived exposure measures quantify direct, transshipment, value-chain, and multi-hop supply chain dependencies independently of the forecasting model. A temporal graph forecasting architecture based on a relational graph neural network then forecasts the evolution of the graph under mass-balance constraints with distribution-free conformal uncertainty quantification. Validation on the global aluminum supply chain shows that the learned graph representations recover economically meaningful supply chain structure, accurately forecast out-of-sample trade relationships, and produce well-calibrated prediction intervals. Subsequent papers apply this computational foundation to exposure assessment, disruption analysis, and scenario-based policy analysis, and extend the framework to multimaterial supply chain modeling and decision support.

36 MATERIALS SCIENCE

Gateway Element and Payload Materials Outgassing Analyses: HALO, HERMES, and ERSA

Gateway was intended to be humanity’s first space station around the Moon, but its development has been paused as the National Aeronautics and Space Administration (NASA) shifts focus to achieving the United States’ National Space Policy goals. Instead of an orbiting lunar outpost, NASA will now pursue the development of a lunar surface base to support a sustained human presence on the Moon. Before the program’s pause, Gateway’s Induced Environments team worked to ensure payloads and elements (i.e., modules) complied with induced environment requirements. Methods developed and insights gained from this work will have applicability to NASA’s Moon Base and the potential repurposing of Gateway elements and payloads, as well as to induced environments modeling for future space stations. The Gateway program’s induced environment included molecular contamination, electric thruster plume sputter and redeposition, and lunar dust transfer from the Human Landing System (HLS). Primary sources of external molecular contamination included materials outgassing, chemical thruster plume contamination, and vacuum venting. The focus of this paper will be on element- and payload-level materials outgassing analyses performed for Gateway Configuration 1, extending the previously-developed framework for Gateway system-level external molecular contamination modeling. Gateway Configuration 1 consisted of the Power and Propulsion Element (PPE) and the Habitation and Logistics Outpost (HALO). It also included payloads like the European Radiation Sensor Array (ERSA) attached to PPE and the Heliophysics Environmental and Radiation Measurement Experiment Suite (HERMES) attached to HALO. The element- and payload-level analyses to be introduced in this paper for HALO, HERMES, and ERSA enabled high-fidelity descriptions of Gateway’s external molecular contamination environment. Approaches to geometric modeling, meshing, outgassing rate assignment, molecular transport modeling, and analysis methodology will be presented. Element and payload contaminant deposition onto sensitive Gateway receiver surfaces will be summarized and results compared to induced environment requirements. While these results incorporate refinements made over the course of the program, they were not intended to be final. Therefore, modeling assumptions and inputs, potential improvements, and lessons-learned will be documented to inform future work on Moon Base, repurposed elements and payloads, and other space stations.

Gateway

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation