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 37 records · Page 2

OH Roaming and Beyond in the Unimolecular Decay of the Methyl-Ethyl-Substituted Criegee Intermediate: Observations and Predictions

Alkene ozonolysis generates short-lived Criegee intermediates that are a significant source of hydroxyl (OH) radicals. Here this study demonstrates that roaming of the separating OH radicals can yield alternate hydroxycarbonyl products, thereby reducing the OH yield. Specifically, hydroxybutanone has been detected as a stable product arising from roaming in the unimolecular decay of the methyl-ethyl-substituted Criegee intermediate (MECI) under thermal flow cell conditions. The dynamical features of this novel multistage dissociation plus a roaming unimolecular decay process have also been examined with ab initio kinetics calculations. Experimentally, hydroxybutanone isomers are distinguished from the isomeric MECI by their higher ionization threshold and distinctive photoionization spectra. Moreover, the exponential rise of the hydroxybutanone kinetic time profile matches that for the unimolecular decay of MECI. A weaker methyl vinyl ketone (MVK) photoionization signal is also attributed to OH roaming. Complementary multireference electronic structure calculations have been utilized to map the unimolecular decay pathways for MECI, starting with 1,4 H atom transfer from a methyl or methylene group to the terminal oxygen, followed by roaming of the separating OH and butanonyl radicals in the long-range region of the potential. Roaming via reorientation and the addition of OH to the vinyl group of butanonyl is shown to yield hydroxybutanone, and subsequent C–O elongation and H-transfer can lead to MVK. A comprehensive theoretical kinetic analysis has been conducted to evaluate rate constants and branching yields (ca. 10–11%) for thermal unimolecular decay of MECI to conventional and roaming products under laboratory and atmospheric conditions, consistent with the estimated experimental yield (ca. 7%).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Extending TOUGH + HYDRATE with a parallel particle transport simulator: numerical investigation of sand production during gas production from hydrate deposits

A new parallel code for simulating particle transport in porous media is integrated with the TOUGH + HYDRATE simulator to investigate sand production associated with gas production from unconsolidated gas hydrate-bearing sediments (HBS). Here, the parallel coupled simulator is named THMPT and uses the integral finite difference method to describe the Darcian and non-Darcian flow of fluids and heat transport, the finite element method to describe the associated geomechanical changes, and the discrete element method to track the trajectory of individual sand particles within the HBS. The THMPT simulator is written in Fortran, incorporates multiple optimized algorithms, and can comprehensively address the coupled flow, thermal, chemical, geomechanical, and particle transport processes that characterize the system behaviors during gas production from HBS. The simulator can capture all processes involved in sand particle transport in porous media, including sand detachment, collision, clogging (i.e., bridging), and migration. A benchmark case study of sand production in the course of depressurization-induced gas production from a representative HBS reveals various distinct microscopic particle migration mechanisms and the adverse impact of sand particle detachment, transport, and clogging. The numerical investigation also examines the effect of bottomhole pressure on mitigating sand production. The simulation results indicate that sand clogging near the wellbore significantly reduces permeability, decreasing gas production by at least 50%. Lastly, the efficiency of gravel packing in mitigating sand production is numerically evaluated, revealing that the structure of the porous media appears to profoundly influence the macroscopic motion behavior of sand particles and sand clogging characteristics.

discrete element method↗

Self-powered Through-wall communication for dry cask storage monitoring

Many nuclear facilities, such as spent fuel storage dry casks and nuclear reactor pressure vessels, are entirely sealed by metal layers to prevent harmful radiation. For safety and security operations, the temperature, pressure, radiation, and humidity inside the vessel needs to be closely monitored. However, no practical technology is currently available to realize the through-wall data communication and monitoring for these vessels due to the inside harsh environment of high temperature and nuclear radiation. In this paper, an innovative self-powered wireless through-wall data communication system for the nuclear environment is presented, which demonstrates a successful solution to such challenges. The presented system is composed of four modules, i.e., energy harvester with power management circuits, ultrasound wireless communication using high-temperature piezoelectric transducers, electronic circuits for sensing and data transmission, and radiation shielding for electronics. Here, constitutive functions of each module were firstly designed and followed by the system integration. Experiments were conducted subsequently to validate the designed functions and evaluate the performance of the integrated system. Results showed that the average power of over 40 mW was harvested from the thermal flow inside the nuclear spent fuel canisters which could provide enough energy to operate the sensing and data communication systems. The gamma radiation test results showed that the thermoelectric energy harvester and ultrasound transceivers can withstand radiation dosing over 100 Mrad. Furthermore, temperature shock tests demonstrated that the entire system including the shielded electronics can survive and maintain their functionalities at temperatures as high as 195°C. Under the in-lab mocked-up high temperature conditions and radiation shielding, the proposed system is foreseen to survive and operate stably for fifty years inside a nuclear spent fuel canister, and send the frequency modulated data out of the canister for 3 s in every 10 min.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Statistical Study of Energy Dissipation in Magnetic Structures During Turbulent Reconnection in the Earth's Magnetotail

