Search NASA⌕ Search

SEARCH · Search NASA

Results for “REACTOR ACCIDENT SIMULATION”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

COMPUTATION FLUID DYNAMICS ANALYSIS FOR GENERIC SMALL MODULAR REACTOR CONTAINMENT SEPARATE EFFECTS TEST

It is desirable for fourth-generation Small Modular Reactors to be passively cooled in standard and accident operations. Passive Containment Cooling Systems can reject heat from the containment structure, without using pumps or blowers. The targeted design containment structure is a large, domed, stainless steel, cylindrical vessel. In a postulated Design Basis Accident, steam will flash inside containment. Steam condensation occurs on the inner containment wall and transfers heat through the steel containment into a large body of water known as the annular reservoir (AR) surrounding the vessel serving as the ultimate heat sink. Natural circulation drives the flow in the AR and heat will be released to the environment by evaporation of water. Unique containment geometry requires a separate effects test (SET) facility for the verification and validation of the computer code and evaluation model development and assessment for reactor licensing efforts. In this study, STAR-CCM+, a computational fluid dynamics (CFD) code was used to inform the decision-making process on the design of the SET. The CFD simulation modeled, a two-phase turbulent flow with fluid film development and heat transfer for different containment geometries. The Reactor Excursion and Leak Analysis Program will also be used in a code-to-code verification against the CFD results.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2023

BISON analysis of FeCrAl and Zircaloy cladding deformation during simulated BWR cyclic dryout conditions

The BISON fuel performance analysis code has been used to study FeCrAl and Zircaloy cladding performance to better understand cladding deformation behavior in boiling water reactors during conditions caused by anticipated operational occurrences (AOOs) or by anticipated transients without scram (ATWS). BISON models of novel separate effects tests simulating ATWS events predicted mechanical responses to thermal cycling boundary conditions. Here, we assessed FeCrAl and Zircaloy deformation and burst models to determine those best suited for cyclic dryout event analysis. In addition, we compared diametral strain profiles of the post-experiment claddings to simulation results. The capability of the BISON code for modeling BWR AOO or ATWS events is promising.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Support for Reactor Operators in Case of Cyber-Security Threats (NEUP Final Report)

Increasing use of digital technology in nuclear power plants (NPPs) poses cyber-security as a crucial threat to public safety and to continuous energy production. Cyber-security risks are comprised of complex known and unknown interactions between various entities, system vulnerabilities, network protocols, human users and malicious attacks. There is little understanding or research geared towards plant operators’ response under cyber-security threats and operational procedures to cope with such threats. This is particularly critical when a cyber-security event masquerades as a safety incident or an evolving accident. Instead of leading operators to the remediation of the accident, the masked cyber event may lead them to circumvent the safety mechanisms of the plant. Given these challenges, the objective of this project is to develop a tool to characterize abnormal NPP events as “cyber” or “safety” incidents, and to further develop this tool as a real-time operator aid to assist in response to the incident. This NEUP Final Report summarizes the research efforts towards the above goals and the results in the following three aspects. First, a methodology for classifying abnormal NPP events as cyber-attacks or safety incidents was developed. This provides the basis for responding to potential cyber-attacks. In this methodology, the causal relationships between various variables (i.e., physical variables, network variables, and control variables) for a nuclear system under consideration are modeled using dynamic Bayesian networks. Based on observations of certain variables (i.e., the evidence), inference is performed to identify the state of hidden variables, which helps indicate the occurrence of cyber-attacks or safety incidents. The proposed methodology was applied to the full scope NPP simulator system at The Ohio State University as a case study. The result shows that the proposed methodology is capable of distinguishing between cyber-attacks and safety incidents. Second, a game-theoretic approach for responding to cyber-attacks on nuclear power plants was developed. This is the second major goal of the project, and the method is used to recommend optimal responses to reactor operators in the case of cyber-attacks. We consider two cases in this research. In the first case, the interaction between the defender (i.e., the reactor operator) and the attacker is modeled as a two-player, nonzero-sum, Markov game with an infinite time horizon. In the second case, we consider a semi-Markov game between the defender and the attacker within a finite time horizon. The model in the second case enables us to consider the time sensitive response of the defender. In each case, the method for deriving the Nash equilibrium is developed. The Nash Equilibrium of the game provides the valid prediction of both players’ actions because no single player can benefit from unilaterally deviating from the equilibrium policy if the other player adheres to his/hers, hence it provides the best response of the defender to cyber-attacks. The proposed method in each case is demonstrated through a simplified benchmark digital feedwater control system. Third, a prototype display is developed to integrate the above information. The prototype display depicts the communication information to the operator and incorporates the cyber-security tool developed by the team into the interface. It consists of the regular display in current control rooms, as well as a panel used to indicate the occurrence of cyber-attacks or safety events and a panel used to guide the operator’s response to detected cyber-attacks.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Flow reversal benchmark of a one-sided heated narrow rectangular channel with CATHARE and RELAP5

