Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational fluid dynamics code”

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 505 records · Page 28

Airfoil Table Generation and Validation for the VR-12 and SSC-A09 Airfoils and Quadrotor Performance Prediction

As part of NASA’s Urban Air Mobility (UAM) mission, the Revolutionary Vertical Lift Technology (RVLT) project is creating rotorcraft reference designs that can be used by the rotorcraft community to develop quiet, efficient, and safer air vehicles. At NASA, this effort is made possible by the use of the RVLT toolchain software. Airfoil tables are often a critical part of the conceptual design process as they inform both comprehensive analysis and Computational Fluid Dynamics (CFD) codes of 2D aerodynamic coefficients, such as lift, drag, and moment coefficients. These airfoil tables, typically generated during experimental testing, are often proprietary and thus not widely distributable; further, there is not always clear or available reference documentation that provides information on the test conditions. The airfoil tables from these experiments are referred to in this work as Legacy airfoil tables. The scope of this study is to generate airfoil tables with the NASA OVERset grid CFD FLOW solver (OVERFLOW) coupled with the AirFoil Table Generator (AFTGen) software and validate these airfoil tables based on performance predictions calculated in Comprehensive Hierarchical Aeromechanics Rotorcraft Model (CHARM). The overall goal is to provide CFD generated airfoil tables, validate them through comprehensive analysis, and publish the airfoil tables for public distribution. The RVLT quadrotor reference model was selected as the basis for all analyses. The quadrotor uses the VR-12 and SSC-A09 airfoils for its rotor blades. Simulations were run in OVERFLOW for the VR-12 and SSC-A09 airfoils for ±20° angles of attack, Mach numbers from 0.3 to 1.0, and a Reynolds-Mach proportionality constant of 1.023x107. The C81 format airfoil tables from AFTGen were then blended with Legacy NACA 0012 airfoil table data for ±20° to ±180° angle of attack range. CHARM was used to compute thrust coefficient, power coefficient, and figure of merit for an isolated rotor in hover using proprietary legacy airfoil tables for the VR-12 and SSC-A09, as well as with the OVERFLOW generated airfoil tables for the VR-12 and SSC-A09. Overall, the OVERFLOW and Legacy performance predictions agreed well, with the best correlation between simulated and Legacy data observed with the VR-12 airfoil tables.

Airfoil↗

Measurement Requirements for Improved Modeling of Arcjet Facility Flows

Current efforts to develop new reusable launch vehicles and to pursue low-cost robotic planetary missions have led to a renewed interest in understanding arc-jet flows. Part of this renewed interest is concerned with improving the understanding of arc-jet test results and the potential use of available computational-fluid- dynamic (CFD) codes to aid in this effort. These CFD codes have been extensively developed and tested for application to nonequilibrium, hypersonic flow modeling. It is envisioned, perhaps naively, that the application of these CFD codes to the simulation of arc-jet flows would serve two purposes: first. the codes would help to characterize the nonequilibrium nature of the arc-jet flows; and second. arc-jet experiments could potentially be used to validate the flow models. These two objectives are, to some extent, mutually exclusive. However, the purpose of the present discussion is to address what role CFD codes can play in the current arc-jet flow characterization effort, and whether or not the simulation of arc-jet facility tests can be used to eva1uate some of the modeling that is used to formu1ate these codes. This presentation is organized into several sections. In the introductory section, the development of large-scale, constricted-arc test facilities within NASA is reviewed, and the current state of flow diagnostics using conventional instrumentation is summarized. The motivation for using CFD to simulate arc-jet flows is addressed in the next section, and the basic requirements for CFD models that would be used for these simulations are briefly discussed. This section is followed by a more detailed description of experimental measurements that are needed to initiate credible simulations and to evaluate their fidelity in the different flow regions of an arc-jet facility. Observations from a recent combined computational and experiment.al investigation of shock-layer flows in a large-scale arc-jet facility are then used to illustrate the current state of development of diagnostic instrumentation, CFD simulations, and general knowledge in the field of arc-jet characterization. Finally, the main points are summarized and recommendations for future efforts are given.

Fletcher, Douglas G.↗

kynema-fmb [SWR-23-07]

