Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermal fluid simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 19 records

Development of Transient Coupled Chemical-Thermal-Fluids Simulation for Unsealed Aluminum-Clad Spent Nuclear Fuel Storage Canister

This study examines a preliminary three-dimensional (3D) multi-physics computational fluid dynamics (CFD) model of the unsealed canisters containing aluminum-clad spent nuclear fuel from INL’s ATR in the irradiated fuel storage facility. This 3D CFD model considers fully resolved fuel plates in an unsealed canister with small mass exchange with surrounding air as time evolves in the simulation and contains a network of chemical reactions to consider due to several of the side reactions occurring in the presence of water vapor, nitrogen, oxygen and carbon dioxide rather than in helium backfilled DOE standard canisters. Prior to the CFD modeling, the full set of chemical reactions under these conditions was analyzed, reduced, and refit to a more manageable number of equations for the computational efficiency of the 3D CFD simulation. To deal with the long time-scale of interest for tracking the possible buildup of trace species such as nitric acid, NOx, and hydrogen, a segregated time stepping algorithm is used, wherein the velocity field is only updated periodically, and assume to be in a pseudo steady state. This study first examines the theory of the model, workflow to establish radiolysis reaction network, and initial simulations of the evolutions of thermal fields, hydrogen gas and nitric acid concentrations within unsealed canisters in conditions relevant to the INL INTEC CPP-603 facility over a long period of time.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Development of Transient Coupled Chemical-Thermal-Fluids Simulation for Unsealed Aluminum-Clad Spent Nuclear Fuel Storage Canister

This study examines a preliminary three-dimensional (3D) multi-physics computational fluid dynamics (CFD) model of the unsealed canisters containing aluminum-clad spent nuclear fuel from INL’s ATR in the irradiated fuel storage facility. This 3D CFD model considers fully resolved fuel plates in an unsealed canister with small mass exchange with surrounding air as time evolves in the simulation and contains a network of chemical reactions to consider due to several of the side reactions occurring in the presence of water vapor, nitrogen, oxygen and carbon dioxide rather than in helium backfilled DOE standard canisters. Prior to the CFD modeling, the full set of chemical reactions under these conditions was analyzed, reduced, and refit to a more manageable number of equations for the computational efficiency of the 3D CFD simulation. To deal with the long time-scale of interest for tracking the possible buildup of trace species such as nitric acid, NOx, and hydrogen, a segregated time stepping algorithm is used, wherein the velocity field is only updated periodically, and assume to be in a pseudo steady state. This study first examines the theory of the model, workflow to establish radiolysis reaction network, and initial simulations of the evolutions of thermal fields, hydrogen gas and nitric acid concentrations within unsealed canisters in conditions relevant to the INL INTEC CPP-603 facility over a long period of time.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Development of Transient Coupled Chemical-Thermal-Fluids Simulation for Unsealed Aluminum-Clad Spent Nuclear Fuel Storage Canister - 19273

This study examines a preliminary three-dimensional (3D) multi-physics computational fluid dynamics (CFD) model of the unsealed canisters containing aluminum-clad spent nuclear fuel from INL's ATR in the irradiated fuel storage facility. This 3D CFD model considers fully resolved fuel plates in an unsealed canister with small mass exchange with surrounding air as time evolves in the simulation and contains a network of chemical reactions to consider due to several of the side reactions occurring in the presence of water vapor, nitrogen, oxygen and carbon dioxide rather than in helium backfilled DOE standard canisters. Prior to the CFD modeling, the full set of chemical reactions under these conditions was analyzed, reduced, and refit to a more manageable number of equations for the computational efficiency of the 3D CFD simulation. To deal with the long time-scale of interest for tracking the possible buildup of trace species such as nitric acid, NOx, and hydrogen, a segregated time stepping algorithm is used, wherein the velocity field is only updated periodically, and assumed to be in a pseudo steady state. This study first examines the theory of the model, workflow to establish radiolysis reaction network, and initial simulations of the evolutions of thermal fields, hydrogen gas and nitric acid concentrations within unsealed canisters in conditions relevant to the INL INTEC CPP-603 facility over a long period of time. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Development of Transient Coupled Chemical-Thermal-Fluid Multiphysics Simulation for Unsealed, Vented Aluminum-clad Spent Nuclear Fuel Storage Canister

