SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS1972

Turbulent properties of a compressible boundary layer.

Turbulent shear-stress, eddy-viscosity, mixing-length, heat-flux and Prandtl number distributions across a hypersonic turbulent boundary layer have been determined from the 'time-averaged' conservation equations using experimental mean profile data obtained at several streamwise locations in a fully developed turbulent boundary layer with negligible pressure gradient. The eddy-viscosity, mixing-length and Prandtl number results show general agreement with previous incompressible and adiabatic compressible correlations. However, when the turbulent Prandtl number is defined using total enthalpy as opposed to static enthalpy no clear correlation of the results in the outer portion of the boundary layer could be obtained.

Horstman, C. C.↗

NASA NTRS1982

Symmetric baroclinic instability of a Hadley cell

A symmetric baroclinic instability is examined in terms of a Boussinesq fluid contained between two horizontal plates to determine the effects of the Ekman and thermal layers. Governing equations are written for a rotating reference frame, taking into account the Rossby, Ekman, and Prandtl numbers. Equations are defined for the perturbation functions, treated as an eigenvalue problem, and a numerical integration of the full eighth order differential system is performed by a shooting technique. An instability is found to occur in the Hadley cell containing both Ekman and thermal boundary layers when the Richardson number is close to unity. If the Prandtl number is fixed the critical Richardson number decreases with an increasing Ekman number until the Ekman number reaches a certain value, at which time the fluid is stable.

Antar, B. N.↗

NASA NTRSDate not supplied

Simulation and Modeling of Hypersonic Turbulent Boundary Layers Subject to Favorable Pressure Gradients due to Streamline Curvature

Direct numerical simulations (DNS) of favorable-pressure-gradient turbulent boundary layers are presented for a nominal freestream Mach number of 5, with the objective of assessing the limitations of the currently available Reynolds-averaged Navier-Stokes (RANS) models. The favorable pressure gradient is induced by the streamwise curvature of the two-dimensional, planar, convex measurement surface used during experiments at the Texas A&M University. The DNS data shows good comparison with the measured velocity profiles, strain rates, and some, but not all, of the Reynolds-stress components. The discrepancies between the predicted and the measured wall-normal as well as shear stress components are primarily attributed to the lower than actual values inferred from typical PIV measurements of turbulent boundary layers. The DNS data shows a zero or slightly negative Reynolds shear stress in the outer part of the boundary-layer, which is indicative of the decaying turbulent motion under a strong favorable pressure gradient. The DNS data is also compared with the results of RANS computations based on commonly used zero, one, and two equation eddy-viscosity models. The RANS models yield reasonable comparisons with the DNS-based skin friction under zero and weak pressure gradients, but significant discrepancies under a strong pressure gradient. The k-w SST model provided the best overall predictions of skin friction, except in the region where the flow transitions from a favorable to an adverse pressure gradient. While the RANS models examined herein also give good predictions of the Reynolds shear stress under a sufficiently weak pressure gradient, none of those models are able to appropriately capture the reduction in the Reynolds stresses when the flow was subjected to a strong pressure gradient. An a priori assessment of the turbulent heat-flux prediction based on the assumption of a constant turbulent Prandtl number with the DNS data shows that while the constant turbulent Prandtl number model is effective in predicting the wall-normal component of turbulent heat flux, it does not capture the turbulent heat transfer in the streamwise direction for all the pressure gradient cases. The failure of the constant turbulent Prandtl number model highlights a requirement for more advanced models of the turbulent heat flux.

turbulent boundary layers↗

NASA NTRS1943

Heat Transfer Through Turbulent Friction Layers

