Search NASA⌕ Search

SEARCH · Search NASA

Results for “electrokinetic flows”

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.

106 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A Full-Induction Magnetohydrodynamics Solver for Liquid Metal Fusion Blankets in Vertex-CFD

Multiphysics modeling of liquid metal fusion blankets, which produce tritium and convert energy of neutrons created via fusion reactions into heat, is crucial for predicting performance, ensuring structural integrity, and optimizing energy production. While traditional blanket modeling of liquid metal flows during normal steady operating conditions commonly employs the inductionless approximation of the magnetohydrodynamics (MHD) equations, transient scenarios, when the plasma-confining magnetic field varies on millisecond time scales, require a full-induction MHD approach that dynamically evolves the magnetic field via the time-dependent induction equation. This paper presents the formulation, implementation, and initial verification of a full-induction MHD solver integrated within the open-source Vertex-CFD framework, which aims to achieve tight multiphysics coupling, a flexible software design enabling easy extension and addition of physics models, and performance portability across computing platforms. The solver utilizes finite element spatial discretization, implicit Runge–Kutta time integration, and an inexact Newton method to solve the resulting discrete nonlinear system, leveraging Trilinos packages for efficient computation. Verification against selected benchmark problems demonstrates accuracy and robustness of the solver. Furthermore, when the solver is applied to an idealized blanket model in 2.5D and full 3D, results obtained with Vertex-CFD are in good agreement with recently published quasi-2D simulations. These findings establish a computational foundation for future simulations of transient MHD phenomena in liquid metal blankets with Vertex-CFD, and open avenues for future extensions and performance optimizations.

Endeve, Eirik [ORNL] (ORCID:0000000312519507)↗
DOE OSTI2026

Two-dimensional non-equilibrium melting of charged colloids

Thermodynamic two-dimensional melting has been extensively studied in experiments and simulations, and is well predicted by theory. For systems in equilibrium, this transition is well described by the Kosterlitz–Thouless–Halperin–Nelson–Young theory, where melting is directly linked to the unbinding of topological defects. For driven, non-equilibrium melting and other non-equilibrium phase transitions, the picture is less clear. Here, in this work, we study the two-dimensional melting of a crystal of charged colloids. By randomly replacing some charged colloids with magnetic colloids, we can melt our system by rotating a fraction of the particles to create non-equilibrium, hydrodynamic random flows and local stresses. We can also melt it thermally by changing the particle number density. We find that an effective temperature approach cannot explain the results of our driven system. Rather, in both experiments and simulations, we observe that plotting the hexatic order parameter and the hexatic correlation’s exponent versus the density of disclinations and dislocations, respectively, yields universal curves. This implies that in our systems, two-dimensional melting depends directly on the density of topological defects and is independent of whether thermal or non-equilibrium forces generate them.

critical phenomena↗
DOE OSTI2026

SimH 2 : an integrated techno-economic modeling framework for hydrogen pipeline infrastructure and network optimization

Large-scale hydrogen (H 2 ) pipeline transport design and network optimization have seldom been reported due to the lack of a cost model accounting for the relationship between transport cost and hydrogen mass flow rate. Here, this work introduced a system-level cost model for hydrogen pipeline transport at supercritical state and integrated it with an existing CO 2 pipeline network tool, SimCCS, for hydrogen-specific pipeline design and optimization. The Intermountain West (I-West) region of the U.S., historically dependent on fossil fuel-based economies, is chosen to demonstrate the capabilities of our H 2 pipeline cost model and transport network optimization platform called SimH 2 . Two scenarios are examined: one where the pipeline is not allowed to pass through disadvantaged communities and the other where it is permitted. The results highlight that incorporating disadvantaged-community constraints lead to longer pipeline routes and increased transport costs, reflecting the trade-offs involved in equitable infrastructure development. It is demonstrated that the newly developed SimH 2 tool not only enables the efficient design of H 2 transportation pipelines but also optimizes the network by accounting for local terrain and the presence of disadvantaged areas.

