Search NASA⌕ Search

SEARCH · Search NASA

Results for “Thermal Flow”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Modeling the Nonlinear Rheology of Polymer Additive Manufacturing

This report summarizes molecular and continuum simulation studies focused on developing physics - based predictive models for the evolution of polymer molecular order during the nonlinear processing flows of additive manufacturing. Our molecular simulations of polymer elongation flows identified novel mechanisms of fluid dissipation for various polymer architectures that might be harnessed to enhance material processability. In order to predict the complex thermal and flow history of polymer realistic additive manufacturing processes, we have developed and deployed a high - performance mesh - free hydrodynamics module in Sandia's LAMMPS software. This module called RHEO – short for Reproducing Hydrodynamics and Elastic Objects – hybridizes an updated - Lagrange reproducing - kernel method for complex fluids with a bonded particle method (BPM) to capture solidification and solid objects in multiphase flows. In combination, our two methods allow rapid, multiscale characterization of the hydrodynamics and molecular evolution of polymers in realistic processing geometries.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Design of experiments to spectroscopically characterize radiation flow in stochastic media

Precise characterization of experimental radiation flow is required to validate the high energy density physics models, numerical methods, and codes that are used to simulate radiation-hydrodynamics phenomena such as thermal radiation transport in stochastic media. The Cassio code is used to simulate thermal radiation flow through inhomogeneous, stochastic-media-foam configurations containing optically thick clumps dispersed within an optically thin background aerogel. Cassio can model small inhomogeneous problems directly, but most problems require approximations to meet computer limitations on run-times and memory usage. Various examples of these approximations are methods that produce, in one calculation, an ensemble-averaged solution and associated standard deviation; reduced spatial dimensionality with approximate geometries; and full material homogenization with no geometric detail. Cassio simulations are used to design experiments at the OMEGA-60 Laser Facility that can measure the radiation flow using the spatially resolved COAX absorption spectroscopy diagnostic. The experimental platforms flow radiation through foam targets ranging from a background-only aerogel, to a single configuration of a specified stochastic medium, to a fully homogenized foam of the background and clump materials. Under constant total clump mass, larger clumps (here, larger than 10 μm diameter) will mix more slowly with the background such that the bulk radiation flow is faster than it would be in a fully homogenized material. The COAX platform can be used to infer temperature and density profiles in both the background material and clumps, simultaneously, and therefore to differentiate radiation flow in a range of stochastic and homogeneous media.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Thermal Modeling and Simulation of the Packed-bed Thermal Energy Storage Combined with INL Thermal Energy Distribution System

Dynamic Energy Transport and Integration Lab (DETAIL) at Idaho National Laboratory is to support experimental demonstration and validation research on Nuclear-Renewable Hybrid Energy System [1]. The Thermal Energy Distribution System (TEDS) is a thermal-hydraulic flow loop in DETAIL with its own dedicated control system 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. Among various TES options, the packed-bed TES is adopted in TEDS because it offers a low-cost single-tank thermal storage option compared to the traditional two-tank TES. However, since the TES tank is filled with granular fillers having different thermophysical properties from those of TES tank wall, there is a potential thermo-mechanical issue to be carefully addressed like thermal ratcheting which may pose a significant design concern for the packed-bed TES tank. Thermal ratcheting is a thermomechanical process caused by the repeated rearrangement of granular filler inside a TES tank during continuous thermal cycling operation of the packed-bed TES system. If the thermally induced stress exceeds the yield strength of a TES tank wall, catastrophic consequences may happen like rupture of the TES tank. Thus, it is important to understand the phenomenon to ensure the robust operation of the packed-bed TES tank. Given that thermal ratcheting is caused by complex interaction of thermal transport in the porous bed and solid mechanics, the accurate prediction of transient thermal behavior of the packed-bed TES tank, which is the focus of this paper, is critical to the reliable thermal ratcheting analysis. This paper discusses the numerical modeling, simulation, and validation studies that are ongoing at INL to investigate the transient thermal behavior of the packed-bed TES. Of particular concern is the transient thermal process occurring in the packed-bed TES unit that is operated in conjunction with the INL TEDS. The main goal of this research is three-fold: (i) provide preliminary insights into the transient thermal behavior of the TEDS TES tank, (ii) support the thermal measurement and validation plan for TEDS experiment, and (iii) provide transient thermal boundary conditions to support the reliable thermal ratcheting analysis of the TEDS TES tank. For the transient thermal modeling and analysis, a CFD model was developed, and the validity of the modeling approach was examined via comparing the numerical simulation results with the experimental data obtained from various design characteristics of packed-bed TES tanks. Then, the present modeling method was applied for the transient thermal analysis of the TEDS TES tank, and the results are discussed along with the potential improvement of data acquisition strategy for the future TEDS experiments for more precise validation study.