As the second step towards modeling three-dimensional (3D) multi-physics computational fluid dynamics (CFD) of the unsealed containers in the IFSF facility, the previous methodology used for developing the sealed DOE standard canister model was utilized here. This 3D CFD model considers unsealed canister with small mass exchange with surrounding air as time evolves in the simulation and contains a larger network of chemical reactions to consider due to several of the side reactions occurring in the presence of water vapor, nitrogen, oxygen and carbon dioxide rather than a helium backfilled canister. Similar to the sealed system, the set of full radiolytic chemical reactions of bulk moist air were established, and then reduced to a more manageable set of equations for the computational efficiency of the 3D CFD simulation. Due to the similar size of the unsealed canisters to the sealed DOE canisters, it is assumed that the chemical and thermal response times to changes in the conditions remains similar such that the same time-stepping algorithm can be utilized. This report documents the theory of the model, workflow to establish radiolysis reaction network, and initial simulations of the evolutions of thermal fields, hydrogen gas and nitric acid concentrations within unsealed canisters in conditions relevant to INL INTEC CPP-603 over a long period of time.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Transient Coupled Chemical-Thermal-Fluid Field Simulation for Sealed Aluminum-clad Spent Nuclear Fuel Storage Canister

As the first step toward developing three-dimensional (3D) multi-physics computational fluid dynamics (CFD) model for unsealed and vented canister storage system, a 3D CFD model coupled with bulk gas radiolysis reactions was developed first for sealed DOE standard canisters filled with inert gas and trace amount of air and water. The workflow for constructing canister-scale 3D CFD models and coupling with gas phase radiolysis reactions were established, which can be readily extended to unsealed, vented canister storage system. This interim milestone report documents the theory of the model, workflow to establish radiolysis reaction network, and initial simulations of the evolutions of thermal fields and hydrogen gas concentrations within sealed DOE standard canisters over long period of time. In addition, a mesh refinement test was done to show that increasing the models mesh refinement had negligible impact upon the temperature profiles.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Modeling and Simulations of Thermal-Fluid Phenomena Related to Pressurized Conduction Cooldown in High Temperature Gas Reactor

Under the support of DOE-NE’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program, an effort is being pursued to support High Temperature Gas Reactors (HTGR) technology development and its modeling and simulation needs. There is a particular need for advanced modeling & simulation tools to predict thermal-fluid behavior in the nuclear reactor primary system—especially the core and lower and upper plena—during safety-related transients.To help accomplish this, we must provide confidence in the reliability and accuracy of the thermal fluid simulation tools being developed under the NEAMS program for HTGR applications. The benchmark simulations of the available experimental tests would help us in this effort.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Thermal Management System of an Outer-Rotor-Motor-Based Traction Drive With Integrated Power Electronics in its Central Cavity

Increasing power densities of electric vehicle traction drive systems necessitates combining the electric motor and the power electronics into one unit. A compact, integrated traction drive unit with fewer components also drives production costs down, enabling wider adoption of electric vehicles. However, the integration of power electronics in the electric machine is associated with challenges of designing an effective thermal management solution for the combined traction drive system. This paper focuses on the thermal management approach selected for Oak Ridge National Laboratory’s outer-rotor-motor-based integrated traction drive and evaluates its potential performance. The outer-rotor-motor configuration provided an opportunity for integration of the six-phase inverter in the available space in the central cavity of the internal stator. A cylindrical inverter enclosure with integrated coolant (water-ethylene glycol) channels in its walls was designed to enable heat removal from the power electronics. Numerical thermal-fluid modeling and initial channel/fin optimization results for the cylindrical heat sink are presented here. As permanent magnets are integrated in a high-speed (20 000 RPM) outer rotor, forced air convection provides cooling for the magnets. The magnets were segmented axially to mitigate the eddy current losses. Heat generated in the stator windings and laminations is removed by water-ethylene glycol coolant circulating in interpolar T-shape ceramic heat exchangers inserted between windings. General design concepts and numerical thermal-fluid simulations illustrating the electric motor thermal management solution are also presented.