Kynema-FMB (FKA: Kynema) is an open-source performance portable flexible multibody (FMB) dynamics solver designed for time-domain simulations. While originally tailored for wind turbine structural dynamics, the formulation and implementation are those of a general flexible-multidbody dynamics solver that can readily be applied to a wide range of systems. Kynema was designed with a narrow focus, namely to provide a lightweight, fast, accurate FMD solver for coupling to computational-fluid-dynamics (CFD) codes, especially the CFD codes in the Kynema suite, for fluid-structure-interaction (FSI) simulations. Kynema-FMB is equipped to model systems that can be represented as a collection of beams and rigid bodies that are connected through constraints. Degrees of freedom are defined in the inertial/global frame of reference and include displacements and rotations (formally as rotation matrices, but stored as quaternions). The underlying formulation is built on a Lie-group time integrator designed for index-3 differential-algebraic equations, which is second-order accurate in time (Bruls et al., 2012). Beam models are based on geometrically exact beam theory and are discretized as high-order spectral finite elements similar to those in BeamDyn (Wang et al., 2017). The governing equations for a FMD system like a wind turbine constitute a highly nonlinear system of constrained partial-differential equations. Kynema-FMB uses analytical Jacobians in the nonlinear-system solves in each time step. Linear systems use sparse storage and several third-party sparse-linear-system solvers are enabled. Ill conditioning of linear systems is mitigated with preconditioning described in Bottasso et al, 2008. Kynema-FMB is integrated with a simple open-source controller (ROSCO). There is an application programming interface (API) for coupling to geometry-resolved CFD (like that in Sharma et al., 2023) and actuator-force CFD (like that in Kuhn et al., 2025). In the latter, for actuator-line models, Kynema-FMB includes an internal blade-element solver that depends on user-provided lookup tables for coefficients of lift and drag, i.e., aerodynamic polars. Kynema-FMB is written in C++ and leverages Kokkos and Kokkos-Kernels (KokkosEcosystem) as its performance portability layer enabling simulations on both CPU and GPU systems. The repository is equipped with extensive automated testing at the unit and regression/system levels. The following describes the high-level development objectives conceived for Kynema: *Kynema will follow modern software development best practices, including test-driven development (TDD), version control, hierarchical automated testing, and continuous integration (CI) for a robust development environment. *The core data structures are memory efficient and enable vectorization and parallelization at multiple levels. *Data structures are data-oriented to exploit methods for accelerated computing including high utilization of chip resources (e.g., single instruction multiple data (SIMD) instruction sets) and parallelization using GP-GPUs. *The computational algorithms incorporate robust open-source libraries for mathematical operations, resource allocation, and data management. *The API design considers multiple stakeholder needs and ensure integration with existing and future ecosystems for data science, machine learning, and AI. *Kynema-FMB is written in modern C++ and leverages Kokkos as its performance-portability library with inspiration from the kynema stack.

Sprague, MichaelA.↗

Progress Towards NQA-1 for Cardinal in FY25

Cardinal is a wrapping of the GPU-oriented spectral element Computational Fluid Dynamics (CFD) code NekRS and the Monte Carlo particle transport code OpenMC within the Multiphysics Object-Oriented Simulation Environment (MOOSE). Cardinal provides high-resolution thermal-hydraulics and/or radiation transport feedback to MOOSE multiphysics simulations. Multiphysics feedback is implemented in a geometry-agnostic manner which eliminates the need for rigid one-to-one mappings. A generic data transfer implementation also allows NekRS and OpenMC to couple to any MOOSE application, enabling a broad set of multiphysics capabilities. Cardinal simulations can also leverage combinations of MPI, OpenMP, and GPU resources. Cardinal continuous development and improvement efforts have led to the software being considered as a high-fidelity design and licensing tool for key areas of nuclear reactor relevant physics, including neutron transport, fluid flow, heat transfer, and mechanical processes. The fast development and expansion of the software from a pure R&D framework towards its application in the nuclear industry and regulation require a focus on developing, enhancing,and maintaining Cardinal’s software quality through strict adherence to a Software Quality Assurance (SQA) framework and SQA program. To facilitate compliance with SQA standards, the Cardinal SQA Program was initiated during Fiscal Year 2023 (FY23). During the development of the Cardinal SQA Program, multiple gaps have been identified. These gaps are primarily related to model verification and code pedigree as they relate to the use of Cardinal as an analysis tool. These gaps were captured in a report published in 2023. A second report highlighted the progress made during Fiscal Year 2024 (FY24) and described Argonne’s effort to document and integrate software verification within Cardinal’s software development process. This report documents the progress made towards NQA-1 for Cardinal in the Fiscal Year 2025 (FY25). All cases in the expanded Continuous Integration (CI) suite of NekRS are included in this report which test the solvers and modules available in NekRS exhaustively. The NekRS tests are integrated with the Cardinal CI suite and made available in publicly accessible Github documentation. Following the CI practice permits integrating of source code changes frequently and ensuring that the integrated codebase clears the verification testing for the software. Also in this report is a brief overview of the development of the Cardinal Software Quality Assurance Plan (SQAP) that was done in FY25, though it should be noted that the rest of the documentation for the SQA program needs to be developed in a future step of this task.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Virtual prototyping of liquid metal blanket performance in fusion pilot plant

