Search NASA⌕ Search

SEARCH · Search NASA

Results for “Liquified pressurized gas electrolyte”

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.

355 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

HTGR Multiphysics Application Drivers FY26 Updates

This report summarizes FY26 progress under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program's high-temperature gas-cooled reactor (HTGR) application driver work, covering a wide range of activities such as code validation and multi-physics code assessment. 1) A detailed SAM model of the High-Temperature Engineering Test Reactor (HTTR) was developed using a unique-block grouping approach, with an extended parallel thermal network method to capture block-to-block conduction and radiation heat transfer, and applied to steady-state simulations of the HTTR 30~MW and 9~MW cases. 2) In another activity, SAM's newly implemented multi-component gas flow model was validated against the Natural convection Shutdown heat removal Test Facility (NSTF) argon ingress experiment, correctly capturing the density-driven suppression and thermal recovery of natural circulation observed when argon is introduced into the air-cooled Reactor Cavity Cooling System (RCCS) loop. 3) For the OECD/NEA High Temperature Test Facility (HTTF) benchmark, we co-led the international benchmark activities as well as the OECD/NEA final benchmark report to be released at the end of this year. 4) Finally, the coupled Griffin-SAM modeling capability for pebble-bed HTGRs was advanced by verifying the Griffin neutronics solution against Serpent Monte Carlo for a realistic non-uniform temperature distribution, resolving several deficiencies in the SAM-to-Griffin temperature transfer scheme, and enabling distinct fuel kernel, moderator, and coolant temperatures for cross section feedback. These new features were demonstrated in a PBR load-following transient.

Lee, Alvin↗
DOE OSTI2026

Performance of Heterostructural TaC/AlGaN Schottky Diodes Based on First Principles Electronic Structure Properties

Advances in ultra-wide bandgap materials, such as high Al-content AlxGa1-xN (AlGaN), are essential for next generation power electronics, but the requirement for lattice matched substrates is currently a significant obstacle. Recently, conductive TaC has emerged as a promising virtual substrate for AlGaN heteroepitaxy, with wurtzite (0001) AlxGa1-xN lattice-matched to rocksalt (111) TaC at x ~ 0.5. Thus, understanding and controlling the electronic properties of the TaC/AlGaN interface is key for developing technological applications based on TaC/AlGaN devices. Using density functional theory and electronic structure calculations, we here investigate TaC/Al0.5Ga0.5N interfaces, where we include explicit alloy models in the slab calculations. We predict the Schottky barrier height and the electric field discontinuity resulting from interface charges. Considering all possible combinations of (Ta, C) substrate termination, (Al/Ga, N) nucleation, and (Al/Ga, N) polarity, we construct a chemical potential phase diagram to identify the stable interfaces that can be accessed through variation of the synthesis conditions. The predicted interface electronic properties are implemented in device performance simulations to demonstrate a practical design for a strain-free, high-efficiency TaC/AlGaN Schottky diode with a low barrier height and without interface charges, underscoring the potential of TaC as a substrate for ultra-wide bandgap devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Iron surface corrosion in supercritical CO2 at atomic scale investigated by molecular dynamics simulations

Understanding the corrosion behavior of steels in supercritical carbon dioxide (S-CO2) is essential for ensuring the safe application of S-CO2 as a heat-transfer fluid in high-temperature energy systems, including advanced nuclear reactors. In this work, molecular dynamics (MD) simulations using ReaxFF potential are performed to explore the atomic-scale corrosion mechanisms of body-centered cubic iron (BCC-Fe) in S-CO2. The results show that CO2 molecules in S-CO2 decompose at the Fe surface, generating free C and O atoms that form Fe-C and Fe-O bonds and subsequently produce oxides and carbides. Concurrently, Fe atoms dissolve from the surface and diffuse into the S-CO2 region, resulting in interdiffusion of Fe, C and O atoms at the interface. The corrosion-layer thickness calculations show that high pressure and temperature induced by S-CO2 have stronger effects than surface orientation on the corrosion process. In addition, surface Fe atoms undergo substantial displacement under S-CO2 exposure, further accelerating corrosion. When a radiation-induced void is introduced near the Fe surface, the corrosion is enhanced. The void-matrix interface expands the reaction surface area and simultaneously induces corrosion reactions inside the bulk, resulting in a deeper penetration of C and O and thicker corrosion layers. All these results indicate that high-temperature, high-pressure and radiation-induced voids can seriously affect the corrosion of Fe in S-CO2, and must be considered to better use S-CO2 in nuclear facilities.

