699 Search Results

This study numerically investigated the effects of various contact strategies on the thermal hydraulic behavior within a structured bed of 100 explicitly modeled pebbles. Four contact strategies and two thermal hydraulic conditions were considered. The strategies to avoid contact singularities include decreasing the pebble diameter, increasing the pebble diameter, bridging the pebble surfaces near the contact region, and capping the pebble surfaces near the contact region. One strategy, Strategy 3a, which involves bridging with a cylinder equal to 10% of the pebble diameter, was selected as the baseline strategy because it addressed the contact singularity while minimizing the geometric changes that affect the bed porosity. The two thermal hydraulic conditions were full-power operation (Case 1) and pressurized loss of forced cooling or PLOFC (Case 2). Simulations of the conjugate heat transfer within the structured bed were performed using the Reynolds-averaged Navier–Stokes approach with the realizable k-ϵ turbulence model and two-layer all y⁺ wall treatment. The thermal-fluid quantities of interest were compared between the contact strategies for each case. In Case 1, the hydraulic behavior was sensitive to the contact strategy, with large differences in the pressure drop (30%) and volume-average velocity (4%). The thermal behavior was not sensitive, with less than a 0.5% difference across the strategies. To better understand the separate effects of each heat transfer mode, Case 2 was divided into the following subcases: conduction (Case 2a); conduction/radiation (Case 2b); and conduction/radiation/convection (Case 2c). Case 2a represents an early phase of the PLOFC transient. Case 2b represents an intermediate phase of the PLOFC transient, with the pebble temperatures sufficiently high for the radiative heat transfer to be non-negligible. Case 2c represents a late phase of the PLOFC transient after the establishment of the natural circulation of the heat transfer fluid. For Case 2, large differences in the contact strategy were observed only in Case 2a with only conduction. The difference in the maximum pebble temperature was 23% in Case 2a, 2% in Case 2b, and 0.3% in Case 2c.

This study presents a comprehensive benchmarking and verification effort of several thermal-hydraulic and multiphysics capabilities for high-temperature gas-cooled reactor (HTGR) applications. The first part of this effort focuses on the running-in verification of Griffin's multiphysics capabilities, specifically for simulating the evolution of Pebble Bed reactor cores from startup to equilibrium. In the absence of validation data, code-to-code comparisons are conducted with Kugelpy, showing good agreement for key quantities like maximum power density and fresh core k-eigenvalue predictions. However, discrepancies in equilibrium core predictions suggest potential issues with cross sections, underscoring the need for further refinement and evaluation. The HTTF system analysis code benchmark involves RELAP5-3D, SAM, and GAMMA+ to assess their predictive capabilities for HTTF behavior under both normal operation and pressurized conduction cooldown (PCC) transient conditions. While there is good agreement in predicting major parameters such as coolant temperature, solid temperature, and flow distribution, discrepancies in transient behavior highlight differences in modeling approaches, nodalizations, and heat transfer models. The HTTF lower plenum CFD benchmark employs nekRS to simulate flow mixing phenomena, successfully capturing relevant flow physics and demonstrating mesh independence in complex geometries. Preliminary results suggest a relatively uniform temperature field but significant unsteadiness in the flow, requiring time-averaging analyses. The GPBR200 system analysis code benchmark uses SAM's core channel and porous media models, incorporating an RCCS loop for decay heat removal. During steady-state and transient conditions, including protected de-pressurized and pressurized loss of forced cooling (DLOFC and PLOFC), both models show good agreement in predicting temperature profiles and key parameters. Notably, while the core channel model underpredicts convective heat transfer effects, both models maintain temperatures well below the TRISO fuel safety limit. These benchmarking efforts collectively enhance the predictive capabilities of the tools used in HTGR design and safety analysis, guiding developments to improve their accuracy and applicability.