The "general Prandtl number" Pr(exp 1) - A(sub q)/A Pr, aside from the Reynolds number determines the ratio of turbulent to molecular heat transfer, and the temperature distribution in turbulent friction layers. A(sub q) = exchange coefficient for heat; A = exchange coefficient for momentum transfer. A formula is derived from the equation defining the general Prandtl number which describes the temperature as a function of the velocity. For fully developed thermal boundary layers all questions relating to heat transfer to and from incompressible fluids can be treated in a simple manner if the ratio of the turbulent shear stress to the total stress T(sub t)/T in the layers near the wall is known, and if the A(sub q)/A can be regarded as independent of the distance from the wall. The velocity distribution across a flat smooth channel and deep into the laminar sublayer was measured for isothermal flow to establish the shear stress ratio T(sub t)/T and to extend the universal wall friction law. The values of T(sub t)/T which resulted from these measurements can be approximately represented by a linear function of the velocity in the laminar-turbulent transition zone. The effect of the temperature relationship of the material values on the flow near the wall is briefly analyzed. It was found that the velocity at the laminar boundary (in contrast to the thickness of the laminar layer) is approximately independent of the temperature distribution. The temperature gradient at the wall and the distribution of temperature and heat flow in the turbulent friction layers were calculated on the basis of the data under two equations. The derived formulas and the figures reveal the effects of the Prandtl number, the Reynolds number, the exchange quantities and the temperature relationship of the material values.

Reichardt, H.↗

NASA NTRS1986

Viscous nongeostrophic baroclinic instability

Calculations have been performed of the (linear) stability of a baroclinic flow to three-dimensional perturbations. Both the simple Eady basic state and the rotating Hadley cell of Antar and Fowlis are considered. The independent influences of the Richardson (Ri), thermal Rossby (baroclinicity), Ekman, and Prandtl numbers are examined, as well as the influences of the angle of orientation of the horizontal wave vector and the wavelength. It is shown that if the wavelength is allowed to vary freely, disturbances of the Eady type are preferred (i.e., have greatest growth rate) unless Ri and Ekman numbers are small enough and the thermal Rossby number is large enough. In the latter case, disturbances whose angles of orientation are almost symmetric and whose wavelengths are mesoscale are preferred. If, on the other hand, the wavelength is fixed at a mesoscale size, only the symmetric and almost symmetric modes have growth. By allowing the wave vector orientation to deviate from purely symmetric, it is noted that the region of instability (i.e., critical Ri) is increased, the extent of which is greater for longer wavelength. For Prandtl number = 1, permitting the angle to be nonsymmetric demonstrates the existence of two maxima in growth rate at opposite angles of orientation and with very different energetics. For Prandtl number far enough from one and for large enough dissipation, only one of these two modes has positive growth rates. Growing oscillatory modes were found for some cases.

Miller, T. L.↗

DOE OSTI2017

Scaling and excitation of combined convection in a rapidly rotating plane layer

The optimum (to my mind) scaling of the combined thermal and compositional convection in a rapidly rotating plane layer is proposed.This scaling follows from self-consistent estimates of typical physical quantities. Similarity coefficients are introduced for the ratio convection dissipation/convection generation (s) and the ratio thermal convection/compositional convection (r). The third new and most important coefficient δ is the ratio of the characteristic size normal to the axis of rotation to the layer thickness. The faster the rotation, the lower δ. In the case of the liquid Earth core, δ ~ 10{sup –3} substitutes for the generally accepted Ekman number (E ~ 10{sup –15}) and s ~ 10{sup –6} substitutes for the inverse Rayleigh number 1/Ra ~ 10{sup –30}. It is found that, at turbulent transport coefficients, number s and the Prandtl number are on the order of unity for any objects and δ is independent of transport coefficients. As a result of expansion in powers of δ, an initially 3D system of six variables is simplified to an almost 2D system of four variables without δ. The problem of convection excitation in the main volume is algebraically solved and this problem for critical values is analytically solved. Dispersion relations and general expressions for critical wavenumbers, numbers s (which determine Rayleigh numbers), other critical parameters, and asymptotic solutions are derived. Numerical estimates are made for the liquid cores in the planets that resemble the Earth. Further possible applications of the results obtained are proposed for the interior of planets, moons, their oceans, stars, and experimental objects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

