SEARCH · Search NASA

Results for “System 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 451 records · Page 25

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Recent progress in the design of the K-DEMO divertor

The preliminary conceptual design of the Korean fusion demonstration reactor (K-DEMO) with a major radius of 6.8 m and the fusion power of 2200 MW has been studied since 2012. The overall configuration of the K-DEMO divetor system based on the ITER-like water-cooled tungsten technology is a double-null type symmetric divertor subdivided into 32 toroidal modules for the vertical maintenance. A detached divertor scenario with impurity seeding was considered as the primary approach for the power exhaust to reduce the peak heat flux lower than the engineering limit of 10 MW/m 2 . The power exhaust performance at the scrap off layer was estimated by using UEDGE-2D code, a two-dimensional fluid transport code for collisional edge plasma and neutral species like N, Ne, and Ar. Particle and heat flux on inboard and outboard divertor targets were calculated for the detached cased depending on parameters such as the impurity seeding rate, pumping rate, and the pedestal density. On the other hand, a magnetic solution like X-divertor, snowflake divertor, and super X-divertor to expand the plasma wet area was considered for K-DEMO since the detached divertor increasing a radiation fraction by impurity seeding might be able to be unstable. However, the extremely high current of poloidal coils was required more than the engineering limit, 20 MA, to form magnetic field lines for the alternative divertors. Based on the physical calculation of the edge plasma, engineering analyses were carried out to find out the thermal and structural reliability. The thermo-hydraulic analysis confirmed thermal stability, whether all comprising materials are operating within their allowable temperature windows when the case of the peak heat flux is set to 10 MW/m 2 on the outboard divertor target. The response surface optimization method derived two optimal design candidates employing two kinds of heat sink materials, respectively: the reduced activation ferritic martensitic (RAFM) steel and CuCrZr alloy. Additionally, the drawbacks and merits of the two materials were definite. The optimal design with applying RAFM steel was vulnerable to withstand thermal and mechanical loads since low thermal conductivity caused too thin thickness of the heat sink. On the other hand, the CuCrZr alloy has critical drawbacks in terms of activation and radioactive waste despite its high thermal conductivity. Meanwhile, preliminary electromagnetic (EM) analysis was carried out to estimate the EM loads caused by the abnormal behaviors of plasma since EM loads are one of the most critical external loads for designing a DEMO divertor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Validation study of sodium pool fire modeling efforts in $\mathrm{MELCOR}$ and $\mathrm{SPHINCS}$ codes

Discharge of sodium coolant into containment from a sodium-cooled fast reactor vessel can occur in the event of a pipe leak or break. In this situation, some of the liquid sodium droplets discharged from the coolant system will react with oxygen in the air before reaching the containment. This phase of the event is normally termed the sodium spray fire phase. Unreacted sodium droplets pool on the containment floor where continued reaction with containment atmospheric oxygen occurs. This phase of the event is normally termed the sodium pool fire phase. Both phases of these sodium-oxygen reactions (or fires) are important to model because of the heat addition and aerosol generation that occur. Any fission products trapped in the sodium coolant may also be released during this progression of events, which if released from containment could pose a health risk to workers and the public. The paper describes progress of an international collaborative research in the area of the sodium fire modeling in the sodium-cooled fast reactors between the United States and Japan under the framework of the Civil Nuclear Energy Research and Development Working Group. In this collaboration between Sandia National Laboratories and Japan Atomic Energy Agency, the validation basis for and modeling capabilities of sodium spray and pool fires in MELCOR of Sandia National Laboratories and SPHINCS of Japan Atomic Energy Agency are being enhanced. Here this study documents MELCOR and SPHINCS sodium pool fire model validation exercises against the JAEA’s sodium pool fire experiments, F7-1 and F7-2. The proposed enhancement of the sodium pool fire models in MELCOR through addition of thermal hydraulic and sodium spreading models that enable a better representation of experimental results is also described.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2026

Techno-Economic Analysis of Geologically Connected Seawater Air Conditioning (GeoSWAC) Concept for District Cooling at the University of Puerto Rico at Rio Piedras: Preprint

