Search NASA⌕ Search

SEARCH · Search NASA

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.

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

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Thermal Hydraulic Experimental Test Article - Fiscal Year 2022 Final Report

The Thermal Hydraulic Experimental Test Article (THETA) is a METL vessel experiment designed for testing and validating sodium fast reactor components and phenomena. THETA has been scaled using a non-dimensional Richardson number approach to represent temperature distributions during nominal and loss of flow conditions in a sodium fast reactor (SFR), this analysis was detailed in the THETA FY19 report. The facility is being constructed with versatility in mind, allowing for the installation of various immersion heaters, heat pipes, and heat exchangers without significant facility modification. THETA was designed in collaboration with systems code experts to inform the geometry and sensor placement to acquire the highest value code validation data.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

MARVEL Thermal-hydraulics and PCAT

Provide a brief overview of MARVEL thermal-hydraulics, focusing on MARVEL integral test faciliy PCAT.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2024

Hierarchical Bayesian modeling for Inverse Uncertainty Quantification of system thermal-hydraulics code using critical flow experimental data

The best estimate plus uncertainty methodology in nuclear system thermal-hydraulic studies necessitates a comprehensive understanding of uncertainties in system code predictions. The forward uncertainty quantification (UQ) process involves the propagation of input uncertainties through the computational models to obtain uncertainties in the outputs. To this end, achieving an accurate estimation of input uncertainties is important, which is the focus of inverse UQ (IUQ). Traditionally, research in Bayesian IUQ within the nuclear engineering domain has largely relied on single-level Bayesian inference. While being effective for relatively small datasets, this approach encounters limitations for cases with large datasets. The use of a single-level model may prove inefficient, as the resultant posterior distributions can significantly differ when distinct subsets of data are employed. To address this issue, we employ an hierarchical Bayesian model for IUQ. Furthermore, this approach involves organizing observations into different groups based on the test conditions, thereby accommodating varying calibration parameters across these distinct groups. In this study, we developed and implemented a hierarchical Bayesian IUQ method to consider the grouping effect of critical flow measurement data from various geometries. Comparing the outcomes of IUQ under different selections of test data using hierarchical Bayesian IUQ against those obtained from single-level Bayesian IUQ, the forward propagation of hierarchical Bayesian IUQ results demonstrates a notably improved agreement with the experimental data.

42 ENGINEERING↗
DOE OSTI2023

Investigation of thermal hydraulic behavior of the High Temperature Test Facility's lower plenum via large eddy simulation

A high-fidelity computational fluid dynamics (CFD) analysis was performed using the Large Eddy Simulation (LES) model for the lower plenum of the High–Temperature Test Facility (HTTF), a ¼ scale test facility of the modular high temperature gas-cooled reactor (MHTGR) managed by Oregon State University. In most next–generation nuclear reactors, thermal stress due to thermal striping is one of the risks to be curiously considered. This is also true for HTGRs, especially since the exhaust helium gas temperature is high. In order to evaluate these risks and performance, organizations in the United States led by the OECD NEA are conducting a thermal hydraulic code benchmark for HTGR, and the test facility used for this benchmark is HTTF. HTTF can perform experiments in both normal and accident situations and provide high-quality experimental data. However, it is difficult to provide sufficient data for benchmarking through experiments, and there is a problem with the reliability of CFD analysis results based on Reynolds–averaged Navier–Stokes to analyze thermal hydraulic behavior without verification. To solve this problem, high-fidelity 3-D CFD analysis was performed using the LES model for HTTF. It was also verified that the LES model can properly simulate this jet mixing phenomenon via a unit cell test that provides experimental information. As a result of CFD analysis, the lower the dependency of the sub-grid scale model, the closer to the actual analysis result. In the case of unit cell test CFD analysis and HTTF CFD analysis, the volume-averaged sub-grid scale model dependency was calculated to be 13.0% and 9.16%, respectively. As a result of HTTF analysis, quantitative data of the fluid inside the HTTF lower plenum was provided in this paper. As a result of qualitative analysis, the temperature was highest at the center of the lower plenum, while the temperature fluctuation was highest near the edge of the lower plenum wall. The power spectral density of temperature was analyzed via fast Fourier transform (FFT) for specific points on the center and side of the lower plenum. FFT results did not reveal specific frequency-dominant temperature fluctuations in the center part. It was confirmed that the temperature power spectral density (PSD) at the top increased from the center to the wake. The vortex was visualized using the well-known scalar Q-criterion, and as a result, the closer to the outlet duct, the greater the influence of the mainstream, so that the inflow jet vortex was dissipated and mixed at the top of the lower plenum. Additionally, FFT analysis was performed on the support structure near the corner of the lower plenum with large temperature fluctuations, and as a result, it was confirmed that the temperature fluctuation of the flow did not have a significant effect near the corner wall. In addition, the vortices generated from the lower plenum to the outlet duct were identified in this paper. It is considered that the quantitative and qualitative results presented in this paper will serve as reference data for the benchmark.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Pronghorn: A Multidimensional Coarse Mesh Application for Advanced Reactor Thermal-Hydraulics

