Search NASA⌕ Search

SEARCH · Search NASA

Results for “PRANDTL NUMBER”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Explicit finite-volume time-marching calculations of total temperature distributions in turbulent flow

A method was developed which calculates two-dimensional, transonic, viscous flow in ducts. The finite volume, time-marching formulation is used to obtain steady flow solutions of the Reynolds-averaged form of the Navier-Stokes equations. The entire calculation is performed in the physical domain. This paper investigates the introduction of a new formulation of the energy equation which gives improved transient behavior as the calculation converges. The effect of variable Prandtl number on the temperature distribution through the boundary layer is also investigated. A turbulent boundary layer in an adverse pressure gradient (M = 0.55) is used to demonstrate the improved transient temperature distribution obtained when the new formulation of the energy equation is used. A flat plate turbulent boundary layer with a supersonic free-stream Mach number of 2.8 is used to investigate the effect of Prandtl number on the distribution of properties through the boundary layer. The computed total temperature distribution and recovery factor agree well with the measurements when a variable Prandtl number is used through the boundary layer.

Nicholson, Stephen↗
DOE OSTI2017

Turbulent Transport by Diffusive Stratified Shear Flows: From Local to Global Models. I. Numerical Simulations of a Stratified Plane Couette Flow

Shear-induced turbulence could play a significant role in mixing momentum and chemical species in stellar radiation zones, as discussed by Zahn. In this paper we analyze the results of direct numerical simulations of stratified plane Couette flows, in the limit of rapid thermal diffusion, to measure the turbulent viscosity and the turbulent diffusivity of a passive tracer as a function of the local shear and the local stratification. We find that the stability criterion proposed by Zahn, namely that the product of the gradient Richardson number and the Prandtl number must be smaller than a critical values (JPr){sub c} for instability, adequately accounts for the transition to turbulence in the flow, with (JPr){sub c}≃0.007. This result recovers and confirms the prior findings of Prat et al. Zahn’s model for the turbulent diffusivity and viscosity, namely that the mixing coefficient should be proportional to the ratio of the thermal diffusivity to the gradient Richardson number, does not satisfactorily match our numerical data. It fails (as expected) in the limit of large stratification where the Richardson number exceeds the aforementioned threshold for instability, but it also fails in the limit of low stratification where the turbulent eddy scale becomes limited by the computational domain size. We propose a revised model for turbulent mixing by diffusive stratified shear instabilities that properly accounts for both limits, fits our data satisfactorily, and recovers Zahn’s model in the limit of large Reynolds numbers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Radiation effects on bifurcation and dual solutions in transient natural convection in a horizontal annulus

Transitions and bifurcations of transient natural convection in a horizontal annulus with radiatively participating medium are numerically investigated using the coupled lattice Boltzmann and direct collocation meshless (LB-DCM) method. As a hybrid approach based on a common multi-scale Boltzmann-type model, the LB-DCM scheme is easy to implement and has an excellent flexibility in dealing with the irregular geometries. Separate particle distribution functions in the LBM are used to calculate the density field, the velocity field and the thermal field. In the radiatively participating medium, the contribution of thermal radiation to natural convection must be taken into account, and it is considered as a radiative term in the energy equation that is solved by the meshless method with moving least-squares (MLS) approximation. The occurrence of various instabilities and bifurcative phenomena is analyzed for different Rayleigh number Ra and Prandtl number Pr with and without radiation. Then, bifurcation diagrams and dual solutions are presented for relevant radiative parameters, such as convection-radiation parameter Rc and optical thickness τ. Numerical results show that the presence of volumetric radiation changes the static temperature gradient of the fluid, and generally results in an increase in the flow critical value. Besides, the existence and development of dual solutions of transient convection in the presence of radiation are greatly affected by radiative parameters. Finally, the advantage of LB-DCM combination is discussed, and the potential benefits of applying the LB-DCM method to multi-field coupling problems are demonstrated.

42 ENGINEERING↗
DOE OSTI2016

Sensitivity of Turbine-Height Wind Speeds to Parameters in Planetary Boundary-Layer and Surface-Layer Schemes in the Weather Research and Forecasting Model

