SEARCH · Engineering Papers

Results for “Thermal fluid physics”

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 253 records · Page 14

All sources NASA NTRS DOE OSTI

DOE OSTI2023

What Controls the Entrainment Rate of Dry Buoyant Thermals with Varying Initial Aspect Ratio?

Here, this study uses theory and numerical simulations to analyze the nondimensional spreading rate a (change in radius with height) of buoyant thermals as they rise and entrain surrounding environmental fluid. A focus is on how a varies with initial thermal aspect ratio A r , defined as height divided by width of the initial buoyancy perturbation. An analytic equation for thermal ascent rate w t that depends on a is derived from the thermal-volume-averaged momentum bud get equation. The thermal top height when w t is maximum, defining a critical height z c , is inversely proportional to a. The height z c also corresponds to the thermal top height when buoyant fluid along the thermal’s vertical axis is fully replaced by entrained nonbuoyant environmental fluid rising from below the thermal. The time scale for this process is controlled by the vertical velocity of parcels rising upward through the thermal’s core. This parcel vertical velocity is approximated from Hill’s analytic spherical vortex, yielding an analytic inverse relation between a and A r . Physically, this a–A r relation is connected to changes in circulation as A r is modified. Numerical simulations of thermals with A r varied from 0.5 to 2 give a values close to the analytic theoretical relation, with a factor of ~3 decrease in a as A r is increased from 0.5 to 2. The theory also explains why a of initially spherical thermals from past laboratory and modeling studies is about 0.15. Overall, this study provides a theoretical underpinning for understanding the entrainment behavior of thermals, relevant to buoyantly driven atmospheric flows.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1982

Systems evaluation of thermal bus concepts

Thermal bus concepts, to provide a centralized thermal utility for large, multihundred kilowatt space platforms, were studied and the results are summarized. Concepts were generated, defined, and screened for inclusion in system level thermal bus trades. Parametric trade studies were conducted in order to define the operational envelope, performance, and physical characteristics of each. Two concepts were selected as offering the most promise for thermal bus development. All of four concepts involved two phase flow in order to meet the required isothermal nature of the thermal bus. Two of the concepts employ a mechanical means to circulate the working fluid, a liquid pump in one case and a vapor compressor in another. Another concept utilizes direct osmosis as the driving force of the thermal bus. The fourth concept was a high capacity monogroove heat pipe. After preliminary sizing and screening, three of these concepts were selected to carry into the trade studies. The monogroove heat pipe concept was deemed unsuitable for further consideration because of its heat transport limitations. One additional concept utilizing capillary forces to drive the working fluid was added. Parametric system level trade studies were performed. Sizing and weight calculations were performed for thermal bus sizes ranging from 5 to 350 kW and operating temperatures in the range of 4 to 120 C. System level considerations such as heat rejection and electrical power penalties and interface temperature losses were included in the weight calculations.

Stalmach, D. D.↗

DOE OSTI2025

Multi-physics Topology OPtimization and Additive Manufacturing for High-temperature Heat Exchangers

This research significantly advances the understanding of high-temperature heat exchanger design through an integrated approach that combines topology optimization (TO), triply periodic minimal surface (TPMS) structures, additive manufacturing (AM) and thermohydraulic testing. Each of these components contributes uniquely to a unified, high-performance design, fabrication and testing workflow. Topology optimization serves as the foundation of the design methodology by providing a systematic way to determine the most effective material layout for separating hot and cold fluids while maximizing thermal performance. The researchers introduced a novel three-material optimization framework using two density fields to represent hot fluid, cold fluid, and solid domains. This approach enables automated discovery of optimal shapes and flow paths that cannot be intuitively designed, especially under constraints imposed by manufacturing technologies. Furthermore, constraints such as minimal wall thickness and overhang angles were embedded into the optimization process, ensuring that resulting designs are not only thermally efficient but also manufacturable using modern additive techniques. In parallel, the study delves into the use of Gyroid-based TPMS geometries for constructing the core of the heat exchanger. TPMS structures are known for their high surface area, excellent fluid mixing capabilities, and minimal pressure drop characteristics. The researchers applied a data-driven modeling framework using Heteroscedastic Sparse Gaussian Process Regression (HSGPR) combined with genetic algorithms. This allowed for the rapid evaluation and optimization of key geometric parameters such as frequency, iso-value, and phase shift. The result was a set of Gyroid structures tailored for high heat transfer and low flow resistance, demonstrating clear improvements over conventional straight-channel designs. After the designing process, additive manufacturing played a critical role by turning these highly complex, optimized geometries into physical components. Utilizing Laser Powder Bed Fusion (LPBF) with Haynes 282, the study demonstrated the feasibility of fabricating these heat exchangers at high precision. Post-processing methods, including dilation-erosion operations, were applied to ensure local features adhered to self-supporting constraints. The fabricated structures were then subjected to thermohydraulic testing under conditions representative of supercritical CO 2 Brayton cycles, validating the predicted performance and confirming the viability of the full design-to-fabrication pipeline. Finally, thermohydraulic testing across the above studies served as a crucial experimental validation of advanced heat exchanger. Under consistent high-temperature and high-pressure conditions using supercritical CO 2 , the testing demonstrated that both TO and Gyroid-based TPMS designs significantly outperformed conventional straight-channel HXs. The TO design achieved a 115% increase in UA and NTU and a 27.6% boost in gravimetric power density, while the data-driven optimized Gyroid design delivered a 166% increase in UA and NTU and improved effectiveness from 68.7% to 86.1%. These results validate the simulation models, confirm the manufacturability of complex geometries under AM constraints, and provide key insights into design-performance trade-offs, thereby advancing the development of high-efficiency, compact heat exchangers for extreme environments.