Flow reversal in narrow coolant channels can be a crucial phenomenon for the safety of research reactors with a downward nominal flow direction. During a loss of forced flow accident, the downward flow stagnates briefly before transitioning into an upward natural circulation flow. The fuel may be damaged if dryout occurs and threshold fuel and/or cladding temperatures are exceeded. A comprehensive study is provided for flow reversal in narrow rectangular channels by examining experimental data and conducting software model analyses. The literature on flow reversal was reviewed, and selected experimental datasets were used to benchmark against CATHARE and RELAP5 models and also compare the code calculations with each other. The experimental data comes from flow reversal tests conducted with a narrow rectangular channel with one-sided heating. The results were compared with experimental data for successful flow reversal tests and predicted dryout power for dryout conditions. Also, the study examined the effects of the pump coastdown period, inlet liquid temperature, system pressure, and localized pressure drops. The experimental results showed that shorter coastdown periods, reduced pressure drops, and lower coolant inlet temperatures increased the dryout power. However, the system pressure did not noticeably affect the results. The simulation results showed that both CATHARE and RELAP5 agreed with experimental data, capturing the trends of the experimental results. Slight differences between each code calculation, as well as the predicted and measured dryout powers, were attributed to experimental uncertainties and the modeling of physical phenomena such as wall nucleation, interfacial heat transfer, drag coefficients, and critical heat flux. Overall, this study provides an understanding of flow reversal and the prediction capabilities of thermal-hydraulics software models. In conclusion, a future study of the flow reversal benchmark of a narrow rectangular channel with two-sided heating may provide additional valuable insights.

CATHARE↗
DOE OSTI2022

Thermal-hydraulic and Fuel Performance Scoping Studies of a Flowing Water Capsule in TREAT

The restart of the Transient Reactor Test (TREAT) facility has provided a much needed capability for integral transient testing of nuclear fuel. This testing is necessary to qualify new fuel concepts such as those developed under the Accident Tolerant Fuel (ATF) program, increasing the burnup limit of light water reactor (LWR) fuels, or a variety of other programs under the Advanced Fuels Campaign (AFC). The ATF campaign has been the main driver behind the development and implementation of a variety of capsules for TREAT experiments. The Separate-Effect Test Holder (SETH) was a dry capsule that enabled testing of ATF concepts under reactivity-initiated accident (RIA) heating conditions. Following SETH, the Static Environment Rodlet Transient Test Apparatus (SERTTA) was developed to enable RIA testing in a static water environment. In efforts to support the need for future Loss-of-Coolant Accident (LOCA) tests, the Transient Water Irradiation System for TREAT (TWIST) capsule has been developed that allows for water to drain from around the fuel rod and lower the pressure to simulate LOCA conditions. All these capsules that have been developed for LWR fuel testing all lack the capability for forced convective cooling which in some applications may limit their ability to test under prototypic conditions. A design and modeling effort has been started to modify the TWIST capsule by adding a flow tube and impeller that can create flowing coolant conditions for the fuel rod. RELAP5-3D and BISON models are being used to study the differences between RIA, LOCA, Anticipated Operational Occurrences (AOO), LWR power cycling, and other transient scenarios under flowing conditions capable in the flowing capsule and the current stagnant water capsules (SERTTA and TWIST). The scoping study will provide guidance on the needed capabilities for the flowing capsule and the limitations of the currently developed capsule for LWR testing in TREAT.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Development of Shielding Benchmarks Using the Godiva IV Assembly