This paper presents an overview of Pronghorn, a multiscale thermal-hydraulic (T/H) application developed by Idaho National Laboratory and the University of California, Berkeley. Pronghorn, built on the open-source finite element Multiphysics Object-Oriented Simulation Environment (MOOSE), leverages state-of-the-art physical models, numerical methods, and nonlinear solvers to deliver fast-running advanced reactor T/H simulation capabilities within a modern software engineering environment. This work summarizes the physical models, multiphysics and multiscale coupling, and numerical discretization in Pronghorn with emphasis on our initial target application to pebble bed reactors (PBRs). A diverse set of applications are shown to depressurized natural circulation in the SANA experiments, forced convection in the Pebble Bed Modular Reactor, three-dimensional (3-D)/one-dimensional coupling of Pronghorn and RELAP-7 systems T/H for loop analysis in the High Temperature Reactor Power Module, and forced convection in the Mark-1 Pebble Bed Fluoride-Salt-Cooled High-Temperature Reactor. A multiphysics coupling of Pronghorn, RELAP-7, and Griffin deterministic neutronics for a gas-cooled PBR demonstrates the capability of the MOOSE framework for reactor design calculations. These applications highlight the verification and validation underlying Pronghorn’s software development while emphasizing features that improve upon capabilities offered by legacy tools in areas such as 3-D unstructured meshing, physics modeling, and multiphysics coupling.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Multi-Physics Investigation of a Natural Circulation Molten Salt Micro-Reactor that Utilizes an Experimental In-Pile Device to Improve Core Physics and System Thermal-Hydraulic Performance.

INL employee PhD Dissertation - The Molten Salt Reactor (MSR) concept is a rapidly evolving Generation IV design that has recently attracted favorable attention due to the potential for reducing waste generation, realizing passive safety features, and seizing on the opportunity for cost effective economics. This thesis investigates the performance benefit of a new device invented by the doctoral candidate. The device is referred to as a Wrapped Helix around an Inclined Plane or WHIP. The WHIP is protected under a provisional patent filed with the USA Patent Office on 28 September 2021 under application number 63/261,776, BEA docket number BA-1254. The WHIP device can be located in-core or near-core to promoted enhanced thermal-hydraulics and neutronics performance. While the WHIP can be employed in a variety of solid or liquid fuel designs, this thesis investigates the device’s benefit in the application of a natural-circulation, micro-molten salt nuclear battery concept (MsNB). This thesis will specifically investigate the temperature coefficient of the MsNB fuel (FLiNaK) using novel temperature sensitivity techniques unique to Serpent particle transport code, evaluate the WHIP’s thermal-hydraulic and neutronic performance effects using both established (STAR-CCM+) and novel (Python code developed to estimate the circulation effective delayed neutron fraction, ßef f ) analytical and numerical methods, evaluate the neutron noise behavior of the MsNB and how the WHIP may alter the character of the MsNB’s transfer function, and how the WHIP affects the autonomous, load-following performance under transient power conditions using Python code developed by the candidate. The sum body of this work, in part, has been published in three journal articles as the timing of the provisional patent process has allowed. Results show that creative utilization of WHIP engineering design and function reduces reactor system volume, fuel loading, control/stability in the buoyant regime