At the University of Puerto Rico at Rio Piedras, a central chilled water plant supplies cooling to several campus buildings, contributing significantly to electricity demand during daytime peak hours, particularly in the summer months. These operational challenges are exacerbated by Puerto Rico's tropical rainforest climate and a power grid vulnerable to frequent disruptions caused by hurricanes and tropical storms. This study presents a techno-economic analysis of the existing chilled water plant serving four representative campus buildings and introduces a conceptual alternative: the Geologically connected Seawater Air Conditioning (GeoSWAC) system. GeoSWAC leverages stable low temperatures of deep ocean water (~1 km depth), hydraulically connected to an inland well, to deliver cooling without the use of vapor-compression refrigeration. Using modeled annual cooling loads and chiller performance data, capital costs, energy consumption, and levelized cost of cooling (LCOC) were evaluated for both systems. While GeoSWAC showed higher capital costs than the chiller-based scenario, operational costs were significantly lower at $26k-$53k annually, resulting in a lower LCOC between $2.3/MWh and $8.2/MWh compared to $30.2/MWh-$33.6/MWh for the chiller scenario. These results suggest that the GeoSWAC system offers a promising, low-energy, and climate-resilient alternative for large-scale cooling in tropical coastal environments, with significant potential to reduce peak electricity demand and improve long-term system reliability.

15 GEOTHERMAL ENERGY↗

DOE OSTI2021

Application of NEAMS Codes to Capture MSR Phenomena

This report documents work completed in FY21 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program’s Molten Salt Reactor (MSR) Application Drivers activity at Argonne. The common focus was on identifying the modeling and simulation functional requirements for designing, licensing, and operating MSRs and applying those capabilities already developed in NEAMS tools to example problems of interest. The four main parts of this report each focuses on a specific area of simulation physics as it relates to MSR phenomena: fuel evolution, chemistry, computational fluid dynamics (CFD), and systems analysis. In terms of fuel evolution, which includes depletion, decay, on-line separations, and transmutation, the current state of computational capabilities for modeling this behavior in liquid-fueled molten salt reactor is discussed. Some of the functional requirements to accomplish the various applications of MSR fuel depletion modeling are highlighted, followed by a summary of recent approaches and code development activities. The chemistry functional requirements were discussed in the context of several applications of high importance for MSRs, such as corrosion, salt chemistry, and species transport. Each of these types of chemistry modeling have considerable impact on various aspects of reactor applications, including informing on reactor designs, improving operational efficiencies, analyzing safety and reactivity concerns, and estimating the mechanistic source term of the reactor. A brief overview is also provided on code development activities ongoing under NEAMS relevant to chemistry modeling of MSRs. In terms of CFD applications, the state-of-the-art spectral element code, Nek5000 was used to model the fluid dynamics within a full core of the Molten Salt Fast Reactor (MSFR) concept designed as part of the Euratom EVOL project. This concept was selected as the challenge problem because of its similar features to several U.S. industry concepts. The goal was to model some of the fast MSR design challenges, including potential large internal re-circulations, the need of accurate tracking of delayed neutron precursors (DNP), and the lack of relevant thermal-hydraulics models/correlations, etc. Therefore, a series of CFD models were created for the MSFR core cavity using a k – τ model two-equation model for the turbulence. These first full core results demonstrated that any potential recirculation zones could be properly identified with the current NEAMS CFD capabilities. The development of these models will also set the stage for future testing of Nek5000’s functionalities to model other MSR thermal-fluid phenomena. Lastly, the validation of SAM against experimental data from the Molten Salt Reactor Experiment, which started in FY19, continues with the inclusion of modeling the reactivity insertion tests. This involved using SAM and its point kinetics model for flowing fuel salt to recreate the time dependent power changes and response after positive reactivity insertions at the 1, 5, and 8 MWt power levels. Through these exercises, several code modifications were suggested to the SAM development team and accommodated to enable closer agreement with solid technical and physical justifications. These include adding a moderator reactivity feedback coefficient as an available input and modifying the solution approach for the point kinetics model. Additional SAM development suggestions for flowing fuel MSRs include adding the capability to allow the moderator power change proportionally with the reactor power and enabling specification of the power and DNP distributions separately.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2022

Thermal feedback coupling in a transient Monte Carlo high-order/low-order scheme

Transient simulations of nuclear systems face the computational challenge of resolving both space and time during reactivity changes. A common strategy for tackling this issue is to split the neutron flux into shape and amplitude functions. This split can be solved with high- order/low-order methods. While this multi-fidelity approach has traditionally been reserved for deterministic methods, it is also possible to implement in Monte Carlo as an efficient alternative to Dynamic Monte Carlo. This work implements the frequency transform method with thermal feedback in high-order/low-order Monte Carlo by blending static coupling methods such as single-batch Monte Carlo, with a simple thermal-fluids calculation. While previous work focused solely on prescribed transients, the addition of time-dependent thermal-fluids allows transients to be self-propagating. Tests were run in a fluids-initiated transient to showcase the basic functionalities of this methodology. Preliminary results behave as expected, paving the way for studying more sophisticated thermal-fluids coupling methods. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

