SEARCH · Search NASA

Results for “FLUID POWER”

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 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Computational fluid dynamic modeling to determine the indoor environment of an electron-ion collider service building

The design for the Electron-Ion Collider (EIC) calls for several service buildings that house various power supplies and control electronics for the collider ring itself. In order to operate within specified conditions, the ambient air entering the power supplies needs to be within a certain temperature range while dissipating the heat from losses. Proper cooling is therefore a necessity in the service buildings to ensure that every aspect of the EIC works as intended. Since the EIC is in the design stage, we are evaluating the indoor environment of the service building using the current design specifications. We have researched multiple forms of literature and performed the necessary calculations to compile a list of boundary conditions that accurately represent the situation at hand. We are using computational fluid dynamics modeling to solve the conservation equations for mass, momentum, and energy (Navier-Stokes). This allows us to perform a finite element analysis which will give us the flow distribution in the room as well as temperature profiles throughout the building. We have obtained a simulation result giving us the temperature profiles for the building and it shows that the placement of the racks and supply vents are essential to obtaining an even temperature distribution. This model will provide a basis for design decisions which will affect the overall cooling of the service building without extending the schedule and avoiding a costly reworking of the cooling system.

43 PARTICLE ACCELERATORS↗

DOE OSTI2021

Investigation of Nonlinear Difference-Frequency Wave Excitation on a Semisubmersible Offshore-Wind Platform with Bichromatic-Wave CFD Simulations: Preprint

The natural surge and pitch frequencies of semisubmersible offshore wind platforms are typically designed to be below the wave frequencies to avoid direct excitation. However, surge or pitch resonance can be excited by the nonlinear low-frequency loads generated by irregular incident waves. Second-order potential-flow models with added Morison drag have been found to underpredict this low-frequency excitation and response. As part of the OC6 project, Computational Fluid Dynamics (CFD) simulations were performed to enable a better understanding of the low-frequency loads and the limitations of lower-fidelity models. The focus of the current paper is to set up a computationally cost-effective CFD simulation of a fixed semisubmersible platform to investigate nonlinear difference-frequency loads and establish the corresponding uncertainty in the results. Due to the high computing cost, CFD simulations of irregular waves can be challenging. Instead, simulations were performed with bichromatic waves having a shorter repeat period. A preliminary comparison with Quadratic Transfer Functions from second-order potential-flow theory shows that CFD models consistently predict higher nonlinear wave loads at the difference frequency, likely due to flow separation and viscous drag not accounted for in potential-flow theory.

2nd order↗

DOE OSTI2020

Steady-state off-design modeling of the supercritical carbon dioxide recompression cycle for concentrating solar power applications with two-tank sensible-heat storage

Concentrating solar power researchers are evaluating the potential of the supercritical carbon dioxide recompression cycle to improve the thermal efficiency and decrease the capital costs of next-generation systems. This analysis investigates the steady-state off-design performance of a recompression cycle integrated with two-tank sensible-heat thermal energy storage as the ambient temperature and heat-transfer fluid (HTF) inlet conditions to the cycle change. This paper presents off-design component models and then cycle convergence and control models to maximize net cycle power output while constraining the high-side pressure and air-cooler fan power to their respective design values and fixing the cycle HTF outlet temperature to its design value. Results show that inventory control and air-cooler fan power are important control parameters that can be optimized to maximize off-design power output as the ambient temperature and HTF mass flow rate diverge from design. The high-side pressure and fan power constraints cause the cycle net power to degrade when the ambient temperature is warmer than the design value, and this study calculates a maximum mass flow rate for hot days above which the cycle cannot achieve the design HTF outlet temperature. Finally, the analysis shows that optimizing compressor shaft speeds can improve cold-day performance by around 1.5 percentage points and part-load performance by up to 2.5 percentage points versus a baseline case with no active compressor control.

14 SOLAR ENERGY↗

DOE OSTI2021

