Engineering Papers⌕ Search

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.

144 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Latest developments in the MOOSE fluid properties module

The fluid properties module in MOOSE serves a variety of fluid simulation applications based on MOOSE, including the MOOSE Navier Stokes module~\cite{moose_ns}, Pronghorn~\cite{pgh}, SAM~\cite{sam}, the MOOSE thermal hydraulics module, RELAP-7~\cite{relap7} and subchannel~\cite{subchannel}. It is used for coarse mesh multi-dimensional thermal-hydraulics~\cite{pgh}, 1D systems analysis~\cite{sam,relap7} in nuclear reactor analysis, and porous flow simulations~\cite{porous} for underground gas storage and water seepage. The use of consistent fluid properties across fluid flow applications facilitates coupled simulations~\cite{anl_sam_pgh}. The module offers a consistent set of interfaces to implement to create a new fluid property. There are numerous fluid properties of interest in the entirety of all fields of fluid flow simulations, and this is exacerbated by the use of different variable sets depending on the compressibility of the fluid. For single-phase fluids, the following variable sets may be used to compute fluid properties: (pressure, temperature) and (specific volume, specific internal energy). Some properties may also be computed using the (pressure, density) or the (specific volume, specific enthalpy) variable sets. In order to reduce the challenge of adding a new fluid property, properties may be implemented partially.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Tabulated Fluid Properties Research Report

The Multiphysics Object-Oriented Simulation Environment (MOOSE) enables a wide range of advanced nuclear reactor simulations.[6] Under the guidance of MOOSE’s Thermal Hydraulics Team,I worked to expand the capabilities of Tabulated Fluid Properties (TFP) in the fluid properties module. The fluid properties module allows the user to determine a variety of fluid properties by interpolating points between tabulated data. I implemented the ability to use bilinear interpolation instead of bicubic interpolation for interpolating tabulated data. I also changed the method of variable set inversions to use a 2-dimensional Newton’s Method utility that I created. Variable set inversions are often done from (v,e) to (p,T), where v is specific volume, e is specific internal energy, p is pressure and T is temperature. New routines have also been added into TFP such that it can be used with more applications, such as the Navier Stokes and Thermal Hydraulics modules in MOOSE for Pronghorn[5] and RELAP-7[1] respectively. This work was spurred by interest from NASA in testing a Nuclear Thermal Propulsion (NTP) engine system. NTP engines have drastically different fluid properties throughout the engine and Tabulated Fluid Properties provides the flexibility needed to properly simulate and test these engines.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Plant-specific Model and Data Analysis using Dynamic Security Modeling and Simulation

The requirements for U.S. nuclear power plants to maintain a large on-site physical security force contribute to their high operational costs. The cost of maintaining the current physical security posture is approximately 10% of the overall operation and maintenance budget for commercial nuclear power plants. The goal of the Light Water Reactor Sustainability (LWRS) program’s physical security pathway is to develop tools, methods, and technologies and provide the technical basis for an optimized physical security posture. The conservatisms built into current security postures may be analyzed and minimized in order to reduce security costs while still ensuring adequate security and operational safety. The research performed at Idaho National Laboratory within LWRS program’s physical security pathway has successfully developed a dynamic force-on-force modeling framework using various computer simulation tools and integrating them with the dynamic assessment Event Modeling Risk Assessment using Linked Diagrams (EMRALD) tool. This document provides an update on the progress in applying a dynamic computational framework that links results from a commercially available force-on-force simulation tool, a commercially available thermal-hydraulic tool, and EMRALD to an operating commercial nuclear power plant. This report is only a summary of the progress and does not contain specific modeling results as those contain sensitive security information. This process of including plant procedures and multiple analysis results is being called Modeling and Analysis for Safety Security using Dynamic EMRALD Framework or MASS-DEF. Previous reports described how a user could integrate their plant-specific force-on-force models with the dynamic simulation tool EMRALD, model operator actions, integrate with probabilistic risk assessment tools, such as CAFTA (Computer Aided Fault Tree Analysis System) or SAPHIRE (Systems Analysis Programs for Hands-on Integrated Reliability Evaluations), and with thermal-hydraulic tools, such as RELAP-5. Previous reports applied various combinations of available simulations codes with EMRALD using generic plant models to demonstrate how to perform the analysis. This report documents the results of applying the dynamic computational framework to an actual nuclear facility using their security scenarios and timelines. This report does not contain any plant's sensitive information and/or Safeguards Information. The purpose of this study was to verify that results achieved using generic models are similar to actual plant results and to refine our guidance on the use of the framework. This assessment enables further analysis, such as what-if scenarios and staff-reduction evaluation, thereby optimizing physical security at plants.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Development of Segregated Thermal-Hydraulics Solvers in MOOSE

