Search NASA⌕ Search

SEARCH · Search NASA

Results for “Prandtl numbers”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

CFD Simulations of Molten Salts Within Plasma-Facing Components with Twisted Tape Inserts Under One-Sided–Heated Conditions

For liquid immersion blankets (LIBs), the heat removal within different plasma-facing components (PFCs) requires substantial heat removal capabilities that could benefit from the usage of heat transfer enhancements. Of consideration are twisted tape insert heat transfer enhancements that are ribbon-like metal pieces inserted within cooling channels that increase fluid mixing, resulting in higher heat transfer and pressure drop. Unlike traditional applications, PFCs such as divertor coolant tubes experience a significantly higher heat flux on one side versus the other or entirely one-sided heating. These conditions are not well studied for molten salts where the higher Prandtl numbers characteristic of the fluids are significantly higher than water, gases, and liquid metals. As an effort to study this knowledge gap, large eddy simulations were done with Prandtl numbers representative of FLiBe at proposed operating conditions. In this study, the effects of applying a surface heat flux to one side of a tube versus the entire tube with a twisted tape insert were investigated for the heat transfer coefficient and effects on the temperature profile. The findings of the study suggest that for Prandtl number fluids above unity, the Nusselt number was largely unaffected by the change in boundary condition. Both boundary conditions showed good agreement with the Manglik and Bergles Nusselt number correlation for tight-fitting twisted tapes. Finally, the outcomes of this study will be useful for future design studies of LIBs to provide fundamental guidance implementing twisted tape inserts and their performance in one-sided–heated conditions.

molten salts↗
DOE OSTI2021

Studying the effect of the upper internal structure design on thermal stratification in sodium fast reactors

Thermal stratification and mixing in large reactor pools can pose significant issues with regard to operation and safety of nuclear plants. While significant work has been done in this area on light water reactors it is more critical for sodium reactors and there is far less data with regard to low Prandtl number fluids. Here this paper focuses on the effect of the upper internal structure on the formation of thermal stratification in a scaled sodium fast reactor upper plenum. Six different upper internal structure configurations were tested to observe the effect on thermal stratification. The degree of stratification was quantified using the stratification number. Each configuration was tested from low to high flow rates at a temperature difference of 50 and 100 °C. Over this broad array of testing, it was observed that higher flow rates, more porosity in the upper internal structure design, and a smaller temperature difference yielded a lower degree of stratification. This information can be used to aid in the design of the upper internal structure in a sodium fast reactor. Furthermore, this high resolved temperature data can be used to verify models under prototypic low Prandtl numbers.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Magnetic fields in the formation of the first stars – II. Results

ABSTRACT Beginning with cosmological initial conditions at z = 100, we simulate the effects of magnetic fields on the formation of Population III stars and compare our results with the predictions of Paper I. We use gadget-2 to follow the evolution of the system while the field is weak. We introduce a new method for treating kinematic fields by tracking the evolution of the deformation tensor. The growth rate in this stage of the simulation is lower than expected for diffuse astrophysical plasmas, which have a very low resistivity (high magnetic Prandtl number); we attribute this to the large numerical resistivity in simulations, corresponding to a magnetic Prandtl number of order unity. When the magnetic field begins to be dynamically significant in the core of the minihalo at z = 27, we map it on to a uniform grid and follow the evolution in an adaptive mesh refinement, MHD simulation in orion2. The non-linear evolution of the field in the orion2 simulation violates flux-freezing and is consistent with the theory proposed by Xu & Lazarian. The fields approach equipartition with kinetic energy at densities ∼1010–1012 cm−3. When the same calculation is carried out in orion2 with no magnetic fields, several protostars form, ranging in mass from ∼1 to 30 M⊙; with magnetic fields, only a single ∼30 M⊙ protostar forms by the end of the simulation. Magnetic fields thus suppress the formation of low-mass Pop III stars, yielding a top-heavy Pop III initial mass function and contributing to the absence of observed Pop III stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Thermovelocimetric characterization of liquid metal convection in a rotating slender cylinder