Design of a 2 MW ZrC/W-based molten-salt-to-sCO 2 PCHE for concentrated solar power

To increase the power cycle efficiency and lower the levelized cost of electricity (LCOE) of concentrated solar power (CSP) plants, printed circuit heat exchangers (PCHEs) capable of operating above 700 °C with molten chloride salt and a sCO 2 -based fluid are needed. In this paper, the design of a high-pressure, high-temperature, 2 MW PCHE comprised of a thermomechanically-robust, zirconium carbide/tungsten (ZrC/W) composite is conducted for CSP plants. The ZrC/W composite is a material with high thermal conductivity, high stiffness, and high failure strength at high temperatures, along with excellent resistance to thermal cycling and thermal shock. In this work, a thermomechanical design analysis was conducted to select appropriate material thicknesses of the ZrC/W plates, and to determine the geometrical dimensions and the thermal performance of the PCHE. The influences of the plate number and heat exchanger length on power density and pressure drop have also been systematically investigated. Economic analyses were conducted to compare the cost of ZrC/W-based PCHEs to those comprised of IN740H (a state-of-the-art, nickel-based superalloy) and 316 stainless steel (316SS). Finally, at a sufficiently high plate production rate, the manufacturing cost of ZrC/W-based PCHEs can be significantly lower, while achieving a much higher power density, compared with state-of-the-art, nickel alloy-based and stainless steel-based PCHEs.

14 SOLAR ENERGY↗

DOE OSTI2022

Modified Lucas-Washburn theory for fluid filling in nanotubes

Ultrafast water transport in carbon nanotubes (CNTs) has drawn a great deal of attention in a number of applications, such as water desalination, power generation, and biomolecule detection. With the recent experimental advances in water filling of isolated CNTs, the Lucas-Washburn theory for capillary rise in tubes needs to be revisited for a better understanding of the physics and dynamics of water filling in CNTs. Here, the Lucas-Washburn theory is corrected for the hydrodynamic entrance effects as well as the variation of capillary pressure and hydrodynamic properties with the radius and length of CNTs. Due to the large slippage in CNTs, inclusion of the entrance effects is important particularly for the initial stages of filling where a L ∝ t scaling, as opposed to L 2 ∝ t , is observed in our molecular dynamics (MD) simulations. The corrected Lucas-Washburn theory is shown to predict the water filling dynamics in CNTs as observed in MD simulations. With the corrected theory, we achieve a better understanding of capillary rise and water filling dynamics in CNTs.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2024

DASSH-F: Subchannel Based Thermal Analysis

The DASSH thermal analysis code is designed to rapidly allow a reactor design engineer to obtain flow rates requirements that satisfy peak temperature constraints in the domain. The advantage of using DASSH over a hand calculation is that it has a more rigorous treatment of the pin power distribution and coolant heat transfer within an assembly and between assemblies. The advantage of using DASSH over a conventional 3D subchannel code or a computational fluid dynamics code (CFD) is that it can obtain the desired solution in a matter of minutes in serial with minor computer memory needs. The DASSH methodology is virtually identical to SUPERENERGY-2 with additional functionalities taken from follow on work to SUPERENERGY-2 done at ANL in the 1980s. DASSH today is an integral component of the Argonne Fast Reactor analysis suite for reactor design work. DASSH obtains the power distribution from a coupled neutron-gamma heating calculation in GAMSOR (including DIF3D) at each time point of a companion fuel cycle analysis calculation with REBUS. The domain in DASSH assumes a hexagonal grid typical for fast reactors with much of the geometry information taken from the DIF3D model. DASSH assumes the assemblies that are loaded into each grid position are ducted to control the coolant flow. Considerable detail is given on the subchannel formulation of DASSH in this document. Much of the formulation and design of the code builds upon research done by previous authors with little new investigation. Thus the decisions made in developing the subchannel model used in DASSH have their origins over 50 years ago. Much of the heat transfer methodology in DASSH is built upon correlations for both the coolant mixing and heat transfer coefficients for pins and ducts. DASSH is thus not a rigorous treatment of a given problem, but a rapid assessment of the temperature field that has known limitations with respect to an experimental measurement or CFD calculation. The DASSH input and output are detailed along with usage of the software. The DASSH output provides tables of evaluated material properties and key coolant and pin temperature results. DASSH can create Python scripts that generate domain summary pictures. DASSH can also generate assembly temperature maps and VTK output files which allow the DASSH solution to be visualized. As the primary purpose of the DASSH software is to compute the coolant and fuel pin temperature distribution for a given model of a reactor, much of the output focus is giving the user quick summary tables needed to assess the performance of a given orifice flow specification. The present version of DASSH has a crude orifice search capability and an efficient orifice flow search capability. The flow search tries to meet user specified constraints for 1) peak 2-sigma clad midwall temperature, 2) peak coolant temperature, and 3) desired bulk outlet temperature. There are many development shortcomings in DASSH detailed in this document, but this version is functional for modern analysis needs. This document serves as the manual for the Fortran based DASSH software that was developed to replace the Python version of DASSH developed as part of the VTR program.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

