SEARCH · Engineering Papers

Results for “simulator study”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 757 records · Page 42

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Simulation of radon-222 with the GEOS-Chem global model: Emissions, seasonality, and convective transport

Radon-222 (222Rn) is a short-lived radioactive gas naturally emitted from land surface, and has long been used to assess convective transport in atmospheric models. In this study, we simulate 222Rn using the GEOS-Chem chemical transport model with aims to improve our understanding of 222Rn emissions and surface concentration seasonality, and to characterize convective transport associated with two Goddard Earth Observing System (GEOS) meteorological products, MERRA and GEOS-FP. We evaluate four available 222Rn emission scenarios by comparing model results with surface observations at 51 global surface sites. The default emission scenario in GEOS-Chem yields a moderate agreement with global surface observations (< 70% data within a factor of 2) and a large underestimate of wintertime surface 222Rn concentrations at Northern Hemisphere mid- and high-latitudes due to an oversimplified formulation of 222Rn emission fluxes (1 atom cm-2 s-1 over land with a reduction by a factor of 3 under freezing conditions). We compose a new global 222Rn emission scenario based on Zhang et al. (2011) and show its potential to improve simulated surface 222Rn concentrations and seasonality. The regional components of this emission scenario include spatially and temporally varying emission fluxes derived from previous measurements of soil radium content and soil exhalation models, which are key factors to 222Rn emission flux rates. However, large model underestimates of surface 222Rn concentrations still exist in Asia, suggesting unusually high regional 222Rn emissions. We propose a conservative up-scaling factor of 1.2 for 222Rn emission fluxes in China, as also constrained by the observed deposition fluxes of 210Pb (decay daughter of 222Rn). With this modification, the model shows better agreement with the observations in Europe and North America (>80% data within a factor of 2), and reasonable agreement in Asia (close to 70%). Further constraints on 222Rn emissions would require additional observations of surface 222Rn concentrations and emission fluxes in central U.S., Canada, Africa, and Asia. We also compare and assess convective transport in model simulations driven by MERRA and GEOS-FP using observed 222Rn vertical profiles during northern mid-latitude summertime and from two short-term airborne campaigns. While the simulations with both GEOS products are able to capture the observed vertical gradient of 222Rn concentrations in the lower troposphere (0-4 km), neither correctly represents the level of convective detrainment, resulting in biases in the middle and upper troposphere. Compared to GEOS-FP, MERRA leads to stronger convective transport of 222Rn, which is partially compensated by its weaker large-scale vertical advection, resulting in similar global vertical distributions of 222Rn concentrations between the two simulations.

Zhang, Bo↗

DOE OSTI2021

A discrete dislocation dynamics study of precipitate bypass mechanisms in nickel-based superalloys

Order strengthening in nickel-based superalloys is associated with the extra stress required for dislocations to bypass the $\gamma'$ precipitates distributed in the $\gamma$ matrix. Depending on the operating conditions and microstructure, a rich variety of bypass mechanism has been identified, with various shearing and Orowan looping processes gradually giving way to climb bypass as the operating conditions change from the low/intermediate temperatures and high stress regime, to the high temperature and low stress regime. When anti phase boundary (APB) shearing and Orowan looping mechanisms operate, the classical picture is that, at for a given volume fraction, the bypass mechanism changes from shearing to looping with increased particle size and within a broad coexistence size window. Another possibility, which is supported by indirect experimental evidence, is that a third “hybrid” transition mechanism may operate. Here, in this paper we use discrete dislocation dynamics (DDD) simulations to study dislocation bypass mechanisms in Ni-based superalloys. We develop a new method to compute generalized stacking fault forces in DDD simulations, based on a concept borrowed from complex analysis and known as the winding number of a closed curve about a point. We use this method to study the mechanisms of bypass of a square lattice of spherical $\gamma'$ precipitates by $a/2\langle110\rangle$ {111} edge dislocations, as a function of the precipitates volume fraction and size. We show that not only the hybrid mechanism is possible, but also that it operates as the transition mechanism between the shearing and looping regimes over a wide range of precipitates volume fraction and radii. Based on our simulation results, we propose a simple model for the strength of this mechanism. We also consider the effects of a $\gamma$/$\gamma'$ lattice misfit on the bypass mechanisms, which we approximate by an additional precipitate stress computed according to Eshelby’s inclusion theory. We show that in the shearing and hybrid looping-shearing regimes, a lattice misfit generally results in an increased bypass stress. For sufficiently high lattice misfit, the critical bypass configuration in attractive dislocation-precipitates interactions changes dramatically, and the bypass stress is controlled by the pinning of the trailing dislocation on the exit side of the precipitates, similar to what has been reported in the high-temperature creep literature.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Extrapolation of thermal conductivity in non-equilibrium molecular dynamics simulations to bulk scale

