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 271 records · Page 15

Modeling Turbulent Combustion for Variable Prandtl and Schmidt Number

This report consists of two abstracts submitted for possible presentation at the AIAA Aerospace Science Meeting to be held in January 2005. Since the submittal of these abstracts we are continuing refinement of the model coefficients derived for the case of a variable Turbulent Prandtl number. The test cases being investigated are a Mach 9.2 flow over a degree ramp and a Mach 8.2 3-D calculation of crossing shocks. We have developed an axisymmetric code for treating axisymmetric flows. In addition the variable Schmidt number formulation was incorporated in the code and we are in the process of determining the model constants.

Hassan, H. A.↗

Stellar convection 2: A multi-mode numerical solution for convection in spheres

The convective flow of a self gravitating sphere of Boussinesq fluid for small Reynolds and Peclet numbers is numerically determined. The decomposition of the equations of motion into modes is reviewed and a relaxation method is developed and presented to compute the solutions to these equations. The stable equilibrium flow for a Rayleigh number of 10 to the 4th power and a Prandtl number of 10 is determined. The 2 and 3 dimensional spectra of the kinetic and thermal energies and the convective flux as a function of wavelengths are calculated in terms of modes. The anisotropy of the flow as a function of wavelength is defined.

Marcus, P. S.↗

Theoretical analyses of Baroclinic flows

A stability analysis of a thin horizontal rotating fluid layer which is subjected to arbitrary horizontal and vertical temperature gradients is presented. The basic state is a nonlinear Hadley cell which contains both Ekman and thermal boundary layers; it is given in closed form. The stability analysis is based on the linearized Navier-Stokes equations, and zonally symmetric perturbations in the form of waves propagating in the meridional direction are considered. Numerical methods were used for the stability problem. It was found that the instability sets in when the Richardson number is close to unity and that the critical Richardson number is a non-monotonic function of the Prandtl number. Further, it was found that the critical Richardson number decreases with increasing Ekman number until a critical value of the Ekman number is reached beyond which the fluid is stable.

Antar, B.↗

Evaporation, Heat Transfer, and Velocity Distribution in Two-Dimensional and Rotationally Symmetrical Laminar Boundary-Layer Flow

The fundamental boundary layer equations for the flow, temperature and concentration fields are presented. Two dimensional symmetrical and unsymmetrical and rotationally symmetrical steady boundary layer flows are treated as well as the transfer boundary layer. Approximation methods for the calculation of the transfer layer are discussed and a brief survey of an investigation into the validity of the law that the Nusselt number is proportional to the cube root of the Prandtl number is presented.

Froessling, Nils↗

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↗

Critical stability of almost adiabatic convection in a rapidly rotating thick spherical shell

In this work, the convection equations in the almost adiabatic approximation is studied for which the choice of physical parameters is primarily based on possible applications to the hydrodynamics of the deep interiors of the Earth and planets and moons of the terrestrial group. The initial system of partial differential equations (PDEs) was simplified to a single second-order ordinary differential equation for the pressure or vertical velocity component to investigate the linear stability of convection. The critical frequencies, modified Rayleigh numbers, and distributions of convection are obtained at various possible Prandtl numbers and in different thick fluid shells. An analytical WKB-type solution was obtained for the case when the inner radius of the shell is much smaller than the outer radius and convective sources are concentrated along the inner boundary.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modeling Scramjet Flows with Variable Turbulent Prandtl and Schmidt Numbers