FSVPy: A python-based package for fluorescent streak velocimetry (FSV)

Predictive constitutive equations that connect easy-to-measure transport properties (e.g., viscosity and conductivity) with system performance variables (e.g., power consumption and efficiency) are needed to design advanced thermal and electrical systems. In this work, we explore the use of fluorescent particle-streak analysis to directly measure the local velocity field of a pressure-driven flow, introducing a new Python package (FSVPy) to perform the analysis. Fluorescent streak velocimetry combines high-speed imaging with highly fluorescent particles to produce images that contain fluorescent streaks, whose length and intensity can be related to the local flow velocity. By capturing images throughout the sample volume, the three-dimensional velocity field can be quantified and reconstructed. We demonstrate this technique by characterizing the channel flow profiles of several non-Newtonian fluids: micellar Cetylpyridinium Chloride solution, Carbopol 940, and Polyethylene Glycol. We then explore more complex flows, where significant acceleration is created due to microscale features encountered within the flow. We demonstrate the ability of FSVPy to process streaks of various shapes and use the variable intensity along the streak to extract position-specific velocity measurements from individual images. Thus, we demonstrate that FSVPy is a flexible tool that can be used to extract local velocimetry measurements from a wide variety of fluids and flow conditions.

Mechanics↗

DOE OSTI2022

Topology optimization of an airfoil fin microchannel heat exchanger using artificial intelligence

High-performance microchannel heat exchangers are needed to supply heat for power conversion for nuclear microreactors. An airfoil fin microchannel design, constructed of Alloy 617 with helium as the working fluid, is analyzed and optimized using a design of experiments with artificial intelligence techniques. The use of airfoil fins offers the potential to reduce pressure drop across the heat exchanger, as compared to other types of channel configurations. A framework for topology optimization of airfoil fin printed circuit heat exchangers (PCHEs) has been developed that can be readily extended to different fin sizes and shapes, as well as different inlet and operating conditions, materials of construction, and working fluids. An optimization procedure is developed that employs computational fluid dynamics for a set of design points identified using Latin hypercube sampling. Computational fluid dynamics is used to analyze a simplified two-channel configuration where five design parameters are varied – inlet angle, fin scale, extent of staggering, transverse and longitudinal pitches. Two methods (a 5D polynomial and a regression neural network) are compared for generating surrogate models and the resulting response surface approximation is input to a genetic algorithm that is used to identify a set of optimal parameters. The optimal geometries are found across six channel Reynolds numbers ranging from 1000 to 5000, since inlet conditions affect flow through the heat exchanger. Additionally, a set of optimal designs that maximizes heat transfer and minimizes pressure drop is identified, and a thermal stress analysis is performed on the optimal design. Correlations for the Nusselt number and Darcy friction factor are developed that can be useful for thermal hydraulic analyses using system codes. Thermal stresses are analyzed and a brief discussion of the status of code cases of PCHEs for nuclear applications is given. Testing and thermomechanical modeling is needed to facilitate future code compliance of PCHEs for high pressure and high temperature applications.