Rotating turbulent convection occurs ubiquitously in natural convective systems encompassing planetary cores, oceans and atmospheres, as well as in many industrial applications. While the global heat and mass transfer of water-like rotating Rayleigh–Bénard convection is well-documented, the dynamics in low-Prandtl-number liquid metals remain less understood. In this study, we experimentally investigate rotating Rayleigh-Bénard convection in liquid gallium (Prandtl number Pr ≈ 0.027) within a slender cylinder (diameter-to-height aspect ratio Γ = D/H =1/2) using novel thermovelocimetric diagnostic techniques that integrate simultaneous multi-point thermometry and ultrasonic Doppler velocity measurements. Our results reveal the formation of a stable, global-scale azimuthal wavenumber m = 2 quadrupolar vortex at low supercriticality. We propose that enhanced wall modes facilitated by the slender cylinder geometry interact with the bulk flow to create these large-scale axialized vortices. Furthermore, our findings imply a distinct scaling behavior for the wall-mode precession frequency in liquid metals, extending previous results obtained for moderate-Pr fluids. In conclusion, this provides new insights into wall-bulk coupling mechanisms of low-Pr rotating convective turbulence.

Magnetohydrodynamics↗
DOE OSTI2023

Data-Driven RANS Turbulence Closures for Forced Convection Flow in Reactor Downcomer Geometry

Recent progress in data-driven turbulence modeling has shown its potential to enhance or replace traditional equation-based Reynolds-averaged Navier-Stokes (RANS) turbulence models. Here, this work utilizes invariant neural network (NN) architectures to model Reynolds stresses and turbulent heat fluxes in forced convection flows (when the models can be decoupled). As the considered flow is statistically one dimensional, the invariant NN architecture for the Reynolds stress model reduces to the linear eddy viscosity model. To develop the data-driven models, direct numerical and RANS simulations in vertical planar channel geometry mimicking a part of the reactor downcomer are performed. Different conditions and fluids relevant to advanced reactors (sodium, lead, unitary-Prandtl-number fluid, and molten salt) constitute the training database. The models enabled accurate predictions of velocity and temperature, and compared to the baseline k–τ turbulence model with the simple gradient diffusion hypothesis, do not require tuning of the turbulent Prandtl number. The data-driven framework is implemented in the open-source graphics processing unit–accelerated spectral element solver nekRS and has shown the potential for future developments and consideration of more complex mixed convection flows.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Validation of Pronghorn for Natural-Circulation Molten Salt Loops

This paper presents the development and validation of a high-fidelity thermal-hydraulic model of a molten salt natural circulation flow loop, designed for integration within a digital twin framework. The study evaluates the performance of Idaho National Laboratory’s Pronghorn against experimental data from Texas A&M University Molten Salt Flow Loop (MSFL) four Hitec-salt test benchmark data. Natural circulation of high-Prandtl-number fluids exhibits complex, counter-intuitive flow patterns that make pointwise thermocouple readings unreliable. Experimental work at TAMU’s MSFL provides benchmark data, including flow visualization at a test-section and centerline steady-state temperature measurements along the loop. Validation includes four single-phase natural circulation test cases with Hitec salt. Key metrics include flow profile agreement and steady-state temperature accuracy. Pronghorn results for two-dimensional single-phase agree qualitatively with the experimental flow profile. This paper illustrates the importance of Computational Fluid Dynamics (CFD) in elucidating the behavior of high-Prandtl-number thermal-hydraulics, along with how misleading centerline temperature measurements can be. Pronghorn reproduces the axial and radial stratification that makes single thermocouple readings unreliable. Future research will focus on reduced-order modeling techniques to enable rapid simulation suitable for real-time digital twin applications. The validated cases provide a basis for developing reduced-order surrogates aimed at real-time digital-twin applications.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
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

Data-driven closure modeling for hypersonic turbulent flows

