SEARCH · Engineering Papers

Results for “Flow Induced Vibration”

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.

SourceResearch connectionFrom yearThrough year

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

532 records · Page 27

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Computational simulation of asymmetric phase transformation in cracked Li₇La₃Zr₂O₁₂: Variant selection and chemo-mechanical implications

Coupling between microcracks and phase transformation in ion-conducting ceramics can jointly affect mechanical responses and ion transport. In this work we investigate the cubic-to-tetragonal phase transformation in Li₇La₃Zr₂O₁₂ in the presence of a microcrack under hydrostatic loading and quantify its implications for crack-tip stress concentration and Li-ion transport using phase-field and molecular dynamics simulations. The phase transformation exhibits a strong asymmetry between hydrostatic tension and compression. Under tension, the crack edge nucleates one tetragonal variant that amplifies the crack-tip stress intensity and promotes crack opening. Under compression, the crack tip nucleates a different tetragonal variant that enhances the stress-induced crack-closure tendency. Effective Li diffusivity analysis shows faster transport degradation under compression due to accelerated transformation kinetics, exposing a trade-off between mechanical stability and ionic conductivity. These results highlight the intertwined nature of cracking, phase transformation, and ionic transport in ion-conducting oxides and provide mechanistic insights into chemo-mechanical degradation of solid electrolytes.

36 MATERIALS SCIENCE↗

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 OSTI2025

Surface deposition of nanometer scale carbon structures on Bi 2 Sr 2 CaCu 2 O 8+δ using a low-cost scanning electron microscope

In this work, we investigate optically active nanostructures on Bi 2 Sr 2 CaCu 2 O 8+δ (BSCCO). The nanostructures are constructed from carbon nanocrystals formed using the electron beam of a scanning electron microscope to locally decompose residual hydrocarbons at the BSCCO surface. Atomic force microscopy is used to characterize the dimensions of these structures as a function of beam exposure time. This technique has been used to induce localized intercalation to form carbon nanoparticles up to tens of nanometers into a variety of 2D materials. Cross-sectional scanning transmission electron microscopy in BSCCO shows the presence of these carbon deposits on the surface, but there is no evidence they penetrate into the BSCCO substrate. Therefore, this technique is not suitable for localized intercalation in BSCCO, and it does have promise as a cost-effective approach to pattern nanometer scale etch stops.

2D materials↗

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 OSTI2026

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↗

DOE OSTI2026

Delamination recycling of multilayer plastic films: Solvent-assisted separation pathways and recovery of solid polyolefins

Multilayer plastic films are excellent packaging materials due to the bound layers of multiple polymers, with each different polymer contributing to the film properties. Desirable properties lead to continuously growing demand for multilayer films. Multilayer plastic films are typically single-use and, as such, their increased production and disposal have led to waste management problems. Multilayer films are not currently recyclable primarily due to the multiple bound polymers. This work advances delamination as a recycling process to separate and sequester valuable polymers from multilayer films, facilitating the incorporation of these polymers into the circular economy. We effect delamination in a physical process that preserves targeted polymers as solid film, hence retaining their embodied energy. This work documents three different pathways of inducing delamination of multilayer films, with appropriate solvents disrupting adhesion between adjacent layers or dissolving a minor component of the film in less than an hour and at temperatures below 90°C, and where all initial polyethylene is maintained in its solid form throughout the process and is recovered at high purity. Solvent-based delamination is exemplified on commercial multilayer films with majority polyethylene (PE), and polyethylene terephthalate (PET) or ethylene vinyl alcohol copolymer (EVOH) or nylon also present. Solvent-based delamination is an energy-efficient and environmentally friendly recycling process that improves upon dissolution-precipitation recycling, since delamination involves very little dissolution and no precipitation, and is superior to pyrolysis which breaks down the polymers, while delamination keeps polymer chains intact. Furthermore, delamination recycling can contribute to the sustainable use of multilayer films as they continue to protect our food and medicine.

Chemical Recycling↗