THERMAL SAFETY MARGIN CALCULATION OF THE MP-2 EXPERIMENT IN THE ADVANCED TEST REACTOR

The Mini-Plate 2 (MP-2) irradiation test is a fueled experiment designed for irradiation in multiple test locations in the Advanced Test Reactor (ATR) at the Idaho National Laboratory (INL). The experiment is a drop-in test where small aluminum-clad fuel plate samples (mini plates) are cooled directly by the ATR Primary Coolant System (PCS) water. The MP-2 fuel plate experiment will be irradiated in several different irradiation locations of the ATR. This fueled experiment contains aluminum-clad fuel mini plates consisting of monolithic U-Mo. Four different types of fuel plates with fuel meat thickness and cladding are part of the MP-2 test. A thermal analysis has been performed on the MP-2 experiment. A method for calculating Departure from Nucleate Boiling Ratio (DNBR) and Flow Instability Ratio (FIR) during a reactivity transient using the commercial finite element and heat transfer code ABAQUS is discussed. At the start of an ATR cycle the heat generation rate of the fueled experiment is high and the heat rate multiplier from the outer shim control cylinders is low, while the reverse is true at the end of the ATR cycle. Thermal analyses at 10 day increments during the cycle calculate the DNBR and FIR during a reactivity transient. This technique calculates DNBR for the fuel plate surfaces and FIR for all water components for each finite element surface and node at various times during the ATR cycle. Heat rates vary with time during the transient calculation that are provided by a detailed physics analysis. Oxide growth on the fuel plates is also incorporated Results from the transient calculations are displayed with the ABAQUS post processor. By calculating these parameters at each location in the finite element model, conservatism is replaced with accuracy. This allows for a greater margin for the thermal hydraulic safety parameters.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

Techno-Economic Analysis of Geologically Connected Seawater Air Conditioning (GeoSWAC) Concept for District Cooling at the University of Puerto Rico at Rio Piedras

15 GEOTHERMAL ENERGY↗

DOE OSTI2024

Towards a NEAMS-based high-fidelity model of the MARVEL reactor

This report outlines the progress of Idaho National Laboratory in developing a high-fidelity and high-resolution model of the Microreactor Applications Research Validation and Evaluation reactor. The model was developed under the Nuclear Energy Advanced Modeling and Simulation microreactor application driver at Idaho National Laboratory. The overarching objective of this activity is the development of a high-fidelity multiphysics MARVEL model using NEAMS tools, and to verify and validate NEAMS tools against MARVEL reference simulation and experimental data, respectively. This is a unique opportunity to conduct multiphysics analysis on a soon-to-be-deployed microreactor. This multiphysics model developed under the NEAMS-funded INL microreactor application driver leverages three single-physics models coupled via the MOOSE’s MultiApp and Transfer systems. The latter systems enable in-memory data transfer between MOOSE-based and MOOSE-wrapped applications. The first single-physics model, that functions as main application, leverages Griffin to model the neutron transport in the core through the discontinuous finite element (DFEM) discrete ordinates solver (SN). Several optimization flags that were developed by the Griffin developer team were beta-tested to enhance the solver’s performance. These include the combined use of using_average_xs and update_averaged_xs_on that enable to avoid expensive on-the-fly cross sections evaluations at each linear iterations in favor of evaluations of the macroscopic cross sections at each Picard iteration. The second single-physics model uses BISON to handle solid heat transfer and asymptotic hydrogen redistribution analysis in the fuel. While the model returns consistent results for the temperature and hydrogen distribution in the fuel, a mismatch was noticed in the calculated temperature in the reflector due to the value of the gap conductance used in our model. Ongoing investigations are being performed to assess the origin of this discrepancy. Finally, the System Analysis Module (SAM) was used to model the flow of the sodium-potassium eutectic in the primary loop. A first verification was also performed showing good agreement in terms of mass flow rate and inlet temperature. All mesh files were generated using the MOOSE Reactor module, removing the need for external meshing tools. Notably, this workscope represents one of the initial applications of the MOOSE Reactor module for modeling highly irregular geometries. The use of the reactor module significantly streamlined the mesh generation process. The full multiphysics mode, that combines all the single physics models, was leveraged to conduct initial steady-state multiphysics simulations to compute power, and temperature distribution in the reactor. Initial testing was performed for transient simulations as well. In this case, the new checkpoint restart capability for eigenvalue calculations was tested showing the capability for streamlined restart of transient calculations. Future work will focus on improving the fidelity of the model by performing comprehensive code-to-code comparisons. For instance, the full-core Griffin neutronics model will be benchmarked against MCNP reference results, that were provided by the MARVEL design team. Additionally, the SAM T/H model will be verified against reference RELAP-5 results for selected accident scenarios. Besides code-to-code verification exercises, the model fidelity will be improved by replacing the single-channel SAM model with a more complex SAM-Pronghorn coupled model, in which the sub-channel capability is deployed to obtain radial temperature resolution in the coolant. This model will be developed in synergy with the NEAMS thermal hydraulics team.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Recent Developments at the Raft River Geothermal Field