Liquid metal blanket is a dominant design option for the next step fusion devices responsible for harvesting energy from fusion reaction, and simultaneously producing fuel for the same reaction through tritium breeding. Liquid metal blankets introduce additional complexity to the design due to fluid motion, fluid structure interaction, and magnetohydrodynamic (MHD) effects arising from the motion of the conducting fluid through the magnetic field. They are also directly affected by the plasma heat flux and neutronic fluence. PPPL is currently developing a virtual prototyping system for numerical analysis of the liquid metal blankets for future fusion devices. The system has a customized 3D computational fluid dynamics (CFD) code in its core, allowing MHD flow and conjugate heat transfer analysis in blankets fluids and solids. The code was successfully used before for dual coolant blanket analysis [A. Khodak et al., Fusion Eng. and Des. 137 (2018)]. Recently the same code was modified to allow verified simulation of MHD flows at high Hartmann numbers of several thousand typical for blanket applications. CFD code receives volumetric heat source distribution from the neutronic analysis based on MCNP code. In addition, direct tritium breeding simulation will be performed allowing optimization of the blanket performance. 2D axisymmetric version of neutronics code will be used for rapid optimization, with 3D version employed for detailed analysis. The surface heat distribution on the plasma facing wall will be defined by the software HEAT allowing 3D modeling of the heat flux based on the magnetic field distribution including gyro-orbit effects. Results of thermal analysis are imported into structural analysis code also included in the system. Finally, direct import of CAD geometry will be used for analyzing all components and as a result design option can be efficiently optimized.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Integrated modeling methodology for ash agglomeration in poly-disperse fluidized beds using particle population framework

This article discusses a unique agglomeration modeling methodology developed based on binary collisions to combine the effects of heterogeneity in ash chemistry and granular physics. A simple population balance is defined to find changes in the particle size distribution (PSD) of a fluidized bed. Thermodynamic equilibrium calculations and a computational fluid dynamics (CFD) code, are used to obtain hydrodynamic parameters. A method to calculate and use a distribution of collision frequencies in fluidized beds, corresponding to a poly-disperse particle size distribution, was developed in order to incorporate the particle level heterogeneities. A distribution of collision frequencies obtained for poly-dispersed particles showed a three orders of magnitude higher collision frequency amongst the smaller particles than the coarser ones, at the initiation of agglomeration. Ash agglomeration occurred when the slag amount (binder) was less than 10 wt% at temperatures <850 °C for Pittsburgh No. 8 coal.

42 ENGINEERING↗

Convolutional Neural Network–Aided Temperature Field Reconstruction: An Innovative Method for Advanced Reactor Monitoring

In this study, the capabilities of a physics-informed convolutional neural network (CNN) for reconstructing the temperature field from a limited set of measurements taken at the boundaries of internal flows are demonstrated. Such an approach enables the development of less invasive monitoring methods for real-time plant diagnostics. As a test case, a Molten Salt Fast Reactor (MSFR) design was selected. This circulating fuel reactor has received interest from both scientific and industrial communities due to its intrinsic safety and sustainability. Molten salt flows in such reactors, however, can present highly localized temperature peaks that can induce significant thermal stresses onto the vessel walls. At these local maxima, the salt temperature may exceed a thousand kelvins, which makes a direct measurement challenging or even unfeasible. The proposed CNN algorithm allows one to detect indirectly such discontinuities through an accurate, albeit indirect, temperature measurement method during reactor operation. The datasets employed to train and test the machine learning models in the present work were generated with Nek5000, a computational fluid dynamics (CFD) code developed at Argonne National Laboratory. The CNN algorithm is trained with CFD results that span a set of MSFR operational power and flow ranges. Here, to demonstrate the efficacy of the algorithm, predictions are made for test cases contained within the training range but for which the CFD data were not used when training. Results demonstrate that the proposed technique properly characterizes temperature peaks and distributions within the domain for a broad range of scenarios.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Integrating Energy-Efficient Computing with Computational Research to Accelerate Energy Technology

