Engineering topics
Publications and source records attributed to Tano Retamales, Mauricio Eduardo.
The MOOSE Navier–Stokes module solves mass, momentum, energy, and passive scalar conservation equations in the context of fluid flow. The module supports solution of these equations in both free flow and porous medium contexts and for a range of fluid compressibility. The conservation equations can be discretized in space using continuous Galerkin finite elements or with cell centered finite volumes.
A Subchannel application (Pronghorn-SC) is developed in MOOSE, which affords the required flow field resolution, while still preserving an engineering-scale approach. This new solver can be natively coupled to Pronghorn and other MOOSE objects to enable full-core, multi-physics, multi-scale engineering studies. This work utilizes the EBR-II SHRT tests to validate the subchannel capabilities. Multi-scale and multi-physics coupling is used to improve the fidelity of the subchannel code calculations.
The National Reactor Innovation Center (NRIC) was created for accelerating the deployment of novel reactor concepts. This will be achieved by providing both physical and virtual spaces for building and testing various components, systems, and complete pilot plants. The Virtual Test Bed (VTB) represents the virtual arm. It is being developed in collaboration with Department of Energy’s (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program.
This report describes the implementation and testing of domain-overlapping coupling of the Pronghorn and SAM codes. Both Pronghorn and SAM are codes developed by the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program for analysis of advanced nuclear reactors. Pronghorn focuses on analyzing multi-dimensional core flow conditions, while SAM focuses on analysis of the entire system including, among other components, piping, pumps, and heat exchangers. It is desirable for many advanced reactor thermal-hydraulics simulations to couple Pronghorn and SAM to obtain self-consistent solutions in the core and the system. The domain-overlapping coupling approach provides a robust and numerically efficient candidate for this coupling. We describe the algorithm for obtaining self-consistent solutions between a SAM network of flow channels and an overlapped multi-dimensional Pronghorn domain. The algorithm is tested for a variety of simple canonical test problems and a model of the Molten Salt Reactor Experiment (MSRE). Results demonstrate that the domain overlapping coupling approach provides consistent pressure, mass flow rate, enthalpy, and passive scalar distributions between Pronghorn and SAM and converges reliably and efficiently (at most 15 iterations, usually much less than 10) for the considered test problems.
Pronghorn is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based thermal-hydraulics code developed at Idaho National Laboratory (INL) for advanced reactor analysis. It has been previously applied to model Pebble-Bed High Temperature Reactors (High Temperature Reactor (HTR)s), Liquid-Metal Cooled Reactors, and Molten Salt Reactors, among others. This work applies the coarse-mesh thermal hydraulics capabilities in Pronghorn to model Prismatic-Core HTRs. In particular, the Oregon State University (OSU)’s High Temperature Test Facility (HTTF) is modeled with Pronghorn. The HTTF is a 1:4 height scaled-down facility of General Atomics’ Modular High Temperature Gas-cooled Reactor.
1.1 OBJECTIVE The United State High Performance Research Reactor (USHPRR) program aims to eliminate more than 200kg of High Enriched Uranium (HEU) from commerce annually by converting five U.S. high-performance research reactors and one associated critical assembly to Low Enriched Uranium (LEU) fuel using a high-density alloy of uranium-10 wt% molybdenum (U-10Mo). The National Bureau of Standards Reactor (NBSR) is one of five research reactors selected for this program. The objective of this report is to provide preliminary thermal-hydraulic and mechanical analyses of the hydrodynamic effects in the NBSR Design Demonstration Element (DDE) under conservative approximations for the plate power distribution. This report provides details on the modeling approach and the simulation results obtained, including pressure, flow velocity, temperature, and oxide layer over the design demonstration experiment for the irradiation cycles in the Belgian Reactor (BR)-2.
This work reports the first step in such an analysis, the verification of isotope removal capability for multi-region MSR depletion in Griffin. Verification studies are performed against SERPENT for chemical species removal from the core to an off-gas system for MSR analysis. A brief review of the Bateman equations including the isotopic removal term is provided along with a description of the example problem used in the verification studies performed in Griffin versus SERPENT. Next, the results of the verification study are presented for isotopic depletion, isotopic depletion with removal, and the specific analysis case of varying the solubility of Iodine in a salt system and how this can affect material buildup in the off-gas system. Lastly, the results are analyzed and future work detailing the fully coupled MSR depletion analysis that incorporates the true multiphysics of the system is discussed.
Earthquakes affect numerous nuclear power plants around the world and safe earthquake design requires that the physics of reactors subjected to earthquakes is well understood. This project attempts to understand the physics of molten salt reactors (MSRs) subjected to earthquakes using coupled multiphysics tools in MOOSE. A beyond-design-basis scenario of an unprotected pump trip (ULOFA) during an earthquake is analyzed.
The VTB hosts 20 distinct models, 8 of which were recently developed during 2022 and will be discussed in more detail here. This article provides an overview of the various models hosted and their capabilities. Part of the intent of the VTB is to showcase a ‘block’ style approach to modeling and simulation: if a capability is not showcased for a specific reactor type, it could, in theory, be easily ported/replicated from another. The breadth of advanced reactor types considered is representative of the large variety of concepts being proposed for demonstration.
Pool-type Molten Salt Reactors (MSRs) are a promising advanced nuclear reactor concept relying on passive physical behavior to provide increased safety characteristics. The high heat capacity, high boiling point, unpressurized liquid salt contains the nuclear fuel and as it passes through the vessel achieves criticality and produces nuclear power. The salt then carries away the heat produced along with the delayed neutron precursors towards heat exchangers and primary pumps before entering the core. This forms a highly coupled multiphysics problem between neutronics and fluid flow, which makes modeling these reactors challenging. The safety/licensing case heavily relies on modeling and simulation to prove their safety, and significant effort has been expanded by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to develop the appropriate simulation tools. Their safety/licensing will require a thorough analysis of their behavior during earthquake-based transients. Earthquakes are a transient condition commonly experienced in many regions, and are safely survived by dozens of reactors every year. Nonetheless, significant conservatisms were introduced with limited predictive modeling and simulation available. The first-of-a-kind capability developed by this research will allow for engineers to analyze ranges of design-basis and beyond-design basis accident scenarios to assess the integrity of the reactor vessel and surrounding system, possibly allowing for limiting the costly conservatisms in the design. This projects uses MASTODON, Griffin, and Pronghorn in the NEAMS ecosystem to model the coupled multiphysics problem posed by an earthquake in a molten salt reactor. It marches through a number of number of coupling schemes, from uncoupled simulations to tightly two-way coupled simulations. The appropriate level of coupling for these simulations will then be determined based on the accuracy in the quantities of interest (QoI) and computational effort involved in each scheme. The QoIs chosen for this paper include the mass flow rate across the heat exchanger (for Pronghorn simulations), power output of the reactor (for Griffin simulations), and the maximum Von Mises stress (for MASTODON simulations). This summary paper reports on the current progress of this project, with 2D tightly coupled multiphysics results.
The United State High Performance Research Reactor (USHPRR) program aims to eliminate more than 200kg of High Enriched Uranium (HEU) from commerce annually by converting five U.S. high-performance research reactors and one associated critical assembly to Low Enriched Uranium (LEU) fuel using a high-density alloy of uranium-10 wt% molybdenum (U-10Mo). Massachusetts Institute of Technology Reactor (MITR) is one of five research reactors selected for this program. The objective of this report is to provide preliminary thermal-hydraulic and mechanical analyses of the hydrodynamic effects in the MITR Design Demonstration Element (DDE) under conservative approximations for the plate power distribution. This report provides details on the modeling approach and the simulation results obtained, including pressure, flow velocity, temperature, and oxide layer over the design demonstration experiment for the irradiation cycles in the Belgian Reactor (BR)-2.
The United State High Performance Research Reactor (USHPRR) program aims to eliminate more than 200kg of High Enriched Uranium (HEU) from commerce annually by converting five U.S. high-performance research reactors and one associated critical assembly to Low Enriched Uranium (LEU) fuel using a high-density alloy of uranium-10 wt% molybdenum (U-10Mo). The Massachusetts Institute of Technology Reactor (MITR) and the National Bureau of Standards Reactor (NBSR) are two of five research reactors selected for this program. Previous studies have addressed the thermal-hydraulics performance of the fuel Design Demonstration Elements (DDE) for these reactors under conservative assumptions in the Belgian Reactor (BR)-2. This report extends these previous studies by analyzing the margins to Onset of Nucleate Boiling (ONB) for these two DDEs under the conservative conditions used in the previous thermal-hydraulics analyses
Earthquakes affect numerous nuclear power plants around the world and safe earthquake design requires that the physics of reactors subjected to earthquakes is well understood. This project attempts to understand the physics of molten salt reactors (MSRs) subjected to earthquakes using coupled multiphysics tools in MOOSE. A beyond-design-basis scenario of an unprotected pump trip (ULOFA) during an earthquake is analyzed.
This report details the improvement of pebble bed reactor modeling capabilities in the Pronghorn code in fiscal year 2022. The following accomplishments are reported: Deployment of weakly compressible finite volume formulation to the HTR- PM reference plan model; Enable modeling of stagnant gas gaps in the finite volume formulation; Enable using all Pronghorn correlations available in the finite element version in the finite volume version; Modeling of decay heat in pebble bed reactors; Simplifying the input for multiphysics equilibrium core calculations and significant reduction of execution time; Implementation of advanced correlations developed by the Center of Excellence for Thermal-Fluids Applications in Nuclear Energy . In addition, this report includes a development plan for Pronghorn and associated NEAMS tools for prismatic gas-cooled reactors.
This report details the development and demonstration of an entirely new capability in Pronghorn, namely the ability to model liquid-metal fast reactor (LMFR) flow conditions on the engineering scale. We developed two modeling approaches for LMFR that can be used separately or be combined into a hybrid simulation: (1) a modern subchannel capability called Pronghorn-Subchannel for square and hexagonal lattices, and (2) a porous flow capability for LMFR geometries. The report emphasizes the novel aspects of the developed subchannel capability and the interoperability of the subchannel capability, porous flow capability, and multiphysics tools within the multiphysics object oriented simulation environment (MOOSE). Here we demonstrate the ability to: (1) Accurately model subchannel flow in hexagonal lattices; (2) Couple the subchannel flow model to multidimensional finite-element method (FEM) or finite-volume method (FVM) heat conduction models; (3) Model LMFRs using Pronghorn’s porous media FVM approach; (4) Couple porous flow FVM and subchannel models in a single simulation; (5) Explicitly model inter-wrapper flows along with conjugate heat transfer from the intra-element flow; and (6) Demonstrate the numerical robustness of the subchannel algorithm by simulating intra-element flow recirculation in a high-buoyancy, low-flow fuel element.
The Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims at developing a simulation tool kit to accelerate the development and deployment of nuclear power technologies. NEAMS multiphysics tools have been designed to provide numerical simulation support for the design and licensing of GEN IV reactors. Pronghorn is NEAMS’s coarse-mesh computational fluid dynamics (CFD) tool, which is designed to run 3D core transients in GEN IV reactors at a reduced computational cost. To increase their accuracy, coarse-mesh CFD simulations require calibrated closure coefficients. One way of computing these coefficients is via the Nek5000, NEAMS’s high-fidelity CFD tool. This article discusses our current research lines in informing Pronghorn closure coefficients via Nek5000 to enable multiphysics simulations of the core cavity of the molten-salt fast reactor. We present an application in which Pronghorn mixing length turbulent viscosity has been calibrated from Nek5000 simulations. The resulting Pronghorn thermal-hydraulics model is then coupled to Griffin, the NEAMS neutron transport solver, to solve for the steady-state configuration of the molten-salt fast reactor.
During the fiscal year of 2021, the finite-volume method (FVM) was deployed in Pronghorn to provide improved efficiency, stability, and accuracy for coarse-mesh, thermal-hydraulics problems. While the main goals of fiscal year 2021 were met, several issues emerged from the early deployment of the finite-volume method in Pronghorn. These issues were: The compressible and incompressible formulations that were implemented are inadequate for many nuclear reactor flow problems; Omission of terms accounting for the porosity and Darcy-Forchheimer body force discontinuities in the Rhie-Chow interpolation lead to oscillations in pressure and velocity at these discontinuities; The FVM lacks a correction for non-orthogonal grids for computing accurate pressure gradient leading to loss of accuracy in regions with skewed elements; The FVM currently uses a monolithic solver. Monolithic solvers have issues dealing with the saddle-point nature of the discretized fluid equations. The result are bad convergence if direct factorization is not used and large memory consumption when direct factorization is used. These four issues are addressed in this report. In particular, we report the completion of the following task: Implementation of a weakly compressible formulation in the MOOSE Navier-Stokes module; Implementation of Moukalled’s method for including body forces in the RCI. Additionally, we identified the need to smooth the porosity using Moukalled’s face-cell smoothing operator; Implementation of a non-orthogonal correction for the Green-Gauss gradient computation; Preliminary implementation of a SIMPLE segregated solver.