The Nuclear Criticality Safety Program (NCSP) is developing a shielding benchmark using the Godiva IV assembly as a source. Its present status is reviewed herein. Even eight decades into the nuclear era, substantial work remains to develop a database of shielding benchmarks to support future nuclear development. A nuclear simulation is only as good as its supporting data, inputs, and validation basis. Uncertainty in these areas is addressed using conservatism, which adds margin and, occasionally, cost. In many shielding situations, high accuracy is not necessary because additional material is not particularly expensive. After all, 1–2 cm of lead often reduces the gamma ray dose substantially. However, in certain areas, conservatism can add unnecessary cost. These areas include mobile shielding applications such as casks, ships, microreactors, and spacecraft, where weight and, thus, margin is expensive. Although these characteristics are side benefits for NCSP shielding benchmark development, the main driver is enabling more reliable placement of criticality accident alarm systems (CAASs) in nuclear material facilities, such as those dedicated to the production of advanced reactor fuels. CAAS placement relies on more than accurate data and code validation. It also relies on sufficiently accurate materials specifications, geometry specifications, and a well-defined, alarm-producing baseline accident. All these things require tacit knowledge and understanding of the problem being evaluated. Benchmarks can help ensure this understanding.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

TRISO Burnup-Dependent Failure Analysis of a HTGR Design-Basis Accident Using BISON

Here, this work assesses the failure behavior of the silicon carbide (SiC) layer in TRIstructural ISOtropic (TRISO) fuel particles with BISON during both steady-state and transient conditions for the MHTGR-350 design by General Atomics. A one-dimensional BISON model for a uranium oxycarbide–bearing TRISO is developed to simulate a power level of 50 mW/particles up to 12.4% fissions per initial metal atom (FIMA). Stress predictions for the materials of interest are presented as a function of burnup, along with the SiC failure probability computed using Weibull statistics under the assumption of realistic SiC quality.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Report on Year-4 of Water NSTF Matrix Testing: Facility Maintenance and Accident Testing