Kekelia, Bidzina↗

Classification of machine learning frameworks for data-driven thermal fluid models

Thermal fluid processes are inherently multi-physics and multi-scale, involving mass-momentum-energy transport phenomena at multiple scales. Thermal fluid simulation (TFS) is based on solving conservative equations, for which – except for “first-principles” direct numerical simulation – closure relations (CRs) are required to provide microscopic interactions or so-called sub-grid-scale physics. In practice, TFS is realized through reduced-order modeling, and its CRs as low-fidelity models can be informed by observations and data from relevant and adequately evaluated experiments and high-fidelity simulations. This paper is focused on data-driven TFS models, specifically on their development using machine learning (ML). Five ML frameworks are introduced including physics-separated ML (PSML or Type I ML), physics-evaluated ML (PEML or Type II ML), physics-integrated ML (PIML or Type III ML), physics-recovered (PRML or Type IV ML), and physics-discovered ML (PDML or Type V ML). The frameworks vary in their performance for different applications depending on the level of knowledge of governing physics, source, type, amount and quality of available data for training. Notably, outlined for the first time in this paper, Type III models present stringent requirements on modeling, substantial computing resources for training, and high potential in extracting value from “big data” in thermal fluid research. Here the current paper demonstrates and investigates ML frameworks in three examples. First, we utilize the heat diffusion equation with a nonlinear conductivity model formulated by convolutional neural networks (CNNs) and feedforward neural networks (FNNs) to illustrate the applications of Type I, Type II, Type III, and Type V ML. The results indicate a preference for Type II ML under deficient data support. Type III ML can effectively utilize field data, potentially generating more robust predictions than Type I and Type II ML. CNN-based closures exhibit more predictability than FNN-based closures, but CNN-based closures require more training data to obtain accurate predictions. Second, we illustrate how to employ Type I ML and Type II ML frameworks for data-driven turbulence modeling using reference works. Third, we demonstrate Type I ML by building a deep FNN-based slip closure for two-phase flow modeling. The results show that deep FNN-based closures exhibit a bounded error in the prediction domain.

42 ENGINEERING↗

Modeling of a Generic Pebble Bed High-temperature Gas-cooled Reactor (PB-HTGR) with SAM

This report presents the modeling of the core of a generic pebble-bed reactor (PBR) at the system level using the System Analysis Module (SAM) code. This work is an extension of a previous work by the authors (Ooi et al. (2021)) that used the so-called 2-D ring model approach to model the PBMR-400. With the new approach, the pebble bed of the reactor is modeled with multiple PBCoreChannel components with spherical heat structures which allows the code to calculate thermal fluid parameters with built-in closure relations. The new core-channel approach is an improvement to the 2-D ring model approach as it does not introduce geometric distortions to the model and thus reduces the uncertainties of the predictions. In addition to thermal fluid simulations, point kinetics (PKE) are included to the model. Simulations are performed under a steady-state normal operation condition and a load-following transient scenario. This particular transient scenario is chosen as it tests both the thermal fluid and neutronics aspects of the model. The predicted results from both the steady-state and transient scenarios are compared with the results by Stew- art et al. (2021) who performed similar simulations with a Griffin-Pronghorn coupled tool. Despite the differences between the codes, with SAM being a system-analysis code and Pronghorn being a porous-medium code, both sets of results compare favorably. The over- all profiles and trends of the predicted temperatures and reactivities from the SAM and Griffin/Pronghorn simulations are similar, with some differences in their predicted values. The first part of the report covers the significance of a relatively fast-running approach that is capable of modeling the pebble bed reactor at the system-level while simultaneously capturing the radial thermal behavior of the core. Then, the modeling approach used in this work is discussed in details. Lastly, the results and comparisons with the Griffin/Pronghorn simulation are presented.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Energy Exascale Computational Fluid Dynamics Simulations With the Spectral Element Method