We evaluate the sensitivity of simulated turbine-height winds to 26 parameters applied in a planetary boundary layer (PBL) scheme and a surface layer scheme of the Weather Research and Forecasting (WRF) model over an area of complex terrain during the Columbia Basin Wind Energy Study. An efficient sampling algorithm and a generalized linear model are used to explore the multiple-dimensional parameter space and quantify the parametric sensitivity of modeled turbine-height winds. The results indicate that most of the variability in the ensemble simulations is contributed by parameters related to the dissipation of the turbulence kinetic energy (TKE), Prandtl number, turbulence length scales, surface roughness, and the von Kármán constant. The relative contributions of individual parameters are found to be dependent on both the terrain slope and atmospheric stability. The parameter associated with the TKE dissipation rate is found to be the most important one, and a larger dissipation rate can produce larger hub-height winds. A larger Prandtl number results in weaker nighttime winds. Increasing surface roughness reduces the frequencies of both extremely weak and strong winds, implying a reduction in the variability of the wind speed. All of the above parameters can significantly affect the vertical profiles of wind speed, the altitude of the low-level jet and the magnitude of the wind shear strength. The wind direction is found to be modulated by the same subset of influential parameters. Remainder of abstract is in attachment.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1996

The structure of variable property, compressible mixing layers in binary gas mixtures

We present the results of a study of the structure of a parallel compressible mixing layer in a binary mixture of gases. The gases included in this study are hydrogen (H2), helium (He), nitrogen (N2), oxygen (02), neon (Ne) and argon (Ar). Profiles of the variation of the Lewis and Prandtl numbers across the mixing layer for all thirty combinations of gases are given. It is shown that the Lewis number can vary by as much as a factor of eight and the Prandtl number by a factor of two across the mixing layer. Thus assuming constant values for the Lewis and Prandtl numbers of a binary gas mixture in the shear layer, as is done in many theoretical studies, is a poor approximation. We also present profiles of the velocity, mass fraction, temperature and density for representative binary gas mixtures at zero and supersonic Mach numbers. We show that the shape of these profiles is strongly dependent on which gases are in the mixture as well as on whether the denser gas is in the fast stream or the slow stream.

Kozusko, F.↗
NASA NTRS1974

Relaxation of an unsteady turbulent boundary layer on a flat plate in an expansion tube

An analysis is presented for the relaxation of a turbulent boundary layer on a semi-infinite flat plate after passage of a shock wave and a trailing driver gas-driven gas interface. The problem has special application to expansion-tube flows. The flow-governing equations have been transformed into the Crocco variables, and a time-similar solution is presented in terms of the dimensionless distance-time variable alpha and the dimensionless velocity variable beta. An eddy-viscosity model, similar to that of time-steady boundary layers, is applied to the inner and outer regions of the boundary layer. A turbulent Prandtl number equal to the molecular Prandtl number is used to relate the turbulent heat flux to the eddy viscosity. The numerical results, obtained by using the Gauss-Seidel line-relaxation method, indicate that a fully turbulent boundary layer relaxes faster to the final steady-state values of heat transfer and skin friction than a laminar boundary layer. The results also give a fairly good estimate of the local skin friction and heat transfer for near steady-flow conditions.

Gurta, R. N.↗
NASA NTRS1993

Supersonic flow calculation using a Reynolds-stress and an eddy thermal diffusivity turbulence model

A second-order model for the velocity field and a two-equation model for the temperature field are used to calculate supersonic boundary layers assuming negligible real gas effects. The modeled equations are formulated on the basis of an incompressible assumption and then extended to supersonic flows by invoking Morkovin's hypothesis, which proposes that compressibility effects are completely accounted for by mean density variations alone. In order to calculate the near-wall flow accurately, correction functions are proposed to render the modeled equations asymptotically consistent with the behavior of the exact equations near a wall and, at the same time, display the proper dependence on the molecular Prandtl number. Thus formulated, the near-wall second order turbulence model for heat transfer is applicable to supersonic flows with different Prandtl numbers. The model is validated against flows with different Prandtl numbers and supersonic flows with free-stream Mach numbers as high as 10 and wall temperature ratios as low as 0.3. Among the flow cases considered, the momentum thickness Reynolds number varies from approximately 4,000 to approximately 21,000. Good correlation with measurements of mean velocity, temperature, and its variance is obtained. Discernible improvements in the law-of-the-wall are observed, especially in the range where the big-law applies.

Sommer, T. P.↗
DOE OSTI2011

STABILITY OF MRI-TURBULENT ACCRETION DISKS