Magnetic reconnection is a ubiquitous plasma phenomenon that plays a critical role in particle heating and energization. During reconnection, the topology of magnetic field rearranges, depositing energy into the surrounding plasma through bulk flow, thermal heating, or non-thermal particle acceleration. While the pathways of this transformation from magnetic energy into kinetic have been studied extensively in recent years through theoretical or case-by-case observations, comprehensive statistical studies remain limited. In this paper, we present a statistical investigation using data from the Magnetospheric Multiscale (MMS) mission, and detail the particle energization mechanisms in magnetic structures found near reconnecting regions in turbulent Earth's magnetotail. We find that electrons with motion perpendicular to the magnetic field dominate $\vec{j}$ ⋅ $\vec{E}$ dissipation. In contrast to the conventional picture of unidirectional energy transfer to particles by laminar two-dimensional (2D) reconnection, we find that energy exchange within magnetic structures during turbulent reconnection tends to be bidirectional with only a small positive bias from electromagnetic fields to particles. Specific electron energization mechanisms are quantified, including those due to parallel electric field, Fermi energization from curvature drift, betatron heating from magnetic field inhomogeneity, and polarization drift.

Wang, Rachel [Princeton University, NJ (United Sta↗

Magnetic reconnection in the era of exascale computing and multiscale experiments

Astrophysical plasmas have the remarkable ability to preserve magnetic topology, which inevitably gives rise to the accumulation of magnetic energy within stressed regions including current sheets. This stored energy is often released explosively through the process of magnetic reconnection, which produces a reconfiguration of the magnetic field, along with high-speed flows, thermal heating and nonthermal particle acceleration. Either collisional or kinetic dissipation mechanisms are required to overcome the topological constraints, both of which have been predicted by theory and validated with in situ spacecraft observations or laboratory experiments. However, major challenges remain in understanding magnetic reconnection in large systems, such as the solar corona, where the collisionality is weak and the kinetic scales are vanishingly small in comparison with macroscopic scales. The plasmoid instability or formation of multiple plasmoids in long, reconnecting current sheets is one possible multiscale solution for bridging this vast range of scales, and new laboratory experiments are poised to study these regimes. In conjunction with these efforts, we anticipate that the coming era of exascale computing, together with the next generation of observational capabilities, will enable new progress on a range of challenging problems, including the energy build-up and onset of reconnection, partially ionized regimes, the influence of magnetic turbulence and particle acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗

ELM-free enhanced D α H-mode with near zero NBI torque injection in DIII-D tokamak

Enhanced D α H-mode (EDA H-mode), an ELM-free H-mode regime, and the concomitant quasi-coherent mode (QCM) are explored in neutral beam heated, lower single null plasmas with near zero external torque injection. This regime exhibits good energy confinement (H 98y2 ~ 1) with β N ~ 2, high density, regime access at low input power, and no ELMs. New time–resolved measurements of electron and ion density, temperature, plasma rotation, and radial electric field as well as turbulence, flows, thermal and particle transport, and linear instability calculations examining the dynamics of the QCM are presented. Measurements using multiple fluctuation diagnostics reveal the QCM to be a separatrix spanning mode, peaking just inside the separatrix, existing in a wide range of k ⊥ ρ s ~ 0.1–1.2 with multiple harmonics, and propagating with a very small phase velocity in the plasma frame, where k ⊥ is the binormal wavenumber and ρ s is the ion sound radius. Linear gyrokinetic simulations of this EDA H-mode discharge with CGYRO indicate that the trapped electron mode (TEM) and electron temperature gradient mode are the dominant instabilities in the radial region where QCM is also unstable. Qualitative analysis indicates that the properties of TEM are consistent with the experimentally observed characteristics of the QCM. These similarities suggest that the QCM is consistent with a TEM instability existing in the edge region of the EDA H-mode plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Two dimensional topology optimization of heat exchangers with the volume fraction method

We perform a comparison study of two topology optimizations methods applied to the design of a two fluids heat exchanger modeled with a coupled thermal-flow problem. The flow follows an isothermal and incompressible Stokes-Brinkman equation and the heat transfer is governed by a convection-diffusion equation without internal generation and high Peclet number. To keep the two fluid phases separated, we solve two Stokes-Brinkman equations, where the Brinkman term models the other phase as a solid. These two velocity fields are then fed to the heat transfer equation. Our goal is to maximize the enthalpy at the cold outlet while constraining the pressure drop. We first solve the design modeling the solid and fluid phases with a volume fraction variable. A SIMP-like penalization in the Brinkman term drives the optimization to a discrete design. The cost and constraint function derivatives are calculated with the library pyadjoint and the optimization is performed by IPOPT. We present optimized designs in two dimensions and discuss the influence of the parameters.

Beck, VictorA.↗

Proton beam power limits for stationary water-cooled tungsten target with different cladding materials

The proton beam power limit for a solid-tungsten spallation target is largely determined by beam induced thermomechanical structural loads and decay heat power deposition, while its lifetime is limited by radiation damage and fatigue life of the target materials. In this paper, we studied the power limits of a stationary water-cooled solid tungsten target concept. Tantalum clad tungsten was considered as a reference case. Being a low activation material, zircaloy 2 cladding option was studied and its decay heat driven power limit was compared with the reference case. Zirconium alloys have proven operations records in spallation target and nuclear fission environments, supported by materials data obtained from post irradiation examinations. Recent study also demonstrated feasibility of diffusion bonding zirconium to tungsten using vanadium foil inter layer. Particle transport simulations code FLUKA was used to calculate energy deposition and decay heat power deposition in the target, based on the beam parameters technically feasible at the Second Target Station of the Spallation Neutron Source at Oak Ridge National Laboratory. The energy deposition data were used for flow, thermal, and structural analyses to determine the beam intensity limit on the target concept studied. The decay heat deposition data were used to calculate the transient temperature evolution in the tungsten volumes in a loss of coolant accident (LOCA) scenario to determine its beam power limit. For a 1.3 GeV proton beam, the power limit on a stationary target was 400 kW for a tantalum clad target model and 800 kW for a zircaloy 2 clad target model.

Lee, Yong Joong↗

Two Dimensional Topology Optimization of Heat Exchangers with the Density and Level-Set Methods

We design heat exchangers using two topology optimization approaches: the density, i.e. volume fraction and level set methods. Our goal is to maximize the heat exchange between two fluids in separate channels while constraining the pressure drop across each channel. The heat exchanger is modeled with a coupled thermal-flow formulation. The flow is governed by an isothermal and incompressible Stokes-Brinkman equation and the heat transfer is governed by a convection-diffusion equation with high Peclet number. We solve one set of Stokes-Brinkman equations per fluid. Each Brinkman term in the flow equation serves to model the other phase as a solid, thereby preventing mixing. We first represent the solid and fluid phases using a volume fraction variable and apply a SIMP-like penalization in the Brinkman term to drive the optimization to a discrete design. The cost and constraint function derivatives are automatically calculated with the library pyadjoint and the optimization is performed by the Method of Moving Asymptotes. In a second optimization formulation, we use the level set approach to define the interface that separates the two fluids. Pyadjoint calculates the shape derivatives of the cost and constraint functions and the Hamilton-Jacobi advects the interface, allowing for topological changes. We present results in two dimensions and discuss the advantages and disadvantages of each approach.

42 ENGINEERING↗

Characterization of Natural Consolidated Halloysite Nanotube Structures

Halloysite is a unique 1:1 clay mineral frequently appearing with nanotubular morphology, and having surfaces of different polarity with interesting and important technological applications. HNTs can be consolidated naturally in the earth by pressure and thermal flows. In this study of natural consolidated HNTs, the strength and hardness of these materials were found to be dependent on the presence of impurities (gibbsite, alunite, quartz, and other silica minerals), which accounted for the increased stability of such samples. In the absence of impurities, the strength of consolidated HNTs was significantly lower. The first 3D mapping of the pore structure of natural consolidated HNT is provided. The contributions of the porosity within the nanotubes and between the nanotubes were delineated using a combination of non-invasive ultra-small and small-angle X-ray scattering (USAXS/SAXS) analyses, BET/BJH pore size analyses, and computed tomography studies. A total porosity of 40%, as determined by X-ray attenuation and He porosimetry, was found for the natural consolidated HNTs, of which about one-third was due to the inter-HNT porosity. Nano-X-ray computed tomography (nano-XCT) analyses also indicated that 76% of the inter-HNT pores were smaller than 150 nm in diameter. The intra-HNT pore size determined by combined USAXS/SAXS and BET/BJH was about 10 nm. This pore network information is essential for the utilization of natural consolidated HNTs as a model geomaterial to investigate the effects of surface characteristics on confined fluid flow.

36 MATERIALS SCIENCE↗

Role of Parallel Solenoidal Electric Field on Energy Conversion in 2.5D Decaying Turbulence with a Guide Magnetic Field

We perform 2.5D particle-in-cell simulations of decaying turbulence in the presence of a guide (out-of-plane) background magnetic field. The fluctuating magnetic field initially consists of Fourier modes at low wavenumbers (long wavelengths). With time, the electromagnetic energy is converted to plasma kinetic energy (bulk flow+thermal energy) at the rate per unit volume of J· E for current density J and electric field E. Such decaying turbulence is well known to evolve toward a state with strongly intermittent plasma current. Here we decompose the electric field into components that are irrotational, E ir , and solenoidal (divergence-free), E so . E ir is associated with charge separation, and J · E ir is a rate of energy transfer between ions and electrons with little net change in plasma kinetic energy. Therefore, the net rate of conversion of electromagnetic energy to plasma kinetic energy is strongly dominated by J · E so , and for a strong guide magnetic field, this mainly involves the component E so,∥ parallel to the total magnetic field B. We examine various indicators of the spatial distribution of the energy transfer rate J ∥ · E so,∥ , which relates to magnetic reconnection, the best of which are (1) the ratio of the out-of-plane electric field to the in-plane magnetic field, (2) the out-of-plane component of the nonideal electric field, and (3) the magnitude of the estimate of current helicity.

79 ASTRONOMY AND ASTROPHYSICS↗

Reduced Order Models Generation for HTGRs Pebble Shuffling Procedure Optimization Studies

This report provides an initial study for producing reduced-order models (ROMs) of pebble-bed high temperature gas reactor (HTGR) models for the purposes of design optimization. As an initial study, this work is meant to be exploratory---identifying useful workflows and methods for ROM generation---and not meant to be a catch-all analysis of HTGR ROM generation and usage for optimization. This report summarizes three tasks performed in Fiscal Year 2022: 1) the creation of HTGR model, 2) the sensitivity analysis of model design parameters, and 3) an introduction to ROM generation techniques. The representative HTGR model created in this work is a multiphysics equilibrium-core using the BlueCRAB (comprehensive reactor analysis bundle) reactor analysis application, coupling four physical phenomena: neutronics, streamline depletion, porous flow thermal hydraulics, and pebble heat conduction. Part of the model creation was identifying some design parameters and quantities of interest that are relevant in an optimization analysis and adjustable in the model. The sensitivity analysis utilized a polynomial chaos expansion methodology to compute global sensitivity metrics. This analysis showed that thermal hydraulics parameters and quantities of interest had a relatively small impact on simulation results. Finally, the ROM generation work involved exploring three different ROM methodologies: polynomial regression, a Gaussian process, and artificial neural networks. Using a cross-validation technique to characterize ROM performance, the Gaussian process and single-layer artificial neural networks showed the most promising results. Overall, this study was insightful and the lessons learned will be invaluable for the eventual development of an HTGR design optimization workflow.