Under support from the Department of Energy (DOE) and the Office of Advanced Reactor Technologies (ART), a large-scale test facility has been constructed at Argonne National Laboratory to generate NQA-1 qualified validation data for passive decay heat removal systems in advanced reactors. The Natural convection Shutdown heat removal Test Facility (NSTF) reflects key features of a ½ scale, water-based, Reactor Cavity Cooling System (RCCS) and is intended to study the behavior, bound performance, and ultimately guide design decisions for passive decay heat removal systems for advanced reactors. In addition to the experimental activities detailed in this report, a supportive computational modeling effort is on-going which has been demonstrated to significantly strengthen the experimental program while also improving accuracy of the computer models. Together these create a mutually beneficial relationship integral to meeting the overall program objective of examining the heat removal performance of the RCCS concept. This report serves as a summary of maintenance and experimental activities during the program’s fourth year of water-based operation. A planned six-month maintenance period began in August 2021, during which major inspections, repairs, cleaning, and installation of new instrumentation and data acquisition hardware were conducted. Most significantly, two heaters that faulted during Year-3 were repaired, allowing the facility to resume use of the full heated section area and full range of available electric power. The remainder of the year consisted of eight months of test operations, during which the facility logged 211 hours of active heating across one bake-out (following the maintenance period) and seven matrix tests; five classifieds as Accepted per NQA-1, one as Trending, and one as Failed. Testing began by performing two repeat cases to confirm expected facility response and behavior during both single- and two-phase flow conditions, ensuring no changes were introduced during the maintenance period that might have altered the thermal-hydraulic characteristics of the facility. In continuation of the power parametric series initiated in previous years, a high-power test case was then performed examining heat removal performance at a decay heat load equivalent to 2.4 MWt, full-scale, a level exceeding maximum design targets. Additional testing then introduced various blockage scenarios along the network piping, examining the effects of partial and complete blockages of the flow paths on the system behavior and heat removal performance. A study of static boiling tests directed at understanding the geysering two-phase instability was also conducted. The loop was filled only to the bottom of the tank outlet, creating an open loop configuration that prevents any natural circulation flow from occurring, and the heaters were powered on until the facility reached saturation conditions. Following, a series of quasi-steady-state conditions were introduced by adjusting the inventory level in the adiabatic chimney piping at decreasingly lower elevations above the heated region. A strong correlation of geysering characteristics and loop level was observed, with flow and temperature excursions decreasing in intensity, but increasing in frequency, as the fill were reduced to lower elevations along the chimney piping. Once the level fell very low in the chimney, at points near the top of the heated section, the system reached a stable state of continuous boiling without any occurrence of geysering eruptions. A final significant testing accomplishment this year was successful completion of an “accident scenario” test, whose operating conditions were based on a prototypic decay heat curve provided by Framatome and scaled for the NSTF. This test began by establishing steady-state, single-phase “normal operation” conditions, before simulating an accident trip where the availability of active cooling systems was lost. Loop temperatures gradually increased until reaching saturation and subsequent two-phase boiling flow. Over the course of an extended operational period along the defined decay heat curve, steam boil-off caused gradual but continued depletion of liquid inventory until reaching a critically low level causing flow stagnation and cessation of natural circulation heat removal. At this point, after nearly 72 hours of continuous operation, a cold refill was performed to replenish the system inventory and allow the facility to re-establish closed loop natural circulation flow and return to a safe operational state.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Development of a BISON validation case for the TRISO transient irradiations in NSRR using effective heat capacity methods

The current tristructural isotropic (TRISO) fuel assessment and validation database in BISON primarily covers steady- state irradiation and high-temperature furnace testing. Transient assessment cases are potentially needed to support U.S. industry efforts in designing and deploying commercial reactors using TRISO fuels. Historical transient tests in- volving TRISO fuels used highly conservative conditions compared to the typical high-temperature gas-cooled reactor accident scenarios. Despite this, modeling historical transient tests is fundamental for evaluating BISON’s predictive capabilities, adapting material properties for high-temperature and high-particle-power regimes, and developing a sys- tematic validation approach for TRISO transient applications. This work developed a 1D model of transient experiments carried out at the Nuclear Safety Research Reactor using BISON. BISON’s predictions of energy deposition, UO 2 melting onset, and molten volume fractions are compared against experimental measurements. Melting was modeled using an effective specific heat capacity model for UO 2 . We found that BISON’s predictions are in reasonable agreement with experimental data for low-energy-deposition cases, and that BISON overpredicts melting at higher energy depositions. We also discuss the potential causes of discrepancies between the simulated and measured results and propose ways to further develop the model. Although these simulations used conservative conditions compared to those expected for actual TRISO-fueled reactors, they extended the range of conditions reflected in the data in the existing BISON database for TRISO fuels.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

BISON High Burnup Structure Modeling Capabilities Validated with a Selection of the Halden IFA-650 Rods