36 MATERIALS SCIENCE↗

DOE OSTI2024

High-Fidelity and High-Performance Computational Simulations for Rapid Design Optimization of Sulfur Thermal Energy Storage (CRADA Report)

NREL and Element 16 collaborated on sulfur thermal energy storage modeling using NREL’s high performance computing (HPC) resources to assist its application in industrial processes. Industrial process heat (IPH) accounts for ~70% of US manufacturing energy use and is primarily produced by fossil fuel combustion. Approximately, 1500 TWht (~60% Terawatt hour thermal) of IPH demand is in the temperature range of 100-300°C. Industrial applications in this temperature range include drying, hydrothermal processing, thermal enhanced oil recovery, food and beverage, bioethanol production, etc. Cost-effective thermal energy storage (TES) that increases the utilization of waste and renewable heat (solar, geothermal, etc.) could provide significant energy savings and reliable heat sources, decrease emissions, and increase US manufacturing competitiveness through reductions in fuel consumption. This HPC4EI project facilitated Element 16’s development of low-cost and high-impact molten sulfur TES for dispatchable IPH. The development of a high-fidelity model validated by experimental data and HPC simulations enabled the successful resolution of the complex interplay between fluid dynamics and heat transfer processes during transient operation of sulfur TES, overcoming the numerical challenges posed by the non-linear temperature-dependent physical properties of sulfur. The project helped accelerate Element 16’s molten sulfur TES product design and support its broad applications.

25 ENERGY STORAGE↗

DOE OSTI2021

Clogging, dynamics, and reentrant fluid for active matter on periodic substrates

Here, we examine the collective states of run-and-tumble active matter disks driven over a periodic obstacle array. When the drive is applied along a symmetry direction of the array, we find a clog-free uniform liquid state for low activity, while at higher activity, the density becomes increasingly heterogeneous and an active clogged state emerges in which the mobility is strongly reduced. For driving along nonsymmetry or incommensurate directions, there are two different clogging behaviors consisting of a drive-dependent clogged state in the low activity thermal limit and a drive-independent clogged state at high activity. These regimes are separated by a uniform flowing liquid at intermediate activity. There is a critical activity level above which the thermal clogged state does not occur, as well as an optimal activity level that maximizes the disk mobility. Thermal clogged states are dependent on the driving direction while active clogged states are not. In the low activity regime, diluting the obstacles produces a monotonic increase in the mobility; however, for large activities, the mobility is more robust against obstacle dilution. We also examine the velocity-force curves for driving along nonsymmetry directions and find that they are linear when the activity is low or intermediate but become nonlinear at high activity and show behavior similar to that found for the plastic depinning of solids. At higher drives, the active clustering is lost. For low activity, we also find a reentrant fluid phase, where the system transitions from a high mobility fluid at low drives to a clogged state at higher drives and then back into another fluid phase at very high drives. We map the regions in which the thermally clogged, partially clogged, active uniform fluid, clustered fluid, active clogged, and directionally locked states occur as a function of disk density, drift force, and activity.

36 MATERIALS SCIENCE↗

NASA NTRS1977

Convection phenomena of importance for materials processing in space