Predictions of the bulk scale thermal conductivity of solids using non-equilibrium molecular dynamics (MD) simulations have relied on the linear extrapolation of the thermal resistivity versus the reciprocal of the system length in the simulations. Several studies have reported deviation of the extrapolation from linearity near the micro-scale, raising a concern of its applicability to large systems. To investigate this issue, present work conducted extensive MD simulations of silicon with two different potentials (EDIP and Tersoff-II) for unprecedented length scales up to 10.3 μm and simulation times up to 530 ns. For large systems ≥0.35 μm in size the non-linearity of the extrapolation of the reciprocal of the thermal conductivity is mostly due to ignoring the dependence of the thermal conductivity on temperature. To account for such dependence, the present analysis fixes the temperature range for determining the gradient for calculating the thermal conductivity values. However, short systems ≤0.23 μm in size show significant non-linearity in the calculated thermal conductivity values using a temperature window of 500 ± 10 K from the simulations results with the EDIP potential. Since these system sizes are shorter than the mean phonon free path in EDIP (~0.22 μm), the nonlinearity may be attributed to phonon transport. For the MD simulations with the Tersoff-II potential there is no significant non-linearity in the calculated thermal conductivity values for systems ranging in size from 0.05 to 5.4 μm.

42 ENGINEERING↗

DOE OSTI2021

Phase evolution and structural modulation during in situ lithiation of MoS 2 , WS 2 and graphite in TEM

Li-ion batteries function by Li intercalating into and through the layered electrode materials. Intercalation is a solid-state interaction resulting in the formation of new phases. The new observations presented here reveal that at the nanoscale the intercalation mechanism is fundamentally different from the existing models and is actually driven by nonuniform phase distributions rather than the localized Li concentration: the lithiation process is a ‘distribution-dependent’ phenomena. Direct structure imaging of 2H and 1T dual-phase microstructures in lithiated MoS 2 and WS 2 along with the localized chemical segregation has been demonstrated in the current study. Li, a perennial challenge for the TEM, is detected and imaged using a low-dose, direct-electron detection camera on an aberration-corrected TEM and confirmed by image simulation. This study shows the presence of fully lithiated nanoscale domains of 2D host matrix in the vicinity of Li-lean regions. This confirms the nanoscale phase formation followed by Oswald ripening, where the less-stable smaller domains dissolves at the expense of the larger and more stable phases.

25 ENERGY STORAGE↗

DOE OSTI2025

Simulations of Heat Transfer Using Tight-Fitting Twisted Tape Inserts for First Wall Cooling in Molten Salt Breeder Blankets

One of the major components in fusion energy systems is the fusion blanket, which has a vacuum vessel to contain the plasma. As part of the fusion blanket/vacuum vessel, the first wall and plasma-facing components require sufficient cooling to prevent material degradation during operation from the superheated plasma. Most fusion blanket concepts involve first wall and divertor coolant channels with heat transfer enhancements (HTEs) that are intended to withstand the incident high heat fluxes of 1 to 5 MW/m 2 . Twisted tape inserts are a proposed HTE that have been investigated previously for first wall cooling and monoblock divertor cooling channels and in other nonfusion heat transfer components. By inserting twisted tapes into straight pipes, the amount of turbulence in the system can be increased at lower Reynolds numbers by swirling the flow. This results in better heat transfer characteristics with marginal increases in frictional pressure losses. In particular, simulations of high-Prandtl-number fluids such as the proposed molten salt FLiBe in twisted tapes, which is prototypic to liquid immersion blankets, have not been previously explored. Here, in this study, we simulate various Prandtl numbers in pipes with twisted tape inserts using large eddy simulations to determine the effects of increasing Prandtl numbers on heat transfer performance. The quantities of particular interest are the Nusselt number and the friction factor, which were recovered using data reduction techniques to determine impacts on heat transfer and pressure losses. This work serves as a starting point for determining the feasibility of twisted tape inserts for liquid immersion blanket concepts.