The U.S. Department of Energy (DOE)’s Nuclear Energy Advanced Modeling and Simulation (NEAMS)program aims to develop predictive capabilities using computational methods for the analysis and design of advanced reactor and fuel cycle systems. This program has been supporting the development of BISON, a high-fidelity and high-resolution fuel performance tool at the engineering scale. This document continues analysis and refinement of capabilities added to BISON early this calendar year in regards to the incorporation of capabilities applicable to extended burnups in response to industry interest. Details are provided on high burnup thermal conductivity models, a refitting of the high-burnup structure (HBS) porosity formation model to include additional data, the coupling of the HBS volume fraction model to thermal conductivity and fine fragmentation models, and validation activities. The IFA-650.4 and IFA-650.9 loss of coolant accident (LOCA) analyses are revisited with the latest developments in this report. A new validation case, IFA-650.14 has also been added to the BISON test suite.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Impact of LWR assembly structural features on cladding burst behavior under LOCA conditions

This provides an initial scoping study on clad balloon and burst behavior for burnup extension of reactor fuel. The associated issues with burnup extension are fuel fragmentation, relocation, and dispersal in the event of cladding failure. The general finding of this work is that the structural features, spacer grids and mixing vanes, locally suppress cladding deformation but have little impact on the overall clad performance during loss-of-coolant accidents. The work detailed in a previous report by Capps et al. focused on core optimization via neutronics, thermal hydraulics and thermomechanical analysis for burnups beyond 62 GWD/tU and enrichments above 5%. Uncertainty of fuel fragmentation relocation and dispersal in high-burnup rods during accident conditions was also investigated. The dispersal aspect of fuel fragmentation depends on cladding rupture. Thus, assessing uncertainties in the rupture behavior is helpful in estimating the dispersal of high-burnup fuel. This study builds on the previous work by assessing the impact of assembly structural features on cladding balloon and burst behavior in a full-length fuel rod. In this work, the BISON fuel performance code was used to generate 2D radial and height meshes containing structural features commonly used in nuclear fuel assemblies. First, meshes were generated with spacer grids. Results were then compared to the cladding burst temperature and balloon strain results from the previous work. A mesh sensitivity study was performed to ensure that mixing vanes and spacer grid effects were appropriately considered, resulting in a more refined mesh than the previous study. The balloon deformation and burst times of the cladding were compared to the original case. Consideration was also given to the effect of rod initial pressure. In conclusion, 3D quarter rod simulations were also performed and found good agreement with the 2D simulations in clad deformation and reasonable agreement in burst times.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Initial development of a generic fluoride salt-cooled reactor model

Fluoride high-temperature reactors (FHRs) are high-temperature, low-pressure reactor concepts that use tri-structural isotropic (TRISO) fuel and molten fluoride salt coolant. These reactors have the potential to provide both electrical power and high-temperature process heat. We used generic FHR parameters for a pebble-bed FHR to develop an initial model with fresh fuel for a generic FHR (gFHR) in MELCOR and SCALE (NEWT and KENO). In this paper, we present the development of our gFHR models, which will serve as the baseline for a sensitivity and uncertainty analysis to quantify the range of possible source terms for FHRs in severe accidents. We present MELCOR results for fuel and coolant temperatures through the core, a nodalization study for the steady-state thermal hydraulic model, and development of reactor physics models in SCALE. As this work progresses, these models will be used to calculate source terms for a loss-of-forced-flow accident and to conduct a sensitivity study on this accident to establish a range of possible source terms. SCALE will provide reactor physics parameters like isotopic inventory, decay heat generation, and temperature coefficients of reactivity. Using the uncertainty quantification tools within SCALE, we will generate distributions for those parameters and will use the uncertainty quantification code RAVEN or DAKOTA to sample those distributions in MELCOR to quantify the impact of reactor physics and thermal hydraulic uncertainties on FHR source terms. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Thermo-Mechanical Phase-Field Modeling of Fracture in High-Burnup UO2 Fuels Under Transient Conditions