The Raft River geothermal field, located in Cassia County in southwestern Idaho, is the site of a Department of Energy Enhanced Geothermal System project. U.S. Geothermal, Inc. currently produces about 11 MWe from Precambrian metamorphic rocks. These lie beneath ~5,000 ft of Quaternary and Tertiary volcaniclastic and volcanic deposits. Maximum temperatures range from 271°F to 300°F. Well RRG-9 ST1, the well targeted for stimulation is located approximately 1 mile south of the main bore field. The open hole section of the well, from 5,551 to 5,900 ft MD, consists of Precambrian Elba Quartzite, the stimulation target, granite and minor diabase. Prior to setting the casing acoustic, gamma ray, and density logs were run. After completing the well, a step rate/step down test was conducted. The maximum injection rate achieved was 18 bpm at a wellhead pressure of 1,150 psig. A borehole televiewer run in the open hole section showed evidence of more than eighty fractures. The majority of these fractures trend from N20⁰W to N20°E and dip from 40° to 60°W. Permeable fractures were encountered in the Elba Quartzite at 5,640-5,660 ft MD. Analysis of the injection test indicates that the minimum in-situ principal stress in this zone is 3,050-3,200 psi, corresponding to a fracture gradient of 0.59-0.62 psi/ft. A discrete fracture network model was developed using measured and inferred fracture orientations, distributions and dimensions. A three-phase stimulation program is proposed for RRG-9 ST-1. During the first two stages, water at 140°F, and later 40°F, will be injected to pre-condition and thermally fracture the reservoir. Here, the third stage will consist of a high rate, large volume conventional hydraulic stimulation.

Bradford, Jacob↗

DOE OSTI2025

Thermo-hydraulic steam pipe models for district heating simulations: Simplifications to balance accuracy and simulation speed

Steam piping networks are essential for optimizing performance in industrial processes and district heating systems. However, dynamic models that balance thermo-hydraulic accuracy with computational efficiency remain limited. In response, this paper presents a new discretized steam pipe model based on the plug flow approach, capturing key thermo-hydraulic behaviors while simplifying steam phase change processes. Implemented in Modelica, the model accurately calculates temperature and pressure distributions along steam pipelines. To improve computational efficiency for district-scale simulations, five model simplifications are introduced: lumped thermo-hydraulic functions, empirical correlations, fluid state approximations, steady-state dynamics and inclusion of flow derivatives. These simplified models achieve 85%-98% accuracy in predicting pressure drop and condensation losses, including dynamic condensate behavior during pipe warm-up—a factor often overlooked in existing models. The models support diverse network configurations, scaling effectively to systems with multiple distribution pipes and connected building loads. Discrete models provide detailed insights but exhibit a cubic increase in simulation time as the network scales by N connected building O(N 2.42 ). In contrast, lumped models simulate 10–28 times faster than discrete, offering quadratic scaling of simulation time O(N 1.73 ). However, they still require 6 times more computation time than a lossless network, highlighting the inherent computational challenges of modeling compressible fluid flow. In conclusion, the steady-state lumped variant, with its near-linear scalability in computational time O(N 1.01 ), emerges as an efficient solution for preliminary design evaluations and extensive parametric studies.

15 GEOTHERMAL ENERGY↗

DOE OSTI2021

Overview of the fusion nuclear science facility, a credible break-in step on the path to fusion energy

