Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reference correlation”

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.

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Revealing rotational symmetry breaking charge density wave order in the kagome superconductor (Rb, K)⁢V 3 ⁢Sb 5 by ultrafast pump-probe experiments

The recently discovered Kagome superconductor 𝐴⁢V 3 ⁢Sb 5 (where 𝐴 refers to K, Rb, Cs) has stimulated widespread research interest due to its interplay of nontrivial topology and unconventional correlated physics including charge-density waves (CDW) and superconductivity. The essential prerequisite to understanding the microscopic mechanisms of this complex electronic landscape is to unveil the configuration and symmetry of the charge-density wave order. As to now, little consensus has been made on what symmetry is broken. Herein, we clarify the microscopic structure and symmetry breaking of the CDW phase in RbV 3 ⁢Sb 5 and KV 3 ⁢Sb 5 by ultrafast time-resolved reflectivity. Our approach is based on extracting coherent phonon spectra induced by three-dimensional CDW and comparing them to calculated phonon frequencies via density-functional theory. The combination of these experimental results and calculations provides compelling evidence that the CDW structure of both compounds prevailing up to 𝑇 CDW is the 2 × 2 × 2 staggered inverse Star-of-David pattern with interlayer 𝜋 phase shift, in which the sixfold rotational symmetry is broken. Finally, these observations thus corroborate sixfold rotational symmetry breaking throughout the CDW phase of RbV 3 ⁢Sb 5 and KV 3⁢ Sb 5 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Contrast Matching Biopolymers: A SANS-Based Approach to Structural Characterization of Chitosan

Immobilizing enzymes in polysaccharide-based matrices has been shown to improve both their stability and their catalytic efficiency. Among available materials, chitosan was selected for its abundance, biocompatibility, and versatility in a range of applications. While the broader aim of this work is to study chitosan as a matrix for immobilizing carbonic anhydrase (CA) for CO₂ capture, this study focused on determining the contrast match point (CMP) of hydrogenated and deuterated chitosan in D₂O:H₂O mixtures to enable future structural investigations using small-angle neutron scattering (SANS). Previous work successfully synthesized deuterated chitosan from Rhizopus oryzae mycelia and assessed the degree of deuteration using Fourier-transform infrared spectroscopy (FTIR) and nuclear magnetic resonance (NMR). Building on this foundation, the current study determined the CMPs of both H-chitosan and D-chitosan by performing SANS measurements across a range of D₂O:H₂O solvent contrasts. The extracted CMPs for both forms correlated well with their respective calculated scattering length densities (SLDs), validating the experimental approach. These results establish a critical reference for selective contrast matching in future SANS studies of carbonic anhydrase immobilized in chitosan matrices, enabling detailed structural analysis of enzyme–matrix interactions.