42 ENGINEERING↗

DOE OSTI2024

Elucidating the Transition of 3D Morphological Evolution of Binary Alloys in Molten Salts with Metal Ion Additives

Molten salts serve as effective high-temperature heat transfer fluids and thermal storage media used in a wide range of energy generation and storage facilities, including concentrated solar power plants, molten salt reactors and high-temperature batteries. However, at the salt–metal interfaces, a complex interplay of charge-transfer reactions involving various metal ions, generated either as fission products or through corrosion of structural materials, takes place. Simultaneously, there is a mass transport of ions or atoms within the molten salt and the parent alloys. The precise physical and chemical mechanisms leading to the diverse morphological changes in these materials remain unclear. Here, to address this knowledge gap, this work employed a combination of synchrotron X-ray nanotomography and electron microscopy to study the morphological and chemical evolution of Ni-20Cr in molten KCl-MgCl 2 , while considering the influence of metal ions (Ni 2+ , Ce 3+ , and Eu 3+ ) and variations in salt composition. Our research suggests that the interplay between interfacial diffusivity and reactivity determines the morphological evolution. The summary of the associated mass transport and reaction processes presented in this work is a step forward toward achieving a fundamental comprehension of the interactions between molten salts and alloys. Overall, the findings offer valuable insights for predicting the diverse chemical and structural alterations experienced by alloys in molten salt environments, thus aiding in the development of protective strategies for future applications involving molten salts.

36 - MATERIALS SCIENCE↗

DOE OSTI2024

MODELING AND METHOD FOR OPTICAL PERFORMANCE OF A RECEIVER COATING AND ITS DEGRADATION DUE TO OPERATION

The optical properties of the coating used on a concentrated solar power (CSP)’s receiver directly impact the performance of the power plant. Receiver coatings are designed to absorb as much heat as possible and transfer it to the circulating fluid. Over time, the absorptivity of the receiver coating degrades due to environmental and operational conditions. This project aims to (1) create a test apparatus that allows coatings to undergo accelerated lifetime testing under conditions that mimic what the coating will experience on-tower; and (2) develop a mathematical model that based on the test data predicts absorptivity over time. To create the mathematical model, samples of BrightSource Energy’s Gen1 and Gen2 coatings were tested first looking at individual failure modes (such as high temperature cycling) and then for combined failure modes in the newly developed apparatus. In addition to taking optical measurements, the samples were destroyed to allow for metalogical evaluation to understand what was happening on a microscopic level. Pre-existing models from the literature that describe accelerated testing as well as individual phenomena observed in the coating were then combined to create the model. Additionally, the Gen2 and Gen3 coatings are being testing at the Plataforma Solar de Almería (PSA) facility in Almeria, Spain to confirm that the apparatus test results are similar to those received on-sun. (PSA’s facility closely mimics on-tower conditions). The goal is to enable industry players to test coatings using the new apparatus, gather absorptivity measurements (without needing to destroy samples), and input the information into the model to determine how the coating will perform over time, thus allowing the optical performance of coatings to be compared even before they are used on-tower. The model will be publicly available by the end of the project (currently projected as December 2024).

14 SOLAR ENERGY↗

DOE OSTI2024

Multilayer Silicon Carbide Composite Material Technology for High-Temperature Concentrated Solar-Thermal Power Components

