SEARCH · Engineering Papers

Results for “RAYLEIGH 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 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS1975

Modal equations for cellular convection

We expand the fluctuating flow variables of Boussinesq convection in the planform functions of linear theory. Our proposal is to consider a drastic truncation of this expansion as a possible useful approximation scheme for studying cellular convection. With just one term included, we obtain a fairly simple set of equations which reproduces some of the qualitative properties of cellular convection and whose steady-state form has already been derived by Roberts (1966). This set of 'modal equations' is analyzed at slightly supercritical and at very high Rayleigh numbers. In the latter regime the Nusselt number varies with Rayleigh number just as in the mean-field approximation with one horizontal scale when the boundaries are rigid. However, the Nusselt number now depends also on the Prandtl number in a way that seems compatible with experiment. The chief difficulty with the approach is the absence of a deductive scheme for deciding which planforms should be retained in the truncated expansion.

Gough, D. O.↗

NASA NTRS1993

A study of the effects of macrosegregation and buoyancy-driven flow in binary mixture solidification

A generalized anisotropic porous medium approach is developed for modelling the flow, heat and mass transport processes during binary mixture solidification. Transient predictions are obtained using FEM, coupled with an implicit time-marching scheme, for solidification inside a two-dimensional rectangular enclosure. A parametric study focusing attention on the effects of solutal buoyancy and thermal buoyancy is presented. It is observed that three parameters, namely the thermal Rayleigh number, the solutal Rayleigh number, and the relative density change parameter, significantly alter the flow fields in the liquid and the mushy regions. Depending upon the nature of these flow fields, the solute enrichment caused by macrosegregation may occur in the top or the bottom region of the enclosure.

Sinha, S. K.↗

NASA NTRS1977

Numerical solutions of single-mode convection equations

In the Boussinesq approximation, single-mode equations describing thermal convection are constructed by expanding the fluctuating velocity and temperature fields in a complete set of functions (or planforms) of the horizontal coordinates and retaining just one term. Numerical solutions of the single-mode equations are investigated, chief consideration being given to hexagonal planforms. Extensive surveys of steady solutions are presented for various Rayleigh numbers, Prandtl numbers, and horizontal wavenumbers. The dependences on Rayleigh number and Prandtl number at very large Rayleigh number are in satisfactory agreement with the results of asymptotic expansions.

Toomre, J.↗

NASA NTRS1995

Rheological transition in mantle convection with a composite temperature-dependent, non-Newtonian and Newtonian rheology

Numerical simulations of mantle convection with a composite temperature-dependent, Newtonian and non-Newtonian creep law have revealed a transition in the dominant creep mechanism with the increasing vigour of convection. Newtonian creep is found to dominate in the low Rayleigh number regime. With sufficiently high effective Rayleigh number, the overall creep mechanism in the convective flow becomes non-Newtonian. The transitional Rayleigh number increases strongly with the activation energy. These results would suggest a scenario that in the early epochs of Earth the flow in the mantle would have been governed by non-Newtonian rheology and would have exhibited both strong spatial and temporal fluctuations. With time the flow mechanism would behave like a Newtonian fluid and would have a different time-dependent character. In time-dependent Newtonian-dominated flows there are still localized features with distinctly non-Newtonian character. Our analysis of the relative contributions to the lateral viscosity field supports the idea that the inference of the nature of lateral viscosity heterogeneities by seismic tomography may be strongly contaminated by the dominant non-Newtonian contributions to the total lateral viscosity field.

Van Den Berg, Arie P.↗

NASA NTRS1972

Finite-amplitude thermal convection within a self-gravitating fluid sphere.

Finite-difference calculations have been carried out to determine the structure of finite-amplitude thermal convection within a self-gravitating fluid sphere with uniform heat release. For a fixed-surface boundary condition, single-cell convection breaks up into double-cell convection at a Rayleigh number of 30,000, at a Rayleigh number of 500,000 four-cell convection is observed. With a free-surface boundary condition only single cell convection is obtained up to a Rayleigh number of 5,000,000.

Hsui, A. T.↗

NASA NTRS1999

Rayleigh Convection in a Rotating Magnetic Field

