Search NASA⌕ Search

SEARCH · Search NASA

Results for “Statistical physics”

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.

531 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Engineering 2D Square Lattice Hubbard Models in 90° Twisted GeX/SnX (X =S, Se) Moiré Superlattices

Because of the large-period superlattices emerging in moiré two-dimensional (2D) materials, electronic states in such systems exhibit low energy flat bands that can be used to simulate strongly correlated physics in a highly tunable setup. While many investigations have thus far focused on moiré flat bands and emergent correlated electron physics in triangular, honeycomb, and quasi-one-dimensional lattices, tunable moiré realizations of square lattices subject to strong correlations remain elusive. Here, in this work, we propose a feasible scheme to construct moiré square lattice systems by twisting two or more layers of 2D materials in a rectangular lattice by 90°. We demonstrate the concept with twisted GeX/SnX (X =S, Se) moiré superlattices and calculate their electronic structures from first principles. We show that the lowest conduction flat band in these systems can be described by a square lattice Hubbard model with parameters which can be controlled by varying the choice of host materials, number of layers, and external electric fields. In particular, twisted double bilayer GeSe realizes a square lattice Hubbard model with strong frustration due to the next-nearest-neighbor hopping that could host unconventional superconductivity, in close analogy to the Hubbard model for copper-oxygen planes of cuprate high-temperature superconductors. The presented scheme uses 90° twisted 2D materials with rectangular unit cells as a promising platform for realizing the physical phenomena of square lattice Hubbard models, establishing a new route for studying its rich phase diagram of magnetism, charge order, and unconventional superconductivity in a highly tunable setting.

2-dimensional systems↗
DOE OSTI2026

Radiative divertor detachment and impurity transport with nitrogen and neon seeding in KSTAR H-mode plasmas

Achieving core-edge compatible divertor detachment is a critical requirement for stable operation in future fusion devices. This study compares nitrogen and neon seeding in KSTAR H-mode plasmas with carbon walls, combining experiments and SOLPS-ITER modelling to evaluate their radiative dissipation and core-edge compatibility. Experimentally, N seeding achieved stronger divertor detachment, with larger reductions in target particle and heat fluxes, a higher divertor radiation fraction, and a lower core radiation fraction than Ne. In contrast, Ne seeding triggered a significant rise in core radiation followed by H–L back transitions, limiting the maximum total radiated power fraction to roughly half that of N. SOLPS-ITER simulations reproduced the experimental trends and revealed that the better core-edge compatibility of N arises from its higher divertor retention in addition to its higher cooling factor. The relative positions of the stagnation points of impurity poloidal velocity and ionization sources did not explain the different divertor compression. Instead, in the present modelling, the higher impurity parallel particle flux, resulting from the higher parallel impurity velocity, explains the stronger nitrogen impurity compression in the divertor region. The parallel temperature distribution with N was more favourable for achieving higher impurity parallel velocity than with Ne, because the impurity velocity is governed by modifications of the main ion flow due to friction and thermal forces, both of which strongly depend on the temperature. Ultimately, this behaviour is attributed to the strongly divertor-localized radiation of N. These results demonstrate that N is more effective than Ne in achieving radiative divertor detachment while maintaining low core contamination in KSTAR, consistent with observations in other present tokamaks.

KSTAR↗
DOE OSTI2026

SODAs: sparse optimization for the discovery of differential and algebraic equations

Differential-algebraic equations (DAEs) integrate ordinary differential equations (ODEs) with algebraic constraints, providing a fundamental framework for developing models of dynamical systems characterized by time-scale separation, conservation laws and physical constraints. While sparse optimization has revolutionized model development by allowing data-driven discovery of parsimonious models from a library of possible equations, existing approaches for dynamical systems assume DAEs can be reduced to ODEs by eliminating variables before model discovery. This assumption limits the applicability of such methods for DAE systems with unknown constraints and time scales. We introduce sparse optimization for differential-algebraic systems (SODAs), a data-driven method for the identification of DAEs in their explicit form. By discovering the algebraic and dynamic components sequentially without prior identification of the algebraic variables, this approach leads to a sequence of convex optimization problems. It has the advantage of discovering interpretable models that preserve the structure of the underlying physical system. To this end, SODAs improves since SODAs is singular numerical stability when handling high correlations between library terms, caused by near-perfect algebraic relationships, by iteratively refining the conditioning of the candidate library. We demonstrate the performance of our method on biological, mechanical and electrical systems, showcasing its robustness to noise in both simulated time series and real-time experimental data.