The Reynolds-averaged Navier–Stokes (RANS) equations remain a workhorse technology for simulating compressible fluid flows of practical interest. Due to model-form errors, however, RANS models can yield erroneous predictions that preclude their use on mission-critical problems. This report summarizes work performed from FY22-FY24 focused on improving RANS models for hypersonic flows using data-driven modeling and scientific machine learning. In this work we: 1. Investigate the current capabilities of RANS models in Sandia’s parallel aerodynamics and re-entry code (SPARC) for hypersonic flows with a focus on shock boundary layer interactions (SBLIs), 2. Assess several established corrections that exist in the literature aimed at improving predictions for SBLIs, 3. Develop improved models for the Reynolds stress tensor using tensor-basis neural networks, 4. Develop a neural-network-based variable turbulent Prandtl number model to reduce errors in wall heating in SBLIs. 5. Begin future investigations including employing the LIFE framework to improve wall heating predictions in SBLIs as well as the ensemble Kalman filter. We find that current RANS models in SPARC are deficient for complex SBLI flows. In particular, no current model jointly predicts wall heat flux, wall shear stress, and wall pressure with reasonable accuracy. Existing corrections help, but do not alleviate this issue altogether. The development of improved models for the Reynolds stress tensor via tensor-basis neural networks results in more predictive RANS models across a suite of low-speed and high-speed cases. For hypersonic boundary layers, the inclusion of the wall-normal Reynolds stress via TBNNs has an appreciable impact on the wall-normal momentum balance and wall quantities. However, we find that improvements to the Reynolds stress tensor do not address the over-prediction in wall heat flux in SBLIs. We find that a neural-network-based variable turbulent Prandtl number model systematically and substantially improves wall heating predictions for a range of SBLI cases.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

General Correlations for Laminar Flow Friction Loss and Heat Transfer in Plain Rectangular Plate-Fin Cores

In this paper, new generalized correlations for predicting the average fanning friction factor f and average Nusselt number Nu for laminar flow in plain plate-fin compact cores of rectangular cross section are presented. These are based on extended experimental data, as well as three-dimensional computational simulations, obtained for a broad range of fin density and geometrical attributes. The results indicate that while the fully developed forced convection scales only with the interfin channel cross-sectional ratio α (fin spacing by fin height), the entrance region hydrodynamic and thermal performance is additionally a function of the fin-core length L, flow Reynolds number Re, and fluid Prandtl number Pr. The developing flow and convection is further shown to scale as: (fRe)~(L/d h Re) 1/2 , and Nu ~(L/d h Re) 1/2 Pr 1/3 Φ(α), where f, Re, and Nu are all based on the hydraulic diameter dh of the interfin flow channel. Generalized correlations for both (fRe) and Nu are developed by the corresponding scaling of the forced convection behavior and asymptotic matching of the entrance or developing flow (short fin-core flow length) and the fully developed flow (large fin-core flow length) region performance. Finally, the predictions from these correlations are found to be within ±15% of all available experimental data for air, water, and glycol (0.71 ≤ Pr ≤ 10), and fin cores with 0 < α ≤ 1.

36 MATERIALS SCIENCE↗
DOE OSTI2021

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↗
DOE OSTI2022

Low-Reynolds-number oscillating boundary layers on adiabatic slopes

We investigate the instabilities and transition mechanisms of Boussinesq stratified boundary layers on sloping boundaries when subjected to oscillatory body forcing parallel to the slope. We examine idealized forms of boundary layers on hydraulically smooth abyssal slopes in tranquil mid- to low-latitude regions, where low-wavenumber internal tides gently heave isopycnals up and down adiabatic slopes in the absence of mean flows, high-wavenumber internal tides, shelf breaks, resonant tide–bathymetry interactions (critical slopes) and other phenomena associated with turbulence ‘hot spots’. In non-rotating low-Reynolds-number flow, increased stratification on the downslope phase has a relaminarizing effect, while on the upslope phase we find transition-to-turbulence pathways arise from shear production triggered by gravitational instabilities. When rotation is significant (low slope Burger numbers) we find that boundary layer turbulence is sustained throughout the oscillation period, resembling stratified Stokes–Ekman layer turbulence. Simulation results suggest that oscillating boundary layers on smooth slopes at low Reynolds number ($Re\leqslant 840$), unity Prandtl number and slope Burger numbers greater than unity do not cause significant irreversible turbulent buoyancy flux (mixing), and that flat-bottom dissipation rate models derived from the tide amplitude are accurate within an order of magnitude.