The onset of steady and time-dependent flows were measured in a cylinder of liquid gallium subjected to a vertical temperature gradient and a horizontal rotating magnetic field (RMF). By varying the applied vertical temperature gradients and the strength of the applied magnetic field, a stability diagram in Rayleigh number - Hartmann number parameter space was constructed. The application of a RMF increased the critical value of the Rayleigh number by up to a factor of 10. However, there was observed a critical value of the Hartmann number where instability set in, regardless of the value of the Rayleigh number. Significant hysteresis in the value of the critical Rayleigh number was observed depending on whether the applied thermal gradient was increasing or decreasing.

Volz, M. P.↗

NASA NTRS2004

Convective Instability in Ice I: Application to Callisto and Ganymede

Laboratory experiments measuring ice rheology suggest that it deforms under the influence of several nonNewtonian creep mechanisms, where the viscosity depends on both strain rate and temperature [Goldsby & Kohlstedt, 2001]. Whether or not a fluid with a purely temperature-dependent viscosity convects can be determined by comparing the Rayleigh number of the system to the critical Rayleigh number (Racr), which depends on rheological, thermal, and physical parameters of the fluid layer. However, in a nonNewtonian fluid where viscosity depends on the strain rate (i.e. velocity), convection can only occur if a temperature or velocity perturbation is issued to the system to lower the viscosity and permit fluid motions. Therefore, whether convection occurs in an ice I layer depends on initial conditions, in addition to rheological, thermal, and physical properties of the layer. We show new results for a scaling between the critical Rayleigh number and perturbation amplitude for grain boundary sliding rheology. This scaling can be used to determine the conditions required to initiate convection in the ice I shell of a generic icy satellite. We use this scaling to judge the convective instability of Ganymede and Callisto's ice shells in the absence of tidal dissipation.

Barr, A. C.↗

NASA NTRS1981

Parametric results for heat transfer across honeycomb sandwich panels

The problem of heat transfer across honeycomb sandwich panels is theoretically investigated. The basic objective is to predict the rear surface temperature as a function of time when the front surface is exposed to a prescribed heat flux; the faces are bounded by planar sheets so that the air circulation in the honeycomb cells is bounded and contained. The influence of the air circulation on convective heat transfer is explicitly included. Drawing upon previous investigators' results that indicated various relations for the Nusselt numbers vs Rayleigh numbers in different regimes of the aspect ratio (defined as the ratio of the cell height to cell width), parametric calculations are performed to predict the rear surface temperatures. Chemical degradation of the material, especially the glue holding the face material on the end, is important but has not yet been included. The results indicate that decreasing the cell width, increasing the cell wall thickness and increasing the cell height all have beneficial effects upon the heat transfer; i.e., all of these variations reduce rear surface temperature for a given time for a prescribed heat flux on the front surface.

Ramohalli, K.↗

NASA NTRS1991

Source function approach to turbulence bulk properties

A turbulence model to compute bulk properties is presented in which a prescribed source function is used to provide the rate of energy input into the turbulent flow, and the EDQNM model is used to treat the nonlinear transfer in the Navier-Stokes equations. The predictions of the model are tested against (1) the measured Nusselt number versus Rayleigh number relation in turbulent laboratory convection, and (2) the measured bulk kinetic energies and dissipation rates in turbulent channel flow for Reynolds numbers Re = 12,300 and 30,800. In addition, a sensitivity study is performed with respect to the choice of Kolmogorov constant. Generally, the predictions of the model are in reasonable accord with the available experimental data.

Canuto, V. M.↗

NASA NTRS1981

Numerical simulation of natural convection in a spherical container due to cooling at the center (idealization of the Lal/Kroes experiment)

Natural convection in a spherical container with cooling at the center was numerically simulated using a numerical fluid dynamics computer program. The numerical analysis was simplified by assuming axisymmetric flow in the spherical container, with the symmetry axis being a sphere diagonal parallel to the gravity vector. This axisymmetric spherical geometry was intended as an idealization of the proposed Lal/Kroes crystal growing experiment to be performed on Spacelab. Results were obtained for a range of Rayleigh numbers from 25 to 10,000. The computed velocities were found to be approximately proportional to the Rayleigh number over the range of Rayleigh numbers investigated.

Robertson, S. J.↗

NASA NTRS1999

Predicted Variations in Flow Patterns in a Horizontal CVD Reactor