The Fusion Nuclear Science Facility (FNSF) is examined here as part of a two step program from ITER to commercial power plants. This first step is considered mandatory to establish the materials and component database in the real fusion in-service environment before proceeding to larger electricity producing facilities. The FNSF can be shown to make tremendous advances beyond ITER, toward a power plant, particularly in plasma duration and fusion nuclear environment. A moderate FNSF is studied in detail, which does not generate net electricity, but does reach the power plant blanket operating temperatures. The full poloidal Dual Coolant Lead Lithium (DCLL) blanket is chosen, with alternates being the Helium Cooled Lead Lithium (HCLL) and Helium Cooled Ceramic Breeder/Pebble Bed (HCCB/PB). Several power plant relevant choices are made in order to follow the philosophy of targeted technologies. Any fusion core component must be qualified by fusion relevant neutron testing and highly integrated non-nuclear testing before it can be installed on the FNSF in order to avoid the high probability of constant failures in a plasma-vacuum system. A range of missions for the FNSF, or any fusion nuclear facility on the path toward fusion power plants, are established and characterized by several metrics. A conservative physics strategy is pursued to accommodate the transition to ultra-long plasma pulses, and parameters are chosen to represent the power plant regime to the extent possible. An operating space is identified, and from this, one point is chosen for further detailed analysis, with R = 4.8 m, a = 1.2 m, IP = 7.9 MA, BT = 7.5 T, βN Gr = 0.9, fBS = 0.52, q95 = 6.0, H98 ∼1.0, and Q = 4.0. The operating space is shown to be robust to parameter variations. A program is established for the FNSF to show how the missions for the facility are met, with a He/H, a DD and 5 DT phases. The facility requires ∼25 years to complete its DT operation, including 7.8 years of neutron production, and the remaining spent on inspections and maintenance. The DD phase is critical to establish the ultra-long plasma pulse lengths. The blanket testing strategy is examined, and shows that many sectors have penetrations for heating and current drive (H/CD), diagnostics, or Test Blanket Modules (TBMs). The hot cell is a critical facility element in order for the FNSF to perform its function of developing the in-service material and component database. The pre-FNSF R&D is laid out in terms of priority topics, with the FNSF phases driving the time-lines for R&D completion. A series of detailed technical assessments of the FNSF operating point are reported in this issue, showing the credibility of such a step, and more detailed emphasis on R&D items to pursue. These include nuclear analysis, thermo-mechanics and thermal-hydraulics, liquid metal thermal hydraulics, transient thermo-mechanics, tritium analysis, maintenance assessment, magnet specification and analysis, materials assessments, core and scrape-off layer (SOL)/divertor plasma examinations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Coupled Investigation of Fracture Permeability Impact on Reservoir Stress and Seismic Slip Behavior (Final Technical Report)

Enhanced Geothermal Systems (EGS) produce clean energy by circulating fluid through hot rock deep underground and bringing that heat to the surface to generate electricity. For this process to work reliably, fluids must be able to move efficiently through networks of natural or engineered fractures in the rock. Enhancing and maintaining subsurface permeability over time is essential for sustainable energy production. However, fluid injection changes the underground temperature, pressure, rock stress, and chemistry, which can alter permeability and sometimes trigger earthquakes. Predicting these interconnected processes remains a key challenge. To address this, we combined high-temperature laboratory experiments with high-fidelity simulations to better understand how fractures in geothermal reservoirs evolve over time. Our experiments measured how fractures respond to stress, slip, slip rate, and chemical reactions under geothermal conditions. These data were integrated into coupled thermal-hydrological-mechanical-chemical and earthquake (THMC+E) models tailored to the Utah FORGE site. The validated modeling framework improves predictions of reservoir performance and seismic response and helps guide operational decisions. This work reduces technical risk and strengthens the scientific foundation needed to make geothermal energy a reliable and scalable clean energy resource.

15 GEOTHERMAL ENERGY↗

DOE OSTI2022

MELCOR Accident Progression and Source Term Demonstration Calculations for a Heat Pipe Reactor

