SEARCH · Search NASA
Results for “Low and high Prandtl fluids”
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.
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.
Taylor–Couette flow for astrophysical purposes
A concise review is given of astrophysically motivated experimental and theoretical research on Taylor–Couette flow. The flows of interest rotate differentially with the inner cylinder faster than the outer, but are linearly stable against Rayleigh’s inviscid centrifugal instability. At shear Reynolds numbers as large as 10 6 , hydrodynamic flows of this type (quasi-Keplerian) appear to be nonlinearly stable: no turbulence is seen that cannot be attributed to interaction with the axial boundaries, rather than the radial shear itself. Direct numerical simulations agree, although they cannot yet reach such high Reynolds numbers. This result indicates that accretion-disc turbulence is not purely hydrodynamic in origin, at least insofar as it is driven by radial shear. Theory, however, predicts linear magnetohydrodynamic (MHD) instabilities in astrophysical discs: in particular, the standard magnetorotational instability (SMRI). MHD Taylor–Couette experiments aimed at SMRI are challenged by the low magnetic Prandtl numbers of liquid metals. High fluid Reynolds numbers and careful control of the axial boundaries are required. The quest for laboratory SMRI has been rewarded with the discovery of some interesting inductionless cousins of SMRI, and with the recently reported success in demonstrating SMRI itself using conducting axial boundaries. Some outstanding questions and near-future prospects are discussed, especially in connection with astrophysics.
Challenge Problem 1: Preliminary Model Development and Assessment of Flexible Heat Transfer Modeling Approaches
This report presents the modeling progress of a university consortium studying flexible modeling for turbulent heat transfer. In the context of advanced nuclear reactor design, where challenges include non-unity Prandtl fluids, free and mixed convection, and transitional flow, computational fluid dynamics methods are underutilized due to low confidence in modeling approaches and high computational cost. This report evaluates the performance of engineering CFD models in diabatic, buoyant turbulent flow. It finds that all turbulence models including several variants of the k-ε and k-ω models struggle to predict accurate turbulent momentum and heat transfer in such flows. The Nusselt numbers have been compared between the models and the DNS data, where calculations have been performed for each case and trends have shown a good agreement between DNS estimated Nusselt numbers and available correlations as well as experimental data. A novel DNS correlation for high Pr cases as they are transitioning from mixed convection to natural convection has been proposed. While Nusselt number errors relative to DNS range from 20% to 50%, the models capture similar trends to DNS with respect to Nusselt suppression and amplification under varying levels of buoyancy effect. This report also highlights model form error as a significant contributor to CFD predictions and proposes a framework for quantifying model error and improving confidence in CFD calculations. To improve the predictive capability of engineering CFD models, data-driven approaches for turbulence models are investigated. Theoretical frameworks based on the invariant tensor / vector basis neural networks for prediction of Reynolds stresses and turbulent heat fluxes are employed. The models are developed using direct numerical simulations data for forced convection flows of different fluids in vertical planar channel domain. The framework is implemented in spectral element solvers Nek5000 / nekRS and has shown a potential for future development and consideration of mixed convection flows.
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.
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.
Unstable turbulent channel flow response to spanwise-heterogeneous heat fluxes: Prandtl's secondary flow of the third kind
Turbulent secondary flows are defined as Prandtl's secondary flow of the first or second kind, the former produced by stretching and/or tilting of vorticity, the latter produced via spatial heterogeneity of Reynolds stresses. Both mechanisms are instantaneously active within inertia-dominated wall turbulence; Reynolds stress spatial heterogeneity is required for Reynolds-averaged secondary flows. Spanwise-variable surface roughness can induce turbulent stress spatial heterogeneity in the spanwise–wall-normal plane and provide sustenance for streamwise-aligned mean secondary flows. Herein, we demonstrate that turbulent secondary flows can also be sustained by spanwise variability in the surface heat flux in unstably stratified turbulent channels, defined hereafter as Prandtl's secondary flow of the third kind. Support for this mechanism is established with scaling arguments, while large-eddy simulation is used to model inertia-dominated channel turbulence responding to a lower boundary with uniform aerodynamic/hydrodynamic roughness but spanwise-variable surface heat flux. Transport equations for streamwise vorticity and turbulent kinetic energy, k, outline the conditions needed for third-kind production: shear and buoyancy production over the elevated heat flux regions necessitates lateral entrainment of low- k fluid, inducing mean counter-rotating secondary cells aligned such that upwelling and downwelling occur over the high and low heat flux regions, respectively. Buoyancy-driven production of k alters aggregate flow response and thus is a distinctly different mechanism responsible for sustenance of secondary flows.
Numerical Simulation and Experimental Comparison of System Analysis Module 1D Mixing Model for Cold Shock Transients in the Gallium Thermal-Hydraulic Mixing Facility
Abstract Liquid metals are being investigated as coolants in many advanced reactor designs because of their high thermal conductivity and effectiveness at high temperatures. However, they often pose challenges to reactor operation and safety because of the complex thermal mixing and stratification in the plenum of pool-type reactor designs. The advanced system analysis code System Analysis Module (SAM) currently under development at Argonne National Laboratory aims to develop and implement thermal mixing models to accurately capture these complex thermal fluid behaviors. In this study, the SAM thermal mixing model was compared against experimental data from the Gallium Thermal-Hydraulic Experiment facility, a scaled liquid metal test facility that uses gallium as a surrogate fluid to investigate the stratification and thermal mixing of low-Prandtl-number fluids in the upper plenum of a liquid metal-cooled reactor. Two cold shock transient cases were used: one with stable stratified flow (Ri = 32) and one with stronger thermal mixing (Ri = 0.5). The resultant temperatures were then compared with the experimental temperatures over the entire plenum to assess the ability of the mixing models to capture the thermal behavior and to better correspond mixing parameters to various flow scenarios. Generally, the zero-dimensional mixing model was more capable of capturing the bulk temperature of the component modeled assuming that an accurate mass flow rate was provided, but it was inherently unable to capture thermal gradients in space. The one-dimensional mixing model was capable of capturing that the thermal gradients provided accurate selection of the mixing coefficients. Further, the temperature at the outlet junction was compared over time for each of the mixing models with the recorded experimental temperature. The implemented mixing models demonstrated the ability to effectively capture the overall thermal behavior for stronger mixing scenarios but struggled with more stably stratified flows. It was found that a system analysis code's covering of the entire range of different operating conditions still remains a challenging task, and it is suggested that further model and closure improvements are necessary to accurately capture complex thermal mixing and stratification phenomena.
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.