Li, Wenhua↗
DOE OSTI2026

Data-Driven Discovery and Experimental Validation of Solvent Polarity Effects on Conjugated Polymer Solution-to-Film Assembly Pathways

Understanding how solvent properties influence the solution-to-film assembly of conjugated polymers remains a critical challenge due to the complex and intertwined nature of polymer–solvent interactions. In this study, we integrate a data-driven framework with experimental validation to identify key parameters influencing the assembly and performance of poly[2,5-(2-octyldodecyl)-3,6-diketopyrrolopyrrole-alt-5,5-(2,5-di(thien-2-yl)thieno[3,2-b]thiophene)] (DPP-DTT) in organic field-effect transistors (OFETs). A machine learning (ML) approach identified the normalized Reichardt polarity parameter (E T N ) as a significant descriptor correlated with DPP-DTT hole mobility (μ). Systematic DPP-DTT devices fabricated using solvents across a wide E T N range revealed that higher E T N solvents yield enhanced μ. To elucidate the structural origins of high μ, we conducted comprehensive analyses using UV–vis–NIR spectroscopy and grazing incidence wide angle X-ray scattering (GIWAXS) measurements. The results revealed that films processed from high E T N solvents exhibit reduced paracrystallinity. By analyzing the solution-state behavior using optical microscopy and solution WAXS, we revealed polymer solubility differences in the various solvents and associated distinct polymer assembly pathways, elucidating why the high E T N solvent produces long-range ordered films. Notably, the high E T N solvent shows a pronounced preference for liquid-crystal (LC)-mediated assembly, providing a mechanistic explanation for the enhanced structural order. Therefore, these results demonstrate that solvent polarity, as evaluated by E T N , serves as an important parameter that plays a significant role in the DPP-DTT assembly pathway and resultant solid-state morphology. This work provides a strategy for integrating data science with experiments to identify critical parameters associated with complex polymer systems and helps guide rational process design for high-performance organic electronics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Propane Activation on Pt Electrodes at Room Temperature: Quantification of Adsorbate Identity and Coverage

Abstract C−H bond activation is the first step in manufacturing chemical products from readily available light alkane feedstock and typically proceeds via carbon‐intensive thermal processes. The ongoing emphasis on decarbonization via electrification motivates low‐temperature electrochemical alternatives that could lead to sustainable chemicals production. Platinum (Pt) electrocatalysts have shown activity towards reacting alkanes; however, little is known about propane electrocatalytic activation and conditions suitable for enabling selective oxidation to valuable products. Herein, we utilize a combination of electrochemical mass spectrometry (ECMS) and density functional theory (DFT) calculations to elucidate the potential dependence of propane activation on Pt electrocatalysts. Results show a strong dependence of adsorption on the applied potential in room‐temperature aqueous acidic electrolyte, with a maximum coverage of propane‐derived adsorbates at 0.30 V vs RHE. Using charge deconvolution and deuterated experiments, the mechanism of adsorption was elucidated, and C 3 H 2 * was determined as the average dehydrogenated propane‐derived adsorbate species. DFT calculations further corroborate these results, showing that the formation of deeply dehydrogenated species is energetically accessible at room temperature. The combined theoretical and experimental findings yield insights for selective activation of paraffinic C−H bonds at room temperature, aqueous conditions—a critical step towards decarbonized chemical manufacturing.

Chemistry↗
DOE OSTI2024

Crack‐Free Single‐Crystalline LiNiO 2 for High Energy Density All‐Solid‐State Batteries

Abstract Single‐crystalline layered oxide (LiNi 1‐ x ‐ y Mn x Co y O 2 ) cathodes have been found to exhibit exceptional electrochemical properties when coupled with various inorganic solid electrolytes (ISEs) in all‐solid‐state batteries (ASSBs). Their advantages stem from the robust morphological integrity with the absence of grain boundaries and the high electrochemical oxidative stability. Here, ASSBs featuring single‐crystalline LiNiO 2 (LNO) with the highest Ni content are reported, offering a high theoretical specific capacity of 275 mAh g ‐1 alongside a high average discharge voltage (3.7 V vs Li + /Li). Through a careful investigation, it is demonstrated that micron‐sized single‐crystalline LNO (µSC‐LNO) composite cathodes with a halide ISE exhibit a high initial discharge capacity of 205 mAh g ‐1 with an outstanding cycle performance over 200 cycles in room‐temperature ASSBs. The significance of engineering parameters is emphasized, such as particle size and specific density, in promoting a homogeneous and fast Li + transport within the composite cathodes. Furthermore, the formation of undesirable interphase between the halide ISE in the cathode and sulfide ISE separator is elucidated, which may be a critical factor impeding long‐term cyclability of ASSBs. This work provides insights into the design of composite cathodes for high‐energy‐density ASSBs.