For code development, Advanced Reactor Technologies - Gas Cooled Reactors Program (ART-GCR) rely on a collaboration with the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, but the cross sections generation and the methodology definition is part of this program area goals. Based on previous studies in FY23, the size of microscopic cross section libraries increases rapidly with the number of tabulations, requiring significant amount of memory and drastically slowing down the Griffin calculations when evaluating cross sections via the multivariate linear interpolation approach. Rising to these challenges, this work investigates constructing Reduced-order Models (ROMs) for the multi-group microscopic cross sections to accelerate the cross section evaluation in Griffin. A database of multigroup cross sections is first collected considering all possible parameters that a designer could change for optimization. Down-selection of the ROM techniques afterward shows Deep Neural Network (DNN) as the best candidate when jointly consider memory efficiency, predictive accuracy, computational cost, scalability, flexibility and ease of implementation of the algorithms in comparison to the multidimensional interpolation. This work develops a specific interface that enables the cross section predictions using pre-trained DNN models into Griffin leveraging the existing ROM capabilities. DNNs have been trained for all isotopes for use in Griffin. Preliminary Griffin testing shows that DNNs exhibit exceptional predictive accuracy and the use of DNNs provides orders of magnitude improvement in memory efficiency compared to conventional interpolation techniques. With such ROM techniques, it holds great promise to further increase the fidelity of the Pebble Bed Reactor (PBR) simulation by increasing the number of tabulations/state variables during cross section evaluation, while maintaining the computational cost affordable in Griffin.

Modeling the run-in of a pebble-bed reactor (PBR) can be challenging as a result of changes in the power, fuel type, and temperatures that occur throughout the run-in period. Previous work utilized high-fidelity neutronics simulations or lower-fidelity coupled neutronics/thermal-hydraulics models to capture the general characteristics of the run-in process. Here, the present work employs high-fidelity neutronics simulations (using Serpent) coupled with thermal-hydraulics simulations (using Griffin–Pronghorn) to capture the thermal feedback present during the run-in and approach to equilibrium for a PBR. Incorporating thermal feedback enables important distinctions to be made about conditions occurring inside the core, as the power distribution, discharge burnup, and isotopic compositions are all affected by the temperature distribution.

Griffin is a reactor multiphysics modeling application based on MOOSE (Multiphysics Object-Oriented Simulation Environment) and specifically targeting transient modeling of advanced reactors. Griffin has been used recently to model pebble-bed reactors for the Nuclear Regulatory Commission (NRC) Office of Nuclear Regulatory Research and the Advanced Reactor Technology program. This modeling work has focused thus far on the direct calculation of equilibrium cores. This report documents an initial demonstration of a new running-in simulation capability. The new running-in capability is verified using the existing direct equilibrium core calculation capability. A simplified pebble-bed reactor model is then used to demonstrate the running-in simulation capability. This demonstration shows that Griffin is able to simulate years of operation during the running-in phase efficiently with each depletion step taking only several seconds. Two new technologies are also presented in this report which have been developed in Griffin that will be essential for improved accuracy both of the direct equilibrium core computation and the new running-in simulation capability. The first technology is an online cross section generation capability specifically targeted for pebble-bed reactors. This will improve the accuracy of the depletion calculation as the cross sections are generated at the exact core status. This also avoids the difficult step of pre-generating a separate standalone multigroup cross section set. Secondly, a newly implemented discretization for discontinuous finite element method (DFEM) SN transport in cylindrical (RZ) coordinates, which can be solved efficiently using the existing SN sweep solver, is discussed and some results are shown demonstrating the usefulness of the additional accuracy transport provides over a diffusion approximation.

As the nuclear industry moves towards licensing and constructing advanced reactors, new attention has been focused on the advanced reactor designs that have past operational experience, such as pebble-bed high-temperature gas-cooled reactors (PB-HTGRs). Pebble-bed reactor designs have many advantages, such as their higher operating temperatures and online refueling capabilities. However, high-fidelity computational modeling of pebble-bed reactor designs, from reactor startup to operation at equilibrium, is more challenging compared to conventionally fueled reactors due to the continuous movement of the fuel pebbles through the reactor during operation. In previous work at Oak Ridge National Laboratory (ORNL), the SCALE Leap-In method for Cores at Equilibrium (SLICE) was developed around tools within the SCALE code system. This iterative method can effectively generate pebble-bed reactor zone-wise fuel inventories at equilibrium core operation within a reasonable computational time. The objective of this work was to further verify the ORNL SLICE method and to investigate the impact of considering temperature profiles during the application of the method. The SLICE method was applied to a modular high-temperature gas-cooled reactor design based upon publicly available design specifications of the Xe-100 pebble-bed reactor. Upon comparing the results from the SLICE method to published literature, the differences in the eigenvalue k_effective were on the order of several hundred pcm (percent millirho). To investigate one possible cause of these differences, a study looking at the sensitivity of the full-core equilibrium k_effective and discharge nuclide inventory to temperature was performed by developing equilibrium cores of two additional temperature profiles. From this temperature study, differences on the order of hundreds of pcm for the full-core equilibrium k_effective, and up to 15% difference for the discharge inventories were found. In conclusion, these results indicated the strong dependence on temperature that needs to be considered for future work in equilibrium modeling of PB-HTGRs.