Based on the characteristics of the magnetorotational instability (MRI) and the MRI-driven turbulence, we construct a steady model for a geometrically thin disk using 'non-standard' {alpha}-prescription. The efficiency of the angular momentum transport depends on the magnetic Prandtl number, Pm = {nu}/{eta}, where {nu} and {eta} are the microscopic viscous and magnetic diffusivities. In our disk model, Shakura-Sunyaev's {alpha}-parameter has a power-law dependence on the magnetic Prandtl number, that is {alpha} {proportional_to} Pm{sup {delta}}, where {delta} is the constant power-law index. Adopting Spitzer's microscopic diffusivities, the magnetic Prandtl number becomes a decreasing function of the disk radius when {delta}>0. The transport efficiency of the angular momentum and the viscous heating rate are thus smaller in the outer part of the disk, while these are impacted by the size of index {delta}. We find that the disk becomes more unstable to the gravitational instability for a larger value of index {delta}. The most remarkable feature of our disk model is that the thermal and secular instabilities can grow in its middle part even if the radiation pressure is negligibly small in the condition {delta}>2/3. In the realistic disk systems, it would be difficult to maintain the steady mass accretion state unless the Pm dependence of MRI-driven turbulence is relatively weak.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

Numerical simulation of supercritical heat transfer under severe axial density gradient in a narrow vertical tube

A number of computational works have been performed so far for the simulation of heat transfer in a supercritical fluid. The simulations, however, faced a lot of difficulties when heat transfer deteriorates due either to buoyancy or by acceleration. When the bulk temperature approaches the pseudo-critical temperature the fluid experiences a severe axial density gradient on top of a severe radial one. Earlier numerical calculations showed, without exception, unrealistic over-predictions, as soon as the bulk temperature exceeded the pseudo-critical temperature. The over-predictions might have been resulted from an inapplicability of widely-used turbulence models. One of the major causes for the difficulties may probably be an assumption of a constant turbulent Prandtl number. Recent research, both numerical and experimental, indicates that the turbulent Prandtl number is never a constant when the gradient of physical properties is significant. This paper describes the applicability of a variable turbulent Prandtl number to the numerical simulation of heat transfer in supercritical fluids flowing in narrow vertical tubes. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2015

Convective heat transfer in a high aspect ratio minichannel heated on one side

Experimental results are presented for single-phase heat transfer in a narrow rectangular minichannel heated on one side. The aspect ratio and gap thickness of the test channel were 29:1 and 1.96 mm, respectively. Friction pressure drop and Nusselt numbers are reported for the transition and fully turbulent flow regimes, with Prandtl numbers ranging from 2.2 to 5.4. Turbulent friction pressure drop for the high aspect ratio channel is well-correlated by the Blasius solution when a modified Reynolds number, based upon a laminar equivalent diameter, is utilized. The critical Reynolds number for the channel falls between 3500 and 4000, with Nusselt numbers in the transition regime being reasonably predicted by Gnielinski's correlation. The dependence of the heat transfer coefficient on the Prandtl number is larger than that predicted by circular tube correlations, and is likely a result of the asymmetric heating. The problem of asymmetric heating condition is approached theoretically using a boundary layer analysis with a two-region wall layer model, similar to that originally proposed by Prandtl. The analysis clarifies the influence of asymmetric heating on the Nusselt number and correctly predicts the experimentally observed trend with Prandtl number. Furthermore, a semi-analytic correlation is derived from the analysis that accounts for the effect of aspect ratio and asymmetric heating, and is shown to predict the experimental results of this study with a mean absolute error (MAE) of less than 5% for 4000 < Re < 70,000.

42 ENGINEERING↗
DOE OSTI2024

Thermal field reconstruction and compressive sensing using proper orthogonal decomposition

Model order reduction allows critical information about sensor placement and experiment design to be distilled from raw fluid mechanics simulation data. In many cases, sensed information in conjunction with reduced order models can also be used to regenerate full field variables. In this paper, a proper orthogonal decomposition (POD) inferencing method is extended to the modeling and compressive sensing of temperature, a scalar field variable. The method is applied to a simulated, critically stable, incompressible flow over a heated cylinder (Re = 1000) with Prandtl number varying between 0.001 and 50. The model is trained on pressure and temperature data from simulations. Field reconstructions are then generated using data from selected sensors and the POD model. Finally, the reconstruction error is evaluated across all Prandtl numbers for different numbers of retained modes and sensors. The predicted trend of increasing reconstruction accuracy with decreasing Prandtl number is confirmed and a Prandtl number/sensor count error matrix is presented.

Matulis, John↗
NASA NTRS1993

Mixing in a stratified shear flow: Energetics and sampling