DAE↗
DOE OSTI2026

A transport-based framework for solidification cracking in Ni-based superalloys processed by laser powder bed fusion

Solidification cracking remains a persistent barrier to laser powder bed fusion (LPBF) processing of solid-solution-strengthened (SSS) Ni-based superalloys and is commonly attributed to carbide formation, liquation-type failure, or elemental segregation. In this work, the cracking behavior of three SSS Ni-based superalloys—Inconel 625 (625), Inconel 617 (617), and Haynes 230 (230)—is examined under comparable LPBF conditions to evaluate the origins of their cracking susceptibility. Carbides were present in both 625 and 230, yet cracking behavior differed significantly. MC-type carbides in 625 remain discrete and preserve liquid connectivity, whereas Cr-rich M23C6 carbides in 230 form at cellular triple junctions that bottleneck the interdendritic liquid network, influencing cracking through liquid connectivity rather than as an independent cause. Alloy 617 exhibited anomalous segregation without resolvable secondary phases yet cracked mildly, indicating that extensive carbide formation is not a prerequisite for cracking. To rationalize these observations, a transport-limited framework is proposed in which cracking occurs when terminal-liquid feeding cannot accommodate solidification-induced strain; among the mechanisms evaluated, this framework is most consistent with the physical constraints imposed by LPBF solidification. Transport-based crack-susceptibility indexes incorporating permeability, viscosity, and solidification kinetics distinguish the alloy hierarchy; permeability ratios Krat increased from approximately 7 (625) to 11.5 (617) and 14.5 (230), with corresponding increases in the μ-weighted mushy zone risk index (CSIMZRM). These results are consistent with transport-limited liquid feeding as a rate-limiting contributor to solidification cracking in LPBF-processed SSS Ni-based superalloys, providing a physically grounded basis for alloy and process design.

Hyer, Holden [ORNL] (ORCID:0000000343915561)↗
DOE OSTI2026

Superconducting phase diagram of multilayer square-planar nickelates

The discovery of superconductivity in square-planar nickelates has offered a rich materials platform to explore the origins of high-temperature superconductivity. However, experimental investigations have largely been limited to the infinite-layer R NiO 2 ( R , rare earth) nickelates. For this work we constructed a phase diagram of multilayer square-planar Nd n+1 Ni n O 2n+2 compounds and found signatures of superconductivity for dimensionality n = 4 to 8. Upon decreasing n , the superconducting anisotropy evolves owing to 4ƒ electron effects, and electronic structure characteristics approach cuprate-like behavior. Magnetic fluctuations persist from within the superconducting regime and into the overdoped, nonsuperconducting regime. The superconducting regime overlaps with that of chemically doped infinite-layer nickelates, demonstrating underlying commonalities as well as differences across varying structural realizations of square-planar nickelates. Our work establishes this layered template for creating new nickel-based superconductors.

36 MATERIALS SCIENCE↗
DOE OSTI1969

Application of Nuclear Technology to Art Identification Problems: First Annual Report

Throughout the centuries, as man has become increasingly affluent, his interest in antiquities and objects of art, and the price he has been willing to pay for them, has increased. Inevitably, one result has been to attract the unscrupulous to the production and sale of counterfeit paintings, sculptures, and other works of art. The skill of forgers varies but often is highly sophisticated. However, dealers, museum curators, and other experts have also become more skillful in the detection of forgeries. Scientific tools of increasing sensitivity and sophistication have gradually been applied to the examination of the materials of art and archaeology. Such tools frequently support and render more reliable the judgment of the experts. The quality of forgeries may improve further as their makers become acquainted with the new methods of detection and in turn learn to circumvent or confound the methods. Therefore, the development of still more advanced methods of examination has great utility in the art world. This is particularly true if the ultimate effect is to make circumvention of the methods of examination so costly that the economic incentives for producing forgeries may be substantially reduced, if not eliminated.

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