Chemistry↗
DOE OSTI2026

Hierarchical low Pt-loading “core-shell” electrocatalysts for the oxygen reduction reaction in fuel cells

The sluggish kinetics of the oxygen reduction reaction (ORR) hinder cost-effective polymer electrolyte fuel cells (PEFCs), which rely on scarce, expensive platinum-based electrocatalysts (ECs). Here, we present a novel synthesis method for ORR ECs achieving exceptional platinum utilization. The design features a hierarchical “multi-carbon” support comprising carbon nanoparticles interacting with graphene nanoplatelets as the “core”, encapsulated by a porous carbon nitride (CN) “shell”. This configuration promotes strong core/shell interactions and a bimodal active site distribution, consisting of chemically dispersed Pt and Ni single-atom complexes and PtNix alloy nanoclusters embedded in the CN shell. These advantages enable high activity and durability, achieving an ORR activity of 1.6 A mgPt−1 at 0.9 V vs. RHE-an order of magnitude higher than Pt/C (0.17 A mgPt−1). A proof-of-concept PEFC demonstrates a specific power of 12.0 kW gPt−1 at 0.60 V. This approach offers a significant step toward more efficient and sustainable PEFC technologies.

Pagot, Gioele [University of Padova, Italy]↗
DOE OSTI2025

CO Adsorption on a Single‐Atom Catalyst Stably Embedded in Graphene

Abstract Confined single metal atoms in graphene‐based materials have proven to be excellent catalysts for several reactions and promising gas sensing systems. However, whether the chemical activity arises from the specific type of metal atom or is a direct consequence of the confinement itself remains unclear.

Chemistry↗
DOE OSTI2026

Prioritizing Uncertainties in Hydrogen Contribution to Risk in Post-Crash Outcomes for Rail

This report presents analysis from Sandia National Laboratories predicting contributions to risk associated with the use of hydrogen technology for rail. Event sequence diagrams are used to describe possible accident scenarios and progressions. Initiating event frequencies and branch event probabilities for each scenario are quantified with uncertainty using distributions fit to Federal Railroad Administration and U.S. Department of Transportation Pipeline and Hazardous Materials Safety Administration data on applicable accidents from 2000 to 2020. Uncertainty is propagated through the event sequence diagram to estimate the frequency and conditional probability of accident end states. The analysis identifies four scenarios with significant contributions to risk from hydrogen that are predicted to occur relatively frequently, which may inform priorities for reducing uncertainty. These scenarios are 1) overpressure events resulting from collisions with hydrogen release due to mechanical damage and delayed ignition, 2) jet fire events resulting from collisions with hydrogen release due to mechanical damage and immediate ignition, 3) jet fires resulting from fire or explosion initiating events involving the hydrogen tank and correct operation of the thermally-activated pressure relief device (TPRD) subsequent to the thermal insult, and 4) pressure burst resulting from fire or explosion initiating events involving the hydrogen tank and failure of the TPRD. Delayed and immediate hydrogen ignition probabilities are identified as being highly uncertain and potential candidates for reducing conservatism in the predicted frequencies for these two scenarios.

08 HYDROGEN↗
DOE OSTI2026

Identifying microstructures susceptible to pulverization in commercially irradiated high burnup UO2 under LOCA conditions

High burnup fuel fragmentation (HBFF) has been a concern in the nuclear industry for many years. When UO2 reaches and exceeds pellet average burnups around 55 GWd/tU, microstructural changes in the fuel pellet begin to occur that render the pellet susceptible to fine fragmentation under loss-of-coolant accident (LOCA) conditions. During a LOCA event, the cladding balloons and bursts, potentially releasing part of these fine fuel fragments into the reactor pressure vessel. This process is known as fuel fragmentation, relocation, and dispersal (FFRD). Presently, the US nuclear industry is developing a safety basis for FFRD with the goal of increasing burnups and pressurized water reactor cycle lengths. Increased cycle lengths would push the burnup limits of the fuel rods past the known threshold for HBFF susceptibility. In this work, the microstructural impact on HBFF was investigated by comparing the microstructures of as-irradiated rod segments to their post-LOCA tested counterparts. This investigation found that varied operational histories result in different microstructural evolutions across the pellet radius, which impacts the fragmentation behavior and radial location of the fragmentation. It is noted that regions with a high bubble density and a high density of grain boundaries fragment under LOCA-relevant conditions. In addition to the fragmentation that has been previously reported at the periphery of the fuel, fragmentation was also noted in the dark zone toward the fuel center. Analysis of the power histories of the fuel samples suggests that the dark zone forms in the central regions of the pellet between temperatures of approximately 740 and 960 °C.