Pebble-bed reactors are an important class of advanced reactors under consideration for various applications where their fuel would give a significant advantage in siting and high-quality heat production. However, the disposal of their fuel is not as thoroughly studied as other fuel forms. This article provides an example of evaluating advanced reactor spent nuclear fuel in MOOSE. In this study, pebble fuel is analyzed in a well-known spent fuel canister design to characterize the behavior of this fuel form over a long time scale. The results indicate that after approximately 100 years, decay heat is significantly reduced and the maximum temperature in the canister equalizes with the external temperature. The simulation goes on to an end time of one million years, demonstrating the capability of efficiently dealing with long time scales efficiently. We conclude that the canister temperatures seem manageable even with very aggressive loading times and while there are several improvements to be implemented in the future, MOOSE is currently capable of simulating the required scenarios.

Griffin, a MOOSE (Multiphysics Object-Oriented Simulation Environment) based application targeting transient modelling of advanced reactors, has been used recently to model pebble-bed reactors (PBRs). The modelling effort has focused thus far on modelling the equilibrium core. A new capability to simulate the running-in phase of PBR operation has been added to Griffin. This work demonstrates the newcapability with a sample multiphysics running-in simulation. The basic features of the new running-in capability were documented previously; however, the sample simulation results presented there did not include multiphysics; the fuel temperatures were assumed to be constant. In this work, Griffin computes power densities in the core at each timestep of the running-in simulation and passes these to Pronghorn which models fluid flow and heat transfer to calculate temperatures that are passed back to Griffin and accounted for with temperature dependent cross-sections.

Pebble-bed reactors are an important class of advanced reactors under consideration for various applications where their fuel would give a significant advantage in siting and high-quality heat production. However, the disposal of their fuel is not as thoroughly studied as other fuel forms. In this study, pebble fuel is analysed in a well known spent fuel canister design to characterize the behavior of this fuel form over a long time scale. The results indicate that after approximately 100 years, decay heat is significantly reduced and the maximum temperature in the canister equalizes with the external temperature. The simulation goes on to an end time of a million years, demonstrating the capability of dealing with long time scales efficiently. We conclude that the canister temperatures seem manageable even with very aggressive loading times and while there are several improvements to be implemented in the future, MOOSE is technically capable of simulating the required scenarios.

Pebble-bed reactors (PBRs) rely on a continual feed of fuel pebbles being cycled through the core. As a result, they require a “run-in” period in order to reach an equilibrium state. The run-in period for a PBR is a complex, time-dependent problem that requires the injection of new fuel, different types of fuel, and power increases. This complexity in the run-in makes it important to capture the physical processes in order to generate an accurate representation. The present work details the creation of a high-fidelity Monte Carlo methodology for analyzing the run-in and subsequent approach to equilibrium for PBRs. The methodology entails a Python module wrapped around Serpent so as to perform neutronics calculations, move pebbles, refuel the core, and discharge pebbles, thereby modeling the explicit behavior of the PBR run-in. Further, three run-in simulations (a constant temperature profile, a linear temperature profile, and a constant temperature profile using control rods) were examined in order to identify the key physical phenomena present in the run-in process. Utilizing kugelpy, we found the inclusion of a temperature profile to be important for accurately capturing a discharge burnup (around 141 MWd/kg), a consistent k-eff (around 1.005), and an average pebble power (around 2.5 kW/pebble) that all fall within acceptable limits.