Microstructure prediction for Ti-22Al-25Nb in laser powder bed fusion

This work presents a physics-informed framework for predicting solidification morphology and defect susceptibility in additively manufactured Ti–22Al–25Nb across a broad processing space. The framework integrates solidification microstructure selection (SMS) analysis with a single-track defect-based printability map to establish a unified methodology linking processing parameters to both interfacial morphology and manufacturability. Thermal gradients G and solidification rates R are first computed using the Thermo-Calc Additive Manufacturing (TC-AM) module, a finite-interface-dissipation (FID) phase-field (PF) model coupled with CALPHAD method is then employed to systematically distinguish planar and dendritic regimes as functions of $G$ and $R$. By superimposing the printability map onto the morphology projections, a comprehensive process–structure framework is obtained. Across most processing conditions, the predicted microstructure is predominantly dendritic, while planar growth emerges only under selected laser power $P$ and scan speed $v$ combinations. In addition to morphology classification, the framework quantifies the dendritic area fraction and introduces a width-based morphology descriptor to characterize the spatial extent of planar/dendritic regions within the melt pool. It provides mechanistic insight into the interplay between solidification physics and defect formation, offering practical guidance for parameter selection and microstructural control in Ti–22Al–25Nb additive manufacturing (AM).

36 MATERIALS SCIENCE↗
DOE OSTI2026

Combined dark matter search towards dwarf spheroidal galaxies with Fermi -LAT, HAWC, H.E.S.S., MAGIC, and VERITAS

Dwarf spheroidal galaxies (dSphs) are excellent targets for indirect dark matter (DM) searches using gamma-ray telescopes because they are thought to have high DM content and a low astrophysical background. The sensitivity of these searches is improved by combining the observations of dSphs made by different gamma-ray telescopes. We present the results of a combined search by the most sensitive currently operating gamma-ray telescopes, namely: the satellite-borne Fermi -LAT telescope; the ground-based imaging atmospheric Cherenkov telescope arrays H.E.S.S., MAGIC, and VERITAS; and the HAWC water Cherenkov detector. Individual datasets were analyzed using a common statistical approach. Results were subsequently combined via a global joint likelihood analysis. We obtain constraints on the velocity-weighted cross section 〈σv〉 for DM self-annihilation as a function of the DM particle mass. This five-instrument combination allows the derivation of up to 2-3 times more constraining upper limits on 〈σv〉 than the individual results over a wide mass range spanning from 5 GeV to 100 TeV. Depending on the DM content modeling, the 95% confidence level observed limits reach 1.5×10 -24 cm 3 s -1 and 3.2×10 -25 cm 3 s -1 , respectively, in the τ + τ - annihilation channel for a DM mass of 2 TeV.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Imaging of a van der Waals spin-orbit torque system using spin ensembles in hBN

Recently, optically active spin defects embedded in two-dimensional (2D) van der Waals (vdW) crystals have emerged as a transformative quantum sensing platform to explore cutting-edge materials science. Taking advantage of excellent solid-state integrability, this new class of spin defects can be readily arranged in nanoscale proximity to target materials, showing great promise for realizing in-situ quantum sensing of microscopic spin and charge behaviors in vdW heterostructures. Here we report hexagonal boron nitride-based quantum imaging of field-free deterministic magnetic switching and electric current distributions in an all-vdW spin-orbit torque (SOT) system. By visualizing variations of nanoscale magnetic stray field profile of room-temperature 2D magnet Fe 3 GaTe 2 under different SOT conditions, we show how the magnetic switching evolves from deterministic to stochastic behavior due to the interplay between spin orientations, anisotropy and Joule heating. Micromagnetic simulations rationalize our results well, revealing the role of field-like SOT in inhibiting thermal fluctuation driven stochastic switching and chaotic multi-domain competition. This understanding, which is otherwise difficult to access by conventional transport measurements, offers valuable insights into material design, testing, and performance evaluation of next-generation vdW spintronic devices.

Imaging techniques↗
DOE OSTI2026

Biochar–polymer composites for 3D printing: a review