MELCOR is an integrated thermal hydraulics, accident progression, and source term code for reactor safety analysis that has been developed at Sandia National Laboratories for the United States Nuclear Regulatory Commission (NRC) since the early 1980s. Though MELCOR originated as a light water reactor (LWR) code, development and modernization efforts have expanded its application scope to includ e non-LWR reactor concepts. Current MELCOR development efforts include providing the NRC with the analytical capabilities to support regulatory readiness for licensing non-LWR techno logies under Strategy 2 of the NRC?s near- term Implementation Action Plans. Beginning with the Next Generation Nuclear Project (NGNP), MELCOR has undergone a range of enha ncements to provide analytical capabilities for modeling the spectrum of advanced non-LWR concepts. This report describes the generic plant model developed to demonstrate MELCOR capabilities to perform heat pipe reactor (HPR) safety evaluations. The generic plant mode l is based on a publicly-available Los Alamos National Laboratory (LANL) Megapower design as modified in the Idaho National Laboratory (INL) Design A description. For plant aspects (e.g., reactor building size and leak rate) that are not described in the LANL and INL references , the analysts made assumptions needed to construct a MELCOR full-plant model. The HP R uses high assay, low-enrichment uranium (HALEU) fuel with steel cladding that uses heat pipes to transfer heat to a secondary Brayton air cycle. The core region is surrounded by a stainless-steel shroud, alumina reflector, core barrel and boron carbide neutron shield. The reactor is secured inside a below-grade cavity, with the operating floor located above the cavity. Example calculations are performed to show the plant response and MELCOR capabilities to characterize a range of accident conditions. The accidents selected for evaluation consider a range of degraded and failed modes of operation for key safety functions providing re activity control, the primary and secondary system heat removal, and the effectiveness of th e confinement natural circulation flow into the reactor cavity (i.e., a flow blockage).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Evaporator Temperature Transient Testing of a High-Performance Sodium Filled Heat Pipe

Heat pipes are two-phase heat transfer devices that enable passive removal of heat from the reactor core to the power conversion system in heat pipe-cooled microreactor designs. Experimental investigations of heat pipe transients are needed for technology demonstration, verification and validation of numerical codes, and the establishment of regulatory requirements. The Single Primary Heat Extraction and Removal Emulator (SPHERE) facility at Idaho National Laboratory (INL) serves as a platform for evaluating the dynamic response of high-temperature heat pipes under a variety of operating conditions. The present work details the experimental investigation of a high-performance, defined as over 2 kW sodium heat pipe subjected to rapid input power fluctuations induced by sudden changes in the evaporator temperature setpoint. In addition, the heat pipe was subjected to an asymmetrical heat load where a subset of heaters operated at 30% and 70% below their nominal power. These experimental conditions were chosen to simulate thermal and operational stresses expected to be encountered in microreactors to provide data on heat pipe behavior during such important transient events. Key data and performance metrics, including time series of temperatures and strains, axial temperature profiles, thermal response times, and heat transfer capabilities, the thermal output over thermal input, were reported and discussed. The results highlight the resilience of heat pipes, revealing their potential to maintain thermal stability and efficiency under varying power loads. Lastly, the paper concludes with a discussion on the significance of the results and their implications for future research.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

Development and Validation of SAM Multi-dimensional Flow Model for Thermal Mixing and Stratification Modeling

Thermal mixing and stratification in large pools or enclosures are very important phenomena that are critical to nuclear reactor safety. Because of the wide ranges of time and length scales associated with such phenomena, accurate modeling and simulation of thermal mixing and stratification remain as the key unresolved, challenging problems for reactor transient analyses. In traditional system analysis codes, simplified zero-dimensional (0-D) models are widely used for their high numerical efficiency, but they generally suffer from very limited prediction accuracies or range of applicability. Like traditional system analysis codes, the current version of the SAM code has implemented such simplified 0-D and one-dimensional mixing models. On the other hand, high-resolution Computational Fluid Dynamics (CFD) tools are often used to model complex thermal mixing and stratification phenomena. They are, however, generally numerically expensive, and they require large amounts of computational resources. It is therefore desirable to implement advanced and efficient thermal mixing and stratification modeling capabilities embedded in a system analysis code. This approach will improve the accuracy of reactor safety analyses when thermal mixing and stratification are involved, and also avoid using the large computational resources needed for high-resolution CFD analysis. Currently, with the support of the U.S. DOE Office of Nuclear Energy’s Nuclear Energy Advanced Modeling and Simulation program, an effort has been launched to develop and implement a multi-dimensional flow model in the system analysis code SAM, and demonstrate its applications to model thermal mixing and stratification phenomena in large enclosures. The main outcomes of this research and development activity are summarized in this report, which presents an attempt to include a built-in advanced multi-dimensional flow model in a system analysis code with the focus on overcoming the simulation challenges of thermal mixing and stratification phenomena. In this report, we start with the introduction of existing SAM code capabilities to simulate thermal mixing and stratification phenomena, which is followed by a short summary of the multi-dimensional model implemented in the SAM code, including both the physical model and the Finite Element Method code implementation. In this study, two options were implemented in the SAM code to model turbulent flows: a relatively simple built-in turbulence model and an interface to accept externally computed turbulent viscosities (e.g., from a high- resolution CFD simulation). Code validation studies on this newly added capability were then carried out to compare SAM simulation results with experimental data from the SUPERCAVNA facility, which was designed to study the complex flow recirculation and thermal stratification phenomena relevant to sodium fast reactor designs. In this study, one transient and two steady-state test cases were used for code validation. Different approaches have been used to model the complex turbulence flow fields in the SUPERCAVNA facility. A highly simplified zero-equation turbulence model was first used, but it was determined that it is too simple to capture the complex turbulence flow fields in these test cases. Subsequently, the code validation continued with the use of turbulent viscosity data from high-resolution STAR-CCM+ CFD simulations to improve the accuracy of the results. Using this approach, the SAM simulation results showed very good agreement with both the SUPERCAVNA experimental data and STAR-CCM+ simulation results. In this report, we demonstrate the development, implementation, and successful validation of a multi-dimensional flow model in the SAM code, which aims to improve the simulation accuracy for complex thermal mixing and stratification phenomena. Lessons have also been learned, including that in cases where the flow fields are not well predicted by the zero-equation model, the iinput of turbulent viscosities from an external source can enhance the overall predictive capabilities needed to accurately capture complex thermal-hydraulic phenomena. Therefore, future research will be needed to further improve the code’s capabilities, such as by developing a more efficient and robust approach to capture the turbulence effects in the SAM code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