42 ENGINEERING↗
DOE OSTI2022

Floquet stability analyses of stratified oscillating boundary layers on adiabatic slopes

The presence of a no-slip, impermeable, adiabatic, sloped boundary in an otherwise quiescent, stably stratified, Boussinesq flow generates baroclinic vorticity within a diffusive boundary layer. Such conditions are typical of the oscillating boundary layers on adiabatic abyssal slopes, sloped lake bathymetry and sloped coastal bathymetry in the absence of high-wavenumber internal waves, mean flows, far-field turbulence on larger scales, and resonant tidal–bathymetric interaction. We investigate the linear stability of the oscillating flow within non-dimensional parameter space typical of the $M_2$ tide and hydraulically smooth, middle-latitude abyssal slopes through Floquet linear stability analysis. The flow dynamics depends on three non-dimensional variables: the Reynolds number for Stokes’ second problem, the Prandtl number, and a frequency ratio that accounts for the resonance conditions ($C$, criticality) of the buoyant restoring force and the tidal forcing. The Floquet analysis results suggest that oscillating laminar boundary layers on adiabatic abyssal slopes are increasingly unstable as Reynolds number, criticality parameter and/or spanwise disturbance wavenumber are increased. We also show that the two-dimensional Floquet linear instability necessarily generates three-dimensional baroclinic vorticity, which suggests that the evolution of the gravitational instabilities may be nonlinear as $t→∞$.

42 ENGINEERING↗
DOE OSTI2021

A first principles framework to predict the transient performance of latent heat thermal energy storage

Thermal energy storage (TES) is increasingly recognized as an essential component of efficient Combined Heat and Power (CHP), Concentrated Solar Power (CSP), Heating Ventilation and Air Conditioning (HVAC), and refrigeration as it reduces peak demand while helping to manage intermittent availability of energy (e.g., from solar or wind). Latent Heat Thermal Energy Storage (LHTES) is a viable option because of its high energy storage density. Parametric analysis of LHTES in terms of dimensionless numbers is highly desired as a tool to model LHTES systems. One approach is to develop a model equation so as to minimize the error between the model and data obtained from experiments or simulations. While this approach can produce an accurate correlation applicable within the range of data used for its creation, it does not provide physical understanding of the rate-limiting process controlling the transient behavior of the device. In this paper we present an alternative approach whereby the potential rate-limiting processes are identified from first principles and then the key process is determined as a function of time as a LHTES device is charged. For example, in a simple geometry, the melt-fraction can be expected to vary linearly in time if the heat transfer rate is limited by natural convection of the phase changing material and we show it scales with the PCM Grashof number as $Gr^1_p$ and PCM Prandtl number as $Pr_p^{(1/3)}$. On the other hand, if surface area of solid PCM limits the heat transfer rate, the melt fraction increases asymptotically to reach full melting. The existence of these linear and asymptotic regions and the $Gr^1_pP r^{1/3}_p$ shape of the melt fraction curve is verified using our database of 64 simulations. Of practical importance in designing LHTES devices is the melt fraction at which the heat transfer rate ceases to be limited by convection, after which the heat storage rate deteriorates. For our geometry, this is found to be about 90%. This test case of our methodology shows the value of our approach, that predicting heat storage rate based on the rate-limiting physical phenomenon as a function of time is an effective approach to modeling LHTES devices.

25 ENERGY STORAGE↗
DOE OSTI2025

Taking control of compressible modes: bulk viscosity and the turbulent dynamo

