Neutronic and Thermal-Hydraulic Design Studies for High Flux Isotope Reactor Conversion to Low-Enriched Uranium High Density Silicide Dispersion Fuel
Explore the source record for details and available documents.
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.
Explore the source record for details and available documents.
The land that is potentially available to produce dedicated cellulosic bioenergy crops, often referred to as 'marginal' land, depends heavily on the underlying assumptions used to classify and identify it. In this study we compare three definitions and types of marginal land to identify the interactions between the bioenergy landscape and the logistics networks needed for the biofuel supply chain. Typical studies of the scale, cost, and greenhouse gas (GHG) mitigation potential of cellulosic biofuel take a land-centered approach which may neglect to account for the trade-offs between establishing bioenergy crops and the supply chain design decisions needed to allow those crops to be converted to liquid fuel. A mathematical programming approach is used to minimize the total annualized cost of a large-scale field-to-product system producing bioethanol in the USA midwest. Results show that a high concentration of marginal land leads to efficient systems and that the bioenergy landscape design becomes increasingly important with a higher emphasis on GHG mitigation. Additionally, targeted landscape design (including fertilization) with a focus on fields with high soil carbon sequestration potential can greatly reduce the system-wide GHG emissions for only a small increase in the unit cost of biofuel.
The limited durability of Pt electrocatalysis toward cathodic oxygen reduction reaction remains challenging, yet crucial for the development of Proton Exchange Membrane Fuel Cell. Here, we present a rational design of a robust catalyst consisting of PtCo nanoparticles supported on Pt-Co-N-graphene nanofiber prepared through electrospun Cobalt-Metal-Organic-Framework. The catalyst delivers unprecedented mass activity of 2.48 A·mgPt -1 , and retains 80% of initial value after 60,000 Accelerated-Stress-Test cycles. Operando X-ray absorption spectroscopies show that the electronic configurations of Pt sites in PtCo and Co sites in Co-N4 in the hybrid catalyst are modified toward high catalytic activities. Density Functional Theory unveils that the enhanced curvature of the substrate induced by the morphology engineering lowers the reaction thermodynamic barrier on Co-N4 sites, favoring the formation of H2O and suppressing that of H2O2. This result along with the strong affinity of PtCo nanoparticles to the Pt-Co-N-graphene fiber endows the catalyst an exceptional durability.
Nuclear Safety is a review journal that covers significant developments in the field of nuclear safety. Its scope includes the analysis and control of hazards associated with nuclear energy, operations involving fissionable materials, and the products of nuclear fission and their effects on the environment. Primary emphasis is on safety in reactor design, construction, and operation; however, the safety aspects of the entire fuel cycle, including fuel fabrication, spent-fuel processing, nuclear waste disposal, handling of radioisotopes, and environmental effects of these operations, are also treated. Table of Contents for this issue follows. GENERAL SAFETY CONSIDERATIONS: 1 Technical Note: A New Approach to Fission Reactor Safety, Yu. V. Petrov; 5 Erratum; ACCIDENT ANALYSIS: 6 Analysis and Modeling of Fission Product Release from Various Uranium-Aluminum Plate-Type Reactor Fuels, R. P. Taleyarkhan; CONTROL AND INSTRUMENTATION: 23 Applications of a Surveillance and Diagnostics Methodology Using Neutron Noise From a Pressurized-Water Reactor, R. T. Wood, L. F. Miller, and R. B. Perez; DESIGN FEATURES: 36 Westinghouse Advanced Passive 600 Plant, B. A. McIntyre and R. K. Beck; 47 System 80+™ PWR Safety Design, C. W. Bagnal, R. A. Matzie, and R. S. Turk; ENVIRONMENTAL EFFECTS: 58 The MATS Experiments—Mesoscale Atmospheric Transport Studies at the Savannah River Site, A. H. Weber, S. Berman, R. J. Kurzeja, and R. P. Addis; 75 Book Review: Environmental Radioactivity in the European Community 1984-1985-1986, C. A. Little; WASTE AND SPENT FUEL MANAGEMENT: 76 Activities Related to Waste and Spent Fuel Management, M. D. Muhlheim and E. G. Silver; OPERATING EXPERIENCES: 87 Aging Assessment of BWR Control Rod Drive Systems, R. H. Greene; 100 Reactor Shutdown Experience, Compiled by J. W. Cletcher; 103 Selected Safety-Related Events, Compiled by G. A. Murphy; 110 Operating U.S. Power Reactors, Compiled by M. D. Muhlheim and E. G. Silver; RECENT DEVELOPMENTS: 129 General Administrative Activities, Compiled by M. D. Muhlheim and E. G. Silver; 139 Reports, Standards, and Safety Guides, D. S. Queener; 143 Proposed Rule Changes as of Sept. 30 1991; ANNOUNCEMENTS: 128 Harvard School of Public Health Announces Short Courses; 128 16th Biennial ANS Topical Meeting on Reactor Operating Experience: Present and Future Technologies—Applying Lessons Learned (Call for Papers); 147 The Authors; 150 Indexes to Nuclear Safety, Volume 32; 154 Errata.
Recently, we studied the phase behavior of ionic fluids under confinement using the classical density functional theory within the framework of the restricted primitive model. Additionally, the theoretical results indicate that narrowing the pore size may lead to a drastic reduction in the electric double layer capacitance, while increasing the surface electrical potential would improve the ionic accessibility of micropores. In this work, we extend the theoretical investigation to systems containing size-asymmetric electrolytes that may exhibit a vapor-liquid like phase transition in the bulk phase. The effects of pore size and surface electric potential on the phase diagram and microscopic structures of the confined electrolytes were studied over a broad range of parameters. We found that decreasing the pore size or increasing the surface potential could destabilize the liquid phase in micropores, and capillary evaporation could occur regardless of the size asymmetry between cations and anions. Compared to that in a symmetric ionic system, the vapor-liquid phase separation is more likely to take place as the size asymmetry becomes more pronounced. The phase transition would alter the “accessibility” of ions to micropores and lead to coexisting micropores with different surface charge densities as identified by Monte Carlo simulation.
Heat pipe microreactor concept involves the use of a solid core block to encapsulate the nuclear fuel and the passive extraction of heat from the core block to the heat exchanger block by advanced heat pipes. The core block is a key component of the heat pipe microreactor design. It integrates the functions of reactor vessel, core structural components and fuel cladding of traditional Gen IV advanced reactor designs into one single structure. The effects of irradiation on the structural failure modes are outside the scope of the rules of the ASME Boiler and Pressure Vessel Code, Section III, Division 5, Subsection HB, Subpart B for metallic Class A components. Methods to account for the irradiation effect were investigated in order to develop design methods that can be applied to the core block design in the heat pipe microreactor. This report presents an initial development of an approach to account for the irradiation effects in the design rules.
One of the most pressing challenges to the continued deployment of nuclear energy systems is in the ultimate management and disposition of discharged fuel assemblies. While reprocessing and recovery of valuable materials from UNF assemblies has been considered as part of an overall strategy for minimization of the volume of reactor-based wastes to be managed, the deployment of commercial-scale reprocessing facilities presents an enormous economic challenge. The MARIE software package has been developed as a means of confronting this challenge. Representing components of a generic fuel reprocessing operation as individual physical processes, MARIE is designed as a modular framework intended to allow for analysis and cost-optimization for a hypothetical reprocessing facility while realistically accounting for the physical characteristics of the used fuel source term, such as decay heat, activity, and radiation dose (informing corresponding shielding requirements). Capabilities supported by MARIE include head-end operations such as fuel shearing, voloxidation, and dissolution; generic solvent extraction operations informed by available open-literature data; a suite of unit operations intended to represent electrochemical processing of used fuel assemblies (i.e., oxide reduction, electrorefining, and electrowinning); and finally, accounting for both costs and physical features of discharged waste streams, which can be used to inform follow-on analyses such as the feasibility of deep-borehole disposal of HLW. This paper presents an overview of the MARIE software capabilities, including how individual unit operations are implemented to enable a larger-scale optimization of a hypothetical reprocessing operation on aspects such as cost and recovery of valuable materials.
Oak Ridge National Laboratory (ORNL) successfully demonstrated the Versatile Remediation Module (VRM), a prototype module designed and built by ORNL for on-site remote repair of welded stainless steel storage containers for spent nuclear fuel and high-level radioactive waste. This paper describes the VRM prototype and its design features and components to support continued long-term storage or off-site transportation of spent nuclear fuel and high-level radioactive waste currently stored in storage containers. A remote (100 ft away from the simulated radiative environment) demonstration of the VRM was successfully performed on a full-scale mock-up welded stainless steel canister. The VRM is designed with features to accommodate remediation techniques beyond those currently selected and described in this paper. Therefore, many of the VRM’s features may benefit other remote nuclear or nonnuclear applications. The VRM is envisioned to serve as a development center to facilitate and enhance further development of new remediation technologies.
Heavy-duty PEM fuel cells could be a low-cost, low emission alternative to combustion turbines for re-electrifying hydrogen if used as part of a long duration grid energy storage systems. Many studies expect heavy duty PEM fuel cell production costs to reduce as manufacturing volumes ramp up and their expected durability of 25,000-30,000 hours aligns well with a 30-year life for hydrogen seasonal energy storage plants that would likely operate less than 10% of the year. The labor, material, and equipment costs associated with installing PEM fuel cells and their required balance of plant for stationary applications have not been thoroughly explored, however. This study performs a detailed design and cost analysis of a 100 MW stationary PEM fuel cell power plant, capturing costs such as cooling, power electronics, pipes, valves, fittings, cabling, conduit, concrete foundations, buildings, and land. It employs methods consistent with NLR's solar PV benchmarking cost analysis and annual technology baseline to derive the total installed costs of stationary PEM fuel cell plants that utilize heavy duty PEM fuel cells.
In this paper, the development of a fission-gas collecting and physical-analysis-enabling instrument was proposed for small-volume determination. Analysis specifications require a design capable of accurately and repeatably determining volumes in the range of 0.07–2.5 mL. This system relies on a series of gas expansions originating from a cylinder with known internal volume. The combined gas law is used to derive the unknown volumes from these expansions. Initial system designs included one of two known volumes, 11.85 ± 0.34 mL and 5.807 ± 0.078 mL, with a manifold volume of 32 mL. Results obtained from modeling this system’s operation showed that 0.07 mL can be determined with a relative expanded uncertainty greater than 300% (k = 2) for a single replicate, which was unacceptable for the proposed experimental design. Initial modeling showed that the volume connecting the known volume and rodlet, i.e., the manifold volume, and the sensitivity of the pressure sensor were key contributors to the expanded uncertainty of the measured rodlet volume. The system’s design limited the available options for pressure sensors, so emphasis was placed on the design of the manifold volume. The final system design reduced the manifold volume to 17 mL. These changes in design, combined with replicate analysis, were able to reduce the relative expanded uncertainty by ±12% (k = 2) for the 0.07 mL volume.
Recent advancements in the characterization of irradiated U-Pu-Zr fuels have revealed complexities that challenge existing understanding of constituent redistribution. Traditionally, models have proposed three concentric regions within the fuel, each characterized by distinctive phases and porosity. However, through detailed analysis of high burnup U-Pu-Zr, we discovered the presence of four distinct constituent redistribution regions. Particularly novel is the observation of significant Pu redistribution, a previously unreported phenomenon that necessitates a reevaluation of current models. This work aims to delve deeper into these findings, seeking to correlate mesoscale measurements of thermal diffusivity and respective thermal conductivity with the phases present in each redistribution region. To achieve this objective, we employed mesoscale thermoreflectance methods using the unique, Idaho National Laboratory (INL) developed, Thermal Conductivity Microscope (TCM) at INL’s Irradiated Materials Characterization Laboratory. The TCM employs two tightly focused lasers: one for heating to generate periodic thermal waves in the substrate, and another spatially separated probe laser to detect changes in the optical reflectivity of the gold-coated substrate resulting from thermal wave diffusion. We conducted several thermal diffusivity measurements within each region of constituent redistribution of a U-19Pu-10Zr fuel pin cross section irradiated to 11 at. % burnup. The TCM measurement positions strategically aligned with transmission electron microscopy (TEM) lift-out locations previously collected from the fuel sample. Complementary microstructural analysis techniques such as optical and scanning electron microscopy (OM/SEM), electron probe microanalysis for chemical compositions, and TEM-based selective area electron diffraction (SAED) analysis for crystallographic insights into each phase were also utilized. This comprehensive approach allowed us to correlate local thermal diffusivity data with microstructural characteristics, enabling the computation of local thermal conductivity at each position. The significance of this contribution lies in its pioneering use of the TCM for ternary fuel mesoscale examination, shedding light on the previously overlooked effects of Pu redistribution on local thermal conductivity. By informing current models capturing constituent redistribution and heat transfer, our findings pave the way for more accurate predictions of metallic fuel performance. Moreover, this work sets the stage for future comparisons with similar TCM examinations on U-19Pu-10Zr fuels at ultra-low burnup, facilitating a comprehensive understanding of thermal property changes across different burnup levels. Ultimately, our study not only enriches our understanding of the thermophysical properties of individual redistribution regions within U-Pu-Zr fuel but also offers valuable insights for the design and operational parameters of proposed next-generation fast reactors.
A preliminary reference core design was previously developed for the Versatile Test Reactor. The main features are that it is a 300 MW th sodium-cooled fast reactor using ternary metallic fuel U-20Pu-10Zr and able to achieve peak fast fluxes (E n > 0.1 MeV) in excess of 4.0x10 15 n/cm 2 -s. The plutonium in the fuel was assumed to be “reactor grade”, with a fissile quality of about 72% (i.e., 239 Pu and 241 Pu comprise 72% of the plutonium isotopes), and the uranium was assumed to be low-enriched uranium with 5% 235 U. This preliminary core design has been used through the CD-0 phase of the VTR project. Progression of the work after CD-0 led to revising, updating and refining the reference VTR core design, based on considerations from the various VTR teams: fuel, experiment, safety, and plant design teams. The objective of this report is to summarize the various changes made and details added to the model, to provide a description of the updated reference VTR core design and of its performance characteristics. This includes the traditional reactor physics characteristics, fuel cycle details, reactivity coefficients, temperature distributions, control rod worths and shutdown requirements. All of these define the updated “reference VTR core” intended to be used through CD-1.
Idaho National Laboratory has unique opportunities to examine fuels, materials, and experiment inside of the available radiological hot cells. Due to the uninhabitable nature of a hot cell environment, everything done inside of the hot cell is operated remotely. Remote operations are optimized by Nuclear Remote System Design (NRSD). NRSD is the designing, altering, or configuring of structures, items, and systems that will be placed into radiological environments.
Problem: Advanced nuclear reactors come in many sizes, materials, and designs; Criticality experiments are necessary to understand, quantify, and validate these reactors. Approach: The National Criticality Experiments Research Center (NCERC) operated by Los Alamos National Laboratory (LANL) conducts criticality experiments; Deimos is a graphite-moderated, beryllium-reflected, high assay low-enriched uranium (HALEU) tri-structural isotropic (TRISO) fueled experiment; Modeling and simulation is needed to design Deimos.
H2FillS: Hydrogen Filling Simulation The Hydrogen Filling Simulation (H2FillS) software is a thermodynamic model designed to track and report on the transient change in hydrogen temperature, pressure, and mass flow when filling a fuel cell electric vehicle (FCEV). H2FillS simulates gas flow from the hydrogen station to the FCEV storage system. Using empirical fueling data sets, the model has been validated over a range of fueling conditions to match common light-duty FCEV fill profiles. Overall, it provides significant benefits to the light-duty fueling market and fill knowledge gaps of the interaction between a hydrogen station and an FCEV. Capabilities Use the comprehensive H2FillS model to: Safely design and operate a hydrogen fueling station Support code refinement with readily available data Develop system and operational improvements to reduce capital or operating costs at hydrogen stations or on-board FCEVs. How It Works H2FillS uses a "drag-and-drop" graphical user interface to simulate station and vehicle systems with preset parameters for common hydrogen station components as default values. When needed, users can define their own parameters to create their own station or vehicle components. Download H2FillS in Two Steps Read the terms of the end user software license agreement. Register to download the software. There are two versions of the model: A full-station model starts the simulation at high-pressure ground storage, runs through a dispenser, and ends at a vehicle storage system. A partial-station model starts at the dispenser breakaway and consists solely of the dispenser components and the vehicle storage system. H2FillS will automatically output fill performance data from the vehicle by tracking pressure and temperature throughout the fill. Users can input their own fill profiles into the model to run a variety of simulations.
Introduction (120 words): Transport of metal ions across the aqueous-organic phase boundary is an essential step in a hydrometallurgical nuclear fuel reprocessing strategy. The study of transport agents for nuclear fuel elements is imperative to guide the design of ligands that boost the separation efficiency of the recovery process from fission products. However, limited studies have been made on the chemistry of these transport agents when complexing with transuranic elements in gas-phase where all surrounding factors are essentially excluded. This work investigates the reagent ligand complexations to transuranic and other metals and their dissociations in the gas phase. Comparisons are made between 4f and 5f elements and between ligands. Methods (120 words): (N,N-diisobutylcarbamoylmethyl)phenyloctylphosphine oxide (CMPO) and N,N,N',N'-tetraoctyldiglycolamide (TODGA) have been selected to complex with metal nitrates. The actinide americium and lanthanides neodymium, samarium, and europium were investigated as part of this work. The lanthanides were selected to act as size and electron configuration analogues of the minor actinides. Metal complexes with two ligands and two nitrates ([M(NO3)2(CMPO)2]+, for example) are studied in Bruker micrOTOF-Q II mass spectrometer equipped with collision-induced dissociation capability. The comparisons of the mass spectra are made in groups of homogenous ligands and mixed TODGA-CMPO ligands clusters. Comparisons are also made based on the complexed metals (Am and lanthanides). Preliminary data (300 words): Collision-induced dissociation mass spectrometry data are collected on two ligands complexed with metal nitrates where the two ligands are homogenous, with (CMPO)2 or (TODGA)2, or heterogeneous, with (TODGA)(CMPO). Several fragmentation patterns are observed among complexes with the CMPO ligand whereas the TODGA ligand commonly dissociates intact from the complex. Most of the metal complexes exhibit similar fragmentation patterns, but there are a few notable deviations in fragmentation patterns between the Am and Ln-bearing complexes. For the [M(CMPO)2(NO3)2]+ complexes, the initial loss of nitrate in the form of nitric acid is observed in all four complexes. However, [Am(CMPO)2(NO3)2]+ exhibits an additional fragmentation not found in the lanthanide complexes. Also, a significantly different ratio of the second nitric acid loss is found in the Am complex. These deviations may indicate the different interaction behaviors between actinides and lanthanides. The [M(TODGA)2(NO3)2]+ complexes exhibit the fragmentation as the loss of one TODGA ligand as an intact form and the loss of nitrate as nitric acid. The Am complex exhibits an additional fragmentation after losing the TOGDA ligand, which is not observed among the Ln complexes. The heterogeneous [M(TODGA)(CMPO)(NO3)2]+ complexes exhibit both similarities and differences between the Am and Ln complexes. For example, the heterogenous Am complex does not exhibit the loss of an intact TODGA ligand while all three Ln complexes do. This indicates that TODGA may bind more strongly to Am than Ln. Additionally, the intensity of the loss of CMPO ligand (as partially or whole) is found to be significantly larger than that of the loss of TODGA (as partially of whole) indicating that TODGA is bound to the metal significantly stronger than CMPO. Planned computational analysis will help understand the deviation in fragmentation behaviors between americium and lanthanide metal centers, or between TODGA and CMPO ligands. Novel aspect (20 words): Gas-phase actinide and lanthanide complex formation and fragmentation provide insight into the coordination environment differences of f-element metals.
Title (20 word): Comparison of Gas Phase Fragmentation Behaviors of Nuclear Fuel Cycle Ligands in Lanthanide and Americium Metal Ligand Nitrate Clusters Introduction (120 words): Transport of metal ions across the aqueous-organic phase boundary is an essential step in a hydrometallurgical nuclear fuel reprocessing strategy. The study of transport agents for nuclear fuel elements is imperative to guide the design of ligands that boost the separation efficiency of the recovery process from fission products. However, limited studies have been made on the chemistry of these transport agents when complexing with transuranic elements in gas-phase where all surrounding factors are essentially excluded. This work investigates the reagent ligand complexations to transuranic and other metals and their dissociations in the gas phase. Comparisons are made between 4f and 5f elements and between ligands. Methods (120 words): (N,N-diisobutylcarbamoylmethyl)phenyloctylphosphine oxide (CMPO) and N,N,N',N'-tetraoctyldiglycolamide (TODGA) have been selected to complex with metal nitrates. The actinide americium and lanthanides neodymium, samarium, and europium were investigated as part of this work. The lanthanides were selected to act as size and electron configuration analogues of the minor actinides. Metal complexes with two ligands and two nitrates ([M(NO3)2(CMPO)2]+, for example) are studied in Bruker micrOTOF-Q II mass spectrometer equipped with collision-induced dissociation capability. The comparisons of the mass spectra are made in groups of homogenous ligands and mixed TODGA-CMPO ligands clusters. Comparisons are also made based on the complexed metals (Am and lanthanides). Preliminary data (300 words): Collision-induced dissociation mass spectrometry data are collected on two ligands complexed with metal nitrates where the two ligands are homogenous, with (CMPO)2 or (TODGA)2, or heterogeneous, with (TODGA)(CMPO). Several fragmentation patterns are observed among complexes with the CMPO ligand whereas the TODGA ligand commonly dissociates intact from the complex. Most of the metal complexes exhibit similar fragmentation patterns, but there are a few notable deviations in fragmentation patterns between the Am and Ln-bearing complexes. For the [M(CMPO)2(NO3)2]+ complexes, the initial loss of nitrate in the form of nitric acid is observed in all four complexes. However, [Am(CMPO)2(NO3)2]+ exhibits an additional fragmentation not found in the lanthanide complexes. Also, a significantly different ratio of the second nitric acid loss is found in the Am complex. These deviations may indicate the different interaction behaviors between actinides and lanthanides. The [M(TODGA)2(NO3)2]+ complexes exhibit the fragmentation as the loss of one TODGA ligand as an intact form and the loss of nitrate as nitric acid. The Am complex exhibits an additional fragmentation after losing the TOGDA ligand, which is not observed among the Ln complexes. The heterogeneous [M(TODGA)(CMPO)(NO3)2]+ complexes exhibit both similarities and differences between the Am and Ln complexes. For example, the heterogenous Am complex does not exhibit the loss of an intact TODGA ligand while all three Ln complexes do. This indicates that TODGA may bind more strongly to Am than Ln. Additionally, the intensity of the loss of CMPO ligand (as partially or whole) is found to be significantly larger than that of the loss of TODGA (as partially of whole) indicating that TODGA is bound to the metal significantly stronger than CMPO. Planned computational analysis will help understand the deviation in fragmentation behaviors between americium and lanthanide metal centers, or between TODGA and CMPO ligands. Novel aspect (20 words): Gas-phase actinide and lanthanide complex formation and fragmentation provide insight into the coordination environment differences of f-element metals.
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.