The basic aspects of convection processes are delineated. It is shown that even in weak gravitational fields buoyancy can induce fluid motion. Furthermore, at reduced gravity, other nongravity forces such as surface or interfacial tensions, g jitter, thermal volume expansions, density differences due to phase changes, and magnetic and electric fields can induce fluid motions. The types of flows possible with these various driving forces are described, and criteria for determining the extent and nature of the resulting flows and heat transfer are presented. The different physical mechanisms that can occur separately and in combination are indicated, and the present state of knowledge of each of the phenomena is outlined. Specific research problems are described for many of the types of convection that are necessary to obtain greater understanding of their implications for space processing.

Ostrach, S.↗

NASA NTRS2013

Long Hole Film Cooling Dataset for CFD Development : Infrared Thermography and Thermocouple Surveys - Part 1

An experiment investigating flow and heat transfer of long (length to diameter ratio of 18) cylindrical film cooling holes has been completed. In this paper, the thermal field in the flow and on the surface of the film cooled flat plate is presented for nominal freestream turbulence intensities of 1.5 and 8 percent. The holes are inclined at 30deg above the downstream direction, injecting chilled air of density ratio 1.0 onto the surface of a flat plate. The diameter of the hole is 0.75 in. (0.01905 m) with center to center spacing (pitch) of 3 hole diameters. Coolant was injected into the mainstream flow at nominal blowing ratios of 0.5, 1.0, 1.5, and 2.0. The Reynolds number of the freestream was approximately 11,000 based on hole diameter. Thermocouple surveys were used to characterize the thermal field. Infrared thermography was used to determine the adiabatic film effectiveness on the plate. Hotwire anemometry was used to provide flowfield physics and turbulence measurements. The results are compared to existing data in the literature. The aim of this work is to produce a benchmark dataset for Computational Fluid Dynamics (CFD) development to eliminate the effects of hole length to diameter ratio and to improve resolution in the near-hole region. In this report, a Time-Filtered Navier Stokes (TFNS), also known as Partially Resolved Navier Stokes (PRNS), method that was implemented in the Glenn-HT code is used to model coolant-mainstream interaction. This method is a high fidelity unsteady method that aims to represent large scale flow features and mixing more accurately.

Shyam, Vikram↗

NASA NTRS2014

Long Hole Film Cooling Dataset for CFD Development - Flow and Film Effectiveness

An experiment investigating flow and heat transfer of long (length to diameter ratio of 18) cylindrical film cooling holes has been completed. In this paper, the thermal field in the flow and on the surface of the film cooled flat plate is presented for nominal freestream turbulence intensities of 1.5 and 8 percent. The holes are inclined at 30 deg above the downstream direction, injecting chilled air of density ratio 1.0 onto the surface of a flat plate. The diameter of the hole is 0.75 in. (approx. 0.02 m) with center to center spacing (pitch) of 3 hole diameters. Coolant was injected into the mainstream flow at nominal blowing ratios of 0.5, 1.0, 1.5, and 2.0. The Reynolds number of the freestream was approximately 11,000 based on hole diameter. Thermocouple surveys were used to characterize the thermal field. Infrared thermography was used to determine the adiabatic film effectiveness on the plate. Hotwire anemometry was used to provide flowfield physics and turbulence measurements. The results are compared to existing data in the literature. The aim of this work is to produce a benchmark dataset for Computational Fluid Dynamics (CFD) development to eliminate the effects of hole length to diameter ratio and to improve resolution in the near-hole region. In this report, a Time Filtered Navier Stokes (TFNS), also known as Partially Resolved Navier Stokes (PRNS), method that was implemented in the Glenn-HT code is used to model coolant-mainstream interaction. This method is a high fidelity unsteady method that aims to represent large scale flow features and mixing more accurately.

Shyam, Vikram↗

DOE OSTI2024

Generative Physics-Informed Neural Network Solving Multi-Scale and Multi-Phase Plasma Chemical Flow Field