NASA NTRSDate not supplied

Performance of a Regenerative Fuel Cell System for the Lunar Surface

Regenerative fuel cells (RFCs) are an attractive energy storage solution for lunar missions as a technology capable of providing a higher specific energy (i.e., W∙h/kg) than state-of-the-art packaged Li-ion battery systems. An RFC consists of the (1 & 2) electrochemical stacks (chemical to electrical energy conversion to supply electricity to an external load, i.e. the fuel cell reaction, and electrical to chemical energy conversion of supplied electrical power to dissociate water into hydrogen and oxygen gases, i.e. water electrolysis), (3) fluidic conditioning, (4) reactant storage, (5) avionics, (6) power management and distribution (PMAD), and (7) thermal management. NASA’s Glenn Research Center has designed, assembled, and tested a breadboard RFC sys-tem capable of operating autonomously for multiple simulated lunar day/night cycles in a laboratory environment. The system is comprised of a non-flow through proton exchange membrane (PEM) fuel cell stack and a liquid-anode feed PEM electrolyzer (EZ) stack designed to electrochemically compress the reactants at balanced pressures up to 12.4 MPa (1800 psia). The fluidic conditioning, avionics, PMAD, and thermal management sub-systems are largely comprised of commercial-off-the-shelf components for this system-level development effort. The hardware is controlled by a CubeSat space processor running an operational program based on core flight architecture that can control the RFC hardware autonomously through a state machine with fault monitoring. The testing results highlighted here were completed with the system in an open-loop configuration such that reactants generated through water electrolysis were vented while gas cylinders supplied fuel cell operation. The breadboard operated autonomously, but there were five unplanned transitions to a safe state that required a manual restart after reviewing the data, determining a root cause, and implementing a solution. Four of the transitions were caused by the thermal management subsystem and the fifth was caused by a water management control issue in the EZ sub-system. The RFC system operated for over 550 hours with the final cycle being slightly abbreviated due to reasons unrelated to system performance.

Kerrigan Cain↗

NASA NTRSDate not supplied

Performance of a Regenerative Fuel Cell System for the Lunar Surface

Kerrigan Cain↗

DOE OSTI2026

Atomic Evolution of Hydrogen Intercalation Wave Dynamics in Palladium Nanocrystals Revealed by Liquid-Phase Transmission Electron Microscopy

Solute-intercalation-induced phase separation creates spatial heterogeneities in host materials, a phenomenon ubiquitous in batteries, hydrogen storage, and other energy devices. Despite many efforts, probing intercalation processes at the atomic scale has been a significant challenge. By utilizing liquid-phase transmission electron microscopy (TEM), we study hydrogen (de)intercalation in palladium nanocrystals as a model system and have achieved unprecedented atomic-resolution imaging of hydrogen intercalation wave dynamics. Our observations reveal that intercalation wave mechanisms, instead of shrinking-core mechanisms, prevail at ambient temperature for palladium nanocubes ranging from ∼60 nm down to ∼10 nm. Systematic image analysis uncovers the atomic evolution of the hydrogen intercalation wave, transitioning from nonplanar and inclined boundaries to those closely aligned with {100} planes. Our kinetic Monte Carlo simulations demonstrate that the observed intercalation wave dynamics correspond to sorption pathways minimizing the lattice mismatch strain at the phase boundary. In conclusion, unveiling the atomic intercalation pathways holds profound implications for engineering intercalation-mediated devices and advancements in energy sciences.