The simulation of fluid flows is an essential part of the design and analysis of nuclear systems. Algorithms able to simulate flows at different fidelity levels are available in the Multiphysics Object-Oriented Simulation Environment (MOOSE) and MOOSE-based applications such as Pronghorn \cite{novak2018pronghorn}, Pronghorn-Subchannel, RELAP-7, and SAM. Currently, significant effort is being invested in the development of coarse-mesh Computational Fluid Dynamics (CFD) capabilities within MOOSE and Pronghorn for the simulation of Generation IV nuclear reactors. Traditionally, the solution algorithms in MOOSE have relied on Newton or quasi-Newton methods (such as the preconditioned Jacobian-free Newton-Krylov method) where residuals and Jacobians (or approximations thereof) are constructed. Both Newton and quasi-Newton methods require the solution of a linear system at each nonlinear Newton iteration with the Jacobian as the system matrix. The Jacobian contains blocks originating from all variables in the problem (i.e., for thermal-hydraulics at least pressure, velocities, and temperature). Due to the formulation of the problem in a general multiphysics setting on unstructured mesh, creating a good preconditioner for the linear system can be challenging, thus many fluid applications have utilized direct solver-based methods such as LU factorization. However, with increasing system size and complexity in multi-dimensional problems, the direct solution of linear systems becomes computationally expensive both in execution time and and memory. For this reason, recent effort has focused on adapting segregated solution algorithms for CFD problems in MOOSE. These algorithms use fixed-point iteration between segregated systems whose assembly and preconditioning are easier those of the monolithic system. Initial results show that the segregated solution algorithm outperforms the monolithic approach in terms of memory usage and for large 3D problems in terms of CPU time as well.

42 ENGINEERING↗
DOE OSTI2024

An Evaluation of The Dynamic Physical Security Risk Assessment Methodology for Fleet-Wide Applications

The requirements for U.S. nuclear power plants to maintain a large onsite physical security force contribute to their high operational costs. The cost of maintaining the current physical security posture is approximately 10% of the overall operation and maintenance budget for commercial nuclear power plants. The goal of the Light Water Reactor Sustainability (LWRS) program’s physical security pathway is to develop tools, methods, and technologies and provide the technical basis for an optimized physical security posture. The conservatisms built into current security postures may be analyzed and minimized to reduce security costs while still ensuring adequate security and operational safety. The research performed at Idaho National Laboratory within LWRS program’s physical security pathway has successfully developed a dynamic force-on-force modeling framework using various computer simulation tools and integrating them with the dynamic assessment Event Modeling Risk Assessment using Linked Diagrams (EMRALD) tool. This integrated process for physical security analysis is named Modeling and Analysis for Safety Security using Dynamic EMRALD Framework (MASS-DEF). This document provides an update on the progress in applying the MASS-DEF process to an operating commercial nuclear power plant as well as additional industry feedback regarding use of the tool for other physical security risk-informed topics. This report is only a summary of the progress and does not contain specific modeling results as those contain sensitive security information. Previous reports described how a user could integrate their plant-specific force-on-force models with the dynamic simulation tool EMRALD, model operator actions, and integrate with probabilistic risk assessment tools, such as CAFTA (Computer Aided Fault Tree Analysis System) or SAPHIRE (Systems Analysis Programs for Hands-on Integrated Reliability Evaluations), and with thermal-hydraulic tools, such as RELAP-5 or MAAP. Previous reports applied various combinations of available simulations codes with EMRALD using generic plant models to demonstrate how to perform the analysis. This report is an update the progress of applying the dynamic computational framework to an actual nuclear facility using their security scenarios and timelines. This report also provides an update to the procedural guidance for the MASS-DEF process and an overview of the generic models available for use by utilities. This report does not contain any plant’s sensitive information and/or safeguards information. This study’s purpose was to verify that the results achieved using generic models are similar to actual plant results and refine our guidance on the use of the framework. This assessment enables further analysis, such as what-if scenarios and staff-reduction evaluation, thereby optimizing physical security at plants.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Development of Segregated Thermal-Hydraulics Solvers in MOOSE