Low-temperature plasmas (LTPs) are non-equilibrium systems with near-room-temperature gas and highly energetic electrons. This makes them ideal for delicate applications in biomedicine and semiconductor manufacturing, enabling processes like wound healing, sterilization, etching, and plasma-enhanced chemical vapor deposition without thermal damage. However, LTPs involve complex chemistries, with hundreds of species and thousands of reactions, complicating their diagnosis, prediction, and control. Conventional diagnostics, such as Fourier-transform infrared spectroscopy (FTIR), laser-induced fluorescence (LIF), and optical emission spectroscopy (OES), offer limited species detection, while mass spectrometry (MS) struggles with low-sensitivity species. Additionally, LTP simulations face multi-scale challenges, as macroscopic fluid dynamics and microscopic particle collisions operate on vastly different timescales. To address these issues, we developed an artificial intelligence (AI) based diagnostic system: a generative physics-informed neural network (PINN-Gen) that can predict spatially resolved species concentrations and temperatures in LTPs by integrating experimental data from planar LIF with microscopic plasma chemical kinetics and macroscopic fluid mechanics, including plasma-liquid interactions at the interface between two phases. PINN-Gen solves no equations but checks the errors of physical laws by substituting the output from neural network, and the comparison with the experimental results. Thus, it naturally avoids the multi-scale difficulty of numerical simulations and predicts the results of conventionally unsolvable multi-scale and multi-phase problems. The real-time prediction will be robust due to the physical information used in the training of such a neural network, and only very limited input of condition required due to its generative feature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Release on the Virtual Test Bed of an MSRE thermal hydraulics model

The nuclear industry is taking leaps in innovations with companies seeking a sustainable energy future through advanced nuclear reactors. The \gls{DOE}’s \gls{neams} program seeks to substantiate and bolster the deployment of advanced reactors through flexible multifidelity, multiphysics simulations of advanced nuclear reactors. Applications like SAM for one-dimensional systems thermal-hydraulics, and Pronghorn for multidimensional coarse mesh thermal-hydraulics, are geared to support innovations in industry by facilitating design, optimization, and licensing of advanced nuclear reactors. Coupling systems thermal-hydraulics and computational fluid dynamics codes can be difficult as the pressure coupling converges slowly; however, it is important to obtain the desired accuracy in each part of the primary loop. The authors of this model created an \gls{odc}~\cite{Mau23} approach to coupling SAM and Pronghorn. Leveraging this coupling technique, a \gls{msre} model was developed and released to the \gls{neams}/\gls{nric} \gls{vtb}. The \gls{msre} was chosen to be modeled because of the wealth of experimental data available and because of the strong physics coupling between the core and primary circuit \cite{doi:10.13182/NT8-2-118}. This document contextualizes the history of the \gls{msre}, describes the thermal hydraulics models used, and detail the implementation of multidimensional thermal-hydraulics and system codes based on the \gls{odc} method~\cite{Penn} for the \gls{msre} model. Finally, this document presents how other modelers could apply the SAM-Pronghorn \gls{odc} for other advanced reactor models. Current item is the set of slides for ANS Winter 23. The release of the model and the ANS summary have already been approved

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

Release on the Virtual Test Bed of an MSRE thermal hydraulics model

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1999

Development of Optical Diagnostic Techniques for Microgravity Materials Processing

Materials processing including crystal growth, either under a gravity environment on ground or a microgravity environment in space, involves complicated phenomena of fluid motions in gas or liquid phases as well as interaction of various species. To obtain important physical insight, it is very necessary to provide gross-field optical diagnostics for monitoring various physical properties. Materials processing inhibits easy access by ordinary instruments and thus characterizing gross-field physical properties is very challenging. Typical properties of importance can be fluid velocity, temperature, and species concentration for fluids, and surface topology and defects for solids. Observing surface grow rate during crystal growth is also important. Material microstructures, i.e., integrity of crystal structures, is strongly influenced by the existence of thermally-induced flow as well as local nucleation of particles during solidification, which may act in many detrimental ways. In both ground-based and microgravity experiments, the nature of product property changes resulting from three-dimensional fluid or particle motions need be characterized. Gross-field diagnostics is thus required to identify their effects on product defects and process deficiencies. The quantitative visualization techniques can also be used for validation of numerical modeling. For optical nonintrusive gross-field diagnostic techniques, two approaches were developed as summer projects. One optical approach allows us to provide information of species concentration and temperature for monitoring in real time. The other approach, that is, the concept which is formulated for detection of surface topography measurement can provide unprecedented spatial resolution during crystal growth.

Cha, Soyoung Stephen↗

NASA NTRS2025

Physics-Based Modeling and Simulation of Self-Reacting Friction Stir Welding Using Computational Fluid Dynamics