McKinney, Casey [ORNL] (ORCID:0000000335383614)↗
DOE OSTI2023

Optimizing Renewable Ammonia Production for a Sustainable Fertilizer Supply Chain Transition

Local renewable ammonia production using electrolytic hydrogen is an emerging approach to alleviate emissions attributed to synthetic nitrogen fertilizer production while also insulating against fluctuations in fertilizer prices and mitigating transportation costs and emissions. However, replacing ammonia currently produced using fossil fuels will not be immediate. To this end, we develop a supply chain transition model, which first optimizes the design and hourly operation of new renewable ammonia facilities to minimize production costs and then optimizes the annual installation timing, production scale, and location of these new renewable facilities along with ammonia transportation to meet county resolution demands. The objective is to augment and eventually replace conventional ammonia market imports in an economically competitive manner. We performed a case study for Minnesota's ammonia supply chain and found that a full transition to in-state renewable production by 2032 is optimal. This is incentivized by the U.S. federal government's clean hydrogen production credits. This transition results in 99 % reduction in carbon intensity along with stable supply costs below $475 per metric tonne. New renewable production facilities are an order of magnitude smaller than existing conventional plants. They use both wind and solar resources and operate dynamically to minimize expensive battery and hydrogen storage capacities.

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

Semi-Dynamic Leach Testing of Densified Silicon-based Iodine Waste Forms

Iodine waste forms (IWF) require a conceptual corrosion release model (CCRM) to provide iodine (I) release rates for performance modeling nuclear waste disposal repositories. To develop a CCRM, an understanding of the corrosion mechanisms of the IWF are required along with data from consistent test methods to parameterize the model. The present study has advanced both areas by providing and assessing the corrosion resistance of IWF types based on their processing history in minor variations of semi-dynamic leach tests. The test included a series of semi-dynamic leach tests using monolithic IWFs in deionized water (leachant). Several experiments were conducted under alternate test conditions with changes to temperature, leachant replacement, leachant pH, leachant volume, masking, and surface finish to elucidate if varying these conditions impacted IWF corrosion behavior. Tests were conducted on two classes of IWFs: (1) I-bearing silver-mordenite (AgZ) materials processed by hot isostatic pressing (HIP) at different temperatures, pressures, sizes, and times; and (2) I-bearing silver-functionalized silica aerogels (SFA) processed by either HIP or spark plasma sintering (SPS). The corrosion susceptibility of AgZ samples was influenced by HIP temperature and pressure. The SPS SFAs retained I far better than HIP SFAs. Additional findings in this study include: (1) The iodine dissolution rate decreased with decreasing temperature, (2) a common ion effect may occur and slow dissolution of the host phase if the leachant is not regularly replaced, (3) pH controls the dissolution rate, and (4) the iodine dissolution rate slows with extended test time (up to 224 days). Based on this work, these parameters should thus be represented when developing a CCRM.

corrosion↗
DOE OSTI2026

Flame stabilization in DME spray flames under engine-relevant conditions characterized by OH* chemiluminescence and formaldehyde laser-induced fluorescence

The transient and quasi-steady flame structures of Dimethyl Ether (DME) fuel sprays, produced by a single-hole injector (Spray D), were investigated using Planar Laser-Induced Fluorescence (PLIF) and chemiluminescence imaging in a constant-volume chamber under Engine Combustion Network (ECN) Spray A conditions (900 K ambient temperature, 60 bar ambient pressure, 1500 bar injection pressure, and 22.8 kg/m 3 ambient density). Low-temperature chemical reaction zones were visualized using formaldehyde (CH 2 O) PLIF with 355 nm excitation, while high-temperature flame regions were captured via chemiluminescence imaging of excited-state hydroxyl radicals (OH*). Both transient and quasi-steady flame structures clearly show the transition from CH 2 O to OH*, highlighting the progression from low- to high-temperature combustion, while the position of the flame is displaced for DME compared to reference hydrocarbon n-dodecane. Homogeneous reactor calculations with detailed chemistry and using adiabatic mixing for initial temperature show that CH 2 O peaks are significantly higher for DME at the same equivalence ratio, with a higher heat-release during the cool-flame regime with respect to the fuel heating value. Thus, the cool-flame dynamic as a precursor to high-temperature combustion and flame stabilization exhibit distinct behavior for DME relative to conventional hydrocarbons, and these phenomena are effectively resolved through the soot-free nature of DME and the high-speed, time-resolved diagnostics.