buoyant flow↗
DOE OSTI2023

Preliminary Monte Carlo and Thermal Hydraulic Analysis using a Hybrid ETF-Corrected-Diffusion Prediction Block

This paper builds upon previous work to accelerate the Picard iteration (PI) method typically applied for coupled Monte Carlo-Thermal hydraulic (MC-TH) solutions. Previously, the use of the generalized transfer functions (GTFs) to predict variation in macroscopic cross sections following a perturbation in TH properties was demonstrated for a subset of simple 3D problems. In addition, the reduced-order transport prediction block relied on the first order perturbation (FOP) method, which was shown to have computational overheads. Recent work replaced the FOP block with a 1-group nodal diffusion solver to eliminate these overheads. While the use of diffusion is desirable for large-scale problems, the new solver introduces significant homogenization error. This work aims to address this issue by using the Jacobian-Free Newton Krylov (JFNK) method to generate a set of super homogenization (SPH) factors to improve the accuracy of the diffusion solution. The SPH factors will be used in conjunction with an improved cross section prediction method – the expanded transfer function (ETF) method – to produce a highly accurate flux prediction for an axial 1D boiling water reactor (BWR) pincell following a large perturbation in moderator density. The ETF-corrected diffusion (ETF-CD) block is shown to be highly accurate for the 1D test case. Future work will investigate the accuracy of the method for a realistic 3D pressurized water reactor core.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Physics-constrained fault diagnosis framework for monitoring a multi-component thermal hydraulic system

A method for diagnosing faults includes receiving a system description of a thermal hydraulic system, the system description indicating a plurality of components and sensors. The method also includes constructing, based on physical conservations laws and using the system description, a plurality of physics-based models for the plurality of components, each of the plurality of physics-based models including unknown parameters. The method further includes receiving historical measurements and calibrating the physics-based models by calculating the unknown parameters of each of the physics-based models using the historical measurements to produce calibrated models. The method also includes receiving sensor measurements of the sensors, and calculating residuals corresponding to differences between measurements predicted by the plurality of calibrated models and the sensor measurements. The method also includes determining, based on the calculated residuals, a fault of a component or a sensor, and generating an alert indicating the fault.

Nguyen, Tat Nghia↗
DOE OSTI2022

Accelerated coupled Monte Carlo-Thermal hydraulic calculations using a hybrid GTF-diffusion-based prediction block: first results

Accurate predictions of spatial power and temperature distributions require the coupling of a neutron transport solver with a thermal-hydraulic (TH) feedback. Nowadays, Monte Carlo (MC) codes are widely coupled to TH solvers, typically via a Picard iteration (PI) method, due to the higher fidelity that such frameworks can produce. To speed up a PI, a prediction step can produce an improved initial guess for a source distribution and feed it to the MC code. Recent work investigated a prediction step that uses generalized transfer functions (GTFs) to predict the macroscopic cross sections' variations following a perturbation in TH properties, such as coolant density. The previous method also relied on first order perturbation (FOP) theory to predict perturbed power profiles, rather than using an expensive MC iterate. The implemented FOP method relied on generating a fission matrix from which the forward and adjoint Eigenmodes were extracted and later used to for power calculations. The generation of the fission matrix can introduce a significant computational overhead, therefore undermining the performance of the proposed hybrid technique when applied to high-dimensional problems, e.g., full core calculations. This work attempts to improve the GTF-FOP prediction step by replacing the FOP solver with a nodal diffusion solver, thus eliminating the need to calculate a fission matrix. The GTF-diffusion step was tested for various moderator density perturbations. In each case, the predicted power distribution showed good agreement with the reference case. The latter is attributed to the generally good prediction of most spatially distributed macroscopic cross sections, except the transport cross section, which will become the focus of future work. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Accelerated Coupled Monte Carlo-Thermal Hydraulic Calculations using a Hybrid GTF-Diffusion-based Prediction Block: First Results