08 HYDROGEN↗
DOE OSTI2023

Host galaxies of ultra-strong Mg ii absorbers at z ∼ 0.7

ABSTRACT We report spectroscopic identification of the host galaxies of 18 ultra-strong Mg ii systems (USMg ii) at 0.6 ≤ z ≤ 0.8. We created the largest sample by merging these with 20 host galaxies from our previous survey within 0.4 ≤ z ≤ 0.6. Using this sample, we confirm that the measured impact parameters ($\rm 6.3\leqslant D[kpc] \leqslant 120$ with a median of 19 kpc) are much larger than expected, and the USMg ii host galaxies do not follow the canonical $\rm {\it W}_{2796}-{\it D}$ anticorrelation. We show that the presence and significance of this anticorrelation may depend on the sample selection. The $\rm {\it W}_{2796}-{\it D}$ anticorrelation seen for the general Mg ii absorbers show a mild evolution at low $\rm W_{2796}$ end over the redshift range 0.4 ≤ z ≤ 1.5 with an increase of the impact parameters. Compared to the host galaxies of normal Mg ii absorbers, USMg ii host galaxies are brighter and more massive for a given impact parameter. While the USMg ii systems preferentially pick star-forming galaxies, they exhibit slightly lower ongoing star-forming rates compared to main sequence galaxies with the same stellar mass, suggesting a transition from star-forming to quiescent states. For a limiting magnitude of mr < 23.6, at least 29 per cent of the USMg ii host galaxies are isolated, and the width of the Mg ii absorption in these cases may originate from gas flows (infall/outflow) in isolated haloes of massive star forming but not starbursting galaxies. We associate more than one galaxy with the absorber in $\ge 21~{{\ \rm per\ cent}}$ cases, where interactions may cause wide velocity spread.

Astronomy & Astrophysics↗
DOE OSTI2022

Efficient electrocatalytic valorization of chlorinated organic water pollutant to ethylene

Electrochemistry can provide an efficient and sustainable way to treat environmental waters polluted by chlorinated organic compounds. However, the electrochemical valorization of 1,2-dichloroethane (DCA) is currently challenged by the lack of a catalyst that can selectively convert DCA in aqueous solutions into ethylene. Here, in this work, we report a catalyst comprising cobalt phthalocyanine molecules assembled on multiwalled carbon nanotubes that can electrochemically decompose aqueous DCA with high current and energy efficiencies. Ethylene is produced at high rates with unprecedented ~100% Faradaic efficiency across wide electrode potential and reactant concentration ranges. Kinetic studies and density functional theory calculations reveal that the rate-determining step is the first C–Cl bond breaking, which does not involve protons—a key mechanistic feature that enables cobalt phthalocyanine/carbon nanotube to efficiently catalyse DCA dechlorination and suppress the hydrogen evolution reaction. The nanotubular structure of the catalyst enables us to shape it into a flow-through electrified membrane, which we have used to demonstrate >95% DCA removal from simulated water samples with environmentally relevant DCA and electrolyte concentrations.

