SEARCH · Engineering Papers

Results for “THERMAL HYDRAULICS”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 469 records · Page 26

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Fluoride-Cooled High-Temperature Pebble-Bed Reactor Reference Plant Model

In this report we present work performed in Fiscal Year 2022 that demonstrates the modeling and simulation of a fully coupled neutronics thermal hydraulics reference plant model for a fluoride-cooled high-temperature pebble-bed reactor. The multiphysics model is developed on the Nuclear Regulatory Commission’s Comprehensive Reactor Analysis Bundle (BlueCRAB) available on the Idaho National Laboratory’s high-performance computer, which natively and seamlessly couples Griffin, Pronghorn, and the BISON Multiphysics Object-Oriented Simulation Environment based applications. Griffin provides reactor physics capabilities, including depletion to the equilibrium core, k-eigenvalue, adjoint, and transient. The unique direct equilibrium core capability in Griffin is based on a streamline methodology to spatially deplete the pebbles into burnup groups. Pronghorn solves the porous medium equations for the fluid regions and conduction in the solid regions and incorporates a fluidic diode model to simulate the transition from forced to natural convection during accident scenarios. MOOSE modules solves thermal conduction problems for the pebbles and tristructural isotropic in the pebble-bed core, thus providing the fuel and moderator spatial fields for each pebble burnup group. The neutronics feedback relies primarily on fuel, moderator, and reflector temperatures as as well as the FLiBe salt density. Here, we present results for the uncoupled equilibrium core and perform comparisons to equivalent Monte Carlo models. The power distributions and kinetic parameters obtained with Griffin are consistent with those computed with Griffin. We demonstrate a noticeable improvement with the use of discrete ordinates method (SN) transport. The coupled steady-state equilibrium core provides the initial condition for two time-dependent problems: a control rod withdrawal event and an unprotected loss of flow event. In both cases, the reactor design is self-stabilizing and the solutions are consistent with the expected physics. Although this model is prototypical regarding BlueCRAB’s capabilities, its results are consistent with published work by Kairos Power and other research entities. Significant improvements to the model are planned in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRSDate not supplied

Nuclear Space System Analysis and Modelling (NSSAM): A Software Tool to Efficiently Analyze the Design Space of Space Reactor Systems

Space reactors have the potential to play a key role in future NASA exploration activities due to their capability to enable sustainable power and advanced propulsion systems. To enable assessment of the space reactor design space, the nuclear space system analysis and modelling (NSSAM) software was developed by Analytical Mechanics Associates. NSSAM leverages a scalable and extensible software architecture which automates reactor analysis to perform coupled engine-reactor and reactor physics-thermal hydraulics calculations. This allows space reactor systems to be evaluated by a wider number of users with a consistent analysis approach to compare designs. NSSAM has been developed with multiple use cases to tailor the analysis to the level of detail desired by the user and computing resources. This summary overviews the NSSAM architecture and development approach, current capabilities (including design variants and use cases) and analysis approach for reactor and system component models.

nuclear thermal propulsion↗

NASA NTRS2008

Gravity Scaling of a Power Reactor Water Shield

Water based reactor shielding is being considered as an affordable option for use on initial lunar surface power systems. Heat dissipation in the shield from nuclear sources must be rejected by an auxiliary thermal hydraulic cooling system. The mechanism for transferring heat through the shield is natural convection between the core surface and an array of thermosyphon radiator elements. Natural convection in a 100 kWt lunar surface reactor shield design has been previously evaluated at lower power levels (Pearson, 2007). The current baseline assumes that 5.5 kW are dissipated in the water shield, the preponderance on the core surface, but with some volumetric heating in the naturally circulating water as well. This power is rejected by a radiator located above the shield with a surface temperature of 370 K. A similarity analysis on a water-based reactor shield is presented examining the effect of gravity on free convection between a radiation shield inner vessel and a radiation shield outer vessel boundaries. Two approaches established similarity: 1) direct scaling of Rayleigh number equates gravity-surface heat flux products, 2) temperature difference between the wall and thermal boundary layer held constant on Earth and the Moon. Nussult number for natural convection (laminar and turbulent) is assumed of form Nu = CRa(sup n). These combined results estimate similarity conditions under Earth and Lunar gravities. The influence of reduced gravity on the performance of thermosyphon heat pipes is also examined.

Reid, Robert S.↗