Nek5000 developments in support of industry and the NRC

This year, the Nuclear Energy Advanced Modeling Simulation program (NEAMS) thermal-hydraulics verification and validation (V&V) work has focused in three areas of Nek5000 V&V-driven development. First, in a close collaborative effort with the U. S. Nuclear Regulatory Commission (NRC) staff, we have continued V&V efforts for the HYMERES-2 project using the OECD/NEA sponsored testing in the PSI PANDA facility. This year’s focus of ANL-NRC collaboration involves Nek5000 setups and validation for a range of problems relevant to and including the HYMERES-2 benchmark from PSI. The primary outcome of this year efforts is a more efficient geometry and inlet modeling simplification after a careful sensitivity study of the inlet profiles and pipe geometries. The resulting modeling choice of a short recycling/fully-developed turbulent inlet is within the experimental uncertainty estimate. This finding simplifies the next step of the cross-V&V HYMERES-2 project. In addition, the ANL team continue to provide assistance to the NRC staff in the form of Nek5000 application support in general and on the use of the HPC platforms of ALCF and INL in particular. This supports the NRC’s assessment of Nek5000 for use with the NRC Blue CRAB code suite. Second, we have implemented and tested more robust model of URANS, namely the k – τ model, a variant of the k-ω model, along with other improvements to RANS Nek5000 modeling in general. Because of its demonstrated robustness and stability, the k – τ model is the only RANS model that has been implemented in the new GPU version of the Nek5000 code, nekRS. Lastly, we report the initial implementation of Jacobian-free Newton Krylov approach to the direct Newton method for steady fluid solvers aimed at acceleration of RANS modeling and at IC improvement for LES campaigns. Also leveraging the Exascale Computing Project (ECP) ANL/CEED & SMR team’s software development effort to support NEAMS problems at large scale of the advanced computing architectures, NekRS, a GPU variant of Nek5000, built on top of kernels from libParanumal using OCCA for portability, has been successfully run on the full system of Summit (4608 nodes, 27648 GPUs).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Multimetallic Layered Composites (MMLCs) for Rapid, Economical Advanced Reactor Deployment (Final Report)