LES↗

DOE OSTI2022

Development and application of ReaxFF methodology for understanding the chemical dynamics of metal carbonates in aqueous solutions

A new ReaxFF reactive force field has been developed for metal carbonate systems including Na + , Ca 2+ , and Mg 2+ cations and the CO 3 2– anion. This force field is fully transferable with previous ReaxFF water and water/electrolyte descriptions. The Me–O–C (Me = metal) three-body valence angle parameters and Me–C non-reactive parameters of the force field have been optimized against quantum mechanical calculations including equations of state, heats of formation, heats of reaction, angle distortions and vibrational frequencies. The new metal carbonate force field has been validated using molecular dynamics simulations to study the solvation and reactivity of metal and carbonate ions in water at 300 K and 700 K. The coordination radius and self-diffusion coefficient show good consistency with existing experimental and simulation results. The angular distribution analysis explains the structural preference of carbonate ions to form carbonates and bicarbonates, where Na + predominantly forms carbonates due to weaker angular strain, while Ca 2+ and Mg 2+ prefer to form bicarbonate monodentate in nature. Residence time distribution analyses on different systems reveal the role of ions in accelerating and decelerating the dynamics of water and carbonate ions under different thermodynamic conditions. The formation and dissolution of bicarbonates and carbonates in solution were explored on the basis of the protonation capability in different systems. Here, the nucleation phenomenon of metal carbonates at ambient and supercritical conditions is explained from the perspective of cluster formation over time: Ca 2+ ions can form prenucleation clusters at ambient temperature but show saturation with increasing temperature, whereas Na + and Mg 2+ ions show a rapid increase in cluster size and amount upon increasing time and temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Linear Solvers for Collector Systems of Generalized Large-scale Inverter-Based Resources

Collector systems for inverter-based resources (IBRs) are typically represented by equivalent circuits for electromagnetic transient (EMT) simulations. Recent studies have revealed that modeling a detailed collector system is essential to accurately represent the behavior of IBRs, especially when dealing with partial tripping during external disturbances. However, there are several challenges in simulating a detailed EMT model of a collector system due to the time required to simulate such systems. Thus, this paper investigates the modeling of a detailed collector system, taking into account its configuration and components as defined in IEEE standard 2800. The configurations include the collector systems of generalized large-scale IBR plants. The components include the main IBR transformer, collector bus, and feeders with lines and/or cables. The EMT model of the collector system is represented by differential algebraic equations (DAEs) that are discretized to form linear equations that are solved using linear solvers. In this paper, linear solvers are proposed based on the Schur complement method, which are utilized for simulation of the EMT model of collector systems of generalized large-scale IBRs to accelerate simulation speed while maintaining the accuracy of the results. The proposed solvers are verified by comparing the performance to that of linear solvers provided in MATLAB.

Choi, Jongchan↗

DOE OSTI2021

Radioactive decay of $\mathrm {{}^{90}Sr}$ in cement: a non-equilibrium first-principles investigation

Cement is an inexpensive and relatively easily manageable material that is used as a last barrier for nuclear waste disposal. Under these conditions, the cement is in contact with low radiation doses, but there is a distinct possibility of being contaminated with radioactive products. Of particular concern is the medium lived half-life product 90 Sr (28.8 years) due to its ability to replace Ca. 90 Sr undergoes β-decay to 90 Y which, in turn, β-decays to stable 90 Zr. In this work, we discuss systematically the chain of non-equilibrium processes that result as a consequence of β-decay events in cement. We first use density functional-based methods to study the consequences of the sudden increase of the nuclear charge from Z to Z+1, a possible induced ionization and the perturbation of the surrounding electronic charge. Secondly, we use molecular dynamics simulations to study the recoil of the daughter nucleus. Finally, we discuss the damage caused by the ionization cascade produced during the propagation of the β-electron and the resulting chemical and structural perturbation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Thermal chains and entrainment in cumulus updrafts, Part 1: Theoretical description