The simulation of fluid flows is an essential part of the design and analysis of nuclear systems. Algorithms able to simulate flows at different fidelity levels are available in the Multiphysics Object-Oriented Simulation Environment (MOOSE) and MOOSE-based applications such as Pronghorn \cite{novak2018pronghorn}, Pronghorn-Subchannel, RELAP-7, and SAM. Currently, significant effort is being invested in the development of coarse-mesh Computational Fluid Dynamics (CFD) capabilities within MOOSE and Pronghorn for the simulation of Generation IV nuclear reactors. Traditionally, the solution algorithms in MOOSE have relied on Newton or quasi-Newton methods (such as the preconditioned Jacobian-free Newton-Krylov method) where residuals and Jacobians (or approximations thereof) are constructed. Both Newton and quasi-Newton methods require the solution of a linear system at each nonlinear Newton iteration with the Jacobian as the system matrix. The Jacobian contains blocks originating from all variables in the problem (i.e., for thermal-hydraulics at least pressure, velocities, and temperature). Due to the formulation of the problem in a general multiphysics setting on unstructured mesh, creating a good preconditioner for the linear system can be challenging, thus many fluid applications have utilized direct solver-based methods such as LU factorization. However, with increasing system size and complexity in multi-dimensional problems, the direct solution of linear systems becomes computationally expensive both in execution time and and memory. For this reason, recent effort has focused on adapting segregated solution algorithms for CFD problems in MOOSE. These algorithms use fixed-point iteration between segregated systems whose assembly and preconditioning are easier those of the monolithic system. Initial results show that the segregated solution algorithm outperforms the monolithic approach in terms of memory usage and for large 3D problems in terms of CPU time as well.

42 ENGINEERING↗
DOE OSTI2018

Narrowing transient testing pulse widths to enhance LWR RIA experiment design in the TREAT facility

The Transient Reactor Test (TREAT) Facility provided thousands of transient irradiations and plays a crucial role in nuclear-heated safety research. The facility’s flexible design and multi-mission nature saw historic experiments for numerous reactor fuels and transient types but was never specifically adapted to address very-brief pulse transients akin to postulated Light Water Reactor (LWR) Reactivity Initiated Accidents (RIA). Since the behaviors of fuel under these conditions depends strongly on energy input duration and resulting cladding time/temperature response under pellet cladding mechanical interactions, this three-year project was conceived to investigate new pulse-narrowing capabilities. Kinetic models showed incremental improvements for minor facility enhancements including increased reactivity step insertions (to initiate the power pulse) and slightly-increased transient rod drive speed for pulse termination (“clipping”). Replacing peak fuel assemblies, repositioning non-transient control rods to hold down the “hot side” of the core, and balancing against required excess reactivity needs, the limiting fuel assembly power can be reduced by ~20 %. Furthermore, this is a remarkable discovery of latent capability in TREAT, not only in enabling the subject capabilities for reduced pulse widths, but also significantly “uprating” the core’s transient energy capacity. Incremental improvements can likely enhance TREAT’s capability into BWR-relevant missions; briefer PWR pulse shapes can only be achieved with an advanced clipping system. Helium-3 ( 3 He) was found to offer the greatest benefits for clipping design and overall performance. A unique concept showed great promise for enabling a 3 He-based system with a reasonable cost; this design concept will likely become the focal point of future work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2014

