Software Requirements Specification for RELAP-7
This document contains the software requirements for the RELAP-7 code.
97 MATHEMATICS AND COMPUTING↗
SEARCH · Engineering Papers
Results for “RELAP”
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.
Modeling a SERTTA-type vessel using RELAP-7
A transient in a SERTTA vessel geometry was modeled using RELAP-7. This transient included a significant amount of boiling in the fuel region, and a resulting transport of vapor into a large expansion tank. A preliminary gas mixture model and transport PDE was used for the presence of non-condensable gases. Results demonstrated that the solution algorithm and/or mixture model needs modification to avoid evaluation of unphysical mixture fluid properties.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Non-Condensable Gases in RELAP-7
This report considers how the model PDEs in RELAP-7 can be amended to account for the presence of non-condensable gases (NCGs) in the flow system. Let the number of NCGs be denoted by $N_{NCG}$. The primary (and condensable) vapor, e.g., steam, will be referred to as the “primary vapor” and will have the subscript “$v$; 0”, and the other constituents are assumed to be non-condensable and will be referred to as “NCGs” and will use the subscript “$v; i$”, with $i$ = 1$,...,N_{NCG}$. Mixture quantities will use the subscript “$v$” (without any numeric index). Before distinguishing between condensable and non-condensable gases, the mixture equations are derived generally; there the mixture components are just referred to as “species” and will just have a subscript “$i$” instead of the corresponding vapor subscript “$v; i$”. The total number of species is denoted by $N_{species}$.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
RELAP-7 Theory Manual
This document summarizes the physical models and mathematical formulations used in the RELAP-7 code.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Sockeye Heat Pipe Code Theory Development: Based on the 7-Equation, Two-Phase Flow Model of RELAP-7
Sockeye is a MOOSE-based computer code designed to simulate the thermal-hydraulic two-phase flows of heat pipes in a nuclear micro-reactor. For this application, a significantly modified version of the fully non-equilibrium, compressible, 7-equation, two-phase flow model of RELAP-7 was developed. This summary report documents the significant points of this model development.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
MULTIPHYSICS STEADY-STATE SIMULATION OF THE HIGH TEMPERATURE TEST REACTOR WITH MAMMOTH, BISON AND RELAP-7
In this work we present a new approach to solving a multiphysics steady-state simulation for the High Temperature Test Reactor (HTTR) using one Super Homogenization (SPH) corrected homogenized 3-D neutronics model in MAMMOTH, multiple heterogeneous 2-D axisymmetric fuel pin models in BISON and several 1-D annular flow thermohydraulics models in RELAP-7. The coupling of the different physics is complicated by the fundamentally disparate models, potentially affecting the stability of the global system. A pseudo-transient to converge the temperature fields is used to circumvent this difficulty. The benefit of doing so as opposed to having fully 3-D heterogenous models for all physics is that the computational burden for reactors as complex as HTTR may quickly become unbearable while coarse models may not capture all the essential phenomena of the problem. The driving purpose of this work is to establish a self-sustaining steady-state solution that can be used to start safety transients.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
MULTIPHYSICS STEADY-STATE SIMULATION OF THE HIGH TEMPERATURE TEST REACTOR WITH MAMMOTH, BISON AND RELAP-7
In this work we present a new approach to solving a multiphysics steady-state simulation for the High Temperature Test Reactor (HTTR) using one Super Homogenization (SPH) corrected homogenized 3-D neutronics model in MAMMOTH, multiple heterogeneous 2-D axisymmetric fuel pin models in BISON and several 1-D annular flow thermohydraulics models in RELAP-7. The coupling of the different physics is complicated by the fundamentally disparate models, potentially affecting the stability of the global system. A pseudo-transient to converge the temperature fields is used to circumvent this difficulty. The benefit of doing so as opposed to having fully 3-D heterogenous models for all physics is that the computational burden for reactors as complex as HTTR may quickly become unbearable while coarse models may not capture all the essential phenomena of the problem. The driving purpose of this work is to establish a self-sustaining steady-state solution that can be used to start safety transients.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
MULTIPHYSICS STEADY-STATE SIMULATION OF THE HIGH TEMPERATURE TEST REACTOR WITH MAMMOTH, BISON AND RELAP-7
In this work we present a new approach to solving a multiphysics steady-state simulation for the High Temperature Test Reactor (HTTR) using one Super Homogenization (SPH) corrected homogenized 3-D neutronics model in MAMMOTH, multiple heterogeneous 2-D axisymmetric fuel pin models in BISON and several 1-D annular flow thermohydraulics models in RELAP-7. The coupling of the different physics is complicated by the fundamentally disparate models, potentially affecting the stability of the global system. A pseudo-transient to converge the temperature fields is used to circumvent this difficulty. The benefit of doing so as opposed to having fully 3-D heterogenous models for all physics is that the computational burden for reactors as complex as HTTR may quickly become unbearable while coarse models may not capture all the essential phenomena of the problem. The driving purpose of this work is to establish a self-sustaining steady-state solution that can be used to start safety transients.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Advanced probabilistic risk analysis using RAVEN and RELAP-7
RAVEN, under the support of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program [1], is advancing its capability to perform statistical analyses of stochastic dynamic systems. This is aligned with its mission to provide the tools needed by the Risk Informed Safety Margin Characterization (RISMC) path-lead [2] under the Department Of Energy (DOE) Light Water Reactor Sustainability program [3]. In particular this task is focused on the synergetic development with the RELAP-7 [4] code to advance the state of the art on the safety analysis of nuclear power plants (NPP). The investigation of the probabilistic evolution of accident scenarios for a complex system such as a nuclear power plant is not a trivial challenge. The complexity of the system to be modeled leads to demanding computational requirements even to simulate one of the many possible evolutions of an accident scenario (tens of CPU/hour). At the same time, the probabilistic analysis requires thousands of runs to investigate outcomes characterized by low probability and severe consequence (tail problem). The milestone reported in June of 2013 [5] described the capability of RAVEN to implement complex control logic and provide an adequate support for the exploration of the probabilistic space using a Monte Carlo sampling strategy. Unfortunately the Monte Carlo approach is ineffective with a problem of this complexity. In the following year of development, the RAVEN code has been extended with more sophisticated sampling strategies (grids, Latin Hypercube, and adaptive sampling). This milestone report illustrates the effectiveness of those methodologies in performing the assessment of the probability of core damage following the onset of a Station Black Out (SBO) situation in a boiling water reactor (BWR). The first part of the report provides an overview of the available probabilistic analysis capabilities, ranging from the different types of distributions available, possible sampling strategies, and post processing analysis capabilities. The first part of the report provides an extensive description of two major developments introduced this year: adaptive sampling for limit surface sampling and multi variate distributions. The document concludes with a description of the demo case (BWR-SBO) and a discussion of the results obtained.
97 MATHEMATICS AND COMPUTING↗
RELAP-7 Progress Report: A Mathematical Model for 1-D Compressible, Single-Phase Flow Through a Branching Junction
In the literature, the abundance of pipe network junction models, as well as inclusion of dissipative losses between connected pipes with loss coefficients, has been treated using the incompressible flow assumption of constant density. This approach is fundamentally, physically wrong for compressible flow with density change. This report introduces a mathematical modeling approach for general junctions in piping network systems for which the transient flows are compressible and single-phase. The junction could be as simple as a 1-pipe input and 1-pipe output with differing pipe cross-sectional areas for which a dissipative loss is necessary, or it could include an active component, between an inlet pipe and an outlet pipe, such as a pump or turbine. In this report, discussion will be limited to the former. A more general branching junction connecting an arbitrary number of pipes with transient, 1-D compressible single-phase flows is also presented. These models will be developed in a manner consistent with the use of a general equation of state like, for example, the recent Spline-Based Table Look-up method [1] for incorporating the IAPWS-95 formulation [2] to give accurate and efficient calculations for properties for water and steam with RELAP-7 [3].
42 ENGINEERING↗
Validation of RELAP-7 for Forced Convective and Natural Circulation Reactor Flows
New reactor designs such as small modular reactors and newly proposed passive safety systems rely on natural circulation two-phase flow cooling of the reactor core under accident conditions. Without a pump to drive the flow, the circulation rate is a strong function of void fraction. Under these conditions where the flow rate, pressure drop, and heat transfer are all closely coupled through the mixture density, the accurate prediction of void fraction is pivotal. Furthermore, while these passive systems boast a more reliable cooling option, the sensitivity of the system to two-phase flow instabilities is more likely than in forced convective cooling and can result in dangerous density wave or flow excursion conditions. Improvement of one-dimensional system codes and further development of three-dimensional code capability are all important efforts in addressing these issues. The proposed work is a critical step towards providing confidence through validation of the capability of thermal-hydraulic tools to meet the emerging challenges in reactor simulation. Historically, system analysis codes develop a stiffness over time due to fine-tuning of adjustable code parameters that result in a reluctance or inability to adopt improved models. Therefore, a rigorous effort should be invested in preparation for RELAP-7 to identify required modeling improvements of constitutive relations and validation of the code applicability to the wide range of possible reactor flow conditions. Part of this validation involves the uncertainty quantification which is a critical step in the Code Scaling, Applicability, and Uncertainty (CSAU) methodology and Best-Estimate Plus Uncertainty (BEPU) licensing method. BEPU allows recovery of safety margins and potentially optimizing plant performance, resulting in economic benefit for the reactor operators. Beginning with the collection of past experimental data covering prototypic LWR conditions, this document reviews three critical datasets available for the validation of system codes. The experimental methods, procedures, facilities, instrumentation, and supplementary information are consolidated to aid the validation efforts of this work and future projects. These three datasets (824 conditions) are chosen due to their wide acceptance and demonstrated application for code validation. Two experimental datasets from the Nuclear Power Engineering Corporation (NUPEC) span prototypic conditions of Boiling Water Reactors (BWR) and Pressurized Water Reactors (PWR). These datasets are labeled as the BFBT (BWR Full-size Fine-mesh Bundle Test) and the PSBT (PWR Subchannel and Bundle Test) datasets respectively. The third dataset is taken in an internally heated annulus that is scaled based on a BWR subassembly and includes two-phase flows important to post-blowdown analysis. Together, the datasets consider the full range of pressure and mass flux expected in LWR forced convective flows from steady-state to long term cooling conditions. Recognizing the deficiency of available data under low-pressure, low-flow, a new natural circulation dataset consisting of 107 conditions is obtained. With system pressures ranging between 145 and 950 kPa, natural circulation is driven by a heat flux up to 275 kW/m2 in a 3 m long heated region of the annulus test section which is followed by a 2 m long unheated region before entering a condenser. Stable natural circulation flow rates in this facility were measured to be between 182 and 590 kg/m2-s. Pressure and temperature are measured at ten locations in the facility. In addition, measurements of void fraction, interfacial area concentration, and gas velocity are obtained at five axial locations in the annulus test section. The RELAP5-3D code demonstrates good accuracy for void fraction predictions in BFBT and PSBT benchmarks, but shows relatively large discrepancy in the annulus experiments at lower pressures. In two-phase natural circulation, large oscillations are observed in many system parameters during the simulation of 19 experiment cases which were conducted at low pressures, and the model shows good predictive performance with the remaining 88 experiment cases. However, some discrepancies in void fraction and temperature predictions are found which may require improved constitutive modeling under these low-pressure, low-flow condition. The quantified uncertainty in void fraction predictions in three forced convection benchmarks due to uncertainties in boundary conditions and geometries are relatively small, except for a few cases in the annulus experiment. The output uncertainty of the natural circulation model is quantified by considering the measurement error in boundary conditions, and it is shown that the output uncertainty ranges are relatively small for the range of conditions considered. Through this work, a clear pathway is provided for validation of RELAP7 once the code is available.
42 ENGINEERING↗
Coupling of Pronghorn and RELAP-7 for a Pebble Bed Reactor
High temperature gas cooled reactors (HTGR) are a candidate for timely Gen-IV reactor technology deployment because of high technology readiness and walk-away safety. Among HTGRs, pebble bed reactors (PBRs) have attractive features such as low excess reactivity and online refueling. Pebble bed reactors pose unique challenges to analysts and reactor designers such as continuous burnup distribution depending on pebble motion and recirculation, radiative heat transfer across a variety of gas-filled gaps, and long design basis transients such as pressurized and depressurized loss of forced circulation. Modeling and simulation is essential for both the PBR’s safety case and design process. In order to verify and validate the new generation codes the Nuclear Energy Agency (NEA) Data bank provide a set of benchmarks data together with solutions calculated by the participants using the state of the art codes of that time. An important milestone to test the new PBR simulation codes is the OECD NEA PBMR-400 benchmark which includes thermal hydraulic and neutron kinetic standalone exercises as well as coupled exercises and transients scenarios. In this work, the reactor multiphysics code MAMMOTH and the thermal hydraulics code Pronghorn, both developed by the Idaho National Laboratory (INL) within the multiphysics object-oriented simulation environment (MOOSE), have been used to solve Phase 1 exercises 1 and 2 of the PBMR-400 benchmark. The steady state results are in agreement with the other participants’ solutions demonstrating the adequacy of MAMMOTH and Pronghorn for simulating PBRs.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Uncertainty analysis of Transportable Fluoride-salt-cooled High-temperature Reactor (TFHR) using coupled DAKOTA with RELAP-3D method
Not Available
Wang, Chenglong (ORCID:000000025351967X)↗
Relap-7: Advanced Nuclear Reactor Systems Analysis Software Package
Advanced Nuclear Reactor Systems Analysis Software Package
Martineau, RichardC.↗
Validation of Thermohydraulic Simulations using RELAP for Critical Dual Purpose Canisters
Simulations using the RELAP5-3D systems analysis code and STAR-CCM+, a high-fidelity CFD code, were performed on a model of a real-scale MPC-32 dual purpose canister undergoing a criticality event. These simulations were performed to provide code-to-code verification of the natural convection capabilities in RELAP5. The STAR-CCM+ simulations revealed various natural convection flow patterns that aided in the cooling of fission heat produced in the fuel rods in the canister. Negligible spatial temperature variations were observed across the canister in the STAR-CCM+ results, supporting the use of RELAP5 in predicting the sub-cooled regime. The surface-averaged rod temperature predicted by RELAP5 compared well with STAR-CCM+, verifying the predictive capability of RELAP5 for sub-cooled conditions.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Reactor core isolation cooling system model development and implementation in RELAP-7
Submitting a poster for the LDRD poster session. Please see the attachment for details.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
THE APPLICATION OF MAMMOTH FOR A DETAILED TIGHTLY COUPLED FUEL PIN SIMULATION WITH A STATION BLACKOUT
Accurate calculation of desired quantities to predict fuel behavior requires the solution of interlinked equations representing different physics. Traditional fuels performance codes often rely on internal empirical models for the pin power density and a simplified boundary condition on the cladding edge. These simplifications are performed because of the difficulty of coupling applications or codes on differing domains and mapping the required data. To demonstrate an approach closer to first principles, the neutronics application Rattlesnake and the thermal hydraulics application RELAP-7 were coupled to the fuels performance application BISON under the master application MAMMOTH. A single fuel pin was modeled based on the dimensions of a Westinghouse 17x17 fuel rod. The simulation consisted of a depletion period of 1343 days, roughly equal to three full operating cycles, followed by a station blackout (SBO) event. The fuel rod was depleted for 1343 days for a near constant total power loading of 65.81 kW. After 1343 days the fission power was reduced to zero (simulating a reactor shut-down). Decay heat calculations provided the time-varying energy source after this time. For this problem, Rattlesnake, BISON, and RELAP-7 are coupled under MAMMOTH in a split operator approach. Each system solves its physics on a separate mesh and, for RELAP-7 and BISON, on only a subset of the full problem domain. Rattlesnake solves the neutronics over the whole domain that includes the fuel, cladding, gaps, water, and top and bottom rod holders. Here BISON is applied to the fuel and cladding with a 2D axi-symmetric domain, and RELAP-7 is applied to the flow of the circular outer water channel with a set of 1D flow equations. The mesh on the Rattlesnake side can either be 3D (for low order transport) or 2D (for diffusion). BISON has a matching ring structure mesh for the fuel so both the power density and local burn up are copied accurately from Rattlesnake. At each depletion time step, Rattlesnake calculates a power density, fission density rate, burn-up distribution and fast flux based on the current water density and fuel temperature. These are then mapped to the BISON mesh for a fuels performance solve. BISON calculates the fuel temperature and cladding surface temperature based upon the current power density and bulk fluid temperature. RELAP-7 then calculates the fluid temperature, water density fraction and water phase velocity based upon the cladding surface temperature. The fuel temperature and the fluid density are then passed back to Rattlesnake for another neutronics calculation. Six Picard or fixed-point style iterations are preformed in this manner to obtain consistent tightly coupled and stable results. For this paper a set of results from the detailed calculation are provided for both during depletion and the SBO event. We demonstrate that a detailed calculation closer to first principles can be done under MAMMOTH between different applications on differing domains.
36 MATERIALS SCIENCE↗
The application of MAMMOTH for a detailed tightly coupled fuel pin simulation with a station blackout
The accurate calculation of desired quantities to predict fuel behavior in a nuclear reactor requires the solution of inter linked equations representing different physics phenomena. Traditional fuels performance codes often rely on internal empirical models for the pin power density and a simplified boundary condition on the cladding outer surface. These simplifications are performed because of the difficulty of coupling applications or codes on differing domains and mapping the required data. In this work the MAMMOTH application is employed in the demonstration of an approach that more consistent with first principles, in a simulation that couples the neutronics application Rattlesnake, the thermal fluids application RELAP-7 and the fuels performance application BISON. A single fuel pin was modeled based on the dimensions of a fuel rod from a Westinghouse 17*17 optimized PWR fuel assembly. The simulation consisted of a depletion period of 1343 days, representing three full operating cycles, followed by a station blackout (SBO) event. The fuel rod was depleted with a total pin power of 65.81 kW. After 1343 days the fission power was reduced to zero (simulating a reactor shutdown). Decay heat calculations provided the time-varying energy source after this time. For this problem, Rattlesnake, BISON, and RELAP-7 are coupled under MAMMOTH in a tightly coupled approach. Each system solves its physics on a separate mesh and, for RELAP-7 and BISON, on only a subset of the full problem domain. Rattlesnake solves the neutronics over the whole domain, which includes the fuel, cladding, gaps, water, and top and bottom rod holders. Here BISON is applied to the fuel and cladding with a 2-D axisymmetric domain, and RELAP-7 is applied to the flow of the circular outer water channel with a set of 1-D flow equations. For this work, the mesh on the Rattlesnake side is either be 3-D (for low order PN transport or diffusion) or 2-D axisymmetric (for diffusion). BISON has a matching ring structure mesh for the fuel so both the power density and local burnup are transferred accurately from Rattlesnake. At each depletion time step, Rattlesnake calculates a power density, fission density rate, burnup distribution and fast flux based on the current water density and fuel temperature. These are then mapped to the BISON mesh to perform a fuels performance solution. BISON calculates the fuel temperature and cladding surface temperature based upon the current power density and bulk fluid temperature. RELAP-7 then calculates the fluid temperature, water density fraction and water phase velocity based upon the cladding surface temperature. The fuel temperature and the fluid density are then passed back to Rattlesnake for another neutronics calculation. For this paper, sets of results from the detailed calculation are provided for both depletion and the SBO event. We demonstrate that a detailed simulation akin to first principles can be achieved using MAMMOTH and the various MOOSE application on differing domains. (authors)
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