SEARCH · Engineering Papers
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.
A drone may be forced down as part of a kinetic Counter-UAS action or as a spoof signal with a command to land immediately. A drone could also be incidentally discovered by an employee or by a random patrol. In each case, the platform should be considered a suspicious package until it can be determined to be nonhazardous by qualified personnel.
The relationship between noble metal phase particles and fission gas bubble production in used nuclear fuel is described. The majority of Te atoms within noble metal phase undergo radioactive decay to form stable Xe within a few hours after particle formation. This results in the production of clusters of Xe atoms contained within the solid metal matrix exhibiting an equivalent gas bubble pressure approaching 1 GPa. These high pressure bubbles are stabilized by the UO2 within the bulk of the fuel. However, when these bubbles form near the fuel / cladding interface, in combination with local and temporal damage caused by fission recoil, they are capable of overcoming the fracture strength of the UO2 and rupturing catastrophically. The force of the resulting bubble rupture is sufficient to eject noble metal phase particles several microns into the cladding. This proposed mechanism explains the observance of noble metal phase in cladding and is consistent with a host of morphological features found near the fuel / cladding interface.
Active metals in used nuclear fuel dissolve into the salt during pyroprocessing and are not removed by electrorefining or drawdown operations. The buildup of 137 Cs over time increases the heat load and ionizing radiation level of the salt such that it must be replaced frequently, resulting in a significant amount of salt waste. An effective means of managing cesium in the molten salt electrolyte would increase the efficiency of pyroprocessing and decrease the volume of salt waste requiring disposal. A previous report summarized issues that must be addressed when developing a removal strategy and assessed the suitability of existing methods and remaining technological gaps to their application (Rose and Thomas 2023). Cesium is extremely stable in molten salt as a chloride—even more stable than the LiCl-KCl eutectic base salt used for pyroprocessing fuel—which makes removing cesium a challenge. However, sufficiently strong atomic interactions occur between active metal species and liquid metals that make the electrodeposition of active metal fission products into liquid metal electrodes energetically favorable. The feasibility of recovering cesium from LiCl-KCl pyroprocessing salt through electrodeposition into liquid metals is eing assessed by identifying potentially effective liquid metals and performing tests to determine the effectiveness of electrodepositing cesium from a LiCl/KCl salt into these liquid metals. Previous studies investigating the electrodeposition of Sr 2+ , and Ba 2+ into zinc, cadmium, bismuth, lead, tin and antimony have shown that alkali and alkaline earth metals can be electrodeposited at liquid metal cathodes (Kim et al., 2018). The removal of Ba 2+ and Sr 2+ was measured to be more efficient than the removal of monovalent cations due to the greater thermochemical driving force for alloying those elements with the liquid metal (Jang et al. 2022). Because the equilibrium potentials are dependent on the interactions of the active metal in the liquid metal, it is likely that the other alkali metals Li + , and K + , will deposit from LiCl/KCl salt with the Cs + . Therefore, application of this method to recover active metals from pyroprocessing salt will benefit from the use of a liquid metal and set of operating conditions that sufficiently increase the reduction potential of cesium to remove cesium from the waste salt with an acceptable amount of co-deposited lithium and potassium.
We investigate some aspects of the connection between mean-field and ab-initio calculations for any-nucleon systems employing the Gogny effective interaction. In particular, we present variational Monte Carlo calculations of 4 He and 16 O, for which we extended the formalism of the linear method to the case of a density-dependent Hamiltonian. Monte Carlo calculations using an uncorrelated basis reduce to Hartree-Fock ones, but Galilean invariance can be preserved, contrary to the usual treatment. An analysis of the results obtained with the Gogny force when used in conjunction with correlated wave functions has also been performed. Here, the calculations show the collapse of the deuteron, leading to a set of unphysical predictions for the structure of heavier nuclei. This behavior remarks the different role of the Gogny effective interaction as compared to bare interactions typically used in fully microscopic calculations.
Measuring the pore size and pore-size distributions and exploring the fluid-exchange dynamics between different types of pores in porous materials remains a significant experimental challenge but is critical to understanding catalysis, chromatography, nutrient cycling, and a whole range of geochemical phenomena, including shale gas and tight gas extraction. Here, we present the results of 1D 13 C NMR and 2D exchange spectroscopy (EXSY) NMR investigations of a porous silica using supercritical methane (scCH 4 ) as a direct probe of pore size and fluid exchange between pore types. The results show that the 13 C chemical shift of scCH 4 adsorbed in nanometer-scale silica pores becomes more negative with increasing pore diameter, in agreement with trends reported for gas hydrates, zeolites, MOFs, and clays and other microporous (<10 nm pores) geochemical materials. These 13 C chemical shifts follow a natural log–linear trend with pore size in the vacuum-dried porous silicas studied here, allowing one to predict pore size based on the 13 C chemical shift in dry silica nanopores. The EXSY and 1D 13 C NMR results both show that CH 4 exchanges between pore and bulk fluid environments over rate scales from 0.01 to 1 kHz in vacuum-dried samples and that exchange dynamics must be considered when interpreting 13 C NMR of pore-adsorbed CH 4 . Introducing H 2 O to the system causes the 13 C chemical shift to become more negative with increasing H 2 O content as a result of H 2 O preferentially filling small pores, forcing CH 4 to occupy larger pores that are better connected to the bulk environment. Likewise, the 1D 13 C data show a decreased distribution of exchange rates between pore-adsorbed and bulk CH 4 in the presence of H 2 O.
We have undertaken an experimental-computational project that addresses a few key technology gaps associated with the use of SiC/SiC composites for light water reactor fuel cladding. Two principal endeavors in this project are to assess the irradiation-induced microstructural changes and swelling in SiC/SiC composites, and to characterize and model the porous oxide surface layer (when exposed to steam) along SiC recession that is detrimental to the clad integrity under accident conditions. The project tasks are: (i) ion irradiation and characterization (University of Tennessee, Knoxville), (ii) mechanical tests and analysis (University of South Carolina), (iii) steam exposure tests (Oak Ridge National Laboratory), (iv) electron microscopy and spectroscopic characterization (NC State University), (v) x-ray microscopy and reconstruction (University of South Carolina and NC State University), (vi) phase field modeling and simulations (NC State University and Idaho National Laboratory). The effects of 10 MeV Au ion irradiation at 350°C on the microstructure evolution in SiC/SiC composites are investigated at doses up to 400 displacements per atom (dpa) at the University of Tennessee, Knoxville. Atomic force microscopy and optical profilometry reveal irradiation induced axial and radial shrinkage of the fibers for doses greater than 10 dpa. Based on detailed electron microscopy characterization, the primary cause of the fiber shrinkage is attributed to irradiation-induced loss of carbon packets. Additionally, the multilayer PyC interface is observed to portray high resistance to irradiation damage. The mechanical response and failure mechanisms of un-irradiated samples is also assessed through loading tests and X-ray imaging at the University of South Carolina. Steam exposure tests are performed at the Oak Ridge National Laboratory in a facility that represents a reactor pressure vessel under a loss-of-coolant-accident scenario. The samples that are analyzed methodically in this report are tested for 32/31 hours at 1200°C with a velocity of 0.25 cm/s for the pressures: 0.1 MPa, 0.45 MPa, 0.92 MPa and 1.38 MPa. Scanning/transmission electron microscopy analysis conducted at the NC State University (NCSU) shows that the oxide layer thickness increases with the steam pressure. While the oxide layer is crystalline (α- cristobalite) for the pressures 0.45 MPa, 0.92 MPa and 1.38 MPa, the layer is amorphous at 0.1 MPa. Results from Raman spectroscopy have also confirmed the formation of α-cristobalite phase of SiO₂. The abrupt increase in the integrated Raman intensity ratio between steam pressures 0.1 and 0.45 MPa suggests the onset of accelerated crystallization. Non-destructive three dimensional X-ray microscopy/tomography (XCT) and computational image processing techniques are employed by the University of South Carolina to probe the porous oxidation features on SiC samples at varying pressures. Interestingly, most pores are observed to be located away from the surface as well as the oxide-SiC interface. The average oxide layer thickness assessed from the XCT analysis is seen to be in excellent agreement with the values determined through scanning electron microscopy analysis for the highest pressures where the oxide layer is relatively more uniform. A phase-field model developed by the NC State University and Idaho National Laboratory for simulating oxidation of SiC by steam captures the paralinear kinetics of SiC oxidation with a high degree of fidelity. Results from quasi-one-dimensional and two dimensional simulations show that the pores with oxidizing species lead to a higher volatilization rate. These simulations indicate that the enhanced apparent volatilization rates at higher pressures observed in experiments can be rationalized by the increased volatilization from the pores. The project team has also successfully developed the capability to import images from experiments for realistic evolution of the oxidizing microstructure.
A new technique for force reconstruction is developed. To estimate the externally applied force, this technique sums the weight-scaled acceleration signals, and is referred to as the Sum of Weighted Accelerations Technique (SWAT). To obtain the scalar weights the inverse of the mode shape matrix is used. Application of this technique is illustrated with both numerical calculations using a mass-spring model and experimental data from a structure impacting a rigid barrier.
The development of algorithms for machine learning and data analysis for the 3013 MIS corrosion surveillance program is a collaborative effort by SRNL, USC and GT. For corrosion detection, LCM image data is extracted from large binary files, with software written to convert the data to physical attributes (i.e. height, color and grayscale values; all as functions of a location in a plane projection). The user interface for the software permits selective downloading of binary data and interrogation of attributes. User input thresholds are used to flag attributes of interest. Machine learning algorithms, developed for this application, are used to determine whether the features are the result of corrosion. To address the fundamental mechanisms of corrosion, machine learning algorithms are being developed to derive interatomic potential force-fields from ab-initio DFT calculations. The goal is to apply molecular modeling on a large enough scale to guide the design of resistant materials.
In the U.S., a key component of the commercialization of advanced reactors is completion of a license application, which must ultimately be approved by the Nuclear Regulatory Commission (NRC). The approval of the license application by the NRC is contingent on satisfactory demonstration of the design basis and the response of the advanced reactor design to transient and accident scenarios using accepted codes and methods. This report describes the qualification and dedication requirements that the advanced reactor safety analysis system software SAS4A/SASSYS-1 are expected to need to fulfill to be used for sodium-cooled pool-type fast reactor licensing. The qualification and dedication requirements are identified through performance critical characteristics and evaluation model acceptance criteria representative of the advanced reactor design considered for licensing. This document captures, additionally, the verification process developed to demonstrate that the software fulfills the qualification and dedication requirements for a generic sodium-cooled pool-type fast reactor as part of the commercial grade dedication process. Like most software that has primarily existed in the research and development space, the most significant challenge facing SAS4A/SASSYS-1 for use in a licensing framework is the availability of a documentation basis describing the code pedigree. SAS4A/SASSYS-1 has been used for licensing of the fast flux test facility (FFTF) and the JOYO sodium-cooled fast reactor in Japan, as well as the design of the CRBR Plant. However, the historical verification and validation (V&V) activities supporting SAS4A/SASSYS-1 development do not align with modern software quality assurance (SQA) and V&V requirements. Two approaches to use of SAS4A/SASSYS-1 in a commercial licensing framework have been identified: commercial-grade dedication (CGD) and software qualification. The methods and requirements prescribed in the ASME NQA-1-2008/2009 Standard and Regulatory Guide 1.203 on the evaluation model development and assessment process (EMDAP) have been used as guidance to define the CGD and qualification processes, respectively. A qualification and dedication requirements matrix has been developed which utilizes fundamental software verification. In this process, software verification is defined as a software quality process aimed at defining software requirement specifications, developing software design documentation, and performing and documenting acceptance testing of the code against requirements. A key element of software qualification and dedication includes determination of software acceptance with respect to critical characteristics relevant to the functional requirements of the software. To assist with identification of cross-cutting transient phenomena and functional requirements, domestic SFR vendor designs have been reviewed to identify a reference SFR design. For this report, the reference design is defined as a pool-type reactor with metal alloy fuel, a liquid-metal intermediate heat transport system, and passive decay heat rejection systems. Given this reference, a series of high-level cross-cutting phenomena was identified for a general class of single-fault undercooling or reactivity insertion transients that scopes the design basis space, with the goal of assisting with prioritization of documentation development efforts for key transient models in SAS: 1) Reactivity feedback response prior to scram; 2) System-wide thermal inertia; 3) Transition in natural circulation flow regime in heat removal systems; 4) Decay heat generation; 5) Steady-state fuel characterization; 5) Clad/fuel behavior at elevated temperatures; 6) Point kinetics and decay heat; 7) Pump coastdown behavior; 8) Core flow redistribution in loss of forced convection; 9) Pool stratification. As a demonstration of CGD of SAS4A/SASSYS-1 for a sodium pool reactor, a software qualification and dedication gap analysis as it relates to code documentation has been performed. This effort leverages the framework established as part of the SAS4A/SASSYS-1 SQA Program. This CGD demonstration provides a framework that vendors can build upon to demonstrate the applicability of the SAS4A/SASSYS-1 software for licensing a sodium-cooled pool-type fast reactor.
Despite surging interest in molten salt reactors and thermal storage systems, knowledge of the physicochemical properties of molten salts are still inadequate due to demanding experiments that require high temperature, impurity control, and corrosion mitigation. Therefore, the ability to predict these properties for molten salts from first-principles computations is urgently needed. Herein, we developed and compared a machine-learned neural network force field (NNFF) and a reparametrized rigid ion model (RIM) for a prototypical molten salt LiF–NaF–KF (FLiNaK). We found that NNFF was able to reproduce both the structural and transport properties of the molten salt with first-principles accuracy and classical-MD computational efficiency. Furthermore, the correlation between the local atomic structures and the dynamics was identified by comparing with RIMs, suggesting the significance of polarization of anions implicitly embedded in the NNFF. Furthermore, this work demonstrated a computational framework that can facilitate the screening of molten salts with different chemical compositions, impurities, and additives, and at different thermodynamic conditions suitable for the next-generation nuclear reactors and thermal energy storage facilities.
Proton conduction in acceptor-doped perovskites is fundamentally a vibronic process: mobile H + and D + do not move independently, but dynamically co-vibrate with the surrounding oxygen–metal framework. Direct experimental evidence for this behavior is presented using in situ 119 𝑆𝑛 nuclear resonance vibrational spectroscopy (NRVS) on hydrated, deuterated, and dry 𝐵𝑎𝑆𝑛 1−𝑥 𝑌 𝑥 𝑂 3−𝛿 . Hydration induces systematic redistributions in the Sn-projected phonon density of states (PDOS), including an upshift of the first spectral moment by about 0.4 meV, indicating a stiffening of the extended Sn–O network. H/D isotopic substitution leaves the Sn-projected PDOS largely unchanged, with only subtle isotope-dependent spectral reweighting, demonstrating that protonic degrees of freedom are not localized oscillators but are embedded in collective lattice modes. These results are rationalized using a classical coupled proton–phonon oscillator model that links the observed PDOS variations to changes in effective force constants and vibrational mass terms. The model captures how H + and D + participate in cooperative lattice dynamics rather than forming isolated OH/OD entities. Overall, NRVS probes proton–lattice coupling in ceramic proton conductors and quantitatively describes how protonic defects modulate host lattice dynamics to enable phonon-assisted long-range proton transport.
This report serves as the first delivery of a four-year applied science effort to transform and advance the error bounds for the yield estimate of an explosion. Each year’s delivery will be in this form, culminating in the submission of this work for peer review to a scientific journal. Importantly, the yearly progress reports can then also be viewed as expanding drafts working towards a formal journal article submission. For the first tranche of funding, we collaborated with Air Force Technical Applications Center (AFTAC) scientists to identify unclassified real-world data that demonstrate and validate our advanced error propagation methods. Collaboration includes visits to AFTAC and telecons. For this development, we illustrate the fusion of seismic, acoustic, optical, and surface effect signatures from an explosion. The mathematics and code being adapted to this specific application (Williams et al., 2021) involves physics models of multiple sensor signatures. We have also identified related physics models and have integrated them into code. Current methods of underground explosion yield estimation for the Threshold Test Ban Treaty (TTBT) have served the US treaty monitoring mission well for decades. A research objective of the Defense Nuclear Nonproliferation Research and Development (DNN R&D) office of the National Nuclear Security Administration (NNSA) has always been to provide new technical capabilities for monitoring lower thresholds. The general model and error propagation code to be developed in this project is based on significant advances in error modeling and propagation needed to analyze data at lower detection thresholds. The second tranche of funding for this project began on May 1, 2022, and planned work for the second tranche includes: i) completing the integration of physical model code into the general error model framework; this code accommodates a wide range of linear/nonlinear source models, fixed/ random effects, and frequentist/Bayesian analyses (the purpose of which is not to dictate to users how to analyze data, but instead to allow users the maximum flexibility in their work); ii) illustrative application of code to identified data, and; iii) initial planning with AFTAC researchers on delivery of code to the Common Development Environment at AFTAC, and continued writing of the planned final journal article submission (year two of this progress report), with particular emphasis on descriptions of data identified for this effort.
A set of subroutines that interact with an existing computer program, Quantum Espresso (QE), that takes a set of electronic orbitals from QE and propagates them forward in time in response to a number of perturbations. These perturbations include the motion of ions and a wide spectrum of applied electromagnetic fields (e.g., static magnetic fields, slowly-varying electric fields, x-rays, etc.). From the time evolution of the electronic orbitals comes time-varying values of physical observables like the electronic charge and current densities and atomic forces. These are further post-processed to yield properties of interest like stopping powers, conductivities, and the dynamic structure factor. The electronic orbitals are described in a plane wave basis and the equations of motion are carried out using a unitary or approximately unitary time propagation algorithm. Related work will include subroutines for interfacing this capability, as well as others that already exist within QE, with a fork of QE that is more extensible and modular.SAND2019-12611 M Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA0003525.
Seismic source characterization is an essential function of global nuclear explosion monitoring (NEM). While large events (roughly with moment magnitude, M w , greater than 5.0) can often be easily detected, located and identified with high signal-to-noise ratios at teleseismic distances (> 20°), trends in NEM research require confident source characterization at much lower magnitudes (say down to 3.0) and exploitation of sparse observations (from only a few stations) at regional distance (< 20°). Regional distance waveform inversion to characterize sources is now widely used and effective (e.g. Ford et al., 2009; Alvizuri and Tape, 2018; Alvizuri et al., 2018; Chiang et al., 2018; Ford et al., 2022). These methods obtain the magnitude, depth and seismic moment tensor, which represents the forces that excited the observed seismic waves (slip on an earthquake fault, explosion, collapse or a combination of various forces). Common to many problems in seismology, the isolation of the source 2 properties requires removal of path propagation effects that waves experience while traveling through the three-dimensional (3D) Earth (the structure exists due to different rock types, material properties, temperature and tectonic processes).
It is important to test the possible existence of fifth forces, as ultralight bosons that would mediate these are predicted to exist in several well-motivated extensions of the Standard Model. Recent work indicated asteroids as promising probes, but applications to real data are lacking so far. Here we use the OSIRIS-REx mission and ground-based tracking data for the asteroid Bennu to derive constraints on fifth forces. Our limits are strongest for mediator masses m ~ (10 -18 -10 -17 ) eV, where we currently achieve the tightest bounds. These can be translated to a wide class of models leading to Yukawa-type fifth forces, and we demonstrate how they apply to U(1) B dark photons and baryon-coupled scalars. Our results demonstrate the potential of asteroid tracking in probing well-motivated extensions of the Standard Model and ultralight bosons near the fuzzy dark matter range.
In this work, we review phenomena relevant to tritium transport in molten-salt reactors, which produce tritium from lithium and beryllium salts at significantly higher levels than other reactor types. Modeling of such phenomena began following MSRE operations, and these early models attempted to predict measured tritium distributions in the MSRE, accounting for turbulent mass-transport processes (using established heat transfer correlations), permeation through a variety of metal structures such as heat exchanger tubes, and transport to and from bubbles introduced into the salt by gas sparging. The models reasonably reproduced the MSRE data, but did so best when the permeability of structures was reduced by about a factor of 1,000. This issue does not appear to have been conclusively resolved, and all of the more recent attempts to model tritium transport in molten salts appear to make use of the same methodology. MSRE remains, however, essentially our only source of integral tritium transport data relevant to MSRs. Here, we generalize the MSRE approach to permeation in order to include potential rate-limiting effects at interfaces, as well the effects of added hydrogen. Appropriately non-dimensionalized, this system of equations identifies two dimensionless numbers whose relative values clearly delineate the conditions under which mass transport, surface effects, permeation, and hydrogen swamping are expected to become rate-limiting. We also describe the preliminary conceptual design of a forced convection FLiBe loop, into which tritium would be introduced for the purpose of providing validation data for such a model. The primary purpose of this is to investigate the coupled transport phenomena described above and identify those that are rate-limiting in MSRs. Additional test section configurations are described that would address other transport phenomena relevant to MSRs, including bubbly flows and graphite interactions.
The desire to improve the economic competitiveness and deployment pace of nuclear energy through modularization, manufacturing, and series production had led to the development of smaller size reactors. As the standard 17x17 fuel technology is mainly maintained in the pressurized water reactors (PWRs) category, this translates into a lower number of fuel assemblies in the core and sometimes a reduced fuel height. To assess the impact of such scale change in core design on fuel cycle cost and spent fuel volume, a scoping analysis tool is developed based on infinite lattice calculations, leakage, fuel management reduced models, and levelized unit cost of electricity (LCOE) estimate. As such, cost dynamics driven by fuel specific power, burnup, core leakage, feed, cycle length, fuel assembly height as well as uranium market data are captured with consistent set of assumptions and analysis methods. A selection of 5 reactor designs representative of leading PWR developers is assessed and compared. Pursuing higher specific powers and optimal burnups are highlighted as the main fuel cost reduction drivers, nevertheless, practical limitations and opportunities must be evaluated to establish the feasibility of such enhanced fuel operation. In consequence, a detailed core design is performed using SIMULATE3 code for 5 PWR variations including natural and forced coolant circulation modes, two reactor scales, power densities of 73, 112, and 123 kW/l and higher discharge burnups. Design and optimization are performed at the lattice level, for the reflector, and at the core loading level. Satisfactory steady-state operation including power distribution, coolant operating limits, and reactivity requirements are analyzed and reported in this paper. The fuel economics of the detailed core designs confirm the scoping analysis findings. Despite the unlocked power uprates in small PWRs, the achievable burnup for a given fuel specific power requires more enrichment and shorter fuel height results in higher fabrication costs per mass of fuel, which makes scaling down core size a more expensive endeavor on the fuel cycle front. Spent fuel volumes are reported for the PWRs designed in this paper. Furthermore, these volumes are driven by the core average discharge burnup regardless of the scale in consideration. Additional cost and core performance aspects related to heavy reflector gains, fuel-reflector substitution, and disposal cost policy in the U.S. are examined.
Kairos Power, LLC, is developing its version of the Fluoride-cooled High-temperature Reactor, the KP-FHR. The design uses a pebble bed core with fluoride salt as a coolant. The pebbles used in the KP-FHR have a diameter of 4 cm, with a shell fuel region where TRISO particles are embedded. A Pebble bed core design is adopted by several Gen IV reactors, They boast many benefits, such as fuel integrity, highly efficient heat transfer, and passive safety. However, it is challenging to accurately predict temperature and flow inside a pebble bed. Traditional approaches use the porous media model, which regards the pebble bed as a continuous medium, but with different temperature fields representing different levels, such as the fluid temperature, pebble surface temperature, and pebble center temperature. Empirical heat transfer correlations are adopted to calculate the heat transfer coefficient between different phases. However, empirical correlations are usually validated with experimental data, which usually lacks detail inside the pebble bed. The available experimental data is also generally at a high Reynolds number, which falls outside of the conditions of KP-FHR. Explicit computational fluid dynamics (CFD) simulations of randomly packed pebble beds have only become feasible recently. This is thanks to the rapid development of computational power and scalable algorithms. In this work, we used the Spectral Element Method (SEM) CFD code NekRS to simulate the randomly packed pebble bed in a cylindrical container. NekRS, which is the GPU variant of Nek5000, but refactored to utilize the computational power of GPUs using the OCCA library to run on hybrid architecture high performance computing systems. It was initially developed with the libParamunal library, but truncated and tuned for large-scale turbulence simulation. As a result, the SEM reaches higher precision with the same degrees of freedom by using a high-order Lagrange polynomial basis distributed on Gauss-Lobatto-Legendre quadrature inside each element, compared to lower-order methods, such the Finite Volume Method and Finite Element Method. The report is divided into five parts. We start with a general discussion of the pebble bed reactor, along with a specific investigation into the KP-FHR. The second part presents the numerical methodology. In the third part, we study a modular pebble bed with 1741 pebbles in a container of 7 pebble-diameter radius. Beyond LES simulations done by NekRS, we also leveraged the thermal radiation model in OpenFOAM to study heat transfer under no-forced-flow scenarios. Then, in the fourth part we simulated a pebble bed similar to the size of the Hermes Test Reactor. The total number of pebbles is in these simulations is 34,374. The container radius is 14 pebble-diameters. Finally, the report concludes in part five, with a discussion of future work.