Kosgallana, Chathurika J [ORNL] (ORCID:00000002851↗
DOE OSTI2026

Quantitative measurement of pebble fuel anisotropy using generalized ellipsometry

The renewed development of high-temperature gas-cooled reactors has created a need to modernize pebble fuel quality control methods. One important specification in these methods is the coefficient of thermal expansion (CTE) anisotropy, which reflects graphite grain alignment and affects fuel performance. This study investigates the feasibility of using two-modulator generalized ellipsometry microscopy (2-MGEM) as a potential surrogate metric for CTE anisotropy. Graphite matrix compacts with controlled levels of anisotropy were fabricated by varying pre-pressures during forming and were subsequently characterized using both optical dilatometry and 2-MGEM. Measured responses were compared to an isotropic reference pebble produced using an industrial manufacturing process to provide a practical baseline. The resulting measurements revealed a strong quantitative correlation between the CTE anisotropy ratio and the corresponding diattenuation ratio. However, propagated measurement uncertainties in the optical method were found to limit definitive pass/fail discrimination, highlighting the need for improved precision before qualification as a quality control tool. Still, the results establish a quantitative relationship between optical and CTE anisotropy and support further development of 2-MGEM as a potential surrogate method for evaluating graphite matrix anisotropy in pebble fuel manufacturing.

Raftery, Alicia [ORNL] (ORCID:0000000259076427)↗
DOE OSTI2026

Electromagnetic Tomography of Radial Flow in the Quark-Gluon Plasma

We present a novel multimessenger approach to extract the effective radial flow of the quark-gluon plasma (QGP) by jointly analyzing thermal photon and dilepton spectra in heavy-ion collisions. A key feature of this method is that it circumvents the need for a directly unmeasurable reference—the photon temperature in the absence of flow—by establishing, within a calibrated model framework, a stable, approximately linear correlation with the dilepton-inferred temperature. This construction defines an experimentally constructible quantity, 𝑣$^{eff}_{𝑟}$, which reflects early-time collectivity and exhibits a strong correlation with the spacetime-averaged radial velocity of the QGP. Together with previous results linking dilepton slopes to the initial QGP temperature, our work establishes a consistent framework for electromagnetic tomography of the QGP. Here, our framework quantifies the experimental precision target, thereby providing a concrete road map for future measurements at RHIC and the LHC and opening a new avenue to probe the early-time dynamics of hot QCD matter.

Collective flow↗
DOE OSTI2025

Convergence of Cloud Droplet Spectral Relative Dispersion During Entrainment‐Mixing Based on Particle‐Resolved Direct Numerical Simulations

Entrainment-mixing processes critically impact cloud microphysical properties, but their effects on the relative dispersion (d) of cloud droplet size distributions (CDSDs) remain elusive. A direct numerical simulation model is initialized with different CDSDs to fill the gap. These results show that d decreases for broad CDSDs and increases for narrow ones, ultimately converging to approximately 0.5 regardless of initial CDSDs during the evaporation-dominated mixing stage. The supersaturation fluctuation and the shape of CDSDs jointly influence the convergence behavior of d. Further sensitivity tests show that the initial microphysical/dynamical/thermodynamical conditions exert negligible effects on the final converged value of d but affect the convergence rate (k). The k generally increases with increasing droplet number concentration and dissipation rate, and increases with decreasing liquid water content, relative humidity of entrained air, and mixing fraction of cloudy air. A conceptual model with two timescales is proposed; k and the timescales are negatively correlated, meaning that slow mixing and/or evaporation process results in slow convergence of d. In conclusion, this finding provides an important reference for improving understanding and parameterization of d during the entrainment-mixing processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Contextual subspace variational quantum eigensolver calculation of the dissociation curve of molecular nitrogen on a superconducting quantum computer

Abstract We present an experimental demonstration of the Contextual Subspace Variational Quantum Eigensolver on superconducting hardware. Calculating the potential energy curve of molecular nitrogen proves challenging for many conventional quantum chemistry techniques, since static correlation dominates in the dissociation limit. Our quantum simulations retain good agreement with the Full Configuration Interaction energy, outperforming all benchmarked single-reference wavefunction techniques in capturing the bond-breaking appropriately. Moreover, our methodology is competitive with multiconfigurational approaches but at a saving of quantum resource, meaning larger active spaces can be treated for a fixed qubit allowance. To achieve this result, we deploy an error mitigation/suppression strategy comprised of Dynamical Decoupling, Measurement-Error Mitigation and Zero-Noise Extrapolation. Circuit parallelization also provides passive noise-averaging and improves the effective shot yield to reduce the measurement overhead. Furthermore, we introduce a modified adaptive ansatz construction algorithm that incorporates hardware awareness into our variational circuits, minimizing the transpilation cost for the target qubit topology.

Physics↗
DOE OSTI2024

Exploring the Strong Interaction of Three-Body Systems at the LHC

Deuterons are atomic nuclei composed of a neutron and a proton held together by the strong interaction. Unbound ensembles composed of a deuteron and a third nucleon have been investigated in the past using scattering experiments, and they constitute a fundamental reference in nuclear physics to constrain nuclear interactions and the properties of nuclei. In this work, K + -d and p-d femtoscopic correlations measured by the ALICE Collaboration in proton-proton (pp) collisions at $\sqrt{s}$ = 13 TeV at the Large Hadron Collider (LHC) are presented. It is demonstrated that correlations in momentum space between deuterons and kaons or protons allow us to study three-hadron systems at distances comparable with the proton radius. The analysis of the K + -d correlation shows that the relative distances at which deuterons and protons or kaons are produced are around 2 fm. The analysis of the p-d correlation shows that only a full three-body calculation that accounts for the internal structure of the deuteron can explain the data. In particular, the sensitivity of the observable to the short-range part of the interaction is demonstrated. These results indicate that correlations involving light nuclei in pp collisions at the LHC will also provide access to any three-body system in the strange and charm sectors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Systematic Improvement of Quantum Monte Carlo Calculations in Transition Metal Oxides: sCI-Driven Wavefunction Optimization for Reliable Band Gap Prediction

Accurate determination of the electronic properties of correlated oxides remains a significant challenge for computational theory. Traditional Hubbard-corrected density functional theory (DFT+U) frequently encounters limitations in precisely capturing electron correlation, particularly in predicting band gaps. We introduce a systematic methodology to enhance the accuracy of diffusion Monte Carlo (DMC) simulations for both ground and excited states, focusing on LiCoO 2 as a case study. By employing a selected configuration interaction (sCI) approach, we demonstrate the capability to optimize wavefunctions beyond the constraints of single-reference DFT+U trial wavefunctions. Here, we show that the sCI framework enables accurate prediction of band gaps in LiCoO 2 , closely aligning with experimental values and substantially improving traditional computational methods. The study uncovers a nuanced mixed state of t 2g and e g orbitals at the band edges that is not captured by conventional single-reference methods, further elucidating the limitations of PBE+U in describing d-d excitations. Our findings advocate for the adoption of beyond-DFT methodologies, such as sCI, to capture the essential physics of excited-state wavefunctions in strongly correlated materials. The improved accuracy in band gap predictions and the ability to generate more reliable trial wavefunctions for DMC calculations underscore the potential of this approach for broader applications in the study of correlated oxides. This work not only provides a pathway for more accurate simulations of electronic structures in complex materials but also suggests a framework for future investigations of the excited states of other challenging systems.

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

From Oxo to Oxyl to Biradical: Systematic Multireference Calculations of Methane Activation at MOF Nodes

Methane C–H activation at transition-metal sites often involves electronic structures that challenge conventional single-reference electronic structure descriptions. Although Kohn–Sham density functional theory (DFT) is widely used to study catalytic trends, its reliability for reactions involving strongly correlated species remains uncertain. Here we present a systematic multireference investigation of methane activation at metal–organic framework (MOF) node catalysts across the 3d transition-metal series. We introduce an automated workflow for active space selection to enable consistent application of multireference methods, including multiconfiguration pair-density functional theory and n-electron valence state perturbation theory, to these catalytic systems. These calculations show substantial static correlation in the C–H activation reaction step and predict activation barriers that differ from DFT by 30–70 kJ mol–1, with DFT often qualitatively disagreeing in barrier height trends across transition metals. Analysis of multireference wave functions shows that reactivity is governed by the electronic structure of the M–O moiety along a continuum from metal–oxo to oxyl radical and O biradical character. Increased oxygen-centered spin density and weakened M–O bonding are identified as descriptors of catalytic activity which correlate with lower activation barriers.

Wardzala, Jacob↗
DOE OSTI2026

Multiresolution Quantum Chemistry: Nonlinear Response Properties at the Basis Set Limit

We benchmark the accuracy of Dunning correlation-consistent Gaussian basis sets for computing frequencydependent second-order hyperpolarizabilities relevant to second-harmonic generation (SHG), using multiresolution analysis (MRA) as a reference. Basis set errors are analyzed using a unit-sphere representation of the effective hyperpolarizability vector, enabling direct assessment of directional error structure. We introduce a relative RMS total error metric that integrates directional deviations over the unit sphere and complement it with signed projection errors that distinguish over- and underestimation. Unsupervised clustering based on these signed directional metrics reveals four distinct convergence behaviors across a set of 68 molecules. Unitsphere visualizations of representative systems show that basis set errors are often highly anisotropic and localized along specific bond directions, even when global error measures appear small. Doubly augmented basis sets consistently outperform singly augmented ones, and core-polarization functions are required for uniform convergence in second-row systems. Overall, this work demonstrates that directional analysis combined with clustering provides a robust framework for understanding basis set convergence in nonlinear optical response properties.

Basis sets↗
DOE OSTI2025

Superadditivity at large charge

The weak gravity conjecture has been invoked to conjecture that the dimensions of charged operators in a CFT should obey a superadditivity relation (sometimes referred to as convexity). In this paper, we study superadditivity of the operator spectrum in theories expanded about the semi-classical saddle point that dominates correlators of large charge operators. We explore this in two contexts. The first is a model with two scalar fields that carry different charges, at a non-trivial Wilson-Fisher fixed point. A careful analysis of the semi-classics for this two field model demonstrates that ‘quantum’ violations of superadditivity (those not forbidden by the conjecture) persist in the large charge regime. We then turn to study the general properties of CFTs at large charge as bottom-up EFTs. By a trial and error procedure we come up with a seemingly consistent family of examples violating the conjecture. In so doing the presence of a genuine dilaton field appears necessary. On the one hand our result demonstrates that the superadditivity conjecture cannot be proven purely on the basis of a bottom-up analysis. On the other hand, the need for a dilaton, with the corresponding infinite fine tuning, indicates the conjecture-violating EFTs are unlikely to be UV completable.

effective field theories↗
DOE OSTI2025

An Ab Initio Electronic Structure Investigation of the Ground and Excited States of ScH + , YH + , and LaH +

Multireference configuration interaction (MRCI), Davidson-corrected MRCI (MRCI+Q), coupled-cluster singles, doubles, and perturbative triples [CCSD(T)], and frozen-core full configuration interaction (fcFCI) calculations were carried out using large, correlation-consistent basis sets to investigate the excited states of the Sc atom and the spin–free and spin–orbit coupled potential energy profiles, energetics, spectroscopic constants, and electron populations of low-lying states of MH + (M = Sc, Y, La). The core electron correlation effects, complete basis set effects, and spin–orbit coupling effects were also evaluated. The first four electronic states of all MH + are 1 2 Δ, 1 2 Σ + , 1 2 Π, and 2 2 Σ + with 1σ 2 1δ 1 , 1σ 2 2σ 1 , 1σ 2 1π 1 , and 1σ 2 3σ 1 single-reference electron configurations, respectively. These states of MH + can be represented by the M 2+ H – ionic structure. The ground states of ScH + , YH + , and LaH + are 1 2 Δ 3/2 , 1 2 Σ + 1/2 , and 1 2 Δ 3/2 with 55.45, 60.54, and 62.34 kcal/mol bond energies, respectively. The core electron correlation was found to be vital for gaining accurate predictions on the ground and excited state properties of MH + . The spin–orbit coupling effects are minor for ScH + but become substantial moving to YH + and LaH + . Overall, the results of this work are in good agreement with the limited set of experimental findings of MH + available in the literature and will be of use for future investigations. Furthermore, the theoretical approaches, findings, and trends reported here are expected to aid studies of similar species.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

A-Site Cation Chemistry in Halide Perovskites

Metal halide perovskites are an important class of semiconductors now being implemented as photovoltaic absorbers and explored for light emission, among other device applications. The semiconducting properties of halide perovskites are deeply intertwined with their composition and structure. Specifically the symmetry, tilting, and distortions of the metal halide octahedra impact the band structure and other optoelectronic properties. In this review, we examine the various compositions of monovalent A-site cations in threedimensional (3D) halide perovskites AMX 3 (M = divalent metal; X = halide). We focus on how the A-site cation templates the inorganic metal-halide perovskite framework, resulting in changes in the crystal structure symmetry, as well as M–X bonding parameters, summarized in a comprehensive table of AMX 3 structures. The A-site cation motion, effects of alloying, and 2D Ruddlesden–Popper perovskite structures with unique A-site cations are further overviewed. Correlations are shown between these A-site cation dominated structural parameters and the resulting optoelectronic properties such as band gap. This review should serve as a reference for the A-site cation structural chemistry of metal halide perovskites and inspire continued research into less explored metal halide perovskite compositions and structures.

14 SOLAR ENERGY↗
DOE OSTI2025

Assessing the Impact of Large Removable Beryllium Reflector Experiments on HFIR Performance and Safety Metrics

This study investigated the effect of various removable beryllium (RB) experiment configurations on High Flux Isotope Reactor (HFIR) metrics to address increasing interest in these facilities for materials and fuels irradiation research. The RB reflector contains eight large and four small irradiation experiment facilities that offer excellent neutron flux conditions to perform fission and fusion reactor materials and fuels irradiation research. This work aims to outline acceptable RB configurations based on safety, performance, and programmatic metrics. This study assessed experiment effects on reactivity, cycle length, fission rate density distributions, and neutron flux distributions using the Shift, HFIRCON, and SCALE ORIGEN codes. Initial evaluations focused on generic experiment materials including aluminum plugs, stainless steel plugs, molybdenum plugs, and aluminum plugs with gadolinium shields. Subsequent analyses of MiniFuel experiments were performed to evaluate a heterogenous experiment, consisting of a more complex geometry and bearing several materials, and to test the correlations developed with the generic materials on a real experiment. Furthermore, the effects of neutron poison concentrations in the standard beryllium plugs were evaluated. When assuming a reference RB configuration with eight large fresh beryllium plugs, perturbed configurations with three aluminum plugs, one stainless steel plug, one molybdenum plug, one gadolinium-shielded aluminum plug, and one MiniFuel experiment resulted in a cycle length reduction of less than 1.3 days, the current threshold before requiring additional approvals. Configurations with five aluminum plugs, one stainless steel plug, one molybdenum plug, one gadolinium-shielded aluminum plug, and two MiniFuel experiments meet the 1.3 day limit if irradiated beryllium plugs are considered. Fuel element fission rate density distributions remained within safety limits, with maximum local increases under 9%. Neutron flux calculations revealed large thermal flux depressions inside and around the perturbed RB facilities, while epithermal and fast neutron fluxes increased because of reduced neutron moderation by the perturbed materials. These findings provide valuable guidance for optimizing RB configurations to balance safety and performance at HFIR.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Noise limits for dc SQUID readout of high- Q resonators below 300 MHz

We present the limits on noise for the readout of cryogenic high-Q resonators using dc Superconducting Quantum Interference Devices (SQUIDs) below 300 MHz. This analysis uses realized first-stage SQUIDs (previously published), whose performance is well described by Tesche–Clarke (TC) theory, coupled directly to the resonators. We also present data from a prototype second-stage dc SQUID array designed to couple to this first-stage SQUID as a follow-on amplifier with high system bandwidth. This analysis is the first full consideration of dc SQUID noise performance referred to a high-Q resonator over this frequency range and is presented relative to the standard quantum limit. We include imprecision, backaction, and backaction–imprecision noise correlations from TC theory, the noise contributed by the second-stage SQUIDs, wiring, and preamplifiers, and optimizations for both on-resonance measurements and off-resonance scan sensitivity. This architecture has modern relevance due to the increased interest in axion searches and the requirements of the DMRadio-m 3 axion search, which uses dc SQUIDs in this frequency range.

Bosons↗
DOE OSTI2024

Optimal Design Approaches for Cost-Effective Manufacturing & Deployment of Chemical Process Families with Economies of Numbers

This work builds on our optimization formulation for process family design and extends it to explicitly include the benefits of economies of numbers. Economies of numbers (sometimes referred to as economies of learning) is a well-documented cost saving phenomenon. It characterizes the manufacturing cost savings due to standardization; in particular, it is capturing the correlation between cost reduction and the number of times a particular product has been manufactured. Following an approach similar to that in Gazzaneo et al. (2022), we develop a costing expression that captures material costs and manufacturing costs as a function of the number of unit modules produced. If the platform has a small number of unit module designs, we will be manufacturing a large number of each of these designs and gaining increased benefits from economies of numbers. However, increasing the number of unit module designs in the platform gives each process variant more choices to consider (at the cost of reducing economies of numbers). The optimization formulation in Stinchfield et al. (2023) pre-specified the number of unit module designs to be included in the platform. Here, by including the economies of numbers explicitly, we allow the mathematical programming formulation to determine the optimal number of unit module designs to include in the platform. We demonstrate this approach on multiple case studies, including MEA-based carbon capture and water desalination.

Stinchfield, Georgia↗
DOE OSTI2026

DEM Modeling and Validation of Pebble Bed Packing Using Chrono::GPU

Accurate prediction of pebble packing structure is important for pebble bed reactors because the spatial distribution of void fraction directly affects coolant flow, pressure drop, heat transfer, and neutronic behavior. However, experimentally validated DEM studies that directly evaluate local void-fraction structure in reactor-relevant pebble beds remain limited. In this work, the pebble bed experiment conducted at Missouri University of Science and Technology is simulated using the graphics processing unit (GPU)-based discrete element method (DEM) code Chrono::GPU. The study focuses on evaluating the ability of Chrono::GPU to reproduce the packing arrangement and void-fraction distribution of a randomly packed spherical pebble bed. The DEM results are first verified against established radial void-fraction correlations, including the Mueller and Vortmeyer-Schuster models, to assess the predicted bulk porosity, near-wall behavior, and oscillatory packing structure. The simulation is then verified against reference DEM data and validated against gamma-ray computed tomography (CT) experimental data at three axial locations. The Chrono::GPU results reproduce the main features of the experimental packing, including the high void fraction near the wall, the first near-wall trough, and the damped oscillatory radial profile caused by wall-induced ordering. Quantitative comparison with DEM data and the CT-based radial profiles shows good agreement, with mean absolute errors on the order of 0.07 and root-mean-square errors below 0.09 for the averaged profiles. These results demonstrate that Chrono::GPU can accurately capture the void-fraction structure of spherical pebble beds and provides a reliable DEM framework for future pebble bed reactor packing, recycling, and thermal-hydraulic studies.

97 - MATHEMATICS AND COMPUTING↗