Numerical modeling of heat transfer in molten silicon during directional solidification process

Numerical investigation is performed for some of the thermal and fluid flow properties of silicon melt during directional solidification by numerical modeling. Dimensionless numbers are extremely useful to understand the heat and mass transfer of fluid flow on Si melt and control the flow patterns during crystal growth processes. The average grain size of whole crystal would increase when the melt flow is laminar. In the silicon growth process, the melt flow is mainly driven by the buoyancy force resulting from the horizontal temperature gradient. The thermal and flow pattern influences the quality of the crystal through the convective heat and mass transport. The computations are carried out in a 2D axisymmetric model using the finite-element technique. The buoyancy effect is observed in the melt domain for a constant Rayleigh number and for different Prandtl numbers. The convective heat flux and Reynolds numbers are studied in the five parallel horizontal cross section of melt silicon region. And also, velocity field is simulated for whole melt domain with limited thermal boundaries. The results indicate that buoyancy forces have a dramatic effect on the most of melt region except central part.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1987

Transport of passive scalars in a turbulent channel flow

A direct numerical simulation of a turbulent channel flow with three passive scalars at different molecular Prandtl numbers is performed. Computed statistics including the turbulent Prandtl numbers are compared with existing experimental data. The computed fields are also examined to investigate the spatial structure of the scalar fields. The scalar fields are highly correlated with the streamwise velocity; the correlation coefficient between the temperature and the streamwise velocity is as high as 0.95 in the wall region. The joint probability distributions between the temperature and velocity fluctuations are also examined; they suggest that it might be possible to model the scalar fluxes in the wall region in a manner similar to the Reynolds stresses.

Kim, John↗

NASA NTRS1987

Transport of passive scalars in a turbulent channel flow

Kim, John↗

DOE OSTI2010

Dynamo Transition

In this article we review the experimental and numerical results related to the dynamo transitions. Recent experiments of Von Karman Sodium (VKS) exhibit various dynamo states including constant, time-periodic, and chaotic magnetic fields. Similarly pseudospectral simulations of dynamo show constant, time-periodic, quasiperiodic, and chaotic magnetic field configurations. One of the windows for the magnetic Prandtl number of unity shows period doubling route to chaos. Quasiperiodic route to chaos has been reported for the Prandtl number of 0.5. The dynamo simulations also reveal coexisting multiple attractors that were obtained for different initial conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2017

EXTENSION OF THE MURAM RADIATIVE MHD CODE FOR CORONAL SIMULATIONS

We present a new version of the MURaM radiative magnetohydrodynamics (MHD) code that allows for simulations spanning from the upper convection zone into the solar corona. We implement the relevant coronal physics in terms of optically thin radiative loss, field aligned heat conduction, and an equilibrium ionization equation of state. We artificially limit the coronal Alfvén and heat conduction speeds to computationally manageable values using an approximation to semi-relativistic MHD with an artificially reduced speed of light (Boris correction). We present example solutions ranging from quiet to active Sun in order to verify the validity of our approach. We quantify the role of numerical diffusivity for the effective coronal heating. We find that the (numerical) magnetic Prandtl number determines the ratio of resistive to viscous heating and that owing to the very large magnetic Prandtl number of the solar corona, heating is expected to happen predominantly through viscous dissipation. We find that reasonable solutions can be obtained with values of the reduced speed of light just marginally larger than the maximum sound speed. Overall this leads to a fully explicit code that can compute the time evolution of the solar corona in response to photospheric driving using numerical time steps not much smaller than 0.1 s. Numerical simulations of the coronal response to flux emergence covering a time span of a few days are well within reach using this approach.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1985

Thermocapillary Flows and Their Stability: Effects of Surface Layers and Contamination