CH2O laser-induced fluorescence↗
DOE OSTI2026

The 10 September 2025 M w 4.1 Earthquake in Northeastern Utah, United States: An Archetypal Continental Mantle Event

The 10 September 2025 M w 4.1 earthquake in northeastern Utah, United States, had a focal depth 68 km beneath sea level, which is ∼20–25 km greater than estimates of local crustal thickness, making it a rare example of a continental mantle earthquake (CME). The focal depth is well resolved from arrival-time inversion (nearest station ∼13 km away) and moment tensor inversion of regional waveforms. Similar to other CMEs in the Intermountain West, there were no obvious aftershocks or foreshocks, and the waveforms were enriched in high-frequency energy. Spectral modeling gives a stress drop of ∼80 MPa and a radiation efficiency of ∼0.08, albeit with large uncertainties. The high stress drop and low radiation efficiency are consistent with a dissipative source process such as thermal runaway. Also similar to previous Intermountain West CMEs, the event occurred along the boundary of the Archean Wyoming craton, where pressure–temperature conditions favor ductile deformation. We hypothesize that edge-driven or regional-scale mantle convection produces increased strain rates near the craton boundary that make either conventional brittle failure or thermal runaway feasible at relatively high pressure–temperature conditions. High conductivity inferred around the edge of the craton may suggest that fluids also contribute to CME occurrence.