A physics-based model was developed to simulate the behavior of material in a self-reacting friction stir welding (SR-FSW) process for the joining of metals. This steady-state model builds upon fundamental computational fluid dynamic (CFD) principles within Ansys Fluent to solve the discretized equations. The effective viscosity is calculated using a viscoplastic model using a Sheppard-Wright formulation of flow stress. Numerous advancements have been made in the incorporated physics including (1) temperature-dependent material properties; (2) locally adaptable flow and thermal boundary conditions; and (3) adapting material properties in nugget in response to microstructural changes. Simulation strategies to accelerate computation and improve numerical stability include adapting the mesh refinement and solver relaxation factors during simulation. The result is a highly robust and computationally efficient model capable of providing the material flow and temperature history across the domain. As material history determines the local microstructure and ultimately weld strength, an accurate and detailed physics-based model has the potential to accelerate SR-FSW process development. The model is highly adaptable to changes in process parameters, tool design, or alloy.

Process Modeling↗

DOE OSTI2022

Sensitivity study of coupled chemical-CFD simulations for analyzing aluminum-clad spent nuclear fuel storage in sealed canisters

We report the United States Department of Energy (DOE) manages over 50 Metric Tons Heavy Metal (MTHM) of aluminum-clad spent nuclear fuel. One main source for DOE’s Aluminum-clad spent nuclear fuel (ASNF) inventory is the advanced test reactor (ATR) at the INL site, which makes this fuel of particular interest for storage scenarios. Road-ready and final disposition packaging configurations for the ATR fuel dictates storage within helium-backfilled, sealed DOE standard canisters. The conditions within these sealed canisters for extended (greater than50 year dry) storage is of interest. To further this goal, a three-dimensional (3D) multi-physics computational fluid dynamics (CFD) model is developed of the sealed DOE standard canisters. This 3D CFD model is one-way coupled with bulk gas radiolysis reactions considering sealed canisters with inert gas and possible trace amounts of air and water vapor. This study looks at the evolution of the thermal history of the canisters over a 50 year time period with a coupling to the chemical reactions occurring from radiolytic breakdown of residual water. A sensitivity study is then carried out over the parameters of the model including the fuel decay heat, residual water content, sealed pressure, canister external temperature, and canister emissivity. In pure helium, hydrogen generation rates are low, under 10 ppm, but hydrogen generation rates are affected greatly by the presence of even 1% residual air, increasing by 50-plus-fold, and nitric acid generation with residual air also occurs ranging from 500 to 4000 ppm after 50 years. The fuel decay heat and the residual water content show the most importance in the generation of hydrogen gas in pure air, and for nitric acid with a residual air condition. External temperature, canister emissivity and sealed pressure all show minor sensitivity effects to the generation of potentially harmful species.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2023

MAS: A versatile Landau-fluid eigenvalue code for plasma stability analysis in general geometry

We have developed a new global eigenvalue code, Multiscale Analysis for plasma Stabilities (MAS), for studying plasma problems with wave toroidal mode number (n) and frequency (ω) in a broad range of interest in general tokamak geometry, based on a five-field Landau-fluid description of thermal plasmas. Beyond keeping the necessary plasma fluid response, we further retain the important kinetic effects including diamagnetic drift, ion finite Larmor radius, finite parallel electric field (E||), ion and electron Landau resonances in a self-consistent and non-perturbative manner without sacrificing the attractive efficiency in computation. The physical capabilities of the code are evaluated and examined in the aspects of both theory and simulation. In theory, the comprehensive Landau-fluid model implemented in MAS can be reduced to the well-known ideal MHD model, electrostatic ion-fluid model, and drift-kinetic model in various limits, which clearly delineates the physics validity regime. In simulation, MAS has been well benchmarked with theory and other gyrokinetic and kinetic-MHD hybrid codes in a manner of adopting the unified physical and numerical framework, which covers the kinetic Alfv\'en wave (KAW), ion sound wave (ISW), low-n kink, high-n ion temperature gradient mode (ITG) and kinetic ballooning mode (KBM). Moreover, MAS is successfully applied to model the Alfv\'en eigenmode (AE) activities in DIII-D discharge #159243, which faithfully captures the frequency sweeping of reversed shear Alfv\'en eigenmode (RSAE), the tunneling damping of toroidal Alfv\'en eigenmode (TAE), as well as the polarization characteristics of kinetic beta-induced Alfv\'en eigenmode (KBAE) and beta-induced Alfv\'en-acoustic eigenmode (BAAE) being consistent with former gyrokinetic theory and simulation. With respect to the key progress contributed to the community, MAS has the advantage of combining rich physics ingredients, realistic global geometry and high computation efficiency together for plasma stability analysis in linear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

A Study on Computational Cost Reduction of Simulations of Phase-Change Material (PCM) Embedded Heat Exchangers