Accurate predictions of spatial power and temperature distributions require the coupling of a neutron transport solver with a thermal-hydraulic (TH) feedback. Nowadays, Monte Carlo (MC) codes are widely coupled to TH solvers, typically via a Picard iteration (PI) method, due to the higher fidelity that such frameworks can produce. To speed up a PI, a prediction step can produce an improved initial guess for a source distribution and feed it to the MC code. Recent work [1, 2] investigated a prediction step that uses generalized transfer functions (GTFs) to predict the macroscopic cross sections’ variations following a perturbation in TH properties, such as coolant density. The previous method also relied on first order perturbation (FOP) theory to predict perturbed power profiles, rather than using an expensive MC iterate. The implemented FOP method relied on generating a fission matrix from which the forward and adjoint eigenmodes were extracted and later used to for power calculations. The generation of the fission matrix can introduce a significant computational overhead, therefore undermining the performance of the proposed hybrid technique when applied to high-dimensional problems, e.g., full core calculations. This work attempts to improve the GTF-FOP prediction step by replacing the FOP solver with a nodal diffusion solver, thus eliminating the need to calculate a fission matrix. The GTF-diffusion step was tested for various moderator density perturbations. In each case, the predicted power distribution showed good agreement with the reference case. The latter is attributed to the generally good prediction of most spatially distributed macroscopic cross sections, except the transport cross section, which will become the focus of future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Development of a calculation framework coupling Monte Carlo neutronics, subchannel thermal-hydraulics and 1.5D fuel thermo-mechanics

A coupled calculation system consisting of the Monte Carlo neutronic solver SERPENT and the fuel behaviour solver TRANSURANUS was extended with the sub-channel thermal-hydraulic solver SUBCHANFLOW to form a coupled calculation framework for multi-physics reactor analyses. The coupling between the 3 solvers was accomplished via a driver program written in Perl that runs the solvers, processes and transfers the output data between the solvers and monitors the convergence. While some features of the coupling are still under development, a steady-state coupled calculation resulting in convergence was performed for a single hexagonal fuel assembly to demonstrate that the base functionalities work as expected. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Physics-constrained fault diagnosis framework for monitoring a standalone component of a thermal hydraulic system

A method for diagnosing faults includes receiving a description of a component of a thermal hydraulic system, where the description also indicates one or more sensors of the component. The method also includes constructing, based on a physical conservation law and using the description, a physics-based model describing operation of the component, the physics-based model including one or more unknown parameters. The method further includes calibrating the physics-based model by calculating the one or more unknown parameters using historical measurements to produce a calibrated model. Further, the method includes receiving sensor measurements captured by the one or more sensors, and calculating residuals corresponding to differences between measurements predicted by the calibrated model and the sensor measurements. The method also includes determining, based on the calculated residuals, a fault of the component or of a sensor of the one or more sensors, and generating an alert indicating the fault.

Nguyen, Tat Nghia↗
DOE OSTI2023

Multiphysics analysis of fuel fragmentation, relocation, and dispersal susceptibility–Part 3: Thermal hydraulic evaluation of large break LOCA under high-burnup conditions

Increasing the peak rod average burnup of pressurized water reactor (PWR) fuel beyond 62 GWd/tU may increase fuel fragmentation, relocation, and dispersal (FFRD) susceptibility during a large break loss of coolant accident (LBLOCA). TRACE thermal hydraulic (TH) LBLOCA analyses were performed for a realistic 24-month high-burnup PWR equilibrium cycle, to inform subsequent transient BISON high-burnup FFRD susceptibility evaluations. Realistic LBLOCA systems behavior was first established by configuring to and comparing with the BEMUSE OECD LBLOCA benchmark. Fuel and operating conditions were then applied from high-burnup VERA depletion calculations. LBLOCA simulations were performed for 281 selected high-burnup rods, for which transient TH boundary conditions were collected for later use in BISON. The TRACE results indicated that rod linear heat rate (rather than burnup) is the main predictor of peak cladding temperature (PCT) during the event. PCT typically occurred at a local burnup lower than the rod-average burnup, especially for twice-burned fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Thermal-Hydraulic-Mechanical Modeling and Simulation of Sodium-Potassium–Cooled MARVEL Microreactor Core