Expressions in terms of common reactor operating parameters were derived for the ratio of the Grashof number to the Reynolds number, Gr/Re, the ratio of the Grashof to the square of 2 the Reynolds number, Gr/Re(exp 2), and the Rayleigh number, Ra. Values for these numbers were computed for an example horizontal CVD reactor and compared to numerical simulations to gauge their effectiveness as predictors of the presence or absence of transverse and longitudinal rolls in the reactor. Comparisons were made for both argon and hydrogen carrier gases over the pressure range 2- 101 kPa. Reasonable agreement was achieved in most cases when using Gr/Re to predict the presence of transverse rolls and Ra to predict the presence of longitudinal rolls. The ratio Gr/Re(exp 2) did not yield useful predictions regarding the presence of transverse rolls. This comparison showed that the ratio of the Grashof number to the Reynolds number, as well as the Rayleigh number, can be used to predict the presence or absence of transverse and longitudinal rolls in a horizontal CVD reactor for a given set of reactor conditions. These predictions are approximate, and care must be exercised when making predictions near transition regions.

Kuczmarski, Maria A.↗

NASA NTRS2001

Three-Dimensional Simulations of Mantle Convection in Io

Io has very high surface heat flow and an abundance of volcanic activity, which are thought to be driven by nonuniform tidal heating in its interior. This nonuniform heat is transported to the base of the lithosphere by very vigorous convection in Io's silicate mantle, the form of which is presumably responsible for the distribution of surface features such as volcanoes and mountains. We here present three-dimensional spherical calculations of mantle convection in Io, in order to ascertain the likely form of this convection and the resulting distribution of heat flow at the surface and core-mantle boundary. Different models of tidal dissipation are considered: the endmember scenarios (identified by M. N. Ross and G. Schubert) of dissipation in the entire mantle, or dissipation in a thin (approximately 100-km-thick) asthenosphere, as well as the 'preferred' distribution of M. N. Ross et al. comprising 1/3 mantle and 2/3 asthenosphere heating. The thermal structure of Io's mantle and asthenosphere is found to be strongly dependent on tidal heating mode, as well as whether the mantle-asthenosphere boundary is permeable or impermeable. Results indicate a large-scale flow pattern dominated by the distribution of tidal heating, with superimposed small-scale asthenospheric instabilities that become more pronounced with increasing Rayleigh number. These small-scale instabilities spread out the surface heat flux, resulting in smaller heat flux variations with increasing Rayleigh number. Scaled to Io's Rayleigh number of O(10(exp 12)) variations of order a few percent are expected. This small but significant variation in surface heat flux may be compatible with the observed distributions of volcanic centers and mountains, which appear fairly uniform at first sight but display a discernible distribution when suitably processed. The observed distribution of volcanic centers is similar to the asthenosphere heating distribution, implying that most of the tidal heating in Io occurs in an asthenosphere.

Tackley, Paul J.↗

NASA NTRS1992

Natural convection between concentric spheres

A finite-difference solution for steady natural convective flow in a concentric spherical annulus with isothermal walls has been obtained. The stream function-vorticity formulation of the equations of motion for the unsteady axisymmetric flow is used; interest lying in the final steady solution. Forward differences are used for the time derivatives and second-order central differences for the space derivatives. The alternating direction implicit method is used for solution of the discretization equations. Local one-dimensional grid adaptation is used to resolve the steep gradients in some regions of the flow at large Rayleigh numbers. The break-up into multi-cellular flow is found at high Rayleigh numbers for air and water, and at significantly low Rayleigh numbers for liquid metals. Excellent agreement with previous experimental and numerical data is obtained.

Garg, Vijay K.↗

NASA NTRS1992

Penetrative convection - Parametrized expression for the growth rates

The equations determining the linear growth rate omega characterizing a convectively unstable fluid with Rayleigh number R(u) bounded below by an impenetrable free boundary and above by a convectively stable fluid with Rayleigh number R(s), are solved numerically. Using the analytical Rayleigh-Benard growth rate omega (RB) as a convenient functional form, it is possible to fit the numerical values for omega if the vertical wave number k(z) = n(pi) and the Rayleigh number R(RB) are taken to be functions of R(s), R(u), and the horizontal wave number k-perpendicular rather than n = integer as in the Rayleigh-Benard case. In addition, contrary to Rayleigh-Benard convection, in which the critical Rayleigh number is fixed, it is found that R super (cr) sub u is variable in the presence of a stable layer, (i.e., it depends on R(s)).