97 MATHEMATICS AND COMPUTING↗

Development of an All-Aqueous Thermally Regenerative Redox Flow Battery to Support Fossil Fuel Assets

Low-temperature thermal energy, a largely untapped resource, holds significant promise for large-scale electrical power generation globally. Various stationary sources, including industrial entities and thermal power plants, emit considerable low-temperature heat that currently remains unutilized. This energy is often overlooked because its low temperature renders it unsuitable for efficient power generation using conventional methods. However, current research is exploring diverse technologies capable of converting low-temperature heat into grid-scale power, aiming to enhance grid efficiency, further decarbonization initiatives, and facilitate a shift toward more decentralized power systems. One such innovative technology is the thermally regenerative battery (TRB), noted for its high power and energy densities compared to similar technologies, positioning it as a potential game-changer in power generation. TRBs integrate two scalable and well-established unit operations: a redox flow battery and a distillation column. This integration suggests that once an effective TRB chemistry is established, the pathway to commercialization could be expedited. The copper-based thermally regenerative ammonia battery (Cu aq -TRAB) stands out as the first TRB that circumvents electrodeposition/dissolution reactions, stabilizing Cu(I) and Cu(II) within the electrolyte and maintaining stability of all electroactive species in an aqueous phase. This stabilization has led to improvements in coulombic efficiency, open circuit potential, and copper solubility, thereby enhancing power density, energy density, and overall energy efficiency. Preliminary tests were conducted to determine the effects of various electrolyte species on the performance metrics of the battery, both theoretically and experimentally. These tests revealed that the solubility of copper in the Cu aq -TRAB electrolyte was constrained by the Cu(I)-NH 3 complex. Adjusting the background electrolyte to 5 M NH4Br and the ligand concentration to 4 M NH 3 enabled the copper concentration to reach a maximum of 0.6 M. This modification led to an estimated theoretical maximum energy density of 9.5 Wh L -1 for the Cu aq -TRAB. Additionally, full cell testing indicated a tradeoff between peak power and energy density with varying copper and ammonia concentrations. Increasing the applied current density during discharge linearly raised the average power output, with a minimal reduction in energy density due to a balance between higher ohmic overpotential and reduced time for undesirable ammonia crossover. Furthermore, a comprehensive numerical sensitivity analysis of the complete Cu aq -TRAB system was performed. This analysis aimed to assess how the battery and the distillation column responded to changes in system input parameters, providing insights into optimal research directions for enhancing system performance. The analysis revealed that at room temperature, battery power was significantly more sensitive to ohmic losses than to mass transfer, with reaction rates having minimal impact. This trend continued even at higher temperatures. Also, the thermal energy required for ammonia separation was studied, showing that increased temperatures generally reduced energy requirements, except in low-pressure scenarios above 65 °C. An investigation into membrane performance in the Cu aq -TRAB was undertaken, given the significant impact of ammonia transport control and ohmic losses on system performance. Various membranes were evaluated to identify key performance metrics. Among the tested membranes, Selemion CMVN exhibited the highest performance, with a peak power density of 84 mW cm -2 and average values of 26 ± 6.8 mW cm -2 for power density and 2.9 Wh L -1 for energy density at an applied current density of 50 mA cm -2 . An economic assessment indicated a levelized cost of storage at $410 per MWh under optimal conditions, highlighting the commercial potential of the Cu aq -TRAB when utilizing cost-effective, readily available materials.

