Search NASA⌕ Search

SEARCH · Search NASA

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.

At least 55 records · Page 3

A SUMMARY OF THE THEORETICAL BASIS OF THE RELAP-7 NUCLEAR REACTOR SYSTEM SAFETY ANALYSIS CODE

The RELAP-7 code is the next generation nuclear reactor system safety analysis code being developed at the Idaho National Laboratory (INL)[1]. The code is based on INL’s scientific software development framework, MOOSE (Multi-physics Object Oriented Simulation Environment)[2]. The primary basis of the RELAP-7 governing theory includes 7-equation two-phase flow, reactor core heat transfer, and reactor kinetics models. This paper will focus upon the compressible, fully nonequilibrium two-phase flow model.

7-equation two-phase model↗

Development and Implementation of Mechanistic Terry Turbine Models in RELAP-7 to Simulate RCIC Normal Operation Conditions

As part of the efforts to understand the unexpected “self-regulating” mode of the RCIC (Reactor Core Isolation Cooling) systems in Fukushima accidents and extend BWR RCIC and PWR AFW (Auxiliary Feed Water) operational range and flexibility, mechanistic models for the Terry turbine, based on Sandia’s original work [1], have been developed and implemented in the RELAP-7 code to simulate the RCIC system. In 2016, our effort has been focused on normal working conditions of the RCIC system. More complex off-design conditions will be pursued in later years when more data are available. In the Sandia model, the turbine stator inlet velocity is provided according to a reduced-order model which was obtained from a large number of CFD (computational fluid dynamics) simulations. In this work, we propose an alternative method, using an under-expanded jet model to obtain the velocity and thermodynamic conditions for the turbine stator inlet. The models include both an adiabatic expansion process inside the nozzle and a free expansion process outside of the nozzle to ambient pressure. The combined models are able to predict the steam mass flow rate and supersonic velocity to the Terry turbine bucket entrance, which are the necessary input information for the Terry turbine rotor model. The analytical models for the nozzle were validated with experimental data and benchmarked with CFD simulations. The analytical models generally agree well with the experimental data and CFD simulations. The analytical models are suitable for implementation into a reactor system analysis code or severe accident code as part of mechanistic and dynamical models to understand the RCIC behaviors. The newly developed nozzle models and modified turbine rotor model according to the Sandia’s original work have been implemented into RELAP-7, along with the original Sandia Terry turbine model. A new pump model has also been developed and implemented to couple with the Terry turbine model. An input model was developed to test the Terry turbine RCIC system, which generates reasonable results. Both the INL RCIC model and the Sandia RCIC model produce results matching major rated parameters such as the rotational speed, pump torque, and the turbine shaft work for the normal operation condition. The Sandia model is more sensitive to the turbine outlet pressure than the INL model. The next step will be further refining the Terry turbine models by including two-phase flow cases so that off-design conditions can be simulated. The pump model could also be enhanced with the use of the homologous curves.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗

Solving implicit multi-mesh flow and conjugate heat transfer problems with RELAP-7

The fully implicit simulation capability of RELAP-7 to solve multi-mesh flow and conjugate heat transfer problems for reactor system safety analysis is presented. Compared to general single-mesh simulations, the reactor system safety analysis-type of code has unique challenges due to its highly simplified, interconnected, one-dimensional, and zero-dimensional flow network describing multiple physics with significantly different time and length scales. To use the Jacobian-free Newton Krylov-type of solver, preconditioning is generally required for the Krylov method. The uniqueness of the reactor safety analysis-type of code in treating the interconnected flow network and conjugate heat transfer also introduces challenges in providing preconditioning matrix. Typical flow and conjugate heat transfer problems involved in reactor safety analysis using RELAP-7, as well as the special treatment on the preconditioning matrix are presented in detail. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Two-zone stratified wetwell model development and implementation for RELAP-7

The Reactor Core Isolation Cooling (RCIC) system consists of a turbine, pump and wetwell and is used in many US Boiling Water Reactors (BWRs) as an important heat removal system. The wetwell is a major heat sink within containment. Experimental investigations have shown that thermal stratification may occur in the wetwell during RCIC operation. Current systems-level analysis codes lack stratified wetwell models. Here, this paper describes the development of a stratified wetwell model for the RELAP-7 code. This paper uses a simple two-zone model for the wetwell water that can capture thermal stratification. This model assumes that a buoyant plume transports heat and mass from the steam injection site to the upper water zone. The buoyant plume assumption requires that steam is injected at a low enough rate so it adds negligible momentum. Model validation against experimental data shows similar trends of temperature field development in the suppression pool.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

RELAP-7 Theory Manual

This document summarizes the physical models and mathematical formulations used in the RELAP-7 code.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

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↗

Developing Fully Coupled Dynamical Reactor Core Isolation System Models in RELAP-7 for Extended Station Black-Out Analysis

The reactor core isolation cooling (RCIC) system in a boiling water reactor (BWR) provides makeup water to the reactor vessel for core cooling when the main steam lines are isolated and the normal supply of water to the reactor vessel is lost. It was one of the very few safety systems still available during the Fukushima Daiichi accidents after the tsunamis hit the plants and the system successfully delayed the core meltdown for a few days for unit 2 & 3. Therefore, detailed models for RCIC system components are indispensable to understand extended station black-out accidents (SBO) for BWRs. As part of the effort to develop the new generation reactor system safety analysis code RELAP-7, major components to simulate the RCIC system have been developed. This paper describes the models for those components such as turbine, pump, and wet well. Selected individual component test simulations and a simplified SBO simulation up to but before core damage is presented. The successful implementation of the simplified RCIC and wet well models paves the way to further improve the models for safety analysis by including more detailed physical processes in the near future.

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↗