Biochar, a bio-based co-product of biofuel production via thermochemical conversion, holds potential as a filler for polymer composites to reduce costs, improve thermomechanical properties, and aid in environmental remediation. 3D-printed biochar composites have received growing interest over the past few years but have experienced difficulties such as poor layer adhesion and nozzle clogging. Currently, no literature review examines 3D-printed biochar composites and related biochar properties in-depth. This work summarizes and discusses recent studies on 3D-printed polymer and biochar composites and examines their mechanical, thermal, and additional properties that result from each study. Technical challenges in printability, such as nozzle clogging from particle size and biochar aggregation, are also discussed. Furthermore, this work discusses the variability of biochar properties resulting from the pyrolysis conditions and feedstock choice in relation to potential 3D printing outcomes. In particular, several studies reported that high lignin feedstocks could be candidates for 3D printing. The post-processing approaches of the biochar via physical and chemical methods are also introduced. Ball milling appears to hold the most promise for physical treatments due to its tunability of particle size, surface area, and functional groups, while chemical treatments with acids or alkalis are used to tailor biochar porosity and wettability. Overall, it was determined that future research needs to be done relating biochar production and post-processing methods to resulting 3D printing parameters as the number of studies is limited.

Day, Rachel [Auburn University, Auburn, AL]↗
DOE OSTI2026

String-Breaking Dynamics in Quantum Adiabatic and Diabatic Processes

Confinement prohibits isolation of color charges, e.g., quarks, in nature via a process called string breaking : the separation of two charges results in an increase in the energy of a color flux, visualized as a string, connecting those charges. Eventually, creating additional charges is energetically favored, hence breaking the string. Such a phenomenon can be probed in simpler models, including quantum spin chains, enabling enhanced understanding of string-breaking dynamics. A challenging task is to understand how string breaking occurs as time elapses, in an out-of-equilibrium setting. This work establishes the phenomenology of dynamical string breaking induced by a gradual increase of string tension over time. It, thus, goes beyond instantaneous quench processes and enables tracking the real-time evolution of strings in a more controlled setting. We focus on domain-wall confinement in a family of quantum Ising chains. Our results indicate that, for sufficiently short strings and slow evolution, string breaking can be described by the transition dynamics of a two-state quantum system akin to a Landau-Zener process. For longer strings, a more intricate spatiotemporal pattern emerges: the string breaks by forming a superposition of bubbles (domains of flipped spins of varying sizes), which involve highly excited states. We finally demonstrate that string breaking driven only by quantum fluctuations can be realized in the presence of sufficiently long-ranged interactions. This work holds immediate relevance for studying string breaking in quantum-simulation experiments.

Ising model↗
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 OSTI2026

Purcell enhancement of directional edge photocurrent in a van der Waals self-cavity

Cavities provide a means to manipulate the optical and electronic responses of quantum materials by selectively enhancing light-matter interaction at specific frequencies and momenta. While cavities typically involve external structures, exfoliated flakes of van der Waals (vdW) materials can form intrinsic self-cavities due to their small finite dimensions, confining electromagnetic fields into plasmonic cavity modes, characterized by standing-wave current distributions. While cavity-enhanced phenomena are well-studied at optical frequencies, the impact of self-cavities on nonlinear electronic responses—such as directional photocurrent—remains largely unexplored, particularly in the terahertz regime, critical for emerging ultrafast optoelectronic technologies. Here, we report a self-cavity-induced Purcell enhancement of directional photocurrents in the vdW semimetal WTe 2 . Using ultrafast optoelectronic circuitry, we measured coherent near-field THz emission resulting from nonlinear photocurrents excited at the sample edges. We observed enhanced emission at finite frequencies, tunable via excitation fluence and sample geometry, which we attribute to plasmonic interference effects controlled by the cavity boundaries. We developed an analytical theory that captures the cavity resonance conditions and spectral response across multiple devices. Our findings establish WTe 2 as a bias-free, geometry-tunable THz emitter and demonstrate the potential of self-cavity engineering for controlling nonlinear, nonequilibrium dynamics in quantum materials.

condensed-matter physics↗
DOE OSTI2024

Observational evidence for Criegee intermediate oligomerization reactions relevant to aerosol formation in the troposphere