In 2012, the U.S Department of Energy defined aggressive targets to achieve lower component costs and higher system efficiencies for concentrated solar-thermal power (CSP), and this, in turn, has led to the exploration of technology options that can operate at higher temperatures [1]. These next-generation CSP options, referred to as Generation 3 (a.k.a. Gen3), are targeting temperatures at or above 700 °C for the energy being delivered to the power cycle, and the more challenging plant conditions have necessitated a review and selection of alternative receiver heat transfer fluids as well as a search for materials that can meet the associated high-temperature component requirements. Nickel-based alloys are currently being considered, but these generally experience a significant drop in strength at temperatures > 775 °C [2] and may not be able to achieve corrosion and other lifetime requirements. Furthermore, these alloys are expensive, frequently have cost and schedule volatility, and offer little potential for lower cost at high production volumes. As an alternative, Ceramic Tubular Products, LLC (CTP) has developed a multilayer silicon carbide composite that can complement or replace alloys currently being considered for these Gen3 CSP applications.

14 SOLAR ENERGY↗

DOE OSTI2021

Design and Optimization of a Gas-Cooled, Airfoil Fin Microchannel Heat Exchanger

High-performance microchannel heat exchangers are needed to supply heat for power conversion for nuclear microreactors. An airfoil fin microchannel design, constructed of Alloy 617 with helium as the working fluid, was analyzed and optimized using a design of experiments with artificial intelligence and machine learning techniques. The use of airfoil fins offers the potential to reduce pressure drop across the heat exchanger, as compared to other types of channel configurations. A framework for topology optimization of airfoil fin PCHEs has been developed that can be readily extended to different fin sizes and shapes, as well as different inlet and operating conditions, materials of construction, and working fluids. An optimization procedure was developed that employs computational fluid dynamics for a set of design points identified using Latin hypercube sampling. STAR-CCM+ was used to analyze a simplified two-channel configuration where five parameters were varied – inlet angle, fin scale, extent of staggering, transverse and longitudinal pitches. Two methods were compared for generating surrogate models – a 5D polynomial and a regression neural network. A response surface approximation was created from the surrogate models and input to a genetic algorithm. The genetic algorithm identified a set of optimal points on the Pareto front. The optimal geometry was found across six channel Reynolds numbers ranging from 1000 to 5000 to analyze how varying inlet conditions affects the optimal design. A set of optimal designs that maximizes heat transfer and minimizes pressure drop was identified, and a thermal stress analysis was performed on the optimal design. This work has developed a digital framework for the expedient topology design and evaluation of PCHE designs for gas-cooled microreactor applications. Correlations for the Nusselt number and Darcy friction factor were developed that can be useful for thermal hydraulic analyses using system codes. A thermal stress analysis was conducted and a brief discussion of the status of code cases of PCHEs for nuclear applications is given. Testing and thermomechanical modeling is needed to facilitate future code compliance of PCHEs for high pressure and high temperature applications.

42 ENGINEERING↗

DOE OSTI2021

Preventing stent failure using adaptive shear responsive endovascular implant

A stent apparatus, system, and method that senses wall shear stress by measuring fluid flow at localized areas within the stent, that processes measured information through an integrated circuit, and selectively sends power to mechanically controllable stent surfaces which results in localized geometric changes. In various embodiments the stent apparatus, system, and method sends data to outside the body in real time.

Mukerjee, Erik V.↗

DOE OSTI2022

NUMERICAL PREDICTIONS OF MEAN PERFORMANCE AND DYNAMIC BEHAVIOR OF A 10 MWe SCO2 COMPRESSOR WITH TEST DATA VALIDATION

High-fidelity aerodynamic analysis has been demonstrated for mean performance and unsteady dynamics in a sCO2 compressor designed by Hanwha Power Systems Americas for a 10 MWe Concentrating Solar Power plant. Simulations were performed with CRUNCH CFD® software tool that was matured to accurately model near critical real fluid effects in sCO2. Pre-test predictions for mean performance were validated with test data collected later. Performance predictions were accurate and captured sensitivity of the efficiency to inlet temperature of CO2 as well as steep drop-off at high flow rates due to condensation in the inlet throat. Detailed analysis was performed to understand the source of these performance losses at near critical conditions. Unsteady dynamic effects in the compressor at off-design conditions were also identified and quantified. In particular a system wide “condensation surge” condition was detected at high flow coefficients that results in large amplitude pulsations with accompanying mass flow fluctuations at low frequencies and has potential to cause damage in closed loop systems.