Many polyatomic astrophysical plasmas are compressible and out of chemical and thermal equilibrium, introducing a bulk viscosity into the plasma via the internal degrees of freedom of the molecular composition, directly impacting the decay of compressible modes, $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$. This is especially important for small-scale, turbulent dynamo processes in the interstellar medium (ISM), which are known to be sensitive to the effects of compression. To control the viscous properties of $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$, we perform trans-sonic, visco-resistive dynamo simulations with additional bulk viscosity $\nu _{\text{bulk}}$, deriving a new $\nu _{\text{bulk}}$ Reynolds number $\text{Re}_{\text{bulk}}$, and viscous Prandtl number $\text{P}\nu \equiv \text{Re}_{\text{bulk}}/ \text{Re}_{\text{shear}}$, where $\text{Re}_{\text{shear}}$ is the shear viscosity Reynolds number. We derive a framework for decomposing $E_{\rm mag}$ growth rates into incompressible and compressible terms via orthogonal tensor decompositions of $\boldsymbol {\nabla }\otimes \mathrm{{\boldsymbol {\mathit {v}}}}$, where $\mathrm{{\boldsymbol {\mathit {v}}}}$ is the fluid velocity. We find that $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$ play a dual role, growing and decaying $E_{\rm mag}$, and that field-line stretching is the main driver of growth, even in compressible dynamos. In the absence of $\nu _{\text{bulk}}$ ($\text{P}\nu \rightarrow \infty$), $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$ pile up on small-scales, creating a spectral bottleneck, which disappears for $\text{P}\nu \approx 1$. As $\text{P}\nu$ decreases, $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$ are dissipated at increasingly larger scales, in turn suppressing incompressible modes through a coupling between high-k modes. We emphasize the importance of further understanding the role of $\nu _{\text{bulk}}$ in compressible astrophysical plasmas, which we estimate could be as strong as the shear viscosity in the cold ISM, and highlight that compressible direct numerical simulations without bulk viscosity have unresolved compressible mode dissipation scales.

MHD↗
DOE OSTI2024

Wall modes and the transition to bulk convection in rotating Rayleigh-Bénard convection

We investigate states of rapidly rotating Rayleigh-Bénard convection in a cylindrical cell over a range of Rayleigh numbers 3 × 10 5 ≤ Ra ≤ 5 × 10 9 and Ekman numbers 10 − 6 ≤ Ek ≤ 10 − 4 for Prandtl number Pr = 0.8 and aspect ratios 1 / 5 ≤ Γ ≤ 5 using direct numerical simulations. We characterize, for perfectly insulating sidewall boundary conditions, the first transition to convection via wall mode instability and the nonlinear growth and instability of the resulting wall mode states, including a secondary transition to time dependence. We show how the radial structure of the vertical velocity u z and the temperature T is captured well by the linear eigenfunctions of the wall mode instability where the radial width of u z is δ u z ∼ Ek 1 / 3 r / H whereas δ T ∼ e − k r ( k is the wave number of a laterally infinite wall mode state). The disparity in spatial scales for Ek = 10 − 6 means that the heat transport is dominated by the radial structure of u z since T varies slowly over the radial scale δ u z . We further describe how the transition to a state of bulk convection is influenced by the presence of the wall mode states. We use temporal and spatial scales as measures of the local state of convection and the Nusselt number Nu as representative of global transport. Our results elucidate the evolution of the wall state of rotating convection and confirm that wall modes are strongly linked with the boundary zonal flow being the robust remnant of nonlinear wall mode states. We also show how the heat transport ( Nu ) contributions of wall modes and bulk modes are related and discuss approaches to disentangling their relative contributions. Published by the American Physical Society 2024

58 GEOSCIENCES↗
DOE OSTI2020

Measurement of the Natural Convection Heat Transfer in a Magnesium Oxide Electrolytic Cell Concept