Criegee intermediates are reactive intermediates that are implicated in transforming the composition of Earth’s troposphere and in the formation of secondary organic aerosol, impacting Earth’s radiation balance, air quality and human health. Yet, direct identification of their signatures in the field remains elusive. Here, from particulate and gas-phase mass-spectrometric measurements in the Amazon rainforest, we identify sequences of masses consistent with the expected signatures of oligomerization of the CH 2 OO Criegee intermediate, a process implicated in ozonolysis-driven aerosol formation. In this study, we assess the potential contributions of oligomerization through laboratory ozonolysis experiments, direct kinetic studies of Criegee intermediate reactions, and high-level theoretical calculations. Global atmospheric models built on these kinetics results indicate that Criegee intermediate chemistry may play a larger role in altering the composition of Earth’s troposphere than is captured in current atmospheric models, especially in areas of high humidity. However, the models still capture only a relatively small fraction of the observed signatures, suggesting considerable underestimates of Criegee intermediate concentrations and reactivity and/or the dominance of other, presently uncharacterized, oxidation mechanisms. Resolving the remaining uncertainties in emission inventories and the effects of atmospheric water vapour on key chemical reactions will be required to definitively assess the role of Criegee intermediate oligomerization reactions.

58 GEOSCIENCES↗
DOE OSTI2026

The gas-phase mass–metallicity relation of dwarf galaxies across large-scale environments using the CAVITY parent sample

Context. The gas-phase mass–metallicity relation (MZR) of galaxies shows a noticeable break in slope and an increased scatter at low stellar masses, suggesting that the physical processes governing chemical enrichment differ between dwarf and high-mass systems. Dwarf galaxies, in particular, are highly susceptible to both internal and environmental mechanisms due to their shallow potential wells. Aims. The primary aim of this work is to assess whether a single, universal MZR can describe dwarf galaxies across diverse large-scale environments, or whether systematic environmental variations emerge. To probe these, we examine the MZR and star formation rate (SFR) of dwarf galaxies with stellar masses in the range of 8.9 < log(M ★ /M ⊙ ) < 9.5. Methods. Using optical spectra from the Sloan Digital Sky Survey, we measured the fluxes of key emission lines via the pyPipe3D full spectral fitting pipeline. Aperture-corrected fluxes, along with multiple metallicity indicators and calibrations, were used to derive the MZR and the SFR for 353, 311, and 22 dwarf galaxies located in voids, filaments, and clusters, respectively. Results. We find a systematic variation in the MZR slope, steeper in voids (0.28 ± 0.03) and progressively flatter in clusters (0.17 ± 0.08), indicating a dependence of the MZR on the large-scale environment in this mass regime. When galaxies are separated by local density, no significant differences are observed between isolated and non-isolated dwarfs in voids. Isolated dwarf galaxies in filaments also exhibit properties similar to those of their counterparts in voids. However, non-isolated filament galaxies exhibit similar MZR slopes comparable to those of cluster dwarfs and flatter slopes than their counterparts in voids. Conclusions. We report both large- and local-scale environmental dependencies in the gas-phase metallicity and in the slope of the MZR for dwarf galaxies. Consistent with the general consensus on the pre-processing of galaxies in filaments, our results indicate that the influence of the local environment becomes increasingly significant within the filamentary regions of the cosmic web, affecting the chemical enrichment and star formation activity of low-mass systems. These findings further suggest that a portion of the scatter commonly observed in the MZR of dwarf galaxies arises from environmental effects.