Koper, Keith D. [Univ. of Utah, Salt Lake City, UT↗
DOE OSTI2026

Seasonal variation of the surface cross-shelf exchange in the northern South China Sea: a Lagrangian perspective

Previous studies on cross-shelf exchange, predominantly adopted an Eulerian perspective, struggled to identify water sources and pathways. Using a high-resolution regional ocean modeling system (ROMS) and Lagrangian particle tracking, this study systematically investigates the seasonal variation and dynamics of surface cross-shelf exchange in the northern South China Sea (NSCS) from a Lagrangian perspective. Based on daily released 30-day drifter trajectories we identify the key pathways, source regions for surface cross-shelf exchange, revealing pronounced seasonal variability. Results show the surface cross-shelf exchange generally following monsoon-driven Ekman transport. However, local dynamics, especially topographic modulation, can outweigh the expected Ekman-driven transport, producing surface exchange patterns opposite to that predicted from the prevailing winds. Topographic effects vary across different scales. In the coastal waters of western Guangdong during winter, despite downwelling-favorable winds, the modulation of alongshore currents by island topography induces an eastward pressure gradient. This gradient, through geostrophic balance, drives offshore flows opposite to wind-driven Ekman onshore transport. Furthermore, the eastern widened shelf exhibits a distinct seasonal variation of cross-shelf exchange, with strong offshore transport (opposite to the direction of Ekman transport) in winter and exceedingly weak exchange in summer. Analysis of the underlying mechanisms reveals that this winter offshore transport is primarily attributed to geostrophic flows driven by surface pressure gradient that is jointly modulated by the Kuroshio intrusion and local widened shelf topography, and enhanced by cumulative submesoscale processes. In summer, a persistent and strong along-isobath jet acts as a dynamic barrier, effectively suppressing the exchange. These findings highlight the important role of topography in regulating surface material transport, and have important implications for predicting the advection and dispersion of plankton or oil spills over the continental shelf influenced by monsoon.

Hao, Dongliang↗
DOE OSTI2026

Phase Stability and Electrochemical Performance of La-Site-Doped Li6La3Zr0.5Nb0.5Ta0.5Hf0.5O12 High-Entropy Garnets

We investigate La-site substitution in the high-entropy garnet Li6La3Zr0.5Nb0.5Ta0.5Hf0.5O12 (LLZNTH) using Ba2+, Sr2+, and Sm3+ to elucidate how dopant governs phase stability, Li-site distribution, and electrochemical behavior. X-ray diffraction shows that Sr2+ is incorporated homogeneously into the garnet lattice, whereas the larger Ba2+ and smaller Sm3+ ions partially exceed the structural tolerance, generating secondary phases. Nevertheless, the Sm-doped composition (x = 0.05) exhibits the highest room-temperature ionic conductivity (2.7 × 10–4 S cm–1). Neutron powder diffraction reveals that Sm substitution drives a redistribution of Li+ from the tetrahedral 24 d sites into the higher-mobility 96 h positions, enhancing the connectivity of the three-dimensional Li-ion migration network. A Sm-doping series (x = 0.01–0.05) further shows that only sufficiently high Sm levels induce this redistribution, whereas lower concentrations retain Li arrangements similar to the undoped garnet. Critical current density measurements demonstrate that La-site dopants also influence interfacial stability against Li metal, underscoring a trade-off between bulk transport enhancement and mechanical robustness. Collectively, these findings reveal that in high-entropy garnets improved ionic conductivity can originate not only from phase-pure structures but also from targeted modification of the Li sublattice, even when accompanied by secondary phases, offering a compositional design principle for garnet electrolytes.

Li, Chang [Mechanical Engineering, School of Scien↗
DOE OSTI2026

Optical transparency of LiF (100) shock compressed to ∼360 GPa

High-purity [100] lithium fluoride (LiF) is the most widely used optical window in dynamic compression experiments due to its wide bandgap and well-characterized mechanical response. Recent plate-impact experiments established the [100] LiF Hugoniot to ∼230 GPa and demonstrated shock-induced melting onset at 182 GPa with complete melting by 195 GPa; theoretical models predict LiF optical transparency to nearly 900 GPa. To experimentally examine the optical transparency of [100] LiF at higher pressures and in the liquid state, laser-driven shock experiments were performed at peak stresses ranging from 223 to 363 GPa. Optical response was examined by measuring the particle velocity histories at the Kapton/LiF interface using laser interferometry at 532 and 1550 nm wavelengths; in-material particle velocities were obtained using established refractive-index corrections. Continuous photonic Doppler velocimetry fringes were observed across the entire stress range, demonstrating that LiF remains transparent to 1550 nm light throughout the multi-megabar regime investigated. At 532 nm, fringe visibility depended on the reflector coating: aluminum mirrors provided signals to ∼235 GPa, while gold mirrors extended this limit to 270.5 GPa, indicating that the shorter-wavelength response is likely sensitive to experimental configuration rather than to the loss of LiF transparency. Continued optical transparency to at least 360 GPa indicates that shock-melted LiF does not display bandgap closure over the stress range explored. Furthermore, these results provide direct experimental constraints on the high-pressure optical response and establish LiF (100) as a robust optical window material for laser-driven dynamic compression experiments approaching 400 GPa.

Renganathan, P. [Argonne National Laboratory (ANL)↗
DOE OSTI2026

Integrated Simulation of Weld Residual Stress Evolution and Crack Propagation Using XFEM

Nuclear power plant components operate in environments that promote multiple degradation mecha- nisms, several of which involve crack initiation and growth. An ongoing effort in the U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program is developing a general capability within the Multiphysics Object Oriented Simulation Environment (MOOSE) framework for simulating three-dimensional crack growth under a range of driving conditions, including fatigue, stress corrosion cracking (SCC), brittle fracture, and stress-relaxation cracking. This report demonstrates an end-to-end workflow that uses this capability to model weld-residual-stress-driven SCC in the J-groove weld of a pressurized-water reactor control rod drive mechanism penetration in the vessel head. The workflow consists of a thermomechanical welding simulation with temperature-dependent plasticity, followed by cooldown to ambient conditions, and a restart of the simulation using the MOOSE extended finite element method (XFEM) module to propagate a three-dimensional crack through the residual stress field. New welding capabilities were developed to properly initialize newly activated elements in the weld region, and robustness improvements were made to the mesh-based algorithm for defining cutting planes in the 3D XFEM algorithm, allowing it to handle complex crack fronts and stress fields. Together these advances allowed the simulated SCC crack to grow from an initial elliptical flaw in the weld, across the weld, through the tube wall, and almost to the triple point (where the weld, tube, and reactor pressure vessel head intersect) over roughly 36 years of simulated service. These results demonstrate a workflow that can be extended to fully three-dimensional welding simulations and more complex crack interaction problems.

42 - ENGINEERING↗