hosangadi, ashvin↗

DOE OSTI2024

Preventing stent failure using adaptive shear responsive endovascular implant

Mukerjee, Erik V.↗

DOE OSTI2024

Multiphysics Demonstration of Temperature-Driven Assembly Bowing in SFRs using MOOSE-Based Codes

Core bowing is an important passive safety mechanism in liquid metal cooled fast reactors. When the core restraint system is properly designed, temperature and flux gradients influence assemblies in the core to bow into less reactive configurations during accident scenarios, resulting in negative reactivity feedback. Prediction of core bowing involves complex interplay of radiation transport, impacts of fluid flow and heat transfer on duct temperature, and mechanical responses to the induced temperature and flux gradients. Under the U.S. Department of Energy Office of Nuclear Energy’s Advanced Modeling and Simulation (NEAMS) Program [1], an integrated multiphysics approach is being developed to model the core bowing phenomena in liquid metal-cooled fast reactors with the Multiphysics Object Oriented Simulation Environment (MOOSE) [2]. In this methodology, the MOOSE-based reactor physics code Griffin [3] will solve the neutron transport equation and determine the power distribution. With the detailed power distribution from Griffin, the subchannel analysis codes MOOSE-Subchannel [4] and Pronghorn [5] are utilized to calculate the assembly temperature distribution. MOOSE’s Solid Mechanics [6] and Contact [7] Modules are leveraged to calculate the thermal expansion and duct bowing displacement with the duct wall temperature from thermal hydraulics calculation. In this work, an initial one-way coupling demonstration of the integrated multiphysics approach has been performed on a seven-assembly problem based on the sodium-cooled fast reactor ABR-1000 design [8]. The neutronics calculation with Griffin is not yet involved in the current simulation. MOOSE-Subchannel and Pronghorn evaluate fluid and solid temperature based on a fixed power distribution. In addition, one-way coupling is utilized in this coupled calculation, via Pronghorn passing the duct temperature data to the MOOSE Solid Mechanics calculation. An assessment of the Solid Mechanics module was performed in parallel to verify duct bowing behavior with duct-to-duct contact phenomenon [9]. The displacement from MOOSE Solid Mechanics is not yet transferred back and utilized in the Pronghorn and MOOSE-Subchannel calculation. This model will be available on the National Reactor Innovation Center (NRIC) Virtual Test Bed (VTB) repository [10]. Future stages of this work will involve solving problems of increasing complexity as well as adding more physics (e.g. reactor physics) to the integrated workflow to reach the end goal of modeling the core bowing phenomenon with an integrated multiphysics workflow.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2026

DASSH-F: Subchannel Based Thermal Analysis