Thermal storage can be implemented using Phase-Change Materials (PCM), which absorb significant latent heat with a relatively small temperature change. PCM phase-change processes are transient and are driven by thermal diffusion and natural convection – the latter, especially for melting process. Modeling and simulation of PCM heat exchangers (HX’s) is typically computationally intensive due to the relatively complex time-dependent physics. Most of the PCM modeling work in the literature uses high-order modeling tools such as Computational Fluid Dynamics (CFD) and Lattice-Boltzmann Method (LBM). For design purposes, the existing PCM modeling approaches are not practical, limiting researchers in their ability to investigate new ideas and different PCM’s with faster turnarounds. This paper presents a study investigating the reduction of computational cost of PCM embedded HX’s CFD models by evaluating the feasibility of spatial reduction without losing accuracy. The analysis consists of comparing full and partial domain under full melting conditions. The subject of this study is a single straight tube with circular transverse fins in the vertical orientation, using PCM’s with 35oC nominal melting temperature. Different tube and fin dimensions are investigated. Results indicate that the reduced domain reproduces - in half the run time - the same behavior as the full domain since the buoyancy effects are localized and patterned. The outputs from the partial domain simulation were used to build a non-general correlation for the PCM heat transfer characteristics and demonstrated how it can be implemented in a Finite Control Volume Reduced Order Model (ROM). The ROM can accurately reproduce the CFD simulations at 4 to 5 orders of magnitude faster.

25 ENERGY STORAGE↗

DOE OSTI2023

Release on the Virtual Test Bed of a Molten Salt Reactor Experiment SAM-Pronghorn Coupled Model using the Domain Overlapping Approach

The nuclear industry is taking leaps in innovations with companies seeking a sustainable energy future through advanced nuclear reactors. The Department of Energy (DOE)’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program seeks to substantiate and bolster the deployment of advanced reactors through flexible multifidelity, multiphysics simulations of advanced nuclear reactors. Applications like SAM for one-dimensional systems thermalhydraulics, and Pronghorn for multidimensional coarse mesh thermal-hydraulics, are geared to support innovations in industry by facilitating design, optimization, and licensing of advanced nuclear reactors. Coupling systems thermal-hydraulics and computational fluid dynamics codes can be difficult as the pressure coupling converges slowly; however, it is important to obtain the desired accuracy in each part of the primary loop. The authors of this model created an Overlapping-Domain Coupling (ODC) approach to coupling SAM and Pronghorn. Leveraging this coupling technique, a Molten Salt Reactor Experiment (MSRE) model was developed and released to the NEAMS/National Reactor Innovation Center (NRIC) Virtual Test Bed (VTB). The MSRE was chosen to be modeled because of the wealth of experimental data available and because of the strong physics coupling between the core and primary circuit. This paper contextualizes the history of the MSRE, describes the thermal hydraulics models used, and detail the implementation of multidimensional thermal-hydraulics and system codes based on the ODC method for the MSRE model. Finally, this paper presents how other modelers could apply the SAM Pronghorn ODC for other advanced reactor models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

Accelerating transients with NekRS: GPU overlapping domain implementation and multi-rate timestepping

The simulation of nuclear transients using Computational Fluid Dynamics (CFD) presents significant computational challenges due to the inherent complexity and the wide separation in temporal scales between various flow physical phenomena. These disparities lead to high computational costs, often making the simulation of transients impractical without advanced techniques. Consequently, multiple research initiatives are being pursued by the NEAMS thermal-hydraulic area, some driven by academic institutions and some by national laboratories. Overall, they are exploring novel methods to make transient simulations more feasible and efficient. This report delves into recent advancements within the CFD code NekRS, specifically those achieved in Fiscal Year 2024 under the CONNECT effort, aimed at improving the performance and feasibility of transient simulations. The first major advancement involves the porting of NekRS to Aurora, one of the Department of Energy’s (DOE) most powerful supercomputers. Additionally, the report discusses the implementation of an overlapping domain capability within NekRS. This novel GPU-accelerated capability allows different spatial regions of the domain to be solved independently, enhancing the code’s efficiency, particularly when running large-scale simulations in complex domains. The scalability of this approach is demonstrated, highlighting its potential to transform how transients are approached in CFD simulations. Lastly, the report focuses on how this overlapping domain capability specifically accelerates transient simulations through multi-rate timestepping. By decoupling different regions and facilitating faster computations, this method offers a promising pathway to making nuclear transient simulations more computationally feasible, addressing one of the critical bottlenecks in the field. Together, these advancements represent a significant leap forward in transient simulation technology, bringing closer the possibility of handling highly complex nuclear scenarios with greater efficiency and accuracy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