25 ENERGY STORAGE↗
DOE OSTI2022

Pressure–strain interaction. III. Particle-in-cell simulations of magnetic reconnection

How energy is converted into thermal energy in weakly collisional and collisionless plasma processes, such as magnetic reconnection and plasma turbulence, has recently been the subject of intense scrutiny. The pressure–strain interaction has emerged as an important piece, as it describes the rate of conversion between bulk flow and thermal energy density. In two companion studies, we presented an alternate decomposition of the pressure–strain interaction to isolate the effects of converging/diverging flow and flow shear instead of compressible and incompressible flow, and we derived the pressure–strain interaction in magnetic field-aligned coordinates. Here, we use these results to study pressure–strain interaction during two-dimensional anti-parallel magnetic reconnection. We perform particle-in-cell simulations and plot the decompositions in both Cartesian and magnetic field-aligned coordinates. We identify the mechanisms contributing to positive and negative pressure–strain interaction during reconnection. Furthermore, this study provides a roadmap for interpreting numerical and observational data of the pressure–strain interaction, which should be important for studies of reconnection, turbulence, and collisionless shocks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

O 2 based resonantly ionized photoemission thermometry analysis of supersonic flows

Characterization of the thermal gradients within supersonic and hypersonic flows is essential for understanding transition, turbulence, and aerodynamic heating. Developments in novel, impactful non-intrusive techniques are key for enabling flow characterizations of sufficient detail that provide experimental validation datasets for computational simulations. In this work, Resonantly Ionized Photoemission Thermometry (RIPT) signals are directly imaged using an ICCD camera to realize the techniques 1D measurement capability for the first time. The direct imaging scheme presented for oxygen-based RIPT (O 2 RIPT) uses the previously established calibration data to direct excite various resonant rotational peaks within the S-branch of the C 3 Π, ( v = 2) ← X 3 Σ( v ′ = 0) absorption band of O 2 . The efficient ionization of O 2 liberates electrons that induce electron avalanche ionization of local N 2 molecules generating N 2 + , which primarily deexcites via photoemissions of the first negative band of N 2 + ( B 2 Σ <#comment/> u + − <#comment/> X 2 Σ <#comment/> g + ) . When sufficient lasing energy is used, the ionization region and subsequent photoemission signal is achieved along a 1D line thus, if directly imaged can allow for gas temperature assignments along said line; demonstrated here of up to five centimeters in length. The temperature gradients present within the ensuing shock train of a supersonic under expanded free jet serves as a basis of characterization for this new RIPT imaging scheme. The O 2 RIPT results are extensively compared and validated against well-known and established techniques (i.e., CARS and CFD). The direct imaging capability fully realizes the technique’s fundamental potential and is expected to be the standard of implementation going forward. The direct imaging capability can play instrumental roles in future scientific studies that rely upon acute characterization of thermal gradients within a medium that cannot be easily resolved by a point. Furthermore, the removal of the spectrometer greatly reduces the cost, complexity, and optical alignment associated with prior RIPT measurements.

42 ENGINEERING↗
DOE OSTI2023

Prediction and Validation of Flow Properties in Porous Lattice Structures

High-porosity metal foams have been extensively studied as an attractive candidate for efficient and compact heat exchanger design. With the advancements in additive manufacturing, such foams can be manufactured with controlled topology to yield highly tailorable mechanical and transport properties. In this study, a lattice Boltzmann method (LBM)-based pore-scale model is implemented to simulate the fluid flow in additively manufactured (AM) metal foams with unit cell topologies of Cube, Face Diagonal (FD)-Cube, Tetrakaidecahedron (TKD), and Octet lattices. The pressure gradient versus average velocity profiles predicted by the LBM model were validated against in-house measurements on the AM lattice samples with the same unit cell topologies. Based on the simulation results, a novel hybrid model is proposed to accurately predict the volume averaged flow properties (permeability and inertial coefficients) of the four structures. Specifically, the linear LBM (neglecting inertial forces) is first implemented to obtain the intrinsic permeability, and then the standard LBM is applied to obtain the inertial coefficient. Convenient correlations for those flow properties as a function of porosity and fiber diameter are constructed. The effects of the AM print qualities on the flow properties are also discussed. The advantages of the hybrid model compared to the polynomial fitting approach for determining flow properties are discussed and compared quantitatively. The hybrid model and presented results are valuable for flow and thermal transport evaluation when designing new metal foams for specific applications and with different materials and topologies. Finally, the presented correlations based on pore-scale simulations can also be conveniently used in volume-averaged models to predict the macroscale flow behavior in such complex structures.

42 ENGINEERING↗
DOE OSTI2021

Boundary zonal flows in rapidly rotating turbulent thermal convection

Recently, in Zhang et al. ( Phys. Rev. Lett. , vol. 124, 2020, 084505), it was found that, in rapidly rotating turbulent Rayleigh–Bénard convection in slender cylindrical containers (with diameter-to-height aspect ratio $\varGamma =1/2$ ) filled with a small-Prandtl-number fluid ( ${Pr}\approx 0.8$ ), the large-scale circulation is suppressed and a boundary zonal flow (BZF) develops near the sidewall, characterized by a bimodal probability density function of the temperature, cyclonic fluid motion and anticyclonic drift of the flow pattern (with respect to the rotating frame). This BZF carries a disproportionate amount ( ${>}60\,\%$ ) of the total heat transport for ${Pr} < 1$ , but decreases rather abruptly for larger ${Pr}$ to approximately $35\,\%$ . In this work, we show that the BZF is robust and appears in rapidly rotating turbulent Rayleigh–Bénard convection in containers of different $\varGamma$ and over a broad range of ${Pr}$ and ${Ra}$ . Furthermore, tirect numerical simulations for Prandtl number $0.1 \leq {\textit {Pr}} \leq 12.3$ , Rayleigh number $10^7 \leq {Ra} \leq 5\times 10^{9}$ , inverse Ekman number $10^{5} \leq 1/{\textit {Ek}} \leq 10^{7}$ and $\varGamma = 1/3$ , 1/2, 3/4, 1 and 2 show that the BZF width $\delta _0$ scales with the Rayleigh number ${Ra}$ and Ekman number ${\textit {Ek}}$ as $\delta _0/H \sim \varGamma ^{0} Pr^{\{-1/4, 0\}} {Ra}^{1/4} {\textit {Ek}}^{2/3}$ ( $\{{\textit {Pr}}<1, {\textit {Pr}}>1\}$ ) and with the drift frequency scales as $\omega /\varOmega \sim \varGamma ^{0} Pr^{-4/3} {Ra}\,{\textit {Ek}}^{5/3}$ , where $H$ is the cell height and $\varOmega$ the angular rotation rate. The mode number of the BZF is 1 for $\varGamma \lesssim 1$ and $2 \varGamma$ for $\varGamma = \{1,2\}$ independent of ${Ra}$ and ${Pr}$ . The BZF is quite reminiscent of wall mode states in rotating convection.

42 ENGINEERING↗
DOE OSTI2024

Optimization of Triply Periodic Minimal Surface Heat Exchanger to Achieve Compactness, High Efficiency, and Low-Pressure Drop

With advancements in additive manufacturing (AM) techniques, high-quality triply periodic minimal surface (TPMS) structures can now be produced. TPMS walled heat exchangers (HX) hold significant potential for industrial applications and are receiving increasing attention. This paper explores the impact of various TPMS design variables on flow and thermal performance to optimize TPMS heat exchangers for compactness, high efficiency, and low pressure drop. The design variables examined include the type of TPMS lattice, unit cell size, wall thickness, aspect ratio, TPMS orientation, and equivalent thickness. The study reveals that the flow and heat transfer performance of TPMS structures are significantly affected by these design variables. For the Gyroid, Diamond, and SplitP lattices, performance is nearly identical when the surface-to-volume ratio is kept constant. The average velocity of the fluid in the TPMS HX should be 0.3 m/s. The corresponding Re is between 300~800. Thin wall thickness, small equivalent thickness, and flat lattice configurations can significantly reduce pressure drop while maintaining the overall heat transfer coefficient. Additionally, the angle between the flow direction and TPMS orientation can increase pressure drop. Three aluminum heat exchangers were successfully printed using an AM machine, and testing results are comparable with theoretical prediction.

42 ENGINEERING↗
DOE OSTI2025

Ab initio investigation of a hypersonic double cone experiment

This article presents a direct molecular simulation (DMS) of a reactive Mach 8.2 oxygen flow over a double cone geometry. The free stream conditions and article configuration generate a flow with thermal and chemical nonequilibrium, which are common attributes of hypersonic flight. This scenario was first studied experimentally at Calspan University of Buffalo Research Center’s test facility. DMS is a particle method that uses quantum mechanically derived interaction potentials to simulate molecular collisions within a flow field. Since these interaction potentials are the only modeling inputs used in the simulation, all flow features can solely be attributed to the ab initio potential energy surfaces. Hence, providing a comparison of a hypersonic ground test and numerical data anchored to quantum mechanics. The fundamental nature of DMS is leveraged to investigate molecular level mechanisms prevalent in the flow, and comparisons with lower fidelity simulations are presented to highlight the role of these first principles calculations as benchmark solutions.

Science & Technology - Other Topics↗
DOE OSTI2020

Behind the Meter Storage for Electric Vehicle Charging, Electrochemical and Thermal Energy Storage, and Solar Photovoltaic

In response to the potentially large and irregular demand from EVs, along with changing load profiles from buildings with on-site generation, utilities are evaluating multiple options for managing dynamic loads, including time-of-use pricing, demand charges, battery storage, and curtailment of variable generation. Buildings, as well as commercial, public, and workplace EV charging operations, can use a combination of electrochemical battery storage and thermal energy storage coupled with on-site generation to manage energy costs as well as provide resiliency and reliability for EV charging and building energy loads. We are completing a behind the meter storage analysis that focuses on determining the optimal system designs and energy flows for thermal and electrochemical behind the meter storage with on-site solar photovoltaic (PV) generation enabling electric vehicle charging in various climates, building types, and utility rate structures. In completing this analysis, we have developed a tool that combines existing battery models via the System Advisor Model (SAM) and building modeling software via EnergyPlus into a single interface. This tool allows us to simulate a building with a detailed battery model to properly size the battery, thermal energy storage, and solar PV systems to maximize profit for the system owner. This also allows us to assess how the battery degrades under various supervisory control dispatch algorithms to control charging/discharging; we can also see how thermal energy storage is created and used to complement the battery to reduce thermal loads in the building. With this project, we can analyze new batteries that are designed specifically for energy storage, rather than designed to be extremely energy dense for electric vehicle applications, using battery lifetime models from other national labs and the existing SAM battery model, which has detailed lifetime and degradation parameters. We can also assess novel thermal storage technologies by integrating them into the whole building energy simulation program EnergyPlus. Because the model calls both SAM and EnergyPlus, required inputs need to be compatible for both models. These inputs include, on a high-level, the following: weather files, building and electric vehicle load profiles, electricity rate tariff information, and system cost information for the stationary battery, solar PV, and thermal storage system. The various buildings we are studying for this analysis are retail big-box grocery store, commercial office building, fleet vehicle depot and operations facility, multi-family residential, and electric vehicle charging station. For these different applications, the battery and thermal storage will be dispatched differently, and the various technologies are sized differently to optimize cost.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2026

SAM Finite Volume Method Development Status Update: GCR Application, Restart, and MultiApp

The System Analysis Module (SAM) is being developed as a modern system analysis code for advanced non-light-water-reactor safety analysis under the U.S. DOE NEAMS program. Previous feasibility studies have demonstrated that a staggered-grid finite volume method (SG-FVM), implemented under the MOOSE framework, can deliver more than an order of magnitude speedup over the existing continuous Galerkin finite element method (CG-FEM) solver for liquid-cooled, incompressible but thermally expandable flow systems. This work extends the previous effort to compressible, gas-cooled reactor applications, where pressure couples directly into the mass equation adding additional nonlinearity into the equation system. New code capabilities are implemented for pebble bed high-temperature gas-cooled reactor (PB-HTGR) analysis, including a pebble bed CoreChannel component, built-in pebble bed effective thermal conductivity model and channel-to-channel crossflow model. The capabilities are tested, benchmarked, and demonstrated for problems with increased level of model and physical complexities, including the HTTU effective thermal conductivity test, the SANA passive cooling test, and a demonstration case using the GPBR200 reactor design covering steady-state operation, DLOFC and PLOFC transients. Across all cases, the SG-FVM solver demonstrated strong robustness and efficiency, and the solutions agree well with reference results and data. The finding of this work proves that SG-FVM is a viable and efficient solver pathway for compressible, gas-cooled reactor system analysis in SAM. In addition, work has been done to successfully support SAM-FVM recover/restart code feature that is essential to reactor safety analysis applications, and MultiApp code feature that is essential to multi-scale and multi-physics simulations. In summary, this work continued from previous feasibility studies, and further demonstrated that the SG-FVM will serve as a strong foundation for SAM’s advanced solver algorithm for future deployment.

Zou, Ling↗
DOE OSTI2021

Crystal plasticity model of residual stress in additive manufacturing using the element elimination and reactivation method

Selective laser melting is receiving increasing interest as an additive manufacturing technique. Residual stresses induced by the large temperature gradients and inhomogeneous cooling process can favour the generation of cracks. In this work, a crystal plasticity finite element model is developed to simulate the formation of residual stresses and to understand the correlation between plastic deformation, grain orientation and residual stresses in the additive manufacturing process. The temperature profile and grain structure from thermal-fluid flow and grain growth simulations are implemented into the crystal plasticity model. An element elimination and reactivation method is proposed to model the melting and solidification and to reinitialize state variables, such as the plastic deformation, in the reactivated elements. The accuracy of this method is judged against previous method based on the stiffness degradation of liquid regions by comparing the plastic deformation as a function of time induced by thermal stresses. The method is used to investigate residual stresses parallel and perpendicular to the laser scan direction, and the correlation with the maximum Schmid factor of the grains along those directions. The magnitude of the residual stress can be predicted as a function of the depth, grain orientation and position with respect to the molten pool. The simulation results are directly comparable to X-ray diffraction experiments and stress–strain curves.

316 stainless steel↗
DOE OSTI2022

Pressure–strain interaction. I. On compression, deformation, and implications for Pi-D

The pressure–strain interaction describes the rate per unit volume that energy is converted between bulk flow and thermal energy in neutral fluids or plasmas. The term has been written as a sum of the pressure dilatation and the collisionless analog of viscous heating referred to as Pi–D , which isolates the power density due to compressible and incompressible effects, respectively. It has been shown that Pi–D can be negative, which makes its identification as collisionless viscous heating troubling. We argue that an alternate decomposition of pressure–strain interaction can be useful for interpreting the underlying physics. Since Pi–D contains both normal deformation and shear deformation, we propose grouping the normal deformation with the pressure dilatation to describe the power density due to converging/diverging flows, with the balance describing the power density purely due to shear deformation. We then develop a kinetic theory interpretation of compression, normal deformation, and shear deformation. We use the results to determine the physical mechanisms that can make Pi–D negative. We argue that both decompositions can be useful for the study of energy conversion in weakly collisional or collisionless fluids and plasmas, and implications are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Exploring thermal runaway propagation in Li-ion batteries through high-speed X-ray imaging and thermal analysis: Impact of cell chemistry and electrical connections

Battery safety design is important to consider from the individual Li-ion cell to the level of the macro-system. On the macro-level, failure in one single cell can lead to propagation of the thermal runaway and rapidly set a whole battery pack on fire. Factors that can impact the propagation outcome, such as cell model/chemistry and electrical connection are here investigated using a combination of measurements. Several abusive tests were conducted, combining two different cell models (Molicel P42A and LG M50, both 21700s) in series and parallel connections (16 tests per configuration). Overall, a propagation outcome of 56% was measured from the 32 conducted tests, a minimum temperature of 150 °C was required to initiate propagation, and the fastest propagation occurred in 123 s. Temperature measurements were higher in series connected cells, initiating the discussion of cell chemistry and internal resistance on this effect. The difference in current-flow during thermal runaway in series and parallel connections, and how this can affect the temperature evolution is further discussed. Spatio-temporal mapping of X-ray radiography allowed us to derive the speed of thermal runaway evolution inside the battery and has shown that series connected cells, in particular P42A, occur faster. It was further observed that deviant sidewall behaviors such as temperature-induced breaches and pressure-induced ruptures occurred in P42As only respective nail-penetrated cells only.

25 ENERGY STORAGE↗
DOE OSTI2022

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles is investigated using a novel transient three-dimensional computational fluid dynamics model. The model couples particle–gas hydrodynamics of solid–gas flows laden with polydisperse particles, radiative heat transfer in non-grey absorbing, emitting and anisotropically-scattering multi-component participating media, conduction heat transfer in the gas phase, and interfacial convection heat transfer. The multiphase particle-in-cell method is used to predict high-fidelity solid–gas flow characteristics, such as the local discrete particle size distribution, with increased computational efficiency by combining the advantages of both Eulerian and Lagrangian methods. The multi-component radiative transfer model is implemented using an advanced collision-based Monte Carlo ray-tracing method. The number of the prescribed discrete particle components is found to be the key parameter affecting the computational accuracy and efficiency, which primarily depends on the size distribution of the particles. For the model particle–gas flow featuring free-falling Gamma-distributed ceramic particles exposed to concentrated solar irradiation, the particle volume fraction, radiative, fluid flow and thermal characteristics appear to converge with the increasing number of the discrete particle components. Five particle components are sufficient to obtain physically meaningful results. A further increase in the number of the particle components only slightly increases the accuracy of the numerical predictions at the expense of a rapidly increasing computational time. For five particle components, the particle vertical velocity at the receiver exit for particles with the diameter of 43.4 μm is 57% of that for the particles with the diameter of 202.8 μm. The temperatures of these two particle components increase from the initial ambient values by factors of 2 and 1.2, respectively, during the simulation time. In conclusion, the model developed allows for increased fidelity of particle–gas flow simulations with significant radiative effects.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

A Fast and Accurate Reduced-Order Model for High-Intensity Transferred Arc Discharges

Arc discharges are widely used in welding, plasma smelting, and other industrial processes, where variations in operating conditions strongly affect arc stability, temperature distribution, and energy transfer. Accurate modeling of these phenomena typically requires computationally expensive high-fidelity simulations. This study presents a hierarchy of three arc discharge models with progressively reduced physical fidelity. The high-fidelity model provides a fully physics-resolved reference, the reduced-order model reproduces the dominant thermal and flow characteristics of the arc with reasonable accuracy, and the Elenbaas-Heller model captures key trends in a simplified, rapid formulation. This hierarchy demonstrates that reduced-order model can effectively balance predictive fidelity and computational efficiency, providing practical tools for arc simulation and parametric studies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Liquid metal compatibility of pre-oxidized FeCrAl in flowing Sn

Here, a flowing Sn thermal convection loop (TCL) experiment was conducted in FeCrAlMo (Kanthal alloy APMT) tubing for 1000 h with a peak temperature of 400 °C. Oxide dispersion strengthened (ODS) Fe-(10–12)Cr-6Al and APMT specimens were exposed and all were pre-oxidized to form an external alumina scale. Unlike static exposures at 400° and 500 °C where pre-oxidation prevented dissolution, large mass losses were observed in the TCL hot leg after the specimens were cleaned in Li. The observed pitting suggests that localized failure of the pre-formed alumina surface oxide caused significant metal loss. Specimen mass gains were not observed on the cold leg implying that mass transfer did not occur and the small levels of Fe and Cr in the Sn after the experiment implies that Sn-rich phases may have formed. Post exposure room temperature tensile testing showed limited impact of the Sn exposure on the 12Cr ODS FeCrAl alloy, comparable to a 1000 h anneal in Ar at 400 °C. These results suggest that pre-oxidation is not sufficient to protect FeCrAl alloys in contact with flowing liquid Sn at 350°-400 °C.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Prediction of PbLi fluid flow and temperature field in a thermal convection loop for qualification of fusion materials

Two modeling tools, a thermohydraulics code and a computational model in COMSOL Multiphysics, have been developed, tested and then applied to the analysis of fluid flow and heat transfer in a thermal convection loop (TCL). Such a device has recently been used to experimentally evaluate corrosion compatibility of APMT (Fe–21Cr–5Al–3Mo) steel with high-temperature molten eutectic lead-lithium (PbLi) alloy at Oak Ridge National Lab, TN, USA. Here, the 1-D thermohydraulics code allows for rapid calculations of the loop parameters as a function of the applied heat flux. The supplemental COMSOL finite-element computations provide detailed 3-D velocity and temperature field data but are much more time consuming. Both modelling tools demonstrate an excellent agreement in the computed circulation velocity as well as maximum and minimum temperatures. The performed TCL analysis focuses on flow development in the “hot” and “cold” legs, formation of Dean vortices in corner regions, blocking effect of the immersed samples, and radiative and convective cooling effects under the experiment-relevant conditions for Prandtl number Pr = 0.015, Grashof numbers Gr ~10 7 , and Reynolds numbers Re ~10 4 .

36 MATERIALS SCIENCE↗