The objective is to study quantitatively thermocapillary flows in two-dimensional and axisymmetric geometries in order to learn the impact of such flows on float-zone refining configurations. The work begins with the analytical study of steady flows in slots and zones to obtain predictions of surface curvature, flow and heat transport. It then examines the instability of such flows. Three new instability mechanisms have been identified and the results show that thermocapillary flows of high Prandtl number fluid become convectively unstable while those of low Prandtl number fluid become hydrodynamically unstable. The work includes numerical simulations of steady, high Marangoni number flows in unit aspect ratio boxes and asymptotic solutions in large boxes. Furthermore, the effect of surface contamination on the steady thermocapillary flows was examined. This fundamental study focused on fluid dynamical systems used to understand events in the flow of melts in float-zone geometries both on Earth and in microgravity environments.

Davis, S. H.↗

NASA NTRS1991

Reynolds shear stress and heat flux calculations in a fully developed turbulent duct flow

The use of a modified form of the Van Driest mixing length for a fully developed turbulent channel flow leads to mean velocity and Reynolds stress distributions that are in close agreement with data obtained either from experiments or direct numerical simulations. The calculations are then extended to a nonisothermal flow by assuming a constant turbulent Prandtl number, the value of which depends on the molecular Prandtl number. Calculated distributions of mean temperature and lateral heat flux are in reasonable agreement with the simulations. The extension of the calculations to higher Reynolds numbers provides some idea of the Reynolds number required for scaling on wall variables to apply in the inner region of the flow.

Antonia, R. A.↗

NASA NTRS1993

Stability and instability of thermocapillary convection in models of the float-zone crystal-growth process

This project was concerned with the determination of conditions of guaranteed stability and instability for thermocapillary convection in a model of the float-zone crystal-growth process. This model, referred to as the half-zone, was studied extensively, both experimentally and theoretically. Our own earlier research determined, using energy-stability theory, sufficient conditions for stability to axisymmetric disturbances. Nearly all results computed were for the case of a liquid with Prandtl Number Pr = 1. Attempts to compute cases for higher Prandtl numbers to allow comparison with the experimental results of other researchers were unsuccessful, but indicated that the condition guaranteeing stability against axisymmetric disturbances would be a value of the Marangoni number (Ma), significantly higher than that at which oscillatory convection was observed experimentally. Thus, additional results were needed to round out the stability picture for this model problem. The research performed under this grant consisted of the following: (1) computation of energy-stability limits for non-axisymmetric disturbances; (2) computation of linear-stability limits for axisymmetric and non-axisymmetric disturbances; (3) numerical simulation of the basic state for half- and full-zones with a deformable free surface; and (4) incorporation of radiation heat transfer into a model energy-stability problem. Each of these is summarized briefly below.

Neitzel, G. P.↗

NASA NTRS1984

Asymmetric instabilities in buoyancy-driven flow in a tall vertical annulus

Linear stability of the one-dimensional flow between infinite vertical coaxial cylinders induced by heating the inner cylinder is considered for various ratios (kappa) of the inner radius to the outer radius, and for Prandtl numbers (P) appropriate to air and water. For air with P = 0.71 the least stable disturbance is nonaxisymmetric for kappa less than 0.44 and is axisymmetric for kappa greater than 0.44, and in either case the instability is due to the action of the shear forces. For P = 3.5, the situation is similar, except that the asymmetric shear mode is superseded by an axisymmetric instability driven by buoyancy forces for kappa = 0.03-0.16. Wave speeds, wavenumbers, and critical Grashof numbers for these cases and for the case of zero Prandtl number are given.

Mcfadden, G. B.↗

DOE OSTI2013

SATURATION OF THE MAGNETO-ROTATIONAL INSTABILITY IN STRONGLY RADIATION-DOMINATED ACCRETION DISKS