This project focused on the development of multi-metallic layered composites (MMLCs) for advanced fission reactor technologies. There are many instances where one alloy or material simply cannot meet all the demands thrown at it by a reactor system, or cannot allow it to perform as strongly as one would like. Instead of focusing all our effort on developing one perfect alloy, we seek to leverage the design principle of “separation of functionality,” used in many other arenas in design, to boost performance beyond single alloys alone. One illustrative example shows the power of this approach for molten salt-cooled reactors: A three meter tall, three meter diameter reactor vessel made of Incoloy 800 was quoted at $\$$500k in 2018. A Hastelloy N vessel was quoted at $\$$5M. An MMLC vessel, in which a layer of Hastelloy N would be weld-overlaid onto Incoloy 800, was quoted at $\$$700k, and it would achieve the same performance. The potential economic gains of leveraging this approach are therefore substantial. At a minimum, each MMLC would contain one core structural layer and one coolant-facing corrosion-resistant layer. Sometimes, MMLCs required buffer layers, as the structural and corrosion-resistant layers were metallurgically incompatible. In other words, they didn’t always play nice, thus separating layers compatible with both functioned as intermediaries to keep the composite together. However, in doing so we inevitably produce new interfaces, where new issues can arise. Therefore, this project focused on what happens at these interfaces from a combination of high temperatures, irradiation, corrosion, and time. After all, a reactor makes money when it is operating, and outages of any kind erode its economic viability. First, we set out to experimentally prove that MMLCs for at least two advanced reactor systems can be made, today, in US domestic facilities. In this respect we were successful – one MMLC (a Ni-201/Incoloy 800H composite) was successfully made and drawn into two-inch coolant piping. Others were attempted, though new issues relating to cracking in vanadium layers for one and radiation damage performance of the corrosion-resistant layer in another prevented us from moving further in those specific arenas – these are engineering problems which deserve continued focus after this project. Additional experimental work focused on long-term corrosion testing of the outermost layers of the salt-cooled and liquid lead-cooled MMLC concepts, which would then be fed into predictions of how long the MMLCs could last. Next, computational (thermodynamics and atomistic) simulation studies studied how much we expect the interfaces to “blend,” due to the mixing action of neutron irradiation. This eats into both the margin for the structural layer of each MMLC, as dilution from the corrosion-resistant layer into the structural layer would decrease the total load-bearing capacity of an MMLC of finite size. On the other hand, dilution of the corrosion-resistant layer into the structural layer further reduced the margin of corrodible material, reducing the lifetime of the MMLC or necessitating extra thickness to be imparted to the MMLC to meet its functional requirements. Work here focused on irradiation-induced segregation to predict new phases which may embrittle the MMLCs, as well as quantifying irradiation-induced mixing at each interface. The results showed that mixing is expected, but it is both steady and therefore predictable, and not lifetime-limiting for most MMLC concepts – it simply has to be accounted for in calculations of reactor performance when utilizing an MMLC. Then, full-core simulations using the experimentally-derived corrosion data, the computationally discovered irradiation-induced mixing data (partially validated by experiment), and existing, benchmarked core designs for large and small sized reactor concepts (one salt-cooled, one lead-cooled) were conducted to quantify any expansion of reactor operating envelopes achieved by utilizing these MMLCs. This new framework, called REX (Reactor Envelope Expansion), incorporates a combination of core neutronics, thermal hydraulics, and the material performance data derived from this project to see how using an MMLC expands advanced fission reactor operating envelopes. It was discovered that in some cases, MMLC utilization does indeed increase the maximum operating temperatures and cycle lengths of reactor concepts, while in other cases it does not. Finally, our tech-to-market (T2M) strategy was not necessarily to create specific embodiments of MMLCs for immediate sale (because getting into the nuclear market is incredibly slow and laden with regulation, this is a long-term goal), but rather immediate stimulation of US industry using the design approach of MMLCs derived from this project. In this respect we were successful, as one of the PhD students funded on this project co-founded Allium Engineering, Inc., which created a stainless steel / low-alloy steel MMLC to function as chloride corrosion-resistant rebar for embedding into concrete structures. Allium Engineering continues to be successful, having recently opened their first factory as of this writing.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Code Benchmark of Depressurized Conduction Cooldown Transient in the High Temperature Test Facility

This paper presents results from modeling of a depressurized conduction cooldown (DCC) transient at the High Temperature Test Facility (HTTF) as part of the OECD-NEA Thermal Hydraulics Code Validation Benchmark for High-Temperature Gas-Cooled Reactors using HTTF Data . This paper briefly describes the benchmark and the models being used. It then presents a comparison of steady state and transient results based on the Problem 2 Exercise 1A and 1B definitions. We compare block and helium temperature distributions, mass flow distribution, and energy balance in steady state. All models show comparable mass flow distributions and energy balances. The temperatures within the core and outer regions are comparable in all models too, but inner reflector temperatures can vary significantly. Despite that, we find that the models are in good agreement for the full-power steady state. In the DCC, we look at block temperature at the core midplane and RCCS water exit temperature. The INL and ANL models are found to be in excellent agreement with one another on block temperature over time, while the agreement when the KAERI and NRG models are added into consideration is good. Differences in the transient heat removal from the RCCS cause the differences in block temperature over time in these models. The CNL models show similar trends to the INL, ANL, KAERI, and NRG models, but the temperatures are high because the volumes used in calculating the average temperature include the heater rods in the CNL models only. The HUN-REN model shows results that suggest significantly lower heat removal in the RCCS which merit further investigation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