25 ENERGY STORAGE↗

In-situ thermal transport measurement of flowing fluid using modulated photothermal radiometry

In situ thermal transport measurement of flowing fluid could be useful for the characterization and diagnosis of practical thermal systems such as fluid heat exchangers and thermal energy storage systems. Despite abundant reports on the ex-situ thermal conductivity measurement of stagnant fluids, a suitable technique for the thermal conductivity measurement of flowing fluid has been rarely reported. Here, this paper presents the thermal conductivity measurement of flowing fluid within a pipe using a non-contact modulated photothermal radiometry (MPR) technique, where the surface of the pipe is heated by an intensity-modulated laser and the heat diffuses into the fluid with suitable modulation frequency. We design a tube section with small wall thickness suitable for the MPR measurements to maximize the sensitivity of the thermal response to the fluid properties while minimizing the lateral heat spreading effect. Intrinsic thermal conductivity of different fluids was obtained within a proper range of frequency and flow velocity where the forced convection effect is negligible. The forced convection effect became prominent at high flowing velocity and at low modulation frequency, leading to higher effective thermal conductivity of the fluid. It is found that the intrinsic thermal conductivity could be obtained when the flow velocity is less than 100 mm/sec and Re D 1/2 Pr 1/3 < 100 for DI water and Xceltherm oil under the specified experimental conditions, where Re D is the Reynolds number and Pr is the Prandtl number.