Development and application of the open-source GPU-based fluid-thermal simulation code, NekRS, are described. Time advancement is based on an efficient kth-order accurate timesplit formulation coupled with scalable iterative solvers. Spatial discretization is based on the high-order spectral element method (SEM), which affords the use of fast, low-memory, matrix-free operator evaluation. Further, recent developments include support for nonconforming meshes using overset grids and for GPU-based Lagrangian particle tracking. Results of large-eddy simulations of atmospheric boundary layers for wind-energy applications as well as extensive nuclear energy applications are presented.

42 ENGINEERING↗

NekRS, a GPU-accelerated spectral element Navier–Stokes solver

The development of NekRS, a GPU-oriented thermal-fluids simulation code based on the spectral element method (SEM) is described. For performance portability, the code is based on the open concurrent compute abstraction and leverages scalable developments in the SEM code Nek5000 and in libParanumal, which is a library of high-performance kernels for high-order discretizations and PDE-based miniapps. Critical performance sections of the Navier–Stokes time advancement are addressed. Performance results on several platforms are presented here, including scaling to 27,648 V100s on OLCF Summit, for calculations of up to 60B gridpoints.

97 MATHEMATICS AND COMPUTING↗

Thermal Hydraulic Computational Fluid Dynamics Simulations and Experimental Investigation of Deformed Fuel Assemblies

The project, Toward a Longer Life Core: Thermal Hydraulic CFD Simulations and Experimental Investigation of Deformed Fuel Assemblies, DOE Project code DE-NE0008321, was a verification and validation project for flow and heat transfer through wire wrapped simulated liquid metal fuel assemblies that included both experiments and computational fluid dynamics simulations of those experiments. This project was a two year collaboration between AREVA, TerraPower, Argonne National Laboratory and Texas A&M University. Experiments were performed by AREVA and Texas A&M University. Numerical simulations of these experiments were performed by TerraPower and Argonne National Lab. Project management was performed by AREVA Federal Services. The first of a kind project resulted in the production of both local point temperature measurements and local flow mixing experiment data paired with numerical simulation benchmarking of the experiments. The project experiments included the largest wire-wrapped pin assembly Mass Index of Refraction (MIR) experiment in the world, the first known wire-wrapped assembly experiment with deformed duct geometries and the largest numerical simulations ever produced for wire-wrapped bundles.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Plant-wide simulation model for transient studies on the XE-100 - 18464

This paper gives an overview of the plant-wide simulation model that is used by X-energy LLC for developing their plant design and control philosophy, as well as assessing plant upset conditions. The simulation model was created in Flownex{sup R} thermal-fluid system simulation software, which solves mass, momentum and energy equations for compressible, incompressible and homogeneous two-phase flow. The plant model encompasses all major components, from reactor to condenser, including main control loop elements. The paper pays particular attention to the validation of the two-phase flow framework used for the secondary steam cycle. Two steady-state validation cases are done using operational data from the THTR. Validation of transient operation is done by simulating density wave oscillations, a type of system-level boiling phenomenon. This is regarded as a thorough test of transient capability because of the challenges this phenomenon poses to several aspects of the calculation and to the consistency of the code as a whole. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Compact Diffusion Bonded Heat Exchanger Fatigue Life Simulations (Final Report)

Compact diffusion bonded heat exchangers are essential for high pressure heat exchange, but they are subject to thermal fatigue and ramp rate limitations. Simulation of these geometries is challenging with a large range of length and time scales from thousands of mm-sized microchannels inside a m-sized heat exchanger. Multi-physics simulations including thermal, fluid, and solid mechanics components are being used to predict stress within the heat exchangers under these conditions. These predictions can then be used to understand thermal ramp rate limitations while keeping maximum stresses low as well as fatigue life predictions from well-known empirical models.

42 ENGINEERING↗

Simulations of multicomponent relativistic thermalization