Savolainen, V.↗

NASA NTRS1988

Investigation of surface tension driven convection as a feasibility study for a micro-gravity experiment

The work performed for the feasibility study of a microgravity surface tension driven convection experiment was reviewed. An experimental investigation of the onset of convection in shallow fluid layers heated uniformly from below and cooled from above by an air layer was made. Results are discussed in relation to the formation of Benard cells. The onset of Rayleigh-Benard convection in thin fluid layers heated uniformly from below were studied experimentally. It was found that in thin fluid layers the onset of Rayleigh-Benard convection is preceded by subcritical convective motions. Secondly, it was found that the onset of Rayleigh-Benard convection in non-Boussinesq fluid layers takes place in the form of hexagonal cells at Rayleigh numbers larger than the critical Rayleigh number R sub C = 1708 which determines the onset of convection in Boussinesq fluid layers.

Koschmieder, E. L.↗

NASA NTRS1977

Studies of heat-source driven natural convection: A numerical investigation

Thermal convection driven by uniform volumetric energy sources was studied in a horizontal fluid layer bounded from above by a rigid, isothermal surface and from below by a rigid, zero heat-flux surface. The side walls of the fluid domain were assumed to be rigid and perfectly insulating. The computations were formally restricted to two-dimensional laminar convection but were carried out for a range of Rayleigh numbers which spans the regimes of laminar and turbulent flow. The results of the computations consists of streamline and isotherm patterns, horizontally averaged temperature distributions, and horizontally averaged Nusselt numbers at the upper surface. Flow and temperature fields do not exhibit a steady state, but horizontally averaged Nusselt numbers reach limiting, quasi-steady values for all Rayleigh numbers considered. Correlations of the Nusselt number in terms of the Rayleigh and Prandtl numbers were determined.

Emara, A. A.↗

NASA NTRS1982

Manufacturing in space: Fluid dynamics numerical analysis

Numerical computations were performed for natural convection in circular enclosures under various conditions of acceleration. It was found that subcritical acceleration vectors applied in the direction of the temperature gradient will lead to an eventual state of rest regardless of the initial state of motion. Supercritical acceleration vectors will lead to the same steady state condition of motion regardless of the initial state of motion. Convection velocities were computed for acceleration vectors at various angles of the initial temperature gradient. The results for Rayleigh numbers of 1000 or less were found to closely follow Weinbaum's first order theory. Higher Rayleigh number results were shown to depart significantly from the first order theory. Supercritical behavior was confirmed for Rayleigh numbers greater than the known supercritical value of 9216. Response times were determined to provide an indication of the time required to change states of motion for the various cases considered.

Robertson, S. J.↗

NASA NTRS1999

Magnetic Control of Convection in Electrically Nonconducting Fluids

Inhomogeneous magnetic fields exert a body force on electrically nonconducting, magnetically permeable fluids. This force can be used to compensate for gravity and to control convection. The effects of uniform and nonuniform magnetic fields on a laterally unbounded fluid layer heated from below or above are studied using a linear stability analysis of the Navier-Stokes equations supplemented by Maxwell's equations and the appropriate magnetic body force. For a uniform oblique field, the analysis shows that longitudinal rolls with axes parallel to the horizontal component of the field are the rolls most unstable to convection. The corresponding critical Rayleigh number and critical wavelength for the onset of such rolls are less than the well-known Rayleigh-Benard values in the absence of magnetic fields. Vertical fields maximize these deviations, which vanish for horizontal fields. Horizontal fields increase the critical Rayleigh number and the critical wavelength for all rolls except longitudinal rolls. For a nonuniform field, our analysis shows that the magnetic effect on convection is represented by a dimensionless vector parameter which measures the relative strength of the induced magnetic buoyancy force due to the applied field gradient. The vertical component of this parameter competes with the gravitational buoyancy effect, and a critical relationship between this component and the Rayleigh number is identified for the onset of convection. Therefore, Rayleigh-Benard convection in such fluids can be enhanced or suppressed by the field. It also shows that magnetothermal convection is possible in both paramagnetic and diamagnetic fluids. Our theoretical predictions for paramagnetic fluids agree with experiments. Magnetically driven convection in diamagnetic fluids should be observable even in pure water using current technology.

Huang, Jie↗