Theory and Implementation of Nuclear Safety System Codes - Part II: System Code Closure Relations, Validation, and Limitations

This is Part II of two articles describing the details of thermal-hydraulic sys- tem codes. In this second part of the article series, the system code closure relationships (used to model thermal and mechanical non-equilibrium and the coupling of the phases) for the governing equations are discussed and evaluated. These include several thermal and hydraulic models, such as heat transfer coefficients for various flow regimes, two phase pressure correlations, two phase friction correlations, drag coefficients and interfacial models be- tween the fields. These models are often developed from experimental data. The experiment conditions should be understood to evaluate the efficacy of the closure models. Code verification and validation, including Separate Effects Tests (SETs) and Integral effects tests (IETs) is also assessed. It can be shown from the assessments that the test cases cover a significant section of the system code capabilities, but some of the more advanced reactor designs will push the limits of validation for the codes. Lastly, the limitations of the codes are discussed by considering next generation power plants, such as Small Modular Reactors (SMRs), analyz- ing not only existing nuclear power plants, but also next generation nuclear power plants. The nuclear industry is developing new, innovative reactor designs, such as Small Modular Reactors (SMRs), High-Temperature Gas-cooled Reactors (HTGRs) and others. Sub-types of these reactor designs utilize pebbles, prismatic graphite moderators, helical steam generators, in- novative fuel types, and many other design features that may not be fully analyzed by current system codes. This second part completes the series on the comparison and evaluation of the selected reactor system codes by discussing the closure relations, val- idation and limitations. These two articles indicate areas where the models can be improved to adequately address issues with new reactor design and development.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2023

Automated power-following control for nuclear thermal propulsion startup and shutdown using MOOSE-based applications

This paper describes an investigation of various automated control strategies applied to a full-core multiphysics Griffin/Bison/RELAP-7 model of a prototypical nuclear thermal propulsion system. In all cases, control is achieved by actuating control drums based on the demanded power and predicted quantities from the numerical model. One key finding is that hybrid proportional integral derivative controllers – a novel type of controller that uses both power and reactivity predicted signals – can demonstrate a level of performance rivaling that of period-generated control. The former requires parameter tuning, while the latter mostly necessitates providing reactivity coefficients and temperature rates of change, which could be very challenging to accurately measure in real-time. In addition, decay heat plays an important role in determining cooling requirements during cooldown phases. A decay heat model, accounting for burn time and throttling, was derived, and is incorporated within the model to simulate the steady-state and shutdown phases and satisfyingly follow the power demand. However, temperature overshoots and non-constant specific impulse during throttling will necessitate further improvements.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2017

Preliminary Two-Phase Terry Turbine Nozzle Models for RCIC Off-Design Operation Conditions