A complete turbulence model, where the turbulent Prandtl and Schmidt numbers are calculated as part of the solution and where averages involving chemical source terms are modeled, is presented. The ability of avoiding the use of assumed or evolution Probability Distribution Functions (PDF's) results in a highly efficient algorithm for reacting flows. The predictions of the model are compared with two sets of experiments involving supersonic mixing and one involving supersonic combustion. The results demonstrate the need for consideration of turbulence/chemistry interactions in supersonic combustion. In general, good agreement with experiment is indicated.

Xiao, X.↗

Simulation of the turbulent Rayleigh-Benard problem using a spectral/finite difference technique

The three-dimensional, incompressible Navier-Stokes and energy equations with the Bousinesq assumption have been directly simulated at a Rayleigh number of 3.8 x 10 to the 5th power and a Prandtl number of 0.76. In the vertical direction, wall boundaries were used and in the horizontal, periodic boundary conditions were used. A spectral/finite difference numerical method was used to simulate the flow. The flow at these conditions is turbulent and a sufficiently fine mesh was used to capture all relevant flow scales. The results of the simulation are compared to experimental data to justify the conclusion that the small scale motion is adequately resolved.

Eidson, T. M.↗

Thermocapillary flow and melt/solid interfaces in floating-zone crystal growth under microgravity

Computer simulation of steady-state axisymmetrical heat transfer and fluid flow was conducted to study thermocapillary flow and melt/solid interfaces in floating-zone crystal growth under microgravity. The effects of key variables on the extent of thermocapillary flow in the melt zone, the shapes of melt/solid interfaces and the length of the melt zone were discussed. These variables are: (1) the temperature coefficient of surface tension (or the Marangoni number), (2) the pulling speed (or the Peclet number), (3) the feed rod radius, (4) the ambient temperature distribution, (5) the heat transfer coefficient (or the Biot number), and (6) the thermal diffusivity of the material (or the Prandtl number).

Lan, C. W.↗

Heat Transmission in the Boundary Layer

In the present paper which deals with the heat transfer between the gas and the wall for large temperature drops and large velocities use is made of the method of Dorodnitsyn of the introduction of a new independent variable, with this difference, however, that the relation between the temperature field (that is, density) and the velocity field in the general case considered is not assumed given but is determined from the solution of the problem. The effect of the compressibility arising from the heat transfer is thus taken into account (at the same time as the effect of the compressibility at the large velocities). A method is given for determining the coefficients of heat transfer and the friction coefficients required in many technical problems for a curved wall in a gas flow at large Mach numbers and temperature drops. The method proposed is applicable both for Prandtl number P = 1 and for P not equal to 1.

Kalikhman, L. E.↗

Penetrative cellular convection in a stratified atmosphere

In the present investigation of penetrative convection within a simple compressible model, the middle one of the three layers of differing stratification prior to the onset of convection is a convectively unstable polytrope bounded above and below by two stably stratified polytropes. One- and two-mode steady solutions with hexagonal planforms have been studied for Rayleigh numbers up to aobut 1000 times critical, and for a range of Prandtl numbers, horizontal wavenumbers, and stratifications. These indicate that the penetration into the lower stable layer by downward plumes is substantially larger in a stratified medium than in a Boussinesq fluid, and produces an extended region of adiabatic stratification. The strong asymmetry between upward and downward penetration in compressible media has major implications for the mixing of stable regions above and below stellar convection zones.

Massaguer, J. M.↗

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 $⩽$Ra$⩽$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 $⩾$ 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↗

Ocean Turbulence, III: New GISS Vertical Mixing Scheme

We have found a new way to express the solutions of the RSM (Reynolds Stress Model) equations that allows us to present the turbulent diffusivities for heat, salt and momentum in a way that is considerably simpler and thus easier to implement than in previous work. The RSM provides the dimensionless mixing efficiencies Gamma-alpha (alpha stands for heat, salt and momentum). However, to compute the diffusivities, one needs additional information, specifically, the dissipation Epsilon. Since a dynamic equation for the latter that includes the physical processes relevant to the ocean is still not available, one must resort to different sources of information outside the RSM to obtain a complete Mixing Scheme usable in OGCMs. As for the RSM results, we show that the Gamma-alpha s are functions of both Ri and Rq (Richardson number and density ratio representing double diffusion, DD); the Gamma-alpha are different for heat, salt and momentum; in the case of heat, the traditional value Gamma-h = 0.2 is valid only in the presence of strong shear (when DD is inoperative) while when shear subsides, NATRE data show that Gamma-h can be three times as large, a result that we reproduce. The salt Gamma-s is given in terms of Gamma-h. The momentum Gamma-m has thus far been guessed with different prescriptions while the RSM provides a well defined expression for Gamma-m(Ri,R-rho). Having tested Gamma-h, we then test the momentum Gamma-m by showing that the turbulent Prandtl number Gamma-m/Gamma-h vs. Ri reproduces the available data quite well. As for the dissipation epsilon, we use different representations, one for the mixed layer (ML), one for the thermocline and one for the ocean;s bottom. For the ML, we adopt a procedure analogous to the one successfully used in PB (planetary boundary layer) studies; for the thermocline, we employ an expression for the variable epsilon/N(exp 2) from studies of the internal gravity waves spectra which includes a latitude dependence; for the ocean bottom, we adopt the enhanced bottom diffusivity expression used by previous authors but with a state of the art internal tidal energy formulation and replace the fixed Gamma-alpha = 0.2 with the RSM result that brings into the problem the Ri, R-rho dependence of the Gamma-alpha; the unresolved bottom drag, which has thus far been either ignored or modeled with heuristic relations, is modeled using a formalism we previously developed and tested in PBL studies. We carried out several tests without an OGCM. Prandtl and flux Richardson numbers vs. Ri. The RSM model reproduces both types of data satisfactorily. DD and Mixing efficiency Gamma-h(Ri,Rq). The RSM model reproduces well the NATRE data. Bimodal epsilon-distribution. NATRE data show that epsilon (Ri < 1) approximately equals 10epsilon(Ri > 1), which our model reproduces. Heat to salt flux ratio. In the Ri much greater than 1 regime, the RSM predictions reproduce the data satisfactorily. NATRE mass diffusivity. The z-profile of the mass diffusivity reproduces well the measurements at NATRE. The local form of the mixing scheme is algebraic with one cubic equation to solve.

Canuto, V. M.↗

Convective Instability of a Gravity Modulated Fluid Layer with Surface Tension Variation

Gravity modulation of an unbounded fluid layer with surface tension variations along its free surface is investigated. In parameter space of (wavenumber, Marangoni number) modulation has a destabilizing effect on the unmodulated neutral stability curve for large Prandtl number, Pr, and small modulation frequency, Omega, while a stabilizing effect is observed for small Pr and large Omega. As Omega yields infinity, the modulated neutral stability curves approach the unmodulated neutral stability curve. At certain values of Pr and L2 multiple minima are observed and the neutral stability curves become highly distorted. Closed regions of subharmonic instability are also observed. Alternating regions of synchronous and subharmonic instability separated by very thin stable regions are observed in (1/Omega,g(sub 1)) space for the singly diffusive cases. Quasiperiodic behavior in addition to the synchronous and subharmonic responses, are observed for the case of a double diffusive fluid layer. Minimum acceleration amplitudes were observed to closely correspond with a subharmonic response, Lambda(sub im) = Omega/2 .

Skarda, J. Raymond Lee↗

Modification of the Two-equation Turbulence Model in NPARC to a Chien Low Reynolds Number K-epsilon Formulation

This report documents the changes that were made to the two-equation k-epsilon turbulence model in the NPARC (National-PARC) code. The previous model based on the low Reynolds number model of Speziale, was replaced with the low Reynolds number k-epsilon model of Chien. The most significant difference was in the turbulent Prandtl numbers appearing in the diffusion terms of the k and epsilon transport equations. A new inflow boundary condition and stability enhancements were also implemented into the turbulence model within NPARC. The report provides the rationale for making the change to the Chien model, code modifications required, and comparisons of the performances of the new model with the previous k-epsilon model and algebraic models used most often in PARC/NPARC. The comparisons show that the Chien k-epsilon model installed here improves the capability of NPARC to calculate turbulent flows.

Georgiadis, Nicholas J.↗

Analysis of turbulent heat transfer, mass transfer, and friction in smooth tubes at high Prandtl and Schmidt numbers

The expression for eddy diffusivity from a previous analysis was modified in order to account for the effect of kinematic viscosity on the turbulence in the region close to a wall. By using the modified expression, good agreement was obtained between predicted and experimental results for heat and mass transfer at Prandtl and Schmidt numbers between 0.5 and 3000. The effects of length-to-diameter ratio and of variable viscosity were also investigated for a wide range of Prandtl numbers.

Deissler, Robert G↗

Unsteady thermocapillary migration of bubbles

Upon the introduction of a gas bubble into a liquid possessing a uniform thermal gradient, an unsteady thermo-capillary flow begins. Ultimately, the bubble attains a constant velocity. This theoretical analysis focuses upon the transient period for a bubble in a microgravity environment and is restricted to situations wherein the flow is sufficiently slow such that inertial terms in the Navier-Stokes equation and convective terms in the energy equation may be safely neglected (i.e., both Reynolds and Marangoni numbers are small). The resulting linear equations were solved analytically in the Laplace domain with the Prandtl number of the liquid as a parameter; inversion was accomplished numerically using a standard IMSL routine. In the asymptotic long-time limit, the theory agrees with the steady-state theory of Young, Goldstein, and Block. The theory predicts that more than 90 percent of the terminal steady velocity is achieved when the smallest dimensionless time, i.e., the one based upon the largest time scale-viscous or thermal-equals unity.

Dill, Loren H.↗

Thermoconvective Instability in a Rotating Magnetic Field

The effect of a rotating magnetic field (RMF) on the stability of a fluid contained in a cylindrical column and heated from below is investigated. The RMF increases the critical Rayleigh number for asymmetric flow modes but does not affect the onset of instability for axisymmetric modes. The critical Rayleigh number is dependent upon the relative penetration of the magnetic field into the cylinder and the Prandtl number of the fluid. Instability first develops in the form of a single asymmetric meridional roll rotating around the axis of the cylinder, driven by the azimuthal component of the magnetic field.

Volz, M. P.↗