The rate of heat transfer by natural convection between the wall and electrolyte of an electrolytic cell that produces magnesium (Mg) from magnesium oxide (MgO) at temperatures near 1000 °C in a molten fluoride salt electrolyte is presented. An experimental model of the cell was developed that enabled measurements of the heat transfer in the absence of electrolysis and at temperatures less than 100 °C over ranges of Rayleigh numbers from 1 × 10 -7 to 7 × 10 -8 and Prandtl numbers from 2 to 6200, ranges that include those anticipated in the operation of the MgO electrolytic cell. The model avoids the substantial experimental challenges associated with the high-temperature, corrosive molten salt to enable a conservative estimate of the heat transfer at a lower cost and greater accuracy than would otherwise be possible. The results are correlated by the expression Nu = 0.412Ra 0.23 Pr 0.15 with Nusselt numbers spanning 30–80. The application of the correlation shows that the heat transfer between the cell wall and the molten fluoride electrolyte at ≈1000 °C is characterized by convection coefficients between 100 and 600 W/m 2 -K and is fast enough to enable heat fluxes up to 10 W/cm 2 without compromising the structural integrity of the steel cell wall.

Engineering↗
DOE OSTI2022

Heat transport and temperature boundary-layer profiles in closed turbulent Rayleigh–Bénard convection with slippery conducting surfaces

Here, we report direct numerical simulations (DNS) of the Nusselt number Nu , the vertical profiles of mean temperature Θ(z) and temperature variance Ω(z) across the thermal boundary layer (BL) in closed turbulent Rayleigh–Bénard convection (RBC) with slippery conducting surfaces ( z is the vertical distance from the bottom surface). The DNS study was conducted in three RBC samples: a three-dimensional cuboid with length L=H and width W=H/4 ( H is the sample height), and two-dimensional rectangles with aspect ratios Γ≡L/H=1 and 10. The slip length b for top and bottom plates varied from 0 to ∞ . The Rayleigh numbers Ra were in the range 10 6 $&LessSlantEqual;$Ra$&LessSlantEqual;$10 10 and the Prandtl number Pr was fixed at 4.3 . As b increases, the normalised Nu/Nu 0 ( Nu 0 is the global heat transport for b =0 ) from the three samples for different Ra and Γ can be well described by the same function Nu/Nu 0 = N 0 tanh( b /λ 0 )+1 , with N 0 =0.8±0.03 . Here λ 0 ≡ L /(2 Nu 0 ) is the thermal boundary layer thickness for b =0 . Considering the BL fluctuations for Pr >1, one can derive solutions of temperature profiles Θ(z) and Ω(z) near the thermal BL for b $&GreaterSlantEqual;$ 0 . When b =0 , the solutions are equivalent to those reported by Shishkina et al. ( Phys. Rev. Lett. , vol. 114, 2015, 114302) and Wang et al. ( Phys. Rev. Fluids , vol. 1, 2016, 082301(R)), respectively, for no-slip plates. For b >0, the derived solutions are in excellent agreement with our DNS data for slippery plates.

42 ENGINEERING↗
DOE OSTI2025

Thermal Performance of Triply Periodic Minimal Surface Lattice Structures in Single-Phase Dielectric Fluid Cooling of Power Electronics

Additive manufacturing has transformed thermal management by enabling the production of complex, optimized geometries that conventional manufacturing methods cannot achieve. This study investigates the single-phase convective heat transfer performance of gyroid triply periodic minimal surface (TPMS) lattice structures with functional porosity. TPMS structures provide high surface area to volume ratios and are amenable to 3D printing. A gyroid numerical model was created and validated against an existing experimental study with a similar feature size to the investigated geometries. The TPMS structure has a periodic width of 1.6 mm, a length of 10 mm, and a height of 4 mm, with a functional porosity ranging from 0.5 to 0.8, decreasing with distance from the heated surface. Three different flow configurations were examined for an inlet fluid temperature of 70 °C. The inlet velocities range from 0.01 to 1.2 m/s, corresponding to a Reynolds number range of 10–900 with a heat flux of 50 W/cm 2 applied at the base. AmpCool ® AC-110 dielectric fluid (Prandtl number 59.5) was used as the coolant. Thermal performance and friction characteristics were studied for the three flow orientations. The parallel flow configuration was identified as the most efficient for heat removal. A detailed analysis of the numerical results highlights the underlying physics behind the thermal performance differences among the flow configurations.

33 ADVANCED PROPULSION SYSTEMS↗