NREL's computational sciences center hosts the largest high performance computing (HPC) capabilities dedicated to energy research while functioning as a living laboratory for energy-efficient computing. NREL's HPC capabilities support the research needs of the Department of Energy's Office of Energy Efficiency and Renewable Energy (EERE). In ten years of operation, HPC use in EERE-sponsored research has grown by a factor of 30, including work in electricity generation, energy efficiency, transportation, and energy system modeling. This paper analyzes this research portfolio, providing examples of individual use cases. The paper documents NREL's history of operating one of the world's most energy-efficient data centers while examining pathways to reduce economic and environmental impact beyond reduction of Power Usage Efficiency (PUE). This paper concludes by examining the unique opportunities created for accelerating improvements in data center efficiency created by combining an HPC system dedicated to energy research and a research program in energy-efficient computing.

97 MATHEMATICS AND COMPUTING↗

Direct Numerical Simulation of the Flow through a Randomly Packed Pebble Bed

The proposition for molten salt and high-temperature gas-cooled reactors has increased the focus on the dynamics and physics in randomly packed pebble beds. Research is being conducted on the validity of these designs as a possible contestant for the fourth-generation nuclear power systems. A detailed understanding of the coolant flow behavior is required in order to ensure proper cooling of the reactor core during normal and accident conditions. In order to increase the understanding of the flow through these complex geometries and enhance the accuracy of lower-fidelity modeling, high-fidelity approaches such as direct numerical simulation (DNS) can be utilized. Nek5000, a spectral-element computational fluid dynamics (CFD) code, was used to develop DNS fluid flow data. The flow domain consisted of 147 pebbles enclosed by a bounding wall. In the work presented, the Reynolds numbers ranged from 430 to 1050 based on the pebble diameter and inlet velocity. Characteristics of the flow domain such as volume averaged porosity, axial porosity, and radial porosity were studied and compared with correlations available in the literature. Friction factors from the DNS results for all Reynolds numbers were compared with correlations in the literature. The first- and second-order statistics show good agreement with the available experimental data. Turbulence length scales were analyzed in the flow. Reynolds stress anisotropy was characterized by utilizing invariant analysis. Overall, the results of the analysis in this study provide deeper understanding of the flow behavior and the effect of the wall in packed beds.

Yildiz, Mustafa Alper↗

Large Eddy Simulation of Convective Heat Transfer in a Random Pebble Bed Using the Spectral Element Method

The development of fluoride-cooled high-temperature reactors has drastically increased the demand for an in-depth understanding of the heat transfer (HT) in packed beds cooled by liquid salts. The complex flow fields and space-dependent porosity found in a pebble bed require a detailed understanding to ensure the proper cooling of the reactor core during normal and accident conditions. As detailed experimental data are complicated to obtain for these configurations, high-fidelity simulation such as large eddy simulation and direct numerical simulation (DNS) can be employed to create a high-resolution heat transfer numerical database that can assist in addressing industrial-driven issues associated with the heat transfer behavior of fluoride-cooled high-temperature reactors. In this paper, we performed a series of large eddy simulation using computational fluid dynamics (CFD) code NekRS to investigate the heat transfer for a bed of 1741 pebbles. Further, the characteristics of the flow, such as average, rms, and time series of velocity and temperature, have been analyzed. Porous media averages have also been performed. The simulation results show a good agreement between non-conjugate heat transfer and conjugate heat transfer. The generated data will be used to benchmark heat transfer modeling methods and local maxima/minima of heat transfer parameters. It will also be used for supporting convective heat transfer quantification for Kairos Power and benchmarking lower fidelity models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