Direct numerical simulations of the time evolution of homogeneous stably stratified shear flows have been performed for Richardson numbers from 0 to 1 and for Prandtl numbers between 0.1 and 2. The results indicate that mixing efficiency R(sub f) varies with turbulent Froude number in a manner consistent with laboratory experiments performed with Prandtl numbers of 0.7 and 700. However, unlike the laboratory results, for a particular Froude number, the simulations do not show a clear dependence on the magnitude of R(sub f) on Pr. The observed maximum value of R(sub f) is 0.25. When averaged over vertical length scales of an order of magnitude greater than either the overturning or Ozmidov scales of the flow, the simulations indicate that the dissipation rate epsilon is only weakly lognormally distributed with an intermittency of about 0.01 whereas estimated values in the ocean are 3 to 7.

Ivey, G. N.↗
NASA NTRS2024

Turbulence Model Validation through Joint Experimental / Computational Studies of Separated Flow Over a Three-Dimensional Tapered Bump: Part II - Computational Investigation

This report serves as documentation of activities related to Computational Fluid Dynamics (CFD) analysis applied to the Smooth Body Separation Experiment designed to produce validation-quality experimental data from wind-tunnel tests performed at the University of Notre Dame. All CFD analysis utilized the Boeing CFD (BCFD) flow solver. BCFD was used to analyze flow over the Boeing Speed Bump with second-order accuracy in space and time. Inviscid fluxes were computed using the HLLE scheme for steady-state simulations, and the Low-Dissipation Roe scheme for time-accurate simulations. The working fluid was modeled as calorically perfect air using Sutherland's law for viscosity. Thermal diffusivity was determined by assuming a constant molecular Prandtl number of 0.72, and a turbulent Prandtl number of 0.9. Reynolds Averaged Navier-Stokes (RANS) simulations were conducted using the one-equation Spalart-Allmaras turbulence model with rotation correction3 and the quadratic constitutive relation (SA-RC-QCR), as well as Menter's two-equation, Shear-Stress Transport (SST) model. Time-accurate computations were conducted using Delayed Detached eddy Simulation (DDES) with the one-equation RANS model, and Large Eddy Stress balance (LESb) with the two-equation RANS model. Implicit dual-time stepping was used to achieve second-order temporal accuracy. The experimental investigation conducted during the course of this study is described in detail in an accompanying report.

Flow separation↗
DOE OSTI2025

Simulations of Heat Transfer Using Tight-Fitting Twisted Tape Inserts for First Wall Cooling in Molten Salt Breeder Blankets

One of the major components in fusion energy systems is the fusion blanket, which has a vacuum vessel to contain the plasma. As part of the fusion blanket/vacuum vessel, the first wall and plasma-facing components require sufficient cooling to prevent material degradation during operation from the superheated plasma. Most fusion blanket concepts involve first wall and divertor coolant channels with heat transfer enhancements (HTEs) that are intended to withstand the incident high heat fluxes of 1 to 5 MW/m 2 . Twisted tape inserts are a proposed HTE that have been investigated previously for first wall cooling and monoblock divertor cooling channels and in other nonfusion heat transfer components. By inserting twisted tapes into straight pipes, the amount of turbulence in the system can be increased at lower Reynolds numbers by swirling the flow. This results in better heat transfer characteristics with marginal increases in frictional pressure losses. In particular, simulations of high-Prandtl-number fluids such as the proposed molten salt FLiBe in twisted tapes, which is prototypic to liquid immersion blankets, have not been previously explored. Here, in this study, we simulate various Prandtl numbers in pipes with twisted tape inserts using large eddy simulations to determine the effects of increasing Prandtl numbers on heat transfer performance. The quantities of particular interest are the Nusselt number and the friction factor, which were recovered using data reduction techniques to determine impacts on heat transfer and pressure losses. This work serves as a starting point for determining the feasibility of twisted tape inserts for liquid immersion blanket concepts.

LES↗
DOE OSTI2019

A General Study of Counterflow Diffusion Flames for Supercritical CO2 Combustion