This report presents the effort to extend the single-phase analytical Terry turbine model to cover two-phase off-design conditions. The work includes: (1) adding well-established two-phase choking models – the Isentropic Homogenous Equilibrium Model (IHEM) and Moody’s model, and (2) theoretical development and implementation of a two-phase nozzle expansion model. The two choking models provide bounding cases for the two-phase choking mass flow rate. The new two-phase Terry turbine model uses the choking models to calculate the mass flow rate, the critical pressure at the nozzle throat, and steam quality. In the divergent stage, we only consider the vapor phase with a similar model for the single-phase case by assuming that the liquid phase would slip along the wall with a much slower speed and will not contribute the impulse on the rotor. We also modify the stagnation conditions according to two-phase choking conditions at the throat and the cross-section areas for steam flow at the nozzle throat and at the nozzle exit. The new two-phase Terry turbine model was benchmarked with the same steam nozzle test as for the single-phase model. Better agreement with the experimental data is observed than from the single-phase model. We also repeated the Terry turbine nozzle benchmark work against the Sandia CFD simulation results with the two-phase model for the pure steam inlet nozzle case. The RCIC start-up tests were simulated and compared with the single-phase model. Similar results are obtained. Finally, we designed a new RCIC system test case to simulate the self-regulated Terry turbine behavior observed in Fukushima accidents. In this test, a period inlet condition for the steam quality varying from 1 to 0 is applied. For the high quality inlet period, the RCIC system behaves just like the normal operation condition with a high pump injection flow rate and a nominal steam release rate through the turbine, with the net addition of water to the primary system; for the low quality inlet period, the RCIC turbine shaft work dramatically decreases and results in a much reduced pump injection flow rate, and the mixture flow rate through the turbine increases due to the high liquid phase flow rate. The net effect for this period is net removal of coolant from the primary loop. With the periodic addition and removal of coolant to the primary loop, the self-regulation mode of the RCIC system can be maintained for a quite long time. Both the IHEM and Moody’s models generate similar phenomena; however noticeable differences can be observed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

RAVEN User Guide

The goal of this document is to provide a set of detailed examples that can help the user to become familiar with the RAVEN code. RAVEN is capable of investigating system response and explore input space using various sampling schemes such as Monte Carlo, grid, or Latin Hypercube. How- ever, RAVEN strength lies in its system feature discovery capabilities such as: constructing limit surfaces, separating regions of the input space leading to system failure, and using dynamic super- vised learning techniques. New users should consult the RAVEN Tutorial to get started.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2018

Validation of RELAP5-3D Using HE-FUS3 Data

The RELAP5-3D computer code was validated using data from the HE-FUS3 experimental facility. HE-FUS3 is a helium-cooled, electrically-heated experimental facility that was used to generate experimental data for a benchmark exercise aimed at the validation of system transient analysis codes for very-high-temperature gas reactor applications. RELAP5-3D was validated using results from seven steady-state tests and one transient that simulated a loss-of-flow accident.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Automated Control for Nuclear Thermal Propulsion Start-Up using MOOSE-based Applications