OVERVIEW OF PARTICLE DEPOSITION MODELS FOR SPENT NUCLEAR FUEL STORAGE SYSTEMS

Deposition models were built to evaluate contaminant deposition on spent nuclear fuel (SNF) canisters. The primary contaminant of concern is chloride, which is dispersed in the atmosphere and then deposits onto the canisters. During dry storage, the primary degradation process is likely to be Chloride Induced Stress Corrosion Cracking (CISCC) at the heat-affected zones of the canister welds. It is known that stainless steel canisters are susceptible to CISCC; however, the rate of chloride deposition onto the canisters is poorly known, based on sparse field data from a small number of sites. The models presented in this study could be useful for determining the rate of deposition on the canisters and the likelihood of CISCC to help with SNF canister ageing management. The deposition models were developed with the commercial computational fluid dynamics (CFD) code STAR-CCM+. Various deposition mechanisms were considered and incorporated into the models, and a sensitivity study was conducted to determine the most important mechanisms for deposition within a SNF storage system. The models included both a vertical and horizontal configuration storage system: NAC International’s Modular, Advanced Generation, Nuclear All-purpose STORage System (MAGNASTOR®) and a NUHOMS® horizontal storage module respectively. The resulting canister deposition on the horizontal canister is visually compared with inspection data taken onsite at the Calvert Cliffs Nuclear Power Plant. These models are preliminary, and development of the models will continue. Future validation exercises are currently being planned, including the Canister Deposition Field Demonstration (CDFD) effort funded by U.S. Department of Energy office of Nuclear Energy office of Spent Fuel Waste Science and Technology. The goal of the modeling presented is to demonstrate a potential modeling technique that could be used to plan and inform SNF canister ageing management programs with predictive models for the timing and occurrence of canister CISCC.

Suffield, Sarah R.↗

Development of a Comprehensive Two-Phase Flow Database for the Validation of NEK-2P

Three-dimensional (3-D) two-phase Computational Fluid Dynamics (CFD) codes are emerging as a powerful and potentially practical tool for applications in which detailed local flow information is needed. However, two-phase flow models and associated closure relations are not well established for CFD applications, which is partly due to the lack of high-quality validation data. The main objective of this work is to develop a comprehensive database of two-phase flows that can be used to validate two-phase CFD codes such as NEK-2P. In this project, four advanced local measurement systems, including Particle Image Velocimetry and Planar Laser-Induced Fluorescence (PIV-PLIF), high-speed imaging, x-ray densitometry, and multi-sensor conductivity probe are employed to measure the local two-phase flow parameters of both gas and liquid phases. By combining these techniques, the local void fraction, bubble velocity, interfacial area concentration, bubble frequency, liquid velocity, turbulence intensity, etc., in various two-phase flow regimes can be obtained. These local measurement techniques are first used in a 25.4 mm circular pipe test section. Seven air-water two-phase flow conditions spanning the bubbly, slug, churn-turbulent, and annular flow regimes are measured in this facility. The obtained database contains the radial profiles of both gas- and liquid-phase parameters at three axial locations along the test section. The second facility used in this work contains a 30 mm × 10 mm rectangular test section, in which three two-phase flow and two single-phase flow conditions are measured. Two-dimensional distributions of local two-phase flow parameters in the cross-sectional plane are measured in this facility at three axial locations as well. A facility featuring a 3×3 electrically heated rod bundle is also designed and being constructed in this project. This facility is specially designed for optical measurements and is expected to provide high-quality boiling data in the future. Preliminary analyses have been performed for the data taken in the 25.4 mm circular pipe. Both center-peaked and wall-peaked void fraction profiles have been observed in the data depending on the two-phase flow conditions and/or developing lengths. The 1-D drift-flux model was evaluated with the newly obtained datasets, in which both gas- and liquid-phases data were directly measured. The distribution parameter model has been optimized based on a new void-profile classification method proposed in this study. The optimized drift-flux model shows a significant improvement in predicting both gas velocity and void fraction. The measured liquid-phase turbulence was used to benchmark Sato’s turbulence model considering the bubble-induced shear stress for the three tested bubbly flows. The benchmark results showed good agreement between the PIV measurements and model predictions. In the bubbly flows tested that have low void fractions less than 3%, the effect of the bubble-induced turbulence was found not significant. However, the bubble-induced shear stress becomes important with the increase of the void fraction. The Conjugate Heat Transfer (CHT) model was developed and implemented in NEK-2P. This model allows the coupled simulation of the solid domain and two-phase fluid domain, allowing the specification of realistic boundary conditions. The CHT implementation was verified first with non-boiling simulations. The predicted temperatures in the fluid domain were shown to be identical in simulations with or without the CHT model. The CHT model was then validated through simulations of three Becker benchmark CHF tests performed under both Dryout (DO) and Departure from Nucleate Boiling (DNB) conditions. Reasonably good agreement was obtained between calculated wall temperatures and corresponding experimental data.