Lee, Daewon [Lawrence Berkeley National Laboratory↗

NASA NTRSDate not supplied

The Origin of Palimpsests

Palimpsests are the most enigmatic impact features observed on icy worlds of the Solar System. Palimpsests are large impact features which are bright, old, mostly flat, and roughly circular. They range from >150 to ~350 km in diameter and sometimes exhibit recognizable associated radial secondary ejecta crater strings beyond their rims. Curiously, among all the Solar System's icy worlds, palimpsests are exclusively seen on the icy Galilean moons. Palimpsests are obviously different than similarly sized impact features seen on other icy satellites. For example, Evander, a 350 km diameter peak-ring impact basin on Saturn’s moon Dione, has a very low depth/diameter ratio of 0.008, indicating significant viscous relaxation. Nevertheless, it still retains short wavelength topography including a prominent central peak, peak ring, and a continuous rim, none of which are seen at palimpsests. Palimpsests display an approximately circular, and dominating, expanse of undulating plains, which is level with surrounding terrain and shows no recognizable rim. At the center of palimpsests are smooth central plains that are enclosed by a low scarp. Central plains material superposes the surrounding undulating plains and is littered with small blocks reminiscent of “small chaos” on Europa. Palimpsests have yet to be explored using modern impact simulation tools. We have begun to investigate two competing hypotheses: (A) Palimpsests form on targets which have a preexisting near-subsurface layer that is liquid or solid but very weak, and is excavated as proximal ejecta; or (B) They form in a solid ice target underlain by a steep thermal gradient that strongly influences the resulting impact feature, including the production of impact-induced melting at depth that is then a major component of proximal ejecta. Our results will have profound implications for the early thermal history of the icy Galilean satellites. We will test our two hypotheses for the formation of palimpsests on Ganymede (where they are both best observed and preserved) by means of 2D iSALE hydrodynamic simulations of impacts into ice targets that will be structured according to the scenarios for each hypothesis, with select runs in 3D. We have already conducted a pilot study using 2D iSALE that has yielded promising results, which we will report out at the meeting.

palimpsests↗

DOE OSTI2026

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↗

DOE OSTI2026

Robustness of pairwise kinematic Sunyaev–Zel’dovich effect to optical-cluster-selection bias

The pairwise kinematic Sunyaev–Zel’dovich (kSZ) effect measures both the pairwise motion between galaxy groups and clusters and the amount of gas within them, providing a tracer for cosmic growth. To interpret the cosmological information in the kSZ measurements, it is crucial to understand the optical-cluster-selection bias on the kSZ observables. Line-of-sight structures that contribute to both the optical observable (e.g., richness) and the cosmological signal can induce a correlation between these two quantities at a fixed cluster mass. The selection bias arising from this correlation is a key systematic effect for cosmological analyses. For cosmological observables such as cluster abundance and weak lensing, controlling this selection bias may help explain the tension between the DES-Y1 results and the Planck constraints. In order to test for a kSZ effect equivalent of such a bias, we adopted an alternative mock richness based on galaxy counts within cylindrical volumes along the line of sight. We applied the cylindrical count method to hydrodynamical simulations across a wide range of galaxy-selection criteria, assigning richness consistent with DES-Y1 to the mock clusters. When comparing optically selected clusters to mass-selected halos, we find no significant bias on pairwise kSZ signals, pairwise velocities, or optical depth within our uncertainty limits of approximately 16, 10, and 8%, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2026

Scientific Highlights and Status of the Mars Science Laboratory Mission, Years 9 through 12

We present the scientific highlights and status from years 9-12 of NASA’s Mars Science Laboratory mission, with its Curiosity rover. A thermochemolysis-based sample analysis revealed a variety of organic molecules preserved in clay mineral-rich lake/lake margin deposits of the Glen Torridon region. The stratigraphically higher clay-sulfate mineral transition on Aeolis Mons was found to be associated with a lithological change from lacustrine to dry aeolian depositional environments. Water intermittently returned to the surface after the increase in aridity, including during the deposition of wave ripple strata which indicate a climate capable of sustaining an ice-free lake in the early Hesperian. The overlying Mg sulfate-bearing unit has both poly- and monohydrated Mg sulfates and is also enriched in iron carbonate of a purity consistent with chemical sedimentation and therefore atmospheric sequestration. A channel-incised canyon within Gediz Vallis and distal alluvial deposits recorded numerous fluvial and debris flow events after the deposition, lithification, and erosion of the clay- and sulfate-bearing units. Bright clasts found in Gediz Vallis canyon are composed of native sulfur. The investigation of the modern atmosphere and environment focused on the influence of the crater and mountain topography and captured high-altitude noctilucent and iridescent clouds. Radiation monitoring now spans an 11-year solar cycle. Rover systems and science instruments remain capable of addressing mission scientific objectives. All ten scientific instruments continue to return high-value data despite degradation of the Chemistry and Camera instrument’s high-voltage system, loss of multispectral imaging capability, and loss of active neutron spectroscopy. The mission has mitigated the loss of redundant brakes on two arm joints, additional wheel wear, and highly reduced capability of the redundant flight computer. Efficiency improvements in rover energy utilization and operations staffing have allowed the mission team to maintain a high level of productivity despite declining rover power generation and funding.

Mars↗

NASA NTRSDate not supplied

Material Properties and Modeling of Room Temperature Vulcanizing Silicone

Room Temperature Vulcanizing silicone (RTV) is a high-temperature adhesive that has successfully been used as a gap-filler between Thermal Protection System (TPS) tiles for heatshields on numerous missions. It is also used to bond instrumentation plugs such as temperature and pressure sensors into the heatshields. While RTV has been traditionally assumed to be a non-porous and non-ablating material, numerous experiments have shown that RTV pyrolyzes and becomes highly porous as it is heated. Heating RTV has also shown swelling, or intumescence, which can pose unique problems that lead to roughness induced boundary-layer transition, surface oxide formation and contamination of heat shield sensors. Therefore, it is crucial to understand and model the intumescence phenomenon of RTV. As data for RTV material properties is limited, the first step in modeling RTV is to collect material properties such as pyrolysis mass-loss, microstructure change, virgin and char porosity, etc. which was performed in our initial study. Additionally, thermomechanical properties such as Young’s modulus and Poisson ratio are required for modeling the intumescence of RTV, which were taken from literature and the coefficient of thermal expansion was collected using in-situ heating and Micro Computed Tomography (µ-CT) in previous studies. Finally, numerous other properties such as pyrolysis gas properties, virgin and char thermal conductivity and specific heat were compiled from previous experiments and literature into a material database that can be used for simulations. In Porous Material Analysis Toolbox based on OpenFOAM (PATO) [4], structural mechanics coupled with material response was used for simulating the intumescence of RTV as it is heated. However, since the permeability of the material is very low, the pyrolysis gas creates an internal pressure build-up as the material is being heated, significantly contributing to the deformation of the material. To correctly characterize this phenomenon, additional physics models were implemented into PATO's stress analysis solver, and results were compared with RTV dilatometry test data as a preliminary verification case. Future work will include experiments of RTV at the Plasmatron X facility and the in-situ heating cell with µ-CT, and improvement of simulation tools to more accurately model RTV intumescence.

PATO↗

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

Discovery of magnetic-field-tunable density modulations and spin tilting in a layered altermagnet

Altermagnets recently emerged as a new class of magnetic materials, arising from specific spin crystal symmetries. They exhibit a spin-polarized electronic band structure similar to ferromagnets, yet possess zero net magnetization, promising exotic properties. Here we study a layered triangular lattice altermagnet, cobalt-intercalated NbSe 2 using scanning tunneling microscopy and spectroscopy (STM/S). Spectroscopic-imaging STM and spin-polarized STM reveals emergent 2 a 0 tri-directional charge and spin density modulations on the selenium surface. Density functional theory simulations suggest these modulations reflect the underlying cobalt superstructure. We discover that an out-of-plane magnetic field tunes the modulation amplitudes and the electronic density-of-states in a manner dependent on the field direction and strength. This behavior is attributed to the field-induced tilting of cobalt spins, which can have profound implications on the electronic properties of the altermagnet. Our results provide atomic-scale insights to uncover a magnetic-field tunable altermagnetic band structure, highlight the importance of understanding spin canting in altermagnets.

condensed-matter physics↗

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↗