Abstract A counterflow diffusion flame for supercritical CO2 combustion is investigated at various CO2 dilution levels and pressures by accounting for real gas effects into both thermal and transport properties. The UCF 1.1 24-species mechanism is used to account the chemistry. The nature of important nonpremixed combustion characteristics such as Prandtl number, thermal diffusivity, Lewis number, stoichiometric scalar dissipation rate, flame thickness, and Damköhler number are investigated with respect to CO2 dilution and pressure. The results show that the aforementioned parameters are influenced by both dilution and pressure; the dilution effect is more dominant. Further, the result shows that Prandtl number increases with CO2 dilution and at 90% CO2 dilution, the difference between the Prandtl number of the inlet jets and the flame is minimal. Also, the common assumption of unity Lewis number in the theory and modeling of nonpremixed combustion does not hold reasonable for sCO2 applications due to large difference of Lewis number across the flame and the Lewis number on the flame drop significantly with an increase in the CO2 dilution. An interesting relation between Lewis number and CO2 dilution is observed. The Lewis number of species drops by 15% when increasing the CO2 dilution by 30%. Increasing the CO2 dilution increases both the flow and chemical timescales; however, chemical timescale increases faster than the flow timescales. The magnitudes of the Damköhler number signify the need to consider finite rate chemistry for sCO2 applications. Further, the Damköhler numbers at 90% sCO2 dilution are very small; hence, laminar flamelet assumptions in turbulent combustion simulations are not physically correct for this application. Also, it is observed that the Damköhler number drops nonlinearly with increasing CO2 dilution in the oxidizer stream. This is a very important observation for the operation of sCO2 combustors. Further, the flame thickness is found to increase with CO2 dilution and reduce with pressure.

Engineering↗
NASA NTRS1961

Experimental Determination of the Recovery Factor and Analytical Solution of the Conical Flow Field for a 20 deg Included Angle Cone at Mach Numbers of 4.6 and 6.0 and Stagnation Temperatures to 2600 degree R

The local recovery factor was determined experimentally along the surface of a thin-walled 20 deg included angle cone for Mach numbers near 6.0 at stagnation temperatures between 1200 deg R and 2600 deg R. In addition, a similar cone configuration was tested at Mach numbers near 4.5 at stagnation temperatures of approximately 612 deg R. The local Reynolds number based on flow properties at the edge of the boundary layer ranged between 0.1 x 10(exp 4) and 3.5 x 10(exp 4) for tests at temperatures above 1200 deg R and between 6 x 10(exp 4) and 25 x 10(exp 4) for tests at temperatures near 612 deg R. The results indicated, generally, that the recovery factor can be predicted satisfactorily using the square root of the Prandtl number. No conclusion could be made as to the necessity of evaluating the Prandtl number at a reference temperature given by an empirical equation, as opposed to evaluating the Prandtl number at the wall temperature or static temperature of the gas at the cone surface. For the tests at temperatures above 1200 deg R (indicated herein as the tests conducted in the slip-flow region), two definite trends in the recovery data were observed - one of increasing recovery factor with decreasing stagnation pressure, which was associated with slip-flow effects and one of decreasing recovery factor with increasing temperature. The true cause of the latter trend could not be ascertained, but it was shown that this trend was not appreciably altered by the sources of error of the magnitude considered herein. The real-gas equations of state were used to determine accurately the local stream properties at the outer edge of the boundary layer of the cone. Included in the report, therefore, is a general solution for the conical flow of a real gas using the Beattie-Bridgeman equation of state. The largest effect of temperature was seen to be in the terms which were dependent upon the internal energy of the gas. The pressure and hence the pressure drag terms were unaffected.

Pfyl, Frank A.↗
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↗
NASA NTRS1985

Kolmogorov and scalar spectral regimes in numerical turbulence

Velocity and passive-scalar spectra for turbulent fields generated by a forced three-dimensional simulation and Taylormicroscale Reynolds numbers up to 83 are shown to have distinct spectral regimes, including a Kolmogorov inertial subrange. Both one- and three-dimensional spectra are shown for comparison with experiment and theory, respectively. When normalized by the Kolmogorov dissipation scales velocity spectra collapse to a single curve and a high-wavenumber bulge is seen. The bulge leads to an artificially high Kolmogorov constant, but is consistent with recent measurements of the velocity spectrum in the dissipation regime and the velocity-derivative skewness. Scalar spectra, when normalized by the Oboukov-Corrsin scales, collapse to curves which depend only on Prandtl number and show a universal inertial-convective subrange, independent of Prandtl number. When normalized by the Batchelor scales, the scalar spectra show a universal dissipation regime which is independent of Prandtl numbers from 0.1 to 1.0. The time development of velocity spectra is illustrated by energy-transfer spectra in which distinct pulses propagate to high wavenumbers.

Kerr, R. M.↗