Multicomponent relativistic fluids have been studied for decades. However, simulating the dynamics of the particles and fluids in such a mixture has been a challenge due to the fact that such simulations are computationally expensive in three spatial dimensions. Here, we report on the development and application of a multidimensional relativistic Monte Carlo code to explore the thermalization process in a relativistic multicomponent environment in a computationally inexpensive way. As an illustration we simulate the fully relativistic three-dimensional Brownian-motion-like solution to the thermalization of a high-mass particle (proton) in a bath of relativistic low-mass particles (electrons). Here, we follow the thermalization and ultimate equilibrium distribution of the Brownian-like particle as can happen in the cosmic plasma during big-bang nucleosynthesis. We also simulate the thermalization of energetic particles injected into the plasma as can occur, for example, by the decay of massive unstable particles during the big bang.

79 ASTRONOMY AND ASTROPHYSICS↗

Thermal Stress Modeling and Analysis of Packed-bed Thermocline Energy Storage Tank for INL Thermal Energy Distribution System (TEDS)

The Thermal Energy Distribution System (TEDS) at Idaho National Laboratory (INL) is a thermal-hydraulic flow loop to support the integration of co-located multiple experimental systems, where a packed-bed thermal energy storage (TES) is installed as a thermal buffer and storage unit for TEDS. The packed-bed TES is adopted in TEDS because of its benefit as a low-cost single-tank storage option compared to the traditional two-tank storage. However, thermal ratcheting is one potential design concern which is caused by the rearrangement of granular filler inside a packed-bed tank during continuous thermal cycling operation of the packed-bed TES tank. If the thermally induced stress exceeds yield strength of the tank wall, it may cause catastrophic consequences like rupture of the thermal storage tank. Thus, it is crucial to understand the phenomenon to ensure the robust operation. Based on the temperature boundary conditions given by transient thermal analyses with computational fluid dynamics (CFD) simulations, the thermal ratcheting analysis is then conducted to evaluate the hoop stress and resultant thermal ratcheting potential of the TES tanks with two different modeling approaches: (1) infinite rigidity model and (2) Drucker-Prager (DP) model. The validity of each modeling method was examined by comparing the numerical simulation with the experimental data obtained from the packed-bed TES tank for Solar One Plant and evaluate the thermal ratcheting potential of the TEDS TES tank.

25 ENERGY STORAGE↗

Compositional simulation of hydrocarbon recovery from oil shale reservoirs with diverse initial saturations of fluid phases by various thermal processes

We have studied the hydrocarbon production from oil shale reservoirs filled with diverse initial saturations of fluid phases by implementing numerical simulations of various thermal in-situ upgrading processes. We use our in-house fully functional, fully implicit, and non-isothermal simulator, which describes the in-situ upgrading processes and hydrocarbon recovery by multiphase-multicomponent systems. We have conducted two sets of simulation cases—five-spot well pattern problems and Shell In-situ Conversion Process (ICP) problems. In the five-spot well pattern problems, we have analyzed the effects of initial fluid phase that fills the single-phase reservoir and thermal processes by four cases—electrical heating in the single-phase-aqueous reservoir, electrical heating in the single-phase-gaseous reservoir, hot water injection in the single-phase-aqueous reservoir, and hot CO 2 -injection in the single-phase-gaseous reservoir. In the ICP problems, we have analyzed the effects of initial saturations of fluid phases that fill two-phase-aqueous-and-gaseous reservoir by three cases—initial aqueous phase saturations of 0.16, 0.44, and 0.72. Through the simulation cases, system response and production behavior including temperature profile, kerogen fraction profile, evolution of effective porosity and absolute permeability, phase production, and product selectivity are analyzed. In the five-spot well pattern problems, it is found that the hot water injection in the aqueous phase reservoir shows the highest total hydrocarbon production, but also shows the highest water-oil-mass-ratio. Productions of phases and components show very different behavior in the cases of electrical heating in the aqueous phase reservoir and the gaseous phase reservoir. In the ICP problems, it is found that the speed of kerogen decomposition is almost identical in the cases, but the production behavior of phases and components is very different. It is found that more liquid organic phase has been produced in the case with the higher initial saturation of aqueous phase by the less production of gaseous phase.

04 OIL SHALES AND TAR SANDS↗