Recent studies have shown that cumulus updrafts often consist of a succession of discrete rising thermals with spherical vortex-like circulations. Herein, a theory is developed for why this “thermal chain” structure occurs. Theoretical expressions are obtained for a passive tracer, buoyancy, and vertical velocity in axisymmetric moist updrafts. Analysis of these expressions suggests that the thermal chain structure arises from enhanced lateral mixing associated with intrusions of dry environmental air below an updraft’s vertical velocity maximum. This dry air entrainment reduces buoyancy locally. Consequently, the updraft flow above levels of locally reduced buoyancy separates from below, leading to a breakdown of the updraft into successive discrete thermals. The range of conditions in which thermal chains exist is also analyzed from the theoretical expressions. A transition in updraft structure from isolated rising thermal, to thermal chain, to starting plume occurs with increases in updraft width, environmental relative humidity, and/or convective available potential energy. Corresponding expressions for the bulk fractional entrainment rate ε are also obtained. These expressions indicate rather complicated entrainment behavior of ascending updrafts, with local enhancement of ε up to a factor of ~2 associated with the aforementioned environmental air intrusions, consistent with recent large eddy simulation (LES) studies. These locally large entrainment rates contribute significantly to overall updraft dilution in thermal chain-like updrafts, while other regions within the updraft can remain relatively undilute. Part 2 of this study compares results from the theoretical expressions to idealized numerical simulations and LES.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Evaluation of Phase Change Plaster/Paste Composites for Building Envelopes

Thermal energy storage (TES) materials increase a building’s thermal mass and provide the flexibility to shave and time-shift energy use. These materials are needed to enable grid-interactive efficient buildings, allowing building occupants to lower costs by leveraging free ambient heating/cooling as well as shifting consumption to cheaper electricity periods. Phase change materials (PCMs) that store heat through a solid-liquid/solid-solid transformation are of particular interest for TES because of their ability to store large amounts of energy across a constant temperature or very small temperature range. Furthermore, previous studies on PCM building envelopes evaluated products that may only be incorporated as part of new construction and are not easily applicable to existing, energy-intensive buildings. In this work, we leverage laboratory PCM characterization and whole-building energy simulation to study the energy and cost savings of PCM plaster/paste coatings, which can be applied post construction to both new builds and retrofits of existing buildings (about 40% of residential buildings in the United States fall into the retrofit category). We evaluated the energy impact of PCM plaster/paste composites for a single-family residential building in the United States, including layer design (location and thickness), composite properties (phase change temperature, energy density, and shape of the enthalpy curve), and U.S. climate zone. The enthalpy curve shape of PCM shows that phase change temperature and energy density are not the only two parameters that must be considered for PCM selection. Maximizing the energy density amid the effective TES zone is key to boosting energy savings. We also studied PCM price targets with the aim of achieving reasonable payback periods based on both block and time-of-use electricity rates. Moreover, we developed an empirical equation to provide a rough estimate of the energy savings potential of incorporating PCM without the need to perform detailed simulations. These results will help guide the selection and design of PCM plaster/paste composite elements for residential building envelopes.

ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATION↗

DOE OSTI2021

Conduction Heat Transfer through Solid in Porous Materials: A Comparative Study by Finite-Element Simulations and Effective Medium Approximations

The development of insulation materials with low effective thermal conductivity is essential for energy savings in various applications, including buildings, food services, pipe insulation, and refrigeration. Such materials can be developed by using micro- or nanoporous structures, as well as low-thermal conductivity gases and/or reduced pressure. A variety of effective medium approximation (EMA) models have been developed to study thermal transport through solid in porous structures. However, in many cases, the impacts of porosity on solid conduction are not well predicted by EMA models because of the assumptions made for simplification. Furthermore, the results can vary by a factor of up to 1.5, depending on the morphology of the pores. Hence, proper guidance is needed to choose the appropriate EMA model for a given morphology. This work presents a finite element method study using COMSOL Multiphysics software for various geometries, including hollow cubes, overlapping hollow spheres (normally stacked and tightly packed), and Voronoi structures, and compares the results to those obtained from various EMA models. Simulation results for the hollow cubes and Voronoi structures matched well with the Maxwell-Eucken and Russell model. The results for hollow spheres closely matched the Bauer model and the Glicksman model with fitting parameters. This work provides guidance on analyzing and designing insulation materials for energy savings in the future.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Preface to the JPCM special issue on intense radiation sources in condensed matter and materials physics