This study presents a novel multiphysics phase-field fracture model to analyze high-burnup uranium dioxide (UO2) fuel behavior under transient reactor conditions. Fracture is treated as a stochastic phase transition, which inherently accounts for the random microstructural effects that lead to variations in the value of fracture strength. Moreover, the model takes into consideration the effects of temperature and burnup on thermal conductivity. Therefore, the model is able to predict crack initiation, propagation, and complex morphologies in response to thermal gradients and stress distributions. Several simulations were conducted to investigate the effects of operational and transient conditions on fracture behavior and the resulting cracking patterns. High-burnup fuels exhibit reduced thermal conductivity, elevating temperature gradients and resulting in extensive radial and circumferential cracks. Transient heating rates and temperatures significantly affect fracture patterns, with higher heating rates generating steeper gradients and more irregular crack trajectories. This approach provides critical insights into fuel integrity during accident scenarios and supports the safety evaluation of extended burnup limits.

Chemistry↗
DOE OSTI2021

Improvement of the $\mathrm{BISON U_3Si_2}$ modeling capabilities based on multiscale developments to modeling fission gas behavior

Uranium silicide (U 3 Si 2 ) is a concept explored as a potential alternative to UO 2 fuel used in light water reactors (LWRs) since it may improve accident tolerance and economics due to its higher thermal conductivity and increased uranium density. U 3 Si 2 has been previously used in research reactors in the form of dispersion fuel, but operated at lower temperatures than commercial LWRs. The research reactor data illustrated that significant gaseous swelling occurs as the fuel burnup increases. Therefore, it is imperative to understand the fission gas behavior of U 3 Si 2 under higher temperature LWR operating conditions. In this work, molecular dynamics and phase-field modeling techniques are used to reduce the uncertainty in select modeling assumptions made in developing the fission gas behavior model for U 3 Si 2 in the BISON fuel performance code. These lower length scale informed models are then utilized in the validation of BISON U 3 Si 2 modeling capabilities to simulate the ATF-1 experiments irradiated in the Advanced Test Reactor (ATR). Sensitivity analysis (SA) and uncertainty quantification (UQ) are included as part of the validation process to identify where further experiments and lower length scale modeling would be beneficial. Here, the multiscale modeling approach utilized in this work can be applied to new fuel concepts being explored for both LWRs and advanced reactors (e.g., uranium nitride, uranium carbide).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Ring Model Development and Validation for Prismatic HTGR Core Thermal-hydraulics and Safety Analysis

Because of the complex core geometry, prismatic high temperature gas-cooled reactors (prismatic HTGRs) often exhibit complex thermal fluid behaviors during both normal operating and transient conditions. Most HTGR designs rely on passive safety system for decay heat removal, such as the reactor cavity cooling system (RCCS). During postulated accidents like Pressurized Conduction Cooldown (PCC) event, the decay heat is first radially transferred from the core region to the reactor vessel outer surface, then to the RCCS cooling panels. The peak fuel temperature is controlled by heat transfer mechanisms with two distinctive characteristic length scales, i.e., the core-wise effective heat conduction and local heat conduction in the fuel pellet scale. As both length scales are essential to determine the fuel temperature, from the modeling perspective, computer codes must be able to capture heat transfer in both scales. This is challenging for both computational fluid dynamics (CFD) tools because of extremely large amount of computation resources required, and for system analysis codes because of the challenge to model the complex core geometry. Under the support of DOE-NE’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, efforts have been pursued to support HTGR technology development and its modeling and simulation needs. There is a particular need for advanced modeling and simulation tools to predict thermal-fluid behavior during safety-related transients. In our previous studies, the ring model was adopted in SAM and further developed to simulate the normal operating condition and a PCC event using the MHTGR-350 design of General Atomics as the reference design. This current work represents a continuation of these previous efforts, and the focus is to critically review and examine simplifications and assumptions made to develop the ring model, and to perform code validation using experimental data from an integral-effect test facility, the High Temperature Test Facility (HTTF) at the Oregon State University. In this study, the test PG-27 from the HTTF test suite was selected for code validation purpose. The test PG-27 is a transient test designed to simulate the PCC event of the MHTGR design. Very good agreements between SAM prediction and experimental measurements were found in both coolant and solid structure temperatures during the transient.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2019

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↗
DOE OSTI2025