The DASSH thermal analysis code is designed to rapidly allow a reactor design engineer to obtain flow rates requirements that satisfy peak temperature constraints in the domain. The advantage of using DASSH over a hand calculation is that it has a more rigorous treatment of the pin power distribution and coolant heat transfer within an assembly and between assemblies. The advantage of using DASSH over a conventional 3D subchannel code or a computational fluid dynamics code (CFD) is that it can obtain the desired solution in a matter of minutes in serial with minor computer memory needs. The DASSH methodology for pin lattice models is virtually identical to SUPERENERGY-2 with additional functionalities taken from follow on work to SUPERENERGY-2 done at ANL in the 1980s. DASSH today is an integral component of the Argonne Fast Reactor analysis suite for reactor design work. DASSH obtains the power distribution from a coupled neutron-gamma heating calculation in GAMSOR (including DIF3D) at each time point of a companion fuel cycle analysis calculation with REBUS. The domain in DASSH assumes a hexagonal grid typical for fast reactors with much of the geometry information taken from the DIF3D model. DASSH assumes the assemblies that are loaded into each grid position are ducted to control the coolant flow. The user can alternatively provide their own geometry and power profile instead of inheriting it from DIF3D. Considerable detail is given on the subchannel formulation of DASSH in this document. Much of the formulation and design of the code builds upon research done by previous authors with little new investigation. Thus the decisions made in developing the subchannel model used in DASSH have their origins over 50 years ago. Much of the heat transfer methodology in DASSH is built upon correlations for both the coolant mixing and heat transfer coefficients for pins and ducts. DASSH is thus not a rigorous treatment of a given problem, but a rapid assessment of the temperature field that has known limitations with respect to an experimental measurement or CFD calculation. The DASSH input and output are detailed along with usage of the software. The DASSH output provides tables of evaluated material properties and key coolant and pin temperature results. DASSH can create Python scripts that generate domain summary pictures. DASSH can also generate assembly temperature maps and VTK output files which allow the DASSH solution to be visualized. As the primary purpose of the DASSH software is to compute the coolant and fuel pin temperature distribution for a given model of a reactor, much of the output focus is giving the user quick summary tables needed to assess the performance of a given orifice flow specification. The present version of DASSH has a crude orifice search capability and a sufficient orifice flow search capability. The flow search tries to meet user specified constraints for 1) peak 2-sigma clad temperature, 2) peak coolant temperature, and 3) desired bulk outlet temperature. This document serves as the manual for the Fortran based DASSH software that was developed to replace the Python version of DASSH developed as part of the VTR program.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

High Fidelity CFD Simulations Supporting the KP-FHR

Kairos Power, LLC, is developing its version of the Fluoride-cooled High-temperature Reactor, the KP-FHR. The design uses a pebble bed core with fluoride salt as a coolant. The pebbles used in the KP-FHR have a diameter of 4 cm, with a shell fuel region where TRISO particles are embedded. A Pebble bed core design is adopted by several Gen IV reactors, They boast many benefits, such as fuel integrity, highly efficient heat transfer, and passive safety. However, it is challenging to accurately predict temperature and flow inside a pebble bed. Traditional approaches use the porous media model, which regards the pebble bed as a continuous medium, but with different temperature fields representing different levels, such as the fluid temperature, pebble surface temperature, and pebble center temperature. Empirical heat transfer correlations are adopted to calculate the heat transfer coefficient between different phases. However, empirical correlations are usually validated with experimental data, which usually lacks detail inside the pebble bed. The available experimental data is also generally at a high Reynolds number, which falls outside of the conditions of KP-FHR. Explicit computational fluid dynamics (CFD) simulations of randomly packed pebble beds have only become feasible recently. This is thanks to the rapid development of computational power and scalable algorithms. In this work, we used the Spectral Element Method (SEM) CFD code NekRS to simulate the randomly packed pebble bed in a cylindrical container. NekRS, which is the GPU variant of Nek5000, but refactored to utilize the computational power of GPUs using the OCCA library to run on hybrid architecture high performance computing systems. It was initially developed with the libParamunal library, but truncated and tuned for large-scale turbulence simulation. As a result, the SEM reaches higher precision with the same degrees of freedom by using a high-order Lagrange polynomial basis distributed on Gauss-Lobatto-Legendre quadrature inside each element, compared to lower-order methods, such the Finite Volume Method and Finite Element Method. The report is divided into five parts. We start with a general discussion of the pebble bed reactor, along with a specific investigation into the KP-FHR. The second part presents the numerical methodology. In the third part, we study a modular pebble bed with 1741 pebbles in a container of 7 pebble-diameter radius. Beyond LES simulations done by NekRS, we also leveraged the thermal radiation model in OpenFOAM to study heat transfer under no-forced-flow scenarios. Then, in the fourth part we simulated a pebble bed similar to the size of the Hermes Test Reactor. The total number of pebbles is in these simulations is 34,374. The container radius is 14 pebble-diameters. Finally, the report concludes in part five, with a discussion of future work.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