This special issue is designed to honour the contributions of our colleague Neville Greaves who was among the first to apply newly-available synchrotron x-ray techniques to study the synthesis, structures and properties of advanced materials [1, 2]. His wide ranging work on ceramics, glasses and nanomaterials was celebrated in a recent issue of the J. Non-Cryst. Solids [3]. Neville's contributions to synchrotron research, highlighted here, focussed on design and construction of new beamline capabilities and sample environments to probe materials under different physical and chemical conditions, integrating various techniques to obtain complementary data on structure and bonding. The innovative nature of his work is illustrated by his studies of liquids and glasses, where he incorporated aerodynamic levitation into these experimental techniques, using both synchrotron x-ray and neutron beams to probe the structural transformations and dynamics combined with molecular dynamics simulations to study atomic scale organisation. His work has had a major impact on fundamental condensed matter physics and chemistry as well as applied materials science; and in this issue, we present contributions from colleagues using and developing new experimental techniques at current and next-generation synchrotron, neutron and free-electron laser sources and in their own laboratories, that probe the fundamental physics, chemistry and functional properties of novel materials and processes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Simulation Analysis of Mode Hopping Impacts on OFDR Sensing Performance

This article examines the impacts of mode hopping on the sensing performance of optical frequency domain reflectometry (OFDR) and explores the potential for developing economical OFDR interrogators employing low-cost distributed feedback (DFB) lasers. By conducting numerical simulations, this study reveals that mode hopping has minimal effects on distance sensing measurements in free space due to the limited duration of beat interference signal at the incorrect frequency within the coherence length. Additionally, the simulations indicate that mode hopping only slightly affects the distributed strain sensing of OFDR, resulting in an error range of less than ±1µε when 100µε is applied to the sensing fiber. These findings highlight the potential of using low-cost DFB lasers with a 1-nm wavelength sweep range and a 1-MHz linewidth as tunable laser sources in OFDR while maintaining reliable and accurate sensing performance.

Wang, Qirui (ORCID:0000000336532984)↗

DOE OSTI2021

Monte Carlo Analysis of the Performance of the ITER Diagnostic Residual Gas Analyzer

The composition of exhausted gas is a key parameter in long-pulse plasma fusion experiments, and its evolution shall be monitored at timescales relevant to plasma dynamics and plasma-wall interactions. A diagnostic residual gas analyzer (DRGA) is a multisensor instrument particularly suited to these studies, and ITER will adopt DRGAs in the equatorial and in the divertor tokamak regions. In this work, we have revisited the design of the ITER divertor DRGA through simple vacuum analytical considerations supported by simulations conducted with Molflow+, a test particle Monte Carlo (TPMC) simulation code commonly used in the particle accelerator community. Here, starting with recommendations on the manufacturing of the vacuum piping of the DRGA, this work is followed by a complete vacuum characterization of the diagnostic vacuum setup (pressure profiles at base pressure and during sampling, orifice diameter, and length optimization), and finally, the in-vessel residence time of the most important gas species is simulated. These studies have allowed us to give insights into some experimental results recently found on the prototype DRGA installed in the Wendelstein W7-X stellarator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

High-Fidelity CFD Simulation of Mixed Convection and Forced Convection in a Pebble Bed Test Reactor Core

The Hermes low-power [35-MW(thermal)] reactor will be built and operated by Kairos Power LLC (KP) to demonstrate its fluoride salt-cooled high-temperature reactor (FHR) technology. In the KP FHR, the reactor core is composed of randomly packed pebbles with TRISO fuel particles inside with FLiBe flow upward through the core acting as a coolant. Previous numerical and experimental studies have been limited to either a small-size bed or to a lack of detailed measurements for heat transfer. Here, to address the lack of high-fidelity heat transfer data in a real-size FHR core, in this study, we simulated a pebble bed core with 34 374 pebbles randomly packed, similar to the Hermes reactor's size. The core radius was 14 times that of the pebble diameter, while the core height was 45 times. In this work, we were particularly interested in a mixed convection regime, where buoyancy is important. Therefore, we performed several large-eddy simulations at different Reynolds numbers (160 to 1000) with gravitational force included. The spectral element computational fluid dynamics code NekRS with graphics processing unit acceleration was used for this study. The low-Mach number approximation was applied to address property changes in the FLiBe and to account for buoyancy. A pure hexahedral mesh with 60 million elements was generated by the Voronoi cell method. At the polynomial order of 5, the total degrees of freedom was 7.5 billion. The developed case in this work is the first of its kind in terms of size and complexity. The local numerical data across the domain were obtained and compared with empirical correlations. After examining the data, we found the following conclusions. For pressure drop, the Reger correlation predicted less than a 5% error. On the other hand, for heat transfer, the Wakao correlation outperformed the others. Based on our findings, we recommend the use of the Wakao correlation for the Nusselt number calculation, and for pressure drop, the KTA (Kerntechnischer Ausschuss) correclation, among the available experimental correlations. In conclusion, the Reger direct numerical simulation-driven correlation for pressure drops should also be considered, given its best agreement with our calculations.