The U.S. Department of Energy's Microreactor Program, with Idaho National Laboratory's development of a nuclear microreactor applications test bed named MARVEL, aims to support R&D for the deployment of small, transportable reactors across civilian, industrial, and defense sectors. The MARVEL microreactor, an 85-kWth thermal fission reactor, incorporates TRIGA nuclear fuel and a sodium-potassium eutectic as its primary coolant, designed for safety and efficiency, with natural circulation eliminating the risk of critical heat flux conditions. The reliance on natural circulation for primary cooling means the reactor avoids using fuel spacers to minimize core pressure drop, which could disrupt the primary coolant's natural flow. However, the reactor core’s tight P/D ratio of 1.05, in the absence of fuel spacers, could pose a risk of fuel rod contact and increased peak cladding temperatures. To ensure the reactor safety, this study conducted computational modeling and simulations to investigate the reactor's thermal-hydraulic-mechanical characteristics, including the reactor core heat transfer coefficients, the potential for rod-to-rod contact, and assessed its impact on peak cladding temperature and overall reactor safety. The computational analyses of the MARVEL microreactor core revealed that the thermal deformation of fuel rods under worst-case scenario may lead to the fuel rod contact, but the peak cladding temperatures will remain significantly lower than the safety criteria, ensuring the safety operation reactor without fuel spacers under normal operating conditions.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2025

Thermal Hydraulics Analysis of a Divertor Monoblock Using SALAMANDER

Divertors are critical components in magnetic confinement fusion devices. One of the Divertor's crictial roles is absorbing the highest heat flux from the plasma. The Divertor accomplishes this through the use of divertor monoblocks and cooling channels. This work presents a thermal hydraulic analysis of a divertor monoblock using the Software for Advanced Large-scale Analysis of MAgnetic confinement for Numerical Design, Engineering and Research (SALAMANDER). SALAMANDER is an open source tool developed at Idaho National Laboratory for conducting multiphysics and multscale analysis of Plasma Facing Components (PFCs). For the monoblock, solid heat conduction is modeled in the block with a constant heat flux of 1E7 W/m^2 applied to the top of the block to represent peak heat flux from the plasma. The cooling channel (which is pressurized water at a Reynolds number of ~1,000,000) is modeled using a k-epsilon turbulence model with standard wall functions. The heat transfer between the solid and fluid domain is modeled using the Dittus-Boelter correlation for the convective heat transfer coefficient. This work aims to show the importance of a multiphysics modeling approach when performing simulations on PFCs.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Toward integrating high-fidelity CFD approaches in the thermal-hydraulic analysis of turbulent dry cask systems

Nuclear power plants have been supplying resilient and reliable electricity for decades, contributing to energy independence of the U.S.. However, nuclear waste management remains one of the most significant challenges in the industry. The safety of dry cask storage systems relies heavily on their thermal-hydraulic performance. Computational Fluid Dynamics (CFD) simulations are often used to demonstrate this performance and ensure that the system design meets safety standards. This study presents reduced numerical models for various types of dry cask systems. These numerical models can produce efficient and fast results based on the employed modeling strategies. Additionally, the study uses a novel approach to high-fidelity simulations to evaluate modeling assumptions in dry cask modeling. Large Eddy Simulations (LES) are used for this purpose, particularly in regions where fluid velocity is relatively high and the turbulence characteristics become important. Furthermore, the results of these high-fidelity simulations will enhance the interpretation of outcomes produced from a lower-fidelity CFD model.

CFD↗