DOE OSTI2025

CFD modeling of natural circulation in LiCl-KCl molten salt closed loop

Characterizing flow within a molten salt closed-loop system is crucial for assessing system requirements, evaluating performance, and identifying potential flaws. Direct flow measurement using instrumentation is challenging due to extreme environmental conditions and the limitations associated with measuring molten salt flow under natural convection. Here, this study aims to provide comprehensive insights into the thermal-hydraulic behavior of a closed loop, with a particular focus on temperature distribution and velocity prediction. The Computational Fluid Dynamics (CFD) model demonstrated the capability to effectively simulate and predict both temperature distributions and flow velocities within the molten salt loop. The CFD model's predictive capability was validated by its ability to replicate temperature measurements under varying boundary conditions. The analysis revealed that the CFD model tends to underpredict temperatures in the cold leg and overpredict them in the hot leg, highlighting the need for continuous model refinement and acknowledging the limitations of using a steady-state approach. Furthermore, the potential of using external temperature measurements to estimate internal molten salt temperatures and predict flow velocity was explored, revealing that this approach could introduce up to a 5.5% error in flow velocity calculations. Line probes mapping temperature distributions across the tube's cross-section and molten salt provided valuable insights into temperature gradients, emphasizing the need for a thermal conductivity equation for molten salt with lower uncertainty to achieve more accurate temperature predictions of the system.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Transient Modeling and Simulation of a Generic Stable Salt Reactor

A SAM system-level model of a generic stable salt reactor has been developed to investigate thermal-hydraulic behavior and safety performance under steady and transient conditions. The model integrates information generated from a reactor physics analysis using PROTEUS and PERSENT, and a computation fluid dynamics (CFD) analysis using STAR-CCM+. A loose, iterative coupling scheme between PROTEUS and SAM is implemented to calculate the equilibrium power and temperature distributions in the steady-state critical core condition. The converged steady-state model is then used in PERSENT to calculate the four reactivity feedback temperature coefficients (Doppler, fuel density, coolant density, and core radial expansion) and kinetic parameters that are needed in SAM to model the temperature feedback effects in transient simulations. Within the fully enclosed liquid fuel pins, natural convection is the dominant heat transfer mechanism. The STAR-CCM+ model of the fuel pin considers conjugate heat transfer from the liquid fuel salt to the pin cladding and external reactor coolant. The CFD results of the axial and radial temperature profiles are used to empirically determine an effective fuel salt thermal conductivity in the SAM fuel pin model so that the temperatures predicted by the SAM model match as closely as possible the CFD results. In the central region of the fuel pin, the effective thermal conductivity is as high as similar to 60 times the physical fuel salt thermal conductivity. The whole-plant SAM model is then used to simulate an unprotected station blackout transient. The results of this simulation showed that the large negative fuel axial expansion reactivity feedback reduces fission power to similar to 2.4% nominal power. The core is cooled by natural circulation, which removes heat in the core to the emergency heat removal system, and ultimately, to the ambient. However, peak fuel salt and cladding temperatures can potentially reach as high as 1500 K, albeit briefly, if the shutdown mechanism fails to operate.

stable salt reactor; transient simulations; system↗

DOE OSTI2022

Analysis of HolosGen Sub-Scale Simulator with Plant Dynamics Code