Mixed Convection↗

DOE OSTI2023

Determining Stress in Metallic Conducting Layers of Microelectronics Devices Using High Resolution Electron Backscatter Diffraction and Finite Element Analysis

Delayed failure due to stress voiding is a concern with some aging microelectronics, as these voids can grow large enough to cause an open circuit. Local measurements of stress in the metallic layers are crucial to understanding and predicting this failure, but such measurements are complicated by the fact that exposing the aluminum conducting lines will relieve most of their stress. In this study, we instead mechanically thin the device substrate and measure distortions on the thinned surface using high resolution electron backscatter diffraction (HREBSD). These measurements are then related to the stresses in the metallic layers through elastic simulations. This study found that in legacy components that had no obvious voids, the stresses were comparable to the theoretical stresses at the time of manufacture (≈300 MPa). Distortion fields in the substrate were also determined around known voids, which may be directly compared to stress voiding models. In conclusion, the technique presented here for stress determination, HREBSD coupled with finite element analysis to infer subsurface stresses, is a valuable tool for assessing failure in layered microelectronics devices.

HREBSD↗

DOE OSTI2023

Dynamic scaling of order parameter fluctuations in model B

Brookhaven National Laboratory has performed a dedicated set of experiments, called the Beam Energy Scan (BES), with the goal of location a possible critical point in the phase diagram of nuclear matter, This point is analogous to the endpoint of the liquid-gas phase transition in water, where the liquid phase corresponds to a phase of hadrons, and the gas phase to a plasma of quarks and gluons. In order to describe the data generated by these experiments we have studied numerical simulations of the stochastic diffusion equation with a conserved charge. Here, we focus on the dynamics in the vicinity of a critical point in the Ising universality class. The model we consider is expected to describe the critical dynamics near a possible QCD critical point if the coupling of the order parameter to the momentum density of the fluid can be neglected. The simulations are performed on a spatial lattice, and the time evolution is performed using a Metropolis algorithm. We determine the dynamical critical exponent z ≃ 3.972(2), which agrees with predictions of the epsilon expansion. We also study non-equilibrium sweeps of the reduced temperature and observe approximate Kibble-Zurek scaling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

The HTR-Proteus Benchmark: Analysis and Use as a Verification and Validation Case

This presentation outlines the evaluation and application of the HTR-Proteus benchmark as a verification and validation (V&V) case for advanced reactor modeling tools. The work supports the U.S. Department of Energy’s HALEU Availability Program (HAP) and the joint DOE/NRC DNCSH project, which aims to reduce criticality safety uncertainties in commercial-scale HALEU fuel cycle and transportation systems. The HTR-Proteus experiments, conducted at the Paul Scherrer Institute, provide high-fidelity data for TRISO-fueled, graphite-moderated pebble bed reactors with high neutron leakage—conditions relevant to HALEU transport scenarios. This study focuses on Cores 4.2 and 4.3 of the HTR-Proteus benchmark, analyzing key sources of uncertainty including pebble packing, TRISO particle positioning, and core height. Using Project Chrono for realistic pebble geometries and Serpent, SHIFT, and MCNP for neutronics simulations, the study quantifies the impact of these uncertainties on the effective multiplication factor (keff). Results show that a sample size of 110 pebble configurations is sufficient to converge keff, with ±30 pcm uncertainty due to packing randomness. TRISO positioning and core height variations also significantly influence keff, highlighting the importance of detailed modeling in V&V efforts. The benchmark serves as a valuable test case for validating the Griffin reactor physics code and improving confidence in HALEU system simulations.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