42 ENGINEERING↗

Simulation and Analysis on Reactor Pressure Vessel (RPV) subjected to Pressurized Thermal Shock (PTS) under SBLOCA scenario by using Cardinal to support the fracture mechanics analyses

The structural components that comprise nuclear reactors and their supporting structures are subjected to harsh operating environments that can challenge their integrity, especially after exposure for extended durations or under accident condition. As one of the most significant components of a Reactor, the Reactor Pressure Vessel (RPV) is exposed to an aggressive environment during the operation time (e.g. more than 40 years). Ageing degradation mechanisms (e.g. thermos-fatigue) could grow initial defects up to a critical size, increasing the susceptibility to failure in the RPV. The conventional methods are mostly based on simple crack and structure geometries. Very limited studies consider the real conditions of the RPV subjected to a thermal shock due to a Loss of Coolant Accident (LOCA). During a LOCA event, the most severe conditions take place when the emergency core cooling (ECC) water is injected inside the cold legs filled initially with hotter water and/or steam. The rapid cooling of the down-comer and the internal RPV surface followed probably by re-pressurization of the RPV causes large temperature gradients and variation of pressure which induces thermal-mechanical stresses. In order to develop the model for integrity assessment of a reactor pressure vessel (RPV) subjected to pressurized thermal shock (PTS), a multi-physics simulation, which includes the thermo-hydraulic, thermo-mechanical and fracture mechanics analyses is necessary. The multi-physics simulations are performed using Cardinal, a wrapping of the GPU-oriented spectral element Computational Fluid Dynamics (CFD) code NekRS and other multi-physics sub-modules within the MOOSE framework. Cardinal now fully supports MOOSE stochastic perturbations of NekRS models with varying boundary conditions, initial conditions, material properties, and any other quantity which is defined by a kernel (such as coefficients in a momentum source model). The implementation is designed in a flexible manner so that scalar values are sent from MOOSE into a user scratch space in NekRS, which can then be applied for any purpose within the NekRS case files (both on the host and device). When modeling PTS, several factors can impact the results significantly. In this report, the impacts of the geometry of the model, Reynolds number and buoyancy effect are investigated. Two geometry, i.e., a simplified model and a realistic RPV model, with both laminar and turbulent flow condition are adopted for the PTS simulation with and without buoyancy effect. The purpose of the investigation is to understand the impact of these factors on the prediction of temperature history of RPV. The accurate prediction on the temperature evolution, which will be exported to Grizzly code for further analyses on the progression of aging mechanisms and their effects on the integrity of RPV structures, is very crucial. Based on the understanding of these factors, a more sophisticated model is built to analysis the PTS under SBLOCA scenario. A literature survey is conducted to pick the SBLOCA scenario for the multi-physics simulation. The analysis helps to explain the form and the transformation of the cold plum when the ECC is activated under SBLOCA. This model can be can be applied to study the PTS effect for different RPV configurations. The results can help to assess structural component degradation for advanced reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Report on Initial Sodium Testing on the Thermal Hydraulic Experimental Test Article (THETA) (Fiscal Year 2024 Final Report)