Bidaran, Bahar [Dpto. de Física Teórica y del Cosm↗
DOE OSTI2026

Single-Nucleon Transfer on Unstable 59 Cu Probes the NiCu Cycle in Astrophysical X-Ray Bursts

Recent models of the rapid proton (𝑟⁢𝑝) capture process indicate that a competition between the 59 Cu⁢(𝑝,𝛾) ⁢60 Zn and 59 Cu⁢(𝑝,𝛼) ⁢56 Ni reactions may result in the formation of a nickel-copper (NiCu) cycle that traps the flux of material between 56 Ni and 60 Zn . Here, in this work, we report the identification of 15 proton-unbound levels in 60 Zn , populated via 59 Cu⁢(𝑑,𝑛) transfer, which govern the rate of the 59 Cu⁢(𝑝,𝛾)⁢ 60 Zn reaction in XRBs. Precise excitation energies for levels in 60 Zn were obtained from observed 𝛾 decays, and spectroscopic factors were determined from angle-integrated cross sections. Incorporating these results into stellar-model calculations, we find that with experimentally constrained uncertainties a NiCu cycle in XRBs is indeed possible, though we limit its branching strength to less than 38%. While modest, such a branching has significant impact on the light curve, motivating further studies of the relevant rates. Our calculations also indicate that a significant NiCu cycle leads to an increase in the amount of odd-𝐴 nuclei in the burst ashes, which may affect Urca cooling processes in neutron star crusts.

direct reactions↗
DOE OSTI2024

Energy conversion and transport in molecular-scale junctions

Molecular-scale junctions (MSJs) have been considered the ideal testbed for probing physical and chemical processes at the molecular scale. Due to nanometric confinement, charge and energy transport in MSJs are governed by quantum mechanically dictated energy profiles, which can be tuned chemically or physically with atomic precision, offering rich possibilities beyond conventional semiconductor devices. While charge transport in MSJs has been extensively studied over the past two decades, understanding energy conversion and transport in MSJs has only become experimentally attainable in recent years. As demonstrated recently, by tuning the quantum interplay between the electrodes, the molecular core, and the contact interfaces, energy processes can be manipulated to achieve desired functionalities, opening new avenues for molecular electronics, energy harvesting, and sensing applications. This Review provides a comprehensive overview and critical analysis of various forms of energy conversion and transport processes in MSJs and their associated applications. We elaborate on energy-related processes mediated by the interaction between the core molecular structure in MSJs and different external stimuli, such as light, heat, electric field, magnetic field, force, and other environmental cues. Key topics covered include photovoltaics, electroluminescence, thermoelectricity, heat conduction, catalysis, spin-mediated phenomena, and vibrational effects. Furthermore, the review concludes with a discussion of existing challenges and future opportunities, aiming to facilitate in-depth future investigation of promising experimental platforms, molecular design principles, control strategies, and new application scenarios.

Charge transport↗
DOE OSTI2026

FY26 Progress on Demonstration of a Multiphysics Steady State Capability for Modeling Core Radial Expansion in SFRs

Under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism in liquid metal-cooled fast reactors and involves multiphysics effects including radiation transport, fluid flow, heat transfer, and mechanical response to temperature and flux gradients. This report summarizes recent progress on developing a multiphysics, MOOSE-based workflow to predict core bowing and associated reactivity feedback. Significant new capabilities in the reactor physics code Griffin - sodium backfill and pin power reconstruction for deformed geometries - were applied in this effort. This year’s work included verification, code comparisons, sensitivity studies, and coupled demonstrations that advance the state of MOOSE-based core bowing workflow. Griffin’s sodium backfill capability was verified by demonstrating that its automated treatment of geometry expansion and material-density updates reproduces manual calculations exactly, confirming solid mass conservation and proper coolant backfilling in expanded geometries. Reconstructed pin powers were compared for Griffin’s ductheterogeneous and ring-heterogeneous treatments in single-, seven-, and nineteen-assembly cases, with best agreement observed in lower-leakage configurations and the duct-heterogeneous approach offering substantially lower computational cost. Thermal-hydraulic sensitivity sensitivities showed that MOOSE SCM, SAM, and CFD are expected to produce similar deformation predictions despite variances in their temperature predictions, and that explicit treatment of inter-assembly flow becomes increasingly important as gap flow rate increases. Finally, coupled demonstrations on small multi-assembly configurations using Griffin, MOOSE Solid Mechanics, MOOSE SCM, and Heat Conduction produced physically consistent reactivity feedback from thermal expansion and bowing. The coupled demonstrations simulated grid plate expansion as well as resultant core bowing at full power conditions. Simplifications were made in current workflow, namely the assumption of instantaneous full power conditions following hot zero power, and pre-expanding the Griffin geometry axially due to lack of an axial fuel pin expansion model and temperature feedback to Griffin.

Wozniak, Nicholas↗