42 ENGINEERING↗

Efficient Thermal Energy Storage with Radial Flow in Packed Beds

Thermal energy storage (TES) is a way to store excess heat in order to generate power later. These types of systems are valuable in the solar industry, where power production can still occur when sunlight is not available. This not only increases system flexibility, it drives down the cost of electricity. One option to store thermal energy is with a packed bed where the storage media resides in a cylindrical container. Flow from one end of the cylinder to the other deposits or recovers heat (axial flow). While a promising technology in terms of energy storage, it exhibits a high pressure drop that lowers overall system efficiency. This project seeks to flow the heat transfer fluid through the storage media in the radial direction. This system offers the potential to retain reasonable thermal performance while substantially lowering the pressure drop. The overall goal was to assess the performance of radial flow experimentally and numerically. Three main designs were considered in this work. The first system utilized a central pipe and four receiving pipes near the wall. The second system utilized a central pipe and an annulus with holes near the wall. The last approach considered segmenting the bed so only select zones would receive flow in the radial direction. In a high-aspect ratio system, the hole patterns in the piping are important for ensuring even flow into the bed, and a major design effort was testing holes patterns to promote even flow. The four pipe system was shown to be not feasible, as preferential flow paths occur that do not allow even and full heating of the packed bed. The annular system can lead to even flow, but the results show this is better accomplished through larger piping instead of considering variations in the hole pattern along the axial length of the pipe. Segmenting the bed can lead to similar exergetic performance when compared to axial or radial flow. However, the pressure drop in the segmented design is likely too high for practical implementation. A baseline commercial design showed exergetic efficiencies for axial, radial, and radial segments of 81.6%, 82.8%, and 80.2%, respectively. Pressure drop for the axial and radial results were 2.36 psi and 2.45 psi, respectively, with segments being nearly an order of magnitude higher. At these large scales, the aspect ratio of the system is important. From a baseline of 0.64, an aspect ratio of 0.32 for radial flow showed an exergetic efficiency of 86% but a pressure drop of 5.39 psi. All of these results provide new insights into packed bed thermal energy storage with radial flow. Competing effects must be considered when designing a radial system, and results show a radial design can show strong thermal results at the expense of system efficiencies from pressure drop. The COMSOL models used in this analysis are available to the public and other researchers. Public benefits include a potential increase in the efficiency of packed bed thermal energy storage. Higher efficiency in storage promotes green energy technologies while reducing costs.