Characterization of Fuel-to-Coolant Heat Transfer During Reactivity-Initiated Accidents Using Tightly Coupled Thermal Hydraulics and Fuel Thermomechanics

The reactivity-initiated accident (RIA) is a complex scenario with several tightly interacting physical phenomena. Accurately predicting fuel behavior during these transients is difficult due to limitations in the modeling of fuel-to-coolant heat transfer. Common approaches to simulate RIAs involve standalone calculations using either a fuel performance code or a thermal-hydraulic code. The complex interdependencies of thermal-hydraulic and fuel mechanical behavior suggest that a tight coupling between these codes may provide more accurate predictions of fuel-to-coolant heat transfer and cladding mechanical response. Here, RELAP5-3D and BISON are coupled in this paper to simulate RIAs, and a sensitivity analysis is performed to rank key thermal properties and two-phase heat transfer parameters relevant for fuel-to-coolant heat transfer and cladding failure mechanisms in UO 2 –Zircaloy-4 systems. Gas gap conductance, film boiling heat transfer uncertainty, pulse width, fuel-specific heat capacity, and cladding-specific heat capacity were identified as important parameters. Variations in figures of merit resulting from changes to pulse width and the material thermal properties indicate that time-dependent heat transfer rates are significant for safety-relevant mechanical parameters due to the time dependence of cladding ductility and pellet-cladding mechanical interaction loading. The results suggest that the thermal-hydraulic factors have a nonnegligible influence on the thermomechanical solution and vice versa. Tight coupling of both sets of physics is recommended to improve prediction of fuel behavior during RIAs. Highlights include the following: 1. The RELAP5-3D thermal-hydraulic code and the BISON fuel performance code are tightly coupled for simulation of RIA transients with energy depositions at the Zircaloy-4 cladding failure threshold. 2. Departure from nucleate boiling occurred for all simulated cases. Due to the ductility of fresh fuel, substantial ballooning occurred in most cases. 3. Gas gap conductance, fuel-specific heat capacity, cladding-specific heat capacity, transient pulse width, and film boiling heat transfer were the dominant thermal factors impacting the safety figures of merit at energy depositions.

Critical Heat Flux (CHF)↗
DOE OSTI2025

A review on the state of thermal hydraulics research on air ingress scenarios in High-Temperature Gas-cooled Reactors following a D-LOFC

With the expectation of near-immediate carbon neutrality, widespread implementation of proven High-Temperature Gas-cooled Reactors (HTGRs) embodies a viable solution pathway given their inherent, passive safety features and high thermal efficiency. This study provides an overview of the current state of research involving the thermal hydraulics associated with air ingress from a depressurized loss of forced cooling (D-LOFC) in HTGRs. Accurately characterizing and predicting the physical phenomena underlying air ingress is of paramount concern, as the integrity of the fuel and core graphite support structures are threatened by the presence of oxygen. Broadly speaking, the air ingress scenario can be delineated into three main stages: (1) Depressurization, (2) Density-Driven Flow, and (3) Natural Convection. In tandem with the underlying fundamental theory, this review collates and synthesizes the existing body of contemporary research concerning the air ingress scenario following a D-LOFC. As evinced by this review, our current understanding and predictive abilities have benefited from extensive research, predominantly concentrated on the rate of air ingestion into the core. Here, additional research is necessary to holistically capture the phenomenology of an air ingress scenario following a D-LOFC by considering an additional variable: the oxygen content of the ingressing air. The latter variable requires investigation into the complex interactions of the fully integrated system. Additionally, while numerical tools are evolving domestically through the Nuclear Energy Advanced Modeling and Simulation program, a sufficiently validated code remains absent.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