The saturation level of the magneto-rotational instability (MRI) in a strongly radiation-dominated accretion disk is studied using a new Godunov radiation MHD code in the unstratified shearing box approximation. Since vertical gravity is neglected in this work, our focus is on how the MRI saturates in the optically thick mid-plane of the disk. We confirm that turbulence generated by the MRI is very compressible in the radiation-dominated regime, as found by previous calculations using the flux-limited diffusion approximation. We also find little difference in the saturation properties in calculations that use a larger horizontal domain (up to four times the vertical scale height in the radial direction). However, in strongly radiation pressure dominated disks (one in which the radiation energy density reaches {approx}1% of the rest mass energy density of the gas), we find that Maxwell stress from the MRI turbulence is larger than the value produced when radiation pressure is replaced with the same amount of gas pressure. At the same time, the ratio between Maxwell stress and Reynolds stress is increased by almost a factor of eight compared with the gas pressure dominated case. We suggest that this effect is caused by radiation drag, which acts like bulk viscosity and changes the effective magnetic Prandtl number of the fluid. Radiation viscosity significantly exceeds both the microscopic plasma viscosity and resistivity, ensuring that radiation-dominated systems occupy the high magnetic Prandtl number regime. Nevertheless, we find that radiative shear viscosity is negligible compared to the Maxwell stress and Reynolds stress in the flow. This may have important implications for the structure of radiation-dominated accretion disks.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2018

Prediction of heat transfer coefficients for forced convective boiling of N 2 -hydrocarbon mixtures at cryogenic conditions using artificial neural networks

A key problem faced in the design of heat exchangers, especially for cryogenic applications, is the determination of convective heat transfer coefficients in two-phase flow such as condensation and boiling of non-azeotropic refrigerant mixtures. Here, this paper proposes and evaluates three models for estimating the convective coefficient during boiling. These models are developed using computational intelligence techniques. The performance of the proposed models is evaluated using the mean relative error (mre), and compared to two existing models: the modified Granryd’s correlation and the Silver-Bell-Ghaly method. The three proposed models are distinguished by their architecture. The first is based on directly measured parameters (DMP-ANN), the second is based on equivalent Reynolds and Prandtl numbers (eq-ANN), and the third on effective Reynolds and Prandtl numbers (eff-ANN). In conclusion, the results demonstrate that the proposed artificial neural network (ANN)-based approaches greatly outperform available methodologies. While Granryd's correlation predicts experimental data within a mean relative error mre = 44% and the S-B-G method produces mre = 42%, DMP-ANN has mre = 7.4% and eff-ANN has mre = 3.9%. Considering that eff-ANN has the lowest mean relative error (one tenth of previously available methodologies) and the broadest range of applicability, it is recommended for future calculations. Implementation is straightforward within a variety of platforms and the matrices with the ANN weights are given in the appendix for efficient programming.

42 ENGINEERING↗

NASA NTRS1982

Computer program for solving laminar, transitional, or turbulent compressible boundary-layer equations for two-dimensional and axisymmetric flow

A numerical algorithm and computer program are presented for solving the laminar, transitional, or turbulent two dimensional or axisymmetric compressible boundary-layer equations for perfect-gas flows. The governing equations are solved by an iterative three-point implicit finite-difference procedure. The software, program VGBLP, is a modification of the approach presented in NASA TR R-368 and NASA TM X-2458, respectively. The major modifications are: (1) replacement of the fourth-order Runge-Kutta integration technique with a finite-difference procedure for numerically solving the equations required to initiate the parabolic marching procedure; (2) introduction of the Blottner variable-grid scheme; (3) implementation of an iteration scheme allowing the coupled system of equations to be converged to a specified accuracy level; and (4) inclusion of an iteration scheme for variable-entropy calculations. These modifications to the approach presented in NASA TR R-368 and NASA TM X-2458 yield a software package with high computational efficiency and flexibility. Turbulence-closure options include either two-layer eddy-viscosity or mixing-length models. Eddy conductivity is modeled as a function of eddy viscosity through a static turbulent Prandtl number formulation. Several options are provided for specifying the static turbulent Prandtl number. The transitional boundary layer is treated through a streamwise intermittency function which modifies the turbulence-closure model. This model is based on the probability distribution of turbulent spots and ranges from zero to unity for laminar and turbulent flow, respectively. Several test cases are presented as guides for potential users of the software.

Harris, J. E.↗