14 SOLAR ENERGY↗

Multi-stage preconditioners for thermal–compositional–reactive flow in porous media

We present a family of multi-stage preconditioners for coupled thermal-compositional-reactive reservoir simulation problems. The most common preconditioner used in industrial practice, the Constrained Pressure Residual (CPR) method, was designed for isothermal models and does not offer a specific strategy for the energy equation. For thermal simulations, inadequate treatment of the temperature unknown can cause severe conver-gence degradation. When strong thermal diffusion is present, the energy equation exhibits significant elliptic behavior that cannot be accurately corrected by CPR's second stage. In this work, we use Schur-complement decompositions to extract a temperature subsystem and apply an Algebraic MultiGrid (AMG) approximation as an additional preconditioning stage to improve the treatment of the energy equation. We present results for several two-dimensional hot air injection problems using an extra heavy oil, including challenging reactive In-Situ Combustion (ISC) cases. We show improved performance and robustness across different thermal regimes, from advection dominated (high Péclet number) to diffusion dominated (low Péclet number). The number of linear iterations is reduced by 40–85% compared to standard CPR for both homogeneous and heterogeneous media, and the new methods exhibit almost no sensitivity to the thermal regime.

42 ENGINEERING↗

Development and Validation of SAM Multi-dimensional Flow Model for Thermal Mixing and Stratification Modeling