The Thermal Hydraulic Experimental Test Article (THETA) is a facility that is used to develop sodium components and instrumentation as well as to acquire experimental data for validation of reactor thermal hydraulic and safety analysis codes. The facility simulates nominal thermal hydraulic conditions as well as protected/unprotected loss of flow accidents in a sodium-cooled fast reactor (SFR). High fidelity distributed temperature profiles of the developed flow field may be acquired with Rayleigh backscatter based optical fiber temperature sensors. The facility was designed in partnership with systems code experts to tailor the experiment to ensure the most relevant and highest quality data for code validation. THETA is comprised of a traditional primary coolant and secondary coolant system. The primary system is submerged in the pool of sodium and consists of a pump, electrically heated core, intermediate heat exchanger, and connected piping and thermal barriers (redan). The secondary system, located outside of the sodium pool, consists of a pump, sodium to air heat exchanger, and connected piping and valves. In fiscal year 2023, thermal stratification tests were completed with the primary system online, while the secondary system was being constructed [1]. These tests had shown that the core barrel and intermediate heat exchanger (IHX) outlet required increased thermal insulation. The THETA primary system was removed from METL, cleaned, thermal insulators installed, and then inserted into METL Test Vessel 4. At the time of this writing the THETA primary and secondary system are operational. During this fiscal year 100+ hours of testing was completed to characterize thermal hydraulic phenomena associated with steady state and transient conditions in a pool type liquid metal cooled reactor. A majority of the testing campaign was completed to satisfy the experimental data acquisition requirements for the GAIN Voucher with Oklo, CRADA 2021-21121. THETA is still operational at the time of this publication and future testing is planned for fiscal year 2025. Work is underway to publish existing and future data to an online database to facilitate collaboration with SFR engineers looking to validate their systems code or computational fluid dynamics models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Glovebox Fires Phase 1: Simulations of Open Burner Experiments

Understanding the duration of a fire required to breach a glovebox via glove ports and/or windows within a facility is important in developing mitigation strategies for safety. To facilitate this understanding, both simulation and experiments are utilized. The approach is to first validate the Computation Fluid Dynamics (CFD) code, Fire Dynamics Simulator (FDS), with experimental data collected by New Mexico Tech (NMT) and then perform simulations of full-scale rooms containing gloveboxes to assess numerous scenarios which would otherwise be cost-prohibitive experimentally. This report provides comparison to first-phase experiments involving a fire without a glovebox. A subsequent report will provide comparison to experiments involving a glovebox. The results indicate good agreement with FDS tending to over predict the pre-mixed and diffusion-mode tests by 6% and 10%, respectively.

42 ENGINEERING↗

Thermochemical Nonequilibrium Modeling in a Continuous-Galerkin, Finite-Element Framework

The presented work discusses the implementation, verification, and validation of Park's two-temperature model in a scalable, computational fluid dynamics (CFD) code developed at the US Department of Energy's Oak Ridge National Laboratory (ORNL). The implementation of Park's two-temperature model was verified through 0D test cases involving an adiabatic reactor and a nitrogen thermal bath. The implementation was then validated through comparisons with other validated CFD codes and experimental data on a hypersonic cylinder and double cones. These are standard benchmark test cases for thermochemical non-equilibrium (TCNE) modeling, and all data are shared publicly. The verification and validation results showed that ORNL's in-house CFD code could model complex, high-speed flow problems with and without TCNE modeling. This work is essential for future research involving 3D shock wave/boundary layer interactions (SBLIs).

Nutter, Nicole↗

Aerodynamic design using numerical optimization

The procedure of using numerical optimization methods coupled with computational fluid dynamic (CFD) codes for the development of an aerodynamic design is examined. Several approaches that replace wind tunnel tests, develop pressure distributions and derive designs, or fulfill preset design criteria are presented. The method of Aerodynamic Design by Numerical Optimization (ADNO) is described and illustrated with examples.

Murman, E. M.↗

Study of non-equilibrium transport phenomena

Nonequilibrium phenomena due to real gas effects are very important features of low density hypersonic flows. The shock shape and emitted nonequilibrium radiation are identified as the bulk flow behavior parameters which are very sensitive to the nonequilibrium phenomena. These parameters can be measured in shock tubes, shock tunnels, and ballistic ranges and used to test the accuracy of computational fluid dynamic (CFD) codes. Since the CDF codes, by necessity, are based on multi-temperature models, it is also desirable to measure various temperatures, most importantly, the vibrational temperature. The CFD codes would require high temperature rate constants, which are not available at present. Experiments conducted at the NASA Electric Arc-driven Shock Tube (EAST) facility reveal that radiation from steel contaminants overwhelm the radiation from the test gas. For the measurement of radiation and the chemical parameters, further investigation and then appropriate modifications of the EAST facility are required.

Sharma, Surendra P.↗