This report presents a Griffin/Bison/RELAP-7 numerical model of a prototypical NTP system that features fuel assemblies arranged in rings, and which was designed to simulate rapid startup transients. The physics modeled include full-core neutronics, assembly-wise heat conduction, and conjugate heat transfer, with the balance of plant mainly imposed through boundary conditions. In addition, various forms of automated reactivity control were deployed by using the MOOSE to autonomously drive the model and simulate the reactor transitioning from assumed initial conditions to nominal power in a fraction of a minute. To generate the cross-sections of the neutronics model, and in an effort to simultaneously account for the tremendous axial temperature gradients in the reactor and to limit the number of state points required for cross-section generation, the average component temperatures and hydrogen densities in the cooling channels were correlated to the average fuel and moderator temperatures, and fixed axial profiles were derived for nominal conditions and then used during the transient. With this approximation, a tractable cross-section library tabulated with fuel/moderator temperatures and CD angles was generated using Serpent. The full-core SPH correction procedure and the CD decusping technology in Griffin, respectively, ensure preservation of the multiplication factor and reaction rates at state points, along with a reasonably accurate reactivity worth between tabulated CD angles, despite using a coarse mesh. Feedback from other physics was calculated by modeling one representative fuel assembly per ring, along with the corresponding fuel and moderator cooling channels. To limit power overshoots during startup, another layer of multiphysics coupling was added to the model in order to automatically control the drums. Two different technologies presented herein showed outstanding performance in this regard: (1) a novel hybrid PID controller based on both power and reactivity signals, and (2) a PGC that relies on kinetics parameters and reactivity coefficients to predict future behavior and adjust the desired signal accordingly. A challenging benchmark was devised, featuring a power demand curve that exponentially increases by a factor of 500 within 30 seconds, then levels out after that. Both control approaches create a simulated power curve that closely follows the power demand curve and limits power overshoots to 1% or less. While the former approach requires more tuning of the internal parameters, the latter requires additional knowledge of the reactivity feedback coefficients and rates of change of the corresponding variables, including fuel and moderator temperature, which could be difficult to dynamically measure for a real NTP system. Fortunately, some inaccuracy in these quantities will not drastically degrade the PGC performance. Subsequently, a more realistic startup sequence was considered, in which the mass flow rate and outlet pressures are ramped up to model bootstrap and thrust build-up phases prior to reaching steady-state conditions, demonstrating the ability of the hybrid PID and PGCs to handle such transients, with both types of controllers exhibiting very similar behavior. Nevertheless, a significant chamber temperature overshoot was observed, caused by the demanded power signal and assumed mass flow rate. This issue could be mitigated by deploying a reactor controller that follows the chamber temperature signal and actuates both the control valves and drums (rather than using a power signal based solely on the drums to control reactivity). Enhancement of the hydrogen fluid properties available in MOOSE, as well as a better understanding of prototypical initial conditions, are also needed to further enhance this startup model. Finally, a study was performed to model decay heat post-shutdown, and to prepare for extending this model to predict shutdown behavior and post-shutdown pulsed cooling requirements.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2016

A Demonstration of Advanced Safety Analysis Tools and Methods Applied to Large Break LOCA and Fuel Analysis for PWRs

The U.S. Nuclear Regulatory Commission (NRC) is currently proposing a rulemaking designated as 10 CFR 50.46c to revise the loss-of-coolant accident (LOCA)/emergency core cooling system acceptance criteria to include the effects of higher burnup on fuel/cladding performance. We propose a demonstration problem of a representative four-loop PWR plant to study the impact of this new rule in the US nuclear fleet. Within the scope of evaluation for the 10 CFR 50.46c rule, aspects of safety, operations, and economics are considered in the industry application demonstration presented in this paper. An advanced safety analysis approach is used, by integrating the probabilistic element with deterministic methods for LOCA analysis, a novel approach to solving these types of multi-physics, multi-scale problems.

LOCA RIMM RISMC RELAP-7↗
DOE OSTI2016

An Equation of State for Liquid Pb83Li17

Thermal hydraulic and accident analysis codes such as RELAP5-3D and MELCOR for Fusion rely on an equation of state to specify all the thermodynamic properties of fusion-relevant working fluids such as PbLi. The existing liquid metal fluid properties in both RELAP5-3D and MELCOR are based on a five parameter "soft sphere" equation of state for which parameter sets that approximately reproduce experiment data were available in the literature. The PbLi parameters were based on a linear, mass-weighted average of the parameters for pure Pb and pure Li, which does not always result in adequate agreement between the equation and experimental data for PbLi, much of which was not available at the time the original equations were implemented. In order to address this shortcoming, we describe here a modification and non-linear least squares fitting of the parameters to closely match available experiment data for the density, sound speed, and specific heat capacity of liquid PbLi over a range of temperatures. The new equation has been implemented in both RELAP5-3D and MELCOR, and we make comparisons with both codes to experiment data to verify the implementation. Transport properties such as the thermal conductivity and viscosity, though not derived from the equation of state, have also been updated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Modeling and Simulation to Support Current and Advanced Reactors : Thermal Hydraulics

Presentation on thermohydraulic modeling and simulation capabilities in support of current and advanced reactors. The presentation includes information on the capabilities of RELAP5-3D, MOOSE and the experimental testbed for this purpose. This presentation is intended for a Taiwan delegation that will be visiting INL in August 2024.

42 ENGINEERING↗