Thermal mixing and stratification in large pools or enclosures are very important phenomena that are critical to nuclear reactor safety. Because of the wide ranges of time and length scales associated with such phenomena, accurate modeling and simulation of thermal mixing and stratification remain as the key unresolved, challenging problems for reactor transient analyses. In traditional system analysis codes, simplified zero-dimensional (0-D) models are widely used for their high numerical efficiency, but they generally suffer from very limited prediction accuracies or range of applicability. Like traditional system analysis codes, the current version of the SAM code has implemented such simplified 0-D and one-dimensional mixing models. On the other hand, high-resolution Computational Fluid Dynamics (CFD) tools are often used to model complex thermal mixing and stratification phenomena. They are, however, generally numerically expensive, and they require large amounts of computational resources. It is therefore desirable to implement advanced and efficient thermal mixing and stratification modeling capabilities embedded in a system analysis code. This approach will improve the accuracy of reactor safety analyses when thermal mixing and stratification are involved, and also avoid using the large computational resources needed for high-resolution CFD analysis. Currently, with the support of the U.S. DOE Office of Nuclear Energy’s Nuclear Energy Advanced Modeling and Simulation program, an effort has been launched to develop and implement a multi-dimensional flow model in the system analysis code SAM, and demonstrate its applications to model thermal mixing and stratification phenomena in large enclosures. The main outcomes of this research and development activity are summarized in this report, which presents an attempt to include a built-in advanced multi-dimensional flow model in a system analysis code with the focus on overcoming the simulation challenges of thermal mixing and stratification phenomena. In this report, we start with the introduction of existing SAM code capabilities to simulate thermal mixing and stratification phenomena, which is followed by a short summary of the multi-dimensional model implemented in the SAM code, including both the physical model and the Finite Element Method code implementation. In this study, two options were implemented in the SAM code to model turbulent flows: a relatively simple built-in turbulence model and an interface to accept externally computed turbulent viscosities (e.g., from a high- resolution CFD simulation). Code validation studies on this newly added capability were then carried out to compare SAM simulation results with experimental data from the SUPERCAVNA facility, which was designed to study the complex flow recirculation and thermal stratification phenomena relevant to sodium fast reactor designs. In this study, one transient and two steady-state test cases were used for code validation. Different approaches have been used to model the complex turbulence flow fields in the SUPERCAVNA facility. A highly simplified zero-equation turbulence model was first used, but it was determined that it is too simple to capture the complex turbulence flow fields in these test cases. Subsequently, the code validation continued with the use of turbulent viscosity data from high-resolution STAR-CCM+ CFD simulations to improve the accuracy of the results. Using this approach, the SAM simulation results showed very good agreement with both the SUPERCAVNA experimental data and STAR-CCM+ simulation results. In this report, we demonstrate the development, implementation, and successful validation of a multi-dimensional flow model in the SAM code, which aims to improve the simulation accuracy for complex thermal mixing and stratification phenomena. Lessons have also been learned, including that in cases where the flow fields are not well predicted by the zero-equation model, the iinput of turbulent viscosities from an external source can enhance the overall predictive capabilities needed to accurately capture complex thermal-hydraulic phenomena. Therefore, future research will be needed to further improve the code’s capabilities, such as by developing a more efficient and robust approach to capture the turbulence effects in the SAM code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Review of heat transfer enhancement techniques for single phase flows

The thermal energy exchange between a flowing fluid and its confining channel is a ubiquitous process in modern society. To enhance the fluid-to-wall or wall-to-fluid heat transfer, several techniques have been developed to maximize the contact area between the fluid and the inner wall and/or disrupt the flow to enhance circulation or induce turbulence. Deployment of channels having features capable of enhancing heat transfer enables the reduction of heat exchanger size while maintaining performance. Reduction in equipment size is critical due to the ability to minimize the required volume of costly working fluids and to mitigate potential safety concerns associated with total system fluid volume. Here, a comprehensive review of single-phase heat transfer enhancement techniques is presented. The article provides a thorough comparison by analyzing the heat transfer rate, pressure drop, and other operational aspects. Single-phase heat transfer enhancement methods are divided into active and passive techniques. Active methods such as electrohydrodynamic (EHD), magnetohydrodynamics (MHD), or mechanical motion require external power to create enhancement. Passive methods such as dimples, fins, or tape inserts do not require external input and rely only on surface modification. Although active methods are more expensive and difficult to implement compared to passive techniques, it enables active control of heat transfer augmentation. Finally, this review develops and summarizes key learning data for design optimization enabled by additive manufacturing and machine learning algorithms, helping to inform these next-generation heat exchanger design methodologies for a plethora of modern applications such as electrification of vehicles, computing, and classical industries.

42 ENGINEERING↗