Choi, Chungseok [Yale Univ., New Haven, CT (United↗
DOE OSTI2023

CFD Simulation of the Dosing Behavior within the Atomic Layer Deposition Feeding System

The effective operation of atomic layer deposition (ALD) feeding system is the premise of realizing specific ALD processes. In the present work, a detailed computational fluid dynamics (CFD) model of the feeding system has been developed and validated, which accounts for the roles of ALD valves and manifolds. A numerical simulation of the compressible fluid flow and heat/mass transfer within the feeding system was conducted. The dosing amounts and the spatiotemporal distributions of the precursors can be accurately predicted using the CFD model, as validated by experimental results. Different precursors, operating conditions, and structures of the feeding system were simulated and analyzed to examine the operating flexibility of the feeding system. The simulation results can be adopted as the upstream boundary conditions for simulations of the ALD process in the reaction chamber. The substrate-scale simulation indicates that the effect of the feeding system on the film deposition is highly related to the surface kinetics of ALD. The present work can serve as a guide for the development and optimization of different ALD-based processes via proper operation and even the design of the feeding system.

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

SAM Finite Volume Method Development Status Update: GCR Application, Restart, and MultiApp

The System Analysis Module (SAM) is being developed as a modern system analysis code for advanced non-light-water-reactor safety analysis under the U.S. DOE NEAMS program. Previous feasibility studies have demonstrated that a staggered-grid finite volume method (SG-FVM), implemented under the MOOSE framework, can deliver more than an order of magnitude speedup over the existing continuous Galerkin finite element method (CG-FEM) solver for liquid-cooled, incompressible but thermally expandable flow systems. This work extends the previous effort to compressible, gas-cooled reactor applications, where pressure couples directly into the mass equation adding additional nonlinearity into the equation system. New code capabilities are implemented for pebble bed high-temperature gas-cooled reactor (PB-HTGR) analysis, including a pebble bed CoreChannel component, built-in pebble bed effective thermal conductivity model and channel-to-channel crossflow model. The capabilities are tested, benchmarked, and demonstrated for problems with increased level of model and physical complexities, including the HTTU effective thermal conductivity test, the SANA passive cooling test, and a demonstration case using the GPBR200 reactor design covering steady-state operation, DLOFC and PLOFC transients. Across all cases, the SG-FVM solver demonstrated strong robustness and efficiency, and the solutions agree well with reference results and data. The finding of this work proves that SG-FVM is a viable and efficient solver pathway for compressible, gas-cooled reactor system analysis in SAM. In addition, work has been done to successfully support SAM-FVM recover/restart code feature that is essential to reactor safety analysis applications, and MultiApp code feature that is essential to multi-scale and multi-physics simulations. In summary, this work continued from previous feasibility studies, and further demonstrated that the SG-FVM will serve as a strong foundation for SAM’s advanced solver algorithm for future deployment.

Zou, Ling↗
DOE OSTI2025

Analysis of Liquids Using the Submerged Plasma for Isotopic Detection and Elemental Resolution (SPIDER)

We developed a sensor called the Submerged Plasma for Isotopic Detection and Elemental Resolution (SPIDER) probe, which uses an atmospheric pressure glow discharge below the surface of liquids to excite species in the liquid. Through emission spectroscopy of molten salts, liquid metals, and heavy water, we demonstrated the SPIDER probe’s high resolution, accuracy, and versatility. We successfully identified trace concentrations of transition and rare-earth metals in molten salts and detected the isotopic shift of the H β → D β emission line. Our analysis revealed unconventional spectral alkali line shapes, indicating two competing excitation modes: film explosion and droplet vaporization. The film explosion mode, characterized by dense plasma, exhibited self-reversal and broadband continuum emission, while the droplet vaporization mode, associated with diffusive plasma, produced narrow-line emissions. Furthermore, by analyzing circuit transients alongside individual plasma events, we observed that film explosions generate higher currents, likely due to a shorter plasma length as the current preferentially flows through the thin liquid layer. Altogether, our results highlight the SPIDER probe’s efficacy and flexibility, making it well-suited for online material quantification of liquids in extreme environments.

AES↗
DOE OSTI2026

Compatibility of Type 304H Stainless Steel in Static Molten FLiBe for Inertial Confinement Fusion Reactors: Role of Impurities and Redox Control

Molten fluoride salts, such as FLiBe (LiF-BeF2), are promising candidates for tritium breeding and heat transfer in fusion reactors, but corrosion of structural materials remains a major challenge. This study investigates the corrosion behavior of austenitic stainless steel 304H in purified and NiF2-containing FLiBe at 500°C and 600°C, focusing on the effects of impurities and redox control. Exposure to purified FLiBe resulted in the concurrent depletion of Cr, Mn, and Fe, with corrosion at 500 °C dominated by the combined oxide formation and elemental dissolution, while at 600°C elemental depletion was predominant. The addition of a controlled NiF2 impurity significantly accelerated corrosion at both temperatures, demonstrating the sensitivity of 304H to the salt redox state. Beryllium additions were effective in mitigating corrosion for both baseline and NiF2-containing FLiBe; minimal depletion of Cr, Mn, and Fe occurred with Be additions as low as 2.5 mg (147 wppm), and no NiBe intermetallics formed at 5 mg (294 wppm), indicating that small Be inventories can provide substantial protection without deleterious phase formation. Thermodynamic equilibrium and coupled thermodynamic-kinetic analyses at the salt-alloy interface suggested low corrosion rates in systems with limited hydrogen fluoride (HF) generation, highlighting the importance of salt redox control. Estimates for a Be addition rate were calculated for the HYLIFE-II fusion reactor that can mitigate corrosion-induced degradation, assuming complete conversion of tritium to tritium fluoride (TF). Overall, 304H shows reasonable compatibility with FLiBe under optimized redox conditions. These results provide quantitative guidance for material selection and salt management in fusion blanket and heat exchanger systems and motivate validation under flowing, nonisothermal, and irradiated conditions.

Pillai, Rishi [ORNL] (ORCID:0000000243688197)↗
DOE OSTI2026

Kinetic control of phase evolution and defect-mediated recombination in AACVD-grown copper antimony sulphide thin films

Copper antimony sulphide (CAS) is a multinary chalcogenide semiconductor in which small deviations from stoichiometry can drive phase competition and strong defect-mediated modulation of optoelectronic properties. However, the systematic roles of copper precursor fraction and deposition time in governing phase evolution, off-stoichiometry, and recombination dynamics in aerosol-assisted chemical vapour deposition (AACVD)-grown undoped CAS thin films remain insufficiently understood. In this work, undoped CAS thin films were deposited by AACVD using Cu(dedtc)2 and Sb(dedtc)3 single-source precursors at 550 °C and a carrier gas flow rate of 150 sccm, while the Cu(dedtc)2 mole fraction (x = 0.15 – 0.55) and deposition time (1 – 2 hours) were systematically varied to probe how growth kinetics influence phase composition, microstructure, and defect-mediated optical properties without post-deposition annealing or extrinsic doping. Increasing copper precursor content drives phase evolution toward tetrahedrite-dominant CAS films at intermediate Cu(dedtc)2 mole fractions, with the film deposited at x = 0.35 exhibiting the strongest tetrahedrite character within the parameter space examined. The films are also copper-rich, antimony-poor, and sulphur-deficient, consistent with off-stoichiometric growth and intrinsic defect formation, plausibly including copper interstitials, Cu-on-Sb antisites, and sulphur vacancies. These growth-dependent compositional deviations are accompanied by tunable indirect optical bandgaps of approximately 1.60 – 2.18 eV and weak visible photoluminescence governed by defect-mediated recombination. Time-resolved photoluminescence reveals bi-exponential decay behaviour with lifetimes of approximately 0.1 – 3.6 ns, with emission dominated by slower donor–acceptor pair recombination and a smaller contribution from faster trap-assisted pathways. Collectively, these results establish an explicit kinetic process–structure–defect–property relationship for AACVD-grown CAS thin films and provide a growth–structure–property framework relevant to future optimization of CAS-based optoelectronic and energy materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Statistical Study of Energy Dissipation in Magnetic Structures During Turbulent Reconnection in the Earth's Magnetotail

Magnetic reconnection is a ubiquitous plasma phenomenon that plays a critical role in particle heating and energization. During reconnection, the topology of magnetic field rearranges, depositing energy into the surrounding plasma through bulk flow, thermal heating, or non-thermal particle acceleration. While the pathways of this transformation from magnetic energy into kinetic have been studied extensively in recent years through theoretical or case-by-case observations, comprehensive statistical studies remain limited. In this paper, we present a statistical investigation using data from the Magnetospheric Multiscale (MMS) mission, and detail the particle energization mechanisms in magnetic structures found near reconnecting regions in turbulent Earth's magnetotail. We find that electrons with motion perpendicular to the magnetic field dominate $\vec{j}$ ⋅ $\vec{E}$ dissipation. In contrast to the conventional picture of unidirectional energy transfer to particles by laminar two-dimensional (2D) reconnection, we find that energy exchange within magnetic structures during turbulent reconnection tends to be bidirectional with only a small positive bias from electromagnetic fields to particles. Specific electron energization mechanisms are quantified, including those due to parallel electric field, Fermi energization from curvature drift, betatron heating from magnetic field inhomogeneity, and polarization drift.

Wang, Rachel [Princeton University, NJ (United Sta↗
DOE OSTI2026

Direct observation of diamond formation in a shock-compressed high explosive

Understanding the formation timescale and structure of carbonaceous reaction products is critical for modeling the high-pressure equation-of-state of organic materials. We use the National Ignition Facility to shock-compress polycrystalline TATB (C 6 H 6 N 6 O 6 ) samples to ~70–130 GPa and ~4000–5500 K, employing in situ nanosecond X-ray diffraction to probe reaction products and velocimetry to measure transmitted compression wave profiles. Our diffraction data is consistent with the formation of diamond over timescales less than ~60 ns. This represents carbon condensation from a molecular explosive on timescales three times faster than previously reported and the earliest observation of diamond produced from reacting TATB. Reactive flow simulations with explicit chemistry reproduce the observed temporal structure within wave profiles to inform the distribution of P-T states. These findings provide direct evidence of ultrafast diamond formation in a reactive system at extreme conditions and provide new constraints for models of shock and detonation chemistry.

Clarke, Samantha M. [Lawrence Livermore National L↗
DOE OSTI2026

Thermal property characterization of phase change materials in building applications: A systematic review of fundamentals, recent progress, and future directions

Phase change materials (PCMs) can reduce building peak loads and enable demand-responsive thermal energy storage (TES), but their deployment depends on reliable measurement and interpretation of thermal properties across laboratory, intermediate, and application scales. Here, this review systematically examines characterization methods, testing protocols, and recent advances for neat PCMs and PCM composites, emphasizing thermal conductivity, enthalpy-related properties (phase change temperature, latent heat, specific heat), and cycling stability. For thermal conductivity, we compare steady-state and transient techniques and note limitations when phase transition and contact resistance affect measurements. For enthalpy–temperature characterization, we discuss differential scanning calorimetry together with intermediate- and bulk-scale methods, including T-history, heat flow meter testing, and three-layer calorimetry (3LC), to generate application-relevant enthalpy–temperature profiles. Cycling stability is organized into four experimental families: thermoelectric–air, fully thermoelectric, water-bath, and in situ chamber approaches, with attention to separating reversible supercooling from true degradation such as phase segregation. We highlight emerging noncontact diagnostics, including infrared thermography and embedded sensing, for spatially resolved validation and multiscale interpretation. Finally, we review the growing use of AI and machine learning for property prediction, inverse characterization from experimental signals, and real-time state estimation in building-integrated TES. Key needs include harmonized protocols, interlaboratory benchmarking, uncertainty reporting, and metadata-rich datasets to accelerate reproducible PCM qualification for grid-flexible buildings.

AI↗
DOE OSTI2025

Preliminary Plan to Inform Testing of a Heat Exchanger Test Article

This report presents a preliminary plan to guide the qualification testing of advanced heat exchanger (HX) components for nuclear-to-industrial heat transfer applications. The objective is to establish a defensible, physics-based methodology that integrates computational modeling, targeted experimentation, and in-service inspection considerations to demonstrate component performance and reliability under representative reactor conditions. The analysis identifies Sodium-cooled Fast Reactor (SFR) and High-Temperature Gas-cooled Reactor (HTGR) systems as reference configurations in terms of temperature, pressure, and chemical environment. Within these operating envelopes, dominant degradation mechanisms— including creep–fatigue interaction, flow-induced vibration, corrosion, and diffusion-bond deterioration—were evaluated to define test requirements. A comprehensive computationalexperimental framework is proposed to support life prediction and qualification activities. The framework couples high-fidelity structural-mechanics, thermal-hydraulic, and fluid-structure interaction models with accelerated degradation testing to produce a traceable linkage between microstructural evolution, mechanical performance, and remaining useful life (RUL). The approach adheres to established Verification, Validation, and Uncertainty Quantification (VVUQ) standards (ASME V&V 10/20; NUREG-2152) and incorporates a digital-twin architecture for continuous model refinement through data assimilation. The plan further outlines testing methodologies, including pre-test analyses, test-loop design parameters, and sensor placement strategies that maximize information yield while maintaining mechanistic fidelity. Complementary sections describe in-service inspection (ISI), on-line monitoring (OLM), and structural-health-monitoring (SHM) techniques applicable to compact HX geometries typical of advanced reactors. Collectively, these activities establish the technical foundation for demonstrating 40-60-year equivalent service life of advanced heat exchangers in support of the U.S. Department of Energy’s Advanced Reactor and Integrated Energy Systems programs. The forthcoming phase will execute the defined pre-test analyses, initiate hardware fabrication, and implement the integrated testing campaign to validate the proposed qualification methodology.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
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 OSTI2026

Distributed Mafic Rock Resources for Carbon Mineralization in Arizona

Ex-situ carbon mineralization is a process by which CO2 is reacted with alkaline silicate minerals and rocks to produce stable carbonate materials, which can be used for other industrial processes. Arizona, U.S.A., hosts abundant surficial mafic rocks in three young volcanic fields, Geronimo-San Bernardino, San Francisco, and Springerville, and other distributed locations throughout the state. We created a Mafic Rock Resource Inventory (MRRI) that categorizes geochemical, physical, and textural characteristics of a diverse suite of surficial mafic rock samples and provide a benchmark reaction dataset parameterizing the temperature, pressure, and pH conditions best suited ex-situ mineralization in different rock types. MRRI data is publicly available online via a map-viewer. We establish two reaction condition sets, varied in temperature, pressure, and pH, where crystal-rich and glassy rocks reach maximum reaction extent and different carbonate phases are formed. Systematic ex-situ mineralization experiments on 21 diverse rock types show trends in geochemical, mineralogical, and reactivity behavior and establish maximum effective capture capacity. From this, scoria cones in three Arizona volcanic fields have a ~62 Gt effective CO₂ storage capacity with one of the fields having a ~42 Gt storage capacity in lava flows. Reactivity results have applications to alkaline mafic rock resources exposed globally, including producing additional effective storage capacity estimates and scaled commercialization of mafic rock ex-situ mineralization, should reaction extents be improved through advances in mineralization techniques. MRRI data were used to create a Direct Air Capture to Mineralization (DACM) systems model, technoeconomic analysis, and life-cycle assessment. These documents are presented as three appendices.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Probing the atomic dynamics of ultrafast melting with femtosecond electron diffraction

Melting is an every-day phase transition that is determined by thermodynamic parameters like temperature and pressure. In contrast, ultra-fast melting is governed by the microscopic response to a rapid energy input and, thus, can reveal the strength and dynamics of atomic bonds as well as the energy flow rate to the lattice. Accurately describing these processes remains challenging and requires detailed insights into transient states encountered. Here, we present data from femtosecond electron diffraction measurements that capture the structural evolution of copper during the ultrafast solid-to-liquid phase transformations. At absorbed energy densities 2-4 times the melting threshold, melting begins at the surface slightly below the nominal melting point followed by rapid homogeneous melting throughout the volume. Molecular dynamics simulations reproduce these observations and reveal a weak electron-lattice energy transfer rate for the given experimental conditions. Both simulations and experiments show no indications of rapid lattice collapse when its temperature surpasses proposed limits of superheating, providing evidence that the inherent dynamics limits the speed of disordering in ultrafast melting of metals.

FOS: Physical sciences↗
DOE OSTI2026

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↗