The Subcritical Power Module Sub-scale Simulator (SPM-SS) has been designed and constructed by HolosGen LLC under the ARPA-E MEITNER program to simulate the thermal-hydraulic and heat transfer behavior of the full-scale Holos-Quad Subcritical Power Modules (SPMs). Four coupled SPMs, each rated at 5.5MW, form the Holos-Quad gas-cooled microreactor design. The SPM-SS represents a substantially scaled-down system with a power rating less than 40 kW, equipped with an electrically heated fuel cartridge heat exchanger, an electrically heated compressor heat exchanger, and a valve actuated turbine heat exchanger, in addition to a recuperator and a cooler heat exchanger. The fuel cartridge represents a portion of the full-scale SPM core, the compressor heat exchanger mimics the temperature changes resulting from the compressor’s turbomachinery inefficiencies, the turbine heat exchanger mimics the expansion process normally occurring through the turbine, while the recuperator and cooler heat exchangers complete the subscale simulator loop. The heaters equipping the fuel cartridge and the compressor heat exchangers are electronically controlled to simulate normal and off-normal SPM operating conditions. The full-scale Holos-Quad SPM design eliminates the traditional balance of plant and executes thermal-to-electric energy conversion by means of an intercooled Brayton cycle with decoupled compressor-turbine turbomachinery. The Holos-Quad full-scale design is equipped with a multi-stage axial Low- and High-Pressure compressor, and a multistage axial turbine. The SPM-SS is designed for testing and validation of selected components which are instead coupled by a traditional balance of plant. The SPMSS is not equipped with turbo-machinery (compressor and turbine) as the development of these components were excluded from the scope of work under the ARPA-E MEITNER funding program. The SPM-SS balance of plant enables modifications, replacement and testing of individual components with different working fluids and is designed to include the turbo-machinery components that will be developed in future research . The SPM-SS can be operated with different gases, variable mass-flow-rates, pressures, and temperatures to obtain test data for selected components, whose performance can be scaled to validate the computer model of the full-scale SPM at various conditions (e.g., start-up, transients conditions). The SPM-SS can operate at the maximum Holos-Quad design pressure of 7 MPa, and a maximum temperature limited to 650 °C by the electrical heaters. Several SPM-SS tests have been conducted and analyzed with the Plant Dynamics Code (PDC) developed at the Argonne National Laboratory (ANL). These tests aimed at validating the PDC modeled predictions of the full-scale Holos-Quad design with data from selected SPM-SS components. In order to address SPM-SS specific characteristics, such as components heat losses and absence of turbomachinery components, some modifications to the PDC have been implemented to factor the design differences from the full-scale Holos-Quad SPM to the SPM-SS. As the PDC offers capabilities to analyze systems with different working fluids, air, nitrogen, and helium were utilized as the SPM-SS working fluids. Air was utilized to fine-tune the SPM-SS Systems Structures and Components (SSCs), nitrogen was utilized to pressure test the SPM-SS loop at the SPMs design maximum pressure of 7MPa. Helium was utilized as the working fluid circulating through the SPM-SS SSCs for specific tests to validate the PDC predictions of the fuel cartridge heat exchanger. SPM-SS tests data were also analyzed with both the steady-state and transient analysis capabilities offered by the PDC. This report describes the PDC analysis of the SPM-SS tests data, including the necessary code modifications and comparison of the code results with the experimental data. Based on the results, a discussion is presented on how the analysis supports design and transient calculations of the full-scale Holos-Quad microreactor. Also based on the results of this work, recommendations are made for future optimizations of the SPM-SS components and PDC model development needs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

High Flux Microchannel Receiver Development (Final Report)

The opportunity addressed in this project was using microscale heat exchanger technology for improving solar thermal central receiver efficiency. This will result in smaller, less expensive components and reduce the levelized cost of electricity from concentrated solar thermal power. This micro-pin solar receiver (MSR) directly heats high pressure (25 MPa) supercritical carbon dioxide (sCO 2 ) from 550°C to 720°C with an incident flux of at least 120 W/cm 2 at a receiver efficiency of >90% (assuming availability of 95% absorptivity coating), and receiver cost target of less than $\$$150/kWt. Realization of the MSR required an understanding of thermal hydraulic performance of sCO 2 in micro-pin passages, mechanical behavior of high nickel alloys at high temperatures including joints, advanced fabrication methods of high nickel alloys (etching, machining, brazing, welding, diffusion bonding), and the cost implications of these when the concept is scaled to production volumes. The technology can be used with gases (sCO 2 , helium, etc.) or with liquids such as molten salt or liquid metal. Given the high pressure associated with sCO 2 , its use demonstrates a “worst case” scenario, and the successful development of the MSR for this application will result in the development of a technology also applicable to a range of other heat transfer fluids that will enable next generation CSP plants, including support of the Gen3CSP program.

14 SOLAR ENERGY↗

DOE OSTI2025

An Overview of the Molten Salt Thermal Properties Database–Thermophysical, Version 4.0 (MSTDB-TP V.4.0)

A central repository of thermophysical and thermochemical properties of molten salt compositions of relevance to molten salt reactors (MSRs) is vital in supporting the broad community of MSR developers, who are at various stages of developing and deploying their reactor designs. In general, these MSR designs differ significantly from developer to developer (e.g., with respect to the hardness of the neutron spectra, level of fissile loading, target multicomponent temperatures and power levels, and moderating capabilities). Therefore, the fuel and coolant salts being considered vary greatly: they may be chlorides or fluorides, they utilize different actinides at different ratios, and the cations in the melt are selected based on perceived advantages and disadvantages. Considering the general need for thermal properties, and the vastness of the array of potential candidate salt mixtures, the Molten Salt Thermal Properties Database (MSTDB) was initiated in 2018 with the goal of providing thermophysical and thermochemical characterization of key molten salt compounds and mixtures across their temperature and compositional domains. The MSTDB is thus divided into the thermophysical arm (MSTDB-TP) and the thermochemical arm (MSTDB-TC). The MSTDB is an effort funded by the Department of Energy, Office of Nuclear Energy (DOE-NE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, and the MSR Campaign. This report provides an overview of the MSTDB-TP v4.0 in terms of the data contained within, the state of the tools used to access the data, the availability of predictive models that leverage the raw data in the database, the preliminary status of developmental efforts that are currently underway, and an account of future goals for MSTDB-TP. The primary goal for the update from MSTDB-TP v.3.1 to v4.0 was the incorporation of surface tension data into the database; this property is important for thermal hydraulics modeling and species transport in other tools that have been developed under the NEAMS program. A breakdown of the surface tension data that have been added into MSTDB-TP v4.0 is provided herein, and the manner in which the quality of the data has been assessed is also documented. For MSTDB-TP v4.0, newly published thermophysical property data—primarily from collaborative experimental efforts under the MSR Campaign—have been incorporated into the database, and the resulting expansion is documented here. Because of the size to which MSTDB-TP has grown, the raw data format has now been recast into JavaScript Object Notation (JSON) format for easier connection with the MSTDB-TP application programming interface (API). Saline; the pre-existing comma-separated value (CSV) format has been deprecated but is still maintained, accessible, and up to date. As a final effort in packaging the MSTDB-TP v4.0 update, the graphical user interface (GUI) for MSTDB has been updated to allow full accessibility to the density and viscosity predictive models, which are based on Redlich-Kister expansions of MSTDB-TP raw data. Some other major aspects of this report, in terms of preliminary and future work, include: (1) documentation of the formalism and preliminary testing of a kinetic theory model that may act as a predictive model for thermal conductivity; (2) documentation of the candidate predictive models that may be considered in the future for surface tension, making use of the surface tension data now in MSTDB-TP v4.0; (3) a preliminary account of a data collection process that will enable the filling of additional gaps within MSTDB-TP, namely with data which have been collected computationally (e.g., through ab initio molecular dynamics).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRSDate not supplied

The Use of Immersion Rigs for High Temperature Hydrogen Exposure Testing Within the Nuclear Thermal Rocket Element Environmental Simulator (NTREES)

The Nuclear Thermal Rocket Element Environmental Simulator (NTREES) was designed to perform non-nuclear testing of nuclear thermal rocket (NTR) fuel elements and materials. NTREES can simulate the thermal hydraulic environment within an NTR and has been leveraged to provide critical information on component performance and material compatibility. NTREES has the ability to expose test articles to pressurized Hydrogen gas up to 1000 psig and 225 g/s. NTREES tests are operated in a deoxygenated environment of Nitrogen or Argon pressurized up to 1000 psig. NTREES simulates fission heat through induction heating powered by a 1.2 MW AC power supply, and NTREES has recorded temperatures exceeding 3700K. To perform hot Hydrogen exposure tests on materials more rapidly and earlier in the fuel development process, NTREES must be able to test small material samples on the order of 1 square inch in volume, individual fuel assembly subcomponents, and particle materials. To meet these needs, multiple testing rigs have been developed to immerse test specimens in hot Hydrogen, with each "immersion rig" varying in its capability of sample type and performance range in order to meet the broadest range of testing needs. Herein is a discussion of a selection of these “immersion rigs.”

Nuclear Thermal Propulsion↗

DOE OSTI2024

Effect of Wettability Gradient on the Scale Formation in Falling Film Flow

Mitigation of scale formation and performance degradation remains a vital challenge for falling film evaporators in various industries. In this work, an experimental study of falling film flow on a horizontal tube is conducted to investigate the effects of wettability gradients on thermal, hydraulic, and fouling behavior. It is revealed that certain hydrophobic coating patterns, such as strip, ring, and grid patterns, lead to unwetted heat transfer area, which results in decreased heat transfer compared to fully wetted plain tube. By adjusting the geometry and position of the wettability gradient, the hybrid coating demonstrates improved heat transfer performance. Based on the characteristics of horizontal tube falling film flow, impinging jet, thin film flow, and liquid retention at the tube bottom, a hybrid coating pattern is developed to improve surface wetting and mitigate the scaling coverage. It is revealed that scale deposition is regulated by wettability gradient. Crystals tend to be dense and compact in hydrophilic areas, while they appear scattered or even absent in hydrophobic regions, depending on the dimension of the hydrophobic area. While at the hydrophilic/hydrophobic boundary, a noticeable scale thickness step is observed, which raises the potential for self-cleaning. Furthermore, the balance of minimization of scaling layer coverage and maximization of wetting area requires an optimal design in coating dimensions, for which a systemic study of both flow dynamics and fouling characteristics on the falling film is necessary in the future.

42 ENGINEERING↗

DOE OSTI2025

A Novel Framework for Performance Evaluation and Design Optimization of PCM Embedded Heat Exchangers for the Built Environment

This research sheds light on the performance evaluation and design optimization of PCM-HXs for the built environment, addressing several barriers to practical issues to PCM-HX commercialization such as modeling aspects (i.e., modeling expertise and computational / time investment, etc.), manufacturing aspects (i.e., at-scale manufacturing, cost assessments, etc.) and experimental performance assessment (i.e., reliable experimental data, assessment of multiple PCM-working fluid combinations, etc.). We present a novel, comprehensive, and experimentally-validated design optimization framework for PCM-HXs capable of simulating any PCM-HX geometry with reasonable accuracy and significant computational time savings when compared to traditional CFD-based design practices. The framework was validated for a wide range of PCM-HX configurations, including a design optimization for a domestic hot water heater application where TES partially replaces electrical heating input. The resulting PCM-HXs were found to deliver 34-68% of the total daily hot water supply with only 5-10% package volume increase from the water heater, thus within U.S. DOE targets for TES systems. To identify the most promising HXs for PCM applications, first-order geometry and cost analyses were conducted based on off-the-shelf HX products. As part of this work, 9 PCM-HX prototypes were manufactured using additive and conventional manufacturing methods. Detailed economy-of-scale assessments were conducted for the most promising PCM-HXs and were found to have a good outlook for the next 5-10 years. The PCM-HX design optimization framework was validated through comprehensive in-house experimental testing using newly-developed PCM-to-fluid test facilities. In total,10 total in-house component-level experiments were conducted using these prototypes, including 9 with water and 1 with refrigerant (R410A) as the working fluid. It was found that the framework can successfully predict experimental thermal-hydraulic performance within ±10-20% the first time without manual design changes, eliminating the need for time-consuming and expensive prototyping efforts as part of the design process. As part of this work, a publicly-available PCM web tool was released which includes a PCM property database (531 PCMs) and PCM-HX modeling tool to assist the design community on common PCM-HX use-cases, e.g., single/multiple flow path(s) fluid-to-PCM and air-to-fluid-to-PCM configurations (https://ceeeweb.umd.edu/pcmapp/). This work will accelerate the design and time to market for next generation PCM-HXs.

25 ENERGY STORAGE↗

DOE OSTI2024

A Reduced-Order Model of a Nuclear Power Plant with Thermal Power Dispatch

This paper presents reduced-order modeling of thermal power dispatch (TPD) from a pressurized water reactor (PWR) for providing heat to nearby heat consuming industrial processes that seek to take advantage of nuclear heat to reduce carbon emissions. The reactor model includes the neutronics of the reactor core, thermal–hydraulics of the primary coolant cycle, and a three-lump model of the steam generator (SG). The secondary coolant cycle is represented with quasi-steady state mass and energy balance equations. The secondary cycle consists of a steam extraction system, high-pressure and low-pressure turbines, moisture separator and reheater, high-pressure and low-pressure feedwater heaters, deaerator, feedwater and condensate pumps, and a condenser. The steam produced by the SG is distributed between the turbines and the extraction steam line (XSL) that delivers steam to nearby industrial processes, such as production of clean hydrogen. The reduced-order simulator is verified by comparing predictions with results from separate validated steady-state and transient full-scope PWR simulators for TPD levels between 0% and 70% of the rated reactor power. All simulators indicate that the flow rate of steam in the main steam line and turbine systems decrease with increasing TPD, which causes a reduction in PWR electric power generation. The results are analyzed to assess the impact of TPD on system efficiency and feedwater flow control. Due to the simplicity of the proposed reduced-order model, it can be scaled to represent a PWR of any size with a few parametric changes. In the future, the proposed reduced-order model will be integrated into a power system model in a digital real-time simulator (DRTS) and physical hardware-in-the-loop simulations.

08 HYDROGEN↗

DOE OSTI2026

Dual Purpose – Heating & Cooling – Thermal Battery for Flexible and Energy-Efficient Heat Pump Systems

The integration of heat pumps with thermal energy storage (HP-TES) systems is gaining attention as a viable solution for managing peak building demand driven by immense cooling and heating loads. With growing reliance on renewable energy sources, thermal energy storage offers an excellent opportunity to mitigate mismatches in thermal load between energy supply and demand. The use of phase-change material (PCM) TES is especially promising, as PCMs offer significant latent energy storage capacity with smaller temperature glides in smaller volumes compared to other TES technologies. However, challenges arise because current HP-TES architectures can load-shift only cooling or heating, not both, requiring two systems and thus doubling cost, weight, and footprint. Furthermore, current research efforts lack specific tools and techniques to advance integrated systems from concept design to end-user application, focusing on only discharge performance. To address these challenges, this research proposes a dual-mode (heating and cooling) integrated HP-TES system that uses room-temperature PCM-TES as a high-temperature heat source in heating mode and a low-temperature heat sink in cooling mode, thereby reducing temperature lifts and compressor power. Design criteria for PCM-TES heat exchangers were developed, balancing thermal-hydraulic performance with practical constraints such as available building space and weight requirements along with PCM selection considerations such as shipping conditions, moisture exposure, and number of available cycles. A detailed transient model for HP-TES systems was developed to enable rapid annual performance assessments in any US climate zone.

42 ENGINEERING↗

NASA NTRS2015

Initial Operation of the Nuclear Thermal Rocket Element Environmental Simulator

The Nuclear Thermal Rocket Element Environmental Simulator (NTREES) facility is designed to perform realistic non-nuclear testing of nuclear thermal rocket (NTR) fuel elements and fuel materials. Although the NTREES facility cannot mimic the neutron and gamma environment of an operating NTR, it can simulate the thermal hydraulic environment within an NTR fuel element to provide critical information on material performance and compatibility. The NTREES facility has recently been upgraded such that the power capabilities of the facility have been increased significantly. At its present 1.2 MW power level, more prototypical fuel element temperatures nay now be reached. The new 1.2 MW induction heater consists of three physical units consisting of a transformer, rectifier, and inverter. This multiunit arrangement facilitated increasing the flexibility of the induction heater by more easily allowing variable frequency operation. Frequency ranges between 20 and 60 kHz can accommodated in the new induction heater allowing more representative power distributions to be generated within the test elements. The water cooling system was also upgraded to so as to be capable of removing 100% of the heat generated during testing In this new higher power configuration, NTREES will be capable of testing fuel elements and fuel materials at near-prototypic power densities. As checkout testing progressed and as higher power levels were achieved, several design deficiencies were discovered and fixed. Most of these design deficiencies were related to stray RF energy causing various components to encounter unexpected heating. Copper shielding around these components largely eliminated these problems. Other problems encountered involved unexpected movement in the coil due to electromagnetic forces and electrical arcing between the coil and a dummy test article. The coil movement and arcing which were encountered during the checkout testing effectively destroyed the induction coil in use at the time and resulted in NTREES being out of commission for a couple of months while a new stronger coil was procured. The new coil includes several additional pieces of support structure to prevent coil movement in the future. In addition, new insulating test article support components have been fabricated to prevent unexpected arcing to the test articles. Additional activities are also now underway to address ways in which the radial temperature profiles across test articles may be controlled such that they are more prototypical of what they would encounter in an operating nuclear engine. The causes of the temperature distribution problem are twofold. First, the fuel element test article is isolated in NTREES as opposed to being in the midst of many other mostly identical fuel elements in a nuclear engine. As a result, the fuel element heat flux boundary conditions in NTREES are far from adiabatic as would normally be the case in a reactor. Second, induction heating skews the power distribution such that power is preferentially deposited near the outside of the fuel element. Nuclear heating, conversely, deposits its power much more uniformly throughout the fuel element. Current studies are now looking at various schemes to adjust the amount of thermal radiation emitted from the fuel element surface so as to essentially vary the thermal boundary conditions on the test article. It is hoped that by properly adjusting the thermal boundary conditions on the fuel element test article, it may be possible to substantially correct for the inappropriate radial power distributions resulting from the induction heating so as to yield a more nearly correct temperature distribution throughout the fuel element.

Emrich, William J., Jr.↗

DOE OSTI2020

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↗

DOE OSTI2022

A digital twin approach to system-level fault detection and diagnosis for improved equipment health monitoring

Automating the task of fault detection and diagnosis is crucial in the effort to reduce the operation and maintenance cost in the nuclear industry. This paper describes a physics-based approach for system-level diagnosis in thermal-hydraulic systems in nuclear power plants. The inclusion of physics information allows for the creation of virtual sensors, which provide improved fault diagnosis capability. The physics information also serves to better constrain diagnostic solutions to the physical domain. As a demonstration, various test cases for fault diagnosis in a high-pressure feedwater system were considered. The use of virtual sensors allows constructing performance models for two first-point feedwater heaters which would not have been possible otherwise due to the limited sensor set. Real-time plant data provided by a utility partner were used to assess the diagnostic approach. The detection of an abnormal event immediate after a plant startup pointed to faulty behaviors in the two first-point feedwater heaters. Further, this double-blind fault diagnosis was subsequently confirmed by the plant operator. In addition, several simulated sensor fault events demonstrated the capability of our algorithms in detecting and discriminating sensor faults.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2023

Discrete element simulation of Pebble Bed Reactors on graphics processing units

Prediction of pebble positions in a Pebble Bed Reactor (PBR) is necessary for both reactor physics and thermal hydraulics simulations as the arrangement of pebbles has a significant impact on the resulting core power, coolant flow, and fuel temperature. Knowledge of pebble movement as the fuel is cycled through the core is also critical for predicting the fuel residence time and subsequently, the fuel burnup. Simulation with the Discrete Element Method (DEM) can provide knowledge of both the fuel packing and the fuel movement during cycling. Previous works that have performed 3D full-core DEM simulation of PBRs have used simplified models that neglect reflector wall features. This work employs a graphics processing unit (GPU)-enabled DEM code, Project Chrono, to analyze the differences in pebble packing and pebble velocities between a simplified smooth PBR reflector and a more realistic reflector that includes circular wall features. Additionally, a sensitivity study is performed on the depth of the wall features to ensure that crystallization is prevented. Project Chrono is also validated for PBR cycling applications using experimental data. It is found that wall features with a depth of at least 0.5 pebble diameters significantly reduce crystallization in the near-wall region, leading to discrepancies in both packing fraction and pebble velocity in this region compared to the simplified reflector models. These discrepancies are found to lead to roughly a 5–10% difference in the prediction of the near-wall porosity and a 10% difference in the prediction of the velocity of pebbles near the wall. As a result of these discrepancies, it is suggested that future DEM simulations of PBRs include wall features to reduce modeling errors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

Multiphysics analysis of fuel fragmentation, relocation, and dispersal susceptibility–Part 1: Overview and code coupling strategies

The US nuclear energy industry is investigating strategies to increase the reactor operating cycle to 24 months, resulting in peak rod average burnups exceeding the current limit of 62 GWd/tU. This increase will in turn increase the probability of fuel fragmentation, relocation, and dispersal (FFRD) in the event of a loss-of-coolant accident (LOCA). This effort couples multiple codes to (1) evaluate full-core power histories for high-burnup fuel operated in a Westinghouse 4-loop pressurized water reactor, (2) model a postulated large-break LOCA, and (3) calculate the mass of fuel susceptible to FFRD. This paper, the first of three describing the work, focuses on code coupling strategies and FFRD susceptibility calculations. The other two companion papers focus on code-specific designs and analyses. Three codes were used in this work. VERA was used to calculate steady-state power histories, TRACE was used to model the transient thermal hydraulics, and BISON was used to model steady-state and transient fuel performance and cladding failure. Herein, several fuel pulverization models were used to calculate FFRD susceptibility in failed rods. Depending on the cladding failure/fuel pulverization model combination, the core-wide FFRD susceptibility during the postulated LOCA range from 0 to over 5,000 kg.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