Engineering Papers⌕ Search

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.

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

At least 55 records · Page 3

Manufacturing in space: Fluid dynamics numerical analysis

Natural convection in a spherical container with cooling at the center was numerically simulated using the Lockheed-developed General Interpolants Method (GIM) numerical fluid dynamic 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 growing experiments to be performed on board Spacelab. Results were obtained for a range of Rayleigh numbers from 25 to 10,000. For a temperature difference of 10 C from the cooling sting at the center to the container surface, and a gravitional loading of 0.000001 g a computed maximum fluid velocity of about 2.4 x 0.00001 cm/sec was reached after about 250 sec. The computed velocities were found to be approximately proportional to the Rayleigh number over the range of Rayleigh numbers investigated.

Robertson, S. J.↗

A numerical investigation of thermal convection in a heat-generating fluid layer

Finite difference solutions of the equations governing thermal convection driven by uniform volumetric energy sources are presented for two-dimensional flows in a rectangular domain. The boundary conditions are a rigid (i.e, zero slip), zero heat-flux lower surface, rigid adiabatic sides, and either a rigid or free (i.e., zero shear) isothermal upper surface. Computations are carried out for Prandtl numbers from 0.05 to 20 and Rayleigh numbers from 5 x 10 to the 4th to 5 x 10 to the 8th. Nusselt numbers and average temperature profiles within the layer are in good agreement with experimental data for rigid-rigid boundaries. For rigid-free boundaries, Nusselt numbers are larger than in the former case. The structure of the flow and temperature fields in both cases is dominated by rolls, except at larger Rayleigh numbers where large-scale eddy transport occurs. Generally, low velocity upflows over broad regions of the layer are balanced by higher velocity downflows when the flow exhibits a cellular structure. The hydrodynamic constraint at the upper surface and the Prandtl number are found to influence only the detailed nature of flow and temperature fields. No truly steady velocity and temperature fields are found despite the fact that average Nusselt numbers reach steady values.

Emara, A. A.↗

Experimental investigation of the Marangoni effect on the stability of a double-diffusive layer

Stability experiments were carried out in 4-cm-thick, salt-stratified fluid layer by heating from below and cooling from above. The bottom boundary was rigid while the top was either free or rigid. The initial solute Rayleigh number varied from 2.5 x 10(exp 6) to 4.6 x 10(exp 7). For the rigid-free case, at initial solute Rayleigh numbers R(sub s) greater than 5.4 x 10(exp 6), thermal Marangoni instabilities were observed to onset along the free surface at a relatively low thermal Rayleigh number, R(sub t). The convection was very weak, and it had almost no effect on the concentration and temperature distributions. Double-diffusive instabilities along the top free surface were observed to onset at a higher R(sub t), with much stronger convection causing changes in the concentration and temperature distributions near the top. At a yet higher R(sub t), double-diffusive convection was observed to onset along the bottom boundary. Fluid motion in the layer then evolved into fully developed thermal convection of a homogeneous fluid without any further increase in the imposed Delta T. For layers with R(sub s) less than 5.4 x 10(exp 6), Marangoni and double-diffusive instabilities onset simultaneously along the free surface first, while double-diffusive instabilities along the bottom wall onset at a higher R(sub t).

Tanny, Josef↗

The onset of time-dependent convection in spherical shells as a clue to chaotic convection in the earth's mantle

This work presents a detailed numerical study of the dynamical behavior of convection in a spherical shell, as applied to mantle convection. From both two-dimensional (120 radial and 360 tangential points) and three-dimensional (60 radial levels and spherical harmonics up to order and degree L = 33, m = 33), it is shown that for a spherical shell (with inner-to-outer radii ratio eta = 0.62) convection becomes time-dependent, with l = 2 dominating, at a Rayleigh number of about 31 times supercritical for a constant-viscosity, base-heated configuration. This secondary instability is characterized by oscillatory time dependence, with higher frequencies involved, at slightly higher Rayleigh numbers. In illustrating the onset of time dependence, the analysis is extended to show that the onset of weak turbulence in spherical-shell convection takes place at about 60 times the critical Rayleigh number via a quasi-periodic mode.

Machetel, Philippe↗

Convective motions in a spherical shell

We compute the axisymmetric convective motions that exist in a spherical shell heated from below with inner to outer radius ratio equal to 0.5. The boundaries are stress-free and gravity is directly proportional to radius. Accurate solutions at large Rayleigh numbers, O(100000), are made feasible by a spectral method that employs diagonal-mode truncation. By examining the stability of axisymmetric motions it is inferred that the preferred form of convection varies dramatically according to the value of the Rayleigh number. While axisymmetric motions with different patterns may exist for modestly nonlinear convection, only a single motion persists at sufficiently large values of the Rayleigh number. This circulation is symmetric about the equator and has two meridional cells with rising motion at the poles. Instability of this single axisymmetric motion determines that the preferred pattern of three-dimensional convection has one azimuthal wave.

Zebib, A.↗

Time-dependent convection with non-Newtonian viscosity

A numerical model of bottom-heated, two-dimensional convection in boxes of aspect ratios 2.5 and 4.0 was used to study the differences in time-dependent convection between Newtonian and stress-dependent viscosity. The onset of time dependence due to boundary layer instability is found at approximately the same effective Rayleigh number for both rheologies. However, with increasing Rayleigh number, the temporal and spatial fluctuations in the flow field become much more pronounced with non-Newtonian rheology; also, the tendency for breakup of long cells into smaller ones is stronger. The presence of stress-dependent rheology can cause long-wavelength lateral viscosity variations up to an order of magnitude; this result could have strong implications for the interpretation of mantle viscosity from postglacial rebound.

Christensen, Ulrich R.↗

The Cool Flames Experiment

A space-based experiment is currently under development to study diffusion-controlled, gas-phase, low temperature oxidation reactions, cool flames and auto-ignition in an unstirred, static reactor. At Earth's gravity (1g), natural convection due to self-heating during the course of slow reaction dominates diffusive transport and produces spatio-temporal variations in the thermal and thus species concentration profiles via the Arrhenius temperature dependence of the reaction rates. Natural convection is important in all terrestrial cool flame and auto-ignition studies, except for select low pressure, highly dilute (small temperature excess) studies in small vessels (i.e., small Rayleigh number). On Earth, natural convection occurs when the Rayleigh number (Ra) exceeds a critical value of approximately 600. Typical values of the Ra, associated with cool flames and auto-ignitions, range from 104-105 (or larger), a regime where both natural convection and conduction heat transport are important. When natural convection occurs, it alters the temperature, hydrodynamic, and species concentration fields, thus generating a multi-dimensional field that is extremely difficult, if not impossible, to be modeled analytically. This point has been emphasized recently by Kagan and co-workers who have shown that explosion limits can shift depending on the characteristic length scale associated with the natural convection. Moreover, natural convection in unstirred reactors is never "sufficiently strong to generate a spatially uniform temperature distribution throughout the reacting gas." Thus, an unstirred, nonisothermal reaction on Earth does not reduce to that generated in a mechanically, well-stirred system. Interestingly, however, thermal ignition theories and thermokinetic models neglect natural convection and assume a heat transfer correlation of the form: q=h(S/V)(T(bar) - Tw) where q is the heat loss per unit volume, h is the heat transfer coefficient, S/V is the surface to volume ratio, and (T(bar) - Tw ) is the spatially averaged temperature excess. This Newtonian form has been validated in spatially-uniform, well-stirred reactors, provided the effective heat transfer coefficient associated with the unsteady process is properly evaluated. Unfortunately, it is not a valid assumption for spatially-nonuniform temperature distributions induced by natural convection in unstirred reactors. "This is why the analysis of such a system is so difficult." Historically, the complexities associated with natural convection were perhaps recognized as early as 1938 when thermal ignition theory was first developed. In the 1955 text "Diffusion and Heat Exchange in Chemical Kinetics", Frank-Kamenetskii recognized that "the purely conductive theory can be applied at sufficiently low pressure and small dimensions of the vessel when the influence of natural convection can be disregarded." This was reiterated by Tyler in 1966 and further emphasized by Barnard and Harwood in 1974. Specifically, they state: "It is generally assumed that heat losses are purely conductive. While this may be valid for certain low pressure slow combustion regimes, it is unlikely to be true for the cool flame and ignition regimes." While this statement is true for terrestrial experiments, the purely conductive heat transport assumption is valid at microgravity (mu-g). Specifically, buoyant complexities are suppressed at mu-g and the reaction-diffusion structure associated with low temperature oxidation reactions, cool flames and auto-ignitions can be studied. Without natural convection, the system is simpler, does not require determination of the effective heat transfer coefficient, and is a testbed for analytic and numerical models that assume pure diffusive transport. In addition, mu-g experiments will provide baseline data that will improve our understanding of the effects of natural convection on Earth.

Pearlman, Howard↗

Surface topography due to convection in a variable viscosity fluid - Application to short wavelength gravity anomalies in the central Pacific Ocean

Finite difference calculations of thermal convection in a fluid layer with a viscosity exponentially decreasing with temperature are performed in the context of examining the topography and gravity anomalies due to mantle convection. The surface topography and gravity anomalies are shown to be positive over regions of ascending flow and negative over regions of descending flow; at large Rayleigh numbers the amplitude of surface topography is inferred to depend on Rayleigh number to the power of 7/9. Compositional stratifications of the mantle is proposed as a mechanism for confining small-scale convection to a thin layer. A comparative analysis of the results with other available models is included.

Lin, J.↗

Influence of non steady gravity on natural convection during micro-gravity solidification of semiconductors. I - Time scale analysis. II - Implications for crystal growth experiments

Consideration is given to the influence of temporal variations in the magnitude of gravity on natural convection during unidirectional solidification of semiconductors. It is shown that the response time to step changes in g at low Rayleigh numbers is controlled by the momentum diffusive time scale. At higher Rayleigh numbers, the response time to increases in g is reduced because of inertial effects. The degree of perturbation of flow fields by transients in the gravitational acceleration on the Space Shuttle and the Space Station is determined. The analysis is used to derive the requirements for crystal growth experiments conducted on low duration low-g vehicles. Also, the effectiveness of sounding rockets and KC-135 aircraft for microgravity experiments is examined.

Griffin, P. R.↗

Bifurcation and stability of low-order steady flows in horizontally and vertically forced convection

A nonlinear spectral model of two-dimensional, shallow Boussinesq convection which responds to heating in both the horizontal and vertical directions is examined. The governing partial differential system is converted to an infinite set of ordinary differential equations and truncated to a small set to permit detailed study of the number and types of transitions from one flow configuration to another. The Hadley number and the Rayleigh number are defined as the horizontal and vertical thermal forcing mechanisms, respectively, for inclusion in the nonlinear spectral model, which is composed of three equations. The model is then used to describe steady states, linearly stable solutions, and balancing factors in unstable stratification. The number and the distribution of the steady states are found to be qualitatively independent of the aspect ratio and the Prandtl number.

Yost, D. A.↗

Thresholds for the onset of fluid and magnetofluid turbulence

Linear stability calculations conducted in plasma physics are based on the theory of hydrodynamic stability of neutral fluids. The present investigation is concerned with the validity of procedures based on linear stability analysis in six much-studied situations. All have the property that the fluid goes from laminar to turbulent at critical values of some dimensionless number, such as the Reynolds number or the Rayleigh number. The question is, at what threshold values of the dimensionless number do the unstable motions set in, and can these thresholds be predicted by a linear analysis of the stability of the laminar state as the threshold is approached from the stable side. In three cases linear stability analysis clearly seems to fail. These cases include Plane Poiseuille Flow, Plane Couette Flow, and Cylindrical Pipe Flow (Hagen-Poiseuille Flow). Situations in which predictions provided by linear stability analysis are correct are related to Rotating Couette Flow, Thermally-Driven Convection, and Instability of Laminar Boundary Layers.

Montgomery, D.↗

Numerical simulations of soft and hard turbulence - Preliminary results for two-dimensional convection

Results on the transition from soft to hard turbulence in simulations of two-dimensional Boussinesq convection are reported. The computed probability densities for temperature fluctuations are exponential in form in both soft and hard turbulence, unlike what is observed in experiments. In contrast, a change is obtained in the Nusselt number scaling on Rayleigh number in good agreement with the three-dimensional experiments.

Deluca, E. E.↗

Numerical investigation on Benard instability in a finite liquid layer

A numerical procedure for directly simulating the Benard-Marangoni instabilities (B-M-I) in a bounded liquid layer is presented in this paper. The procedure consists of applying a finite amplitude disturbance to the basic static state, and then integrating the Navier-Stokes equations to determine whether the disturbance will die down or will reach a state of finite-strength steady convection. The critical Marangoni number (Mac) for the onset of B-M-I can thus be determined and can be correlated as a function of the aspect ratio (Ar), Prandtl number (Pr), and Rayleigh number (Ra). The Biot number (B) between the liquid and the air is analyzed to approximate the heat transfer condition along the free surface. A 2D calculation is performed to investigate the effect of various initial disturbances, and the Mac is determined for Ar = 2, Ra = 0, and Pr = 0.7. Current results show that disturbances of different nature and amplitude have little effect on Mac. The Mac determined in this study also clearly demonstrates the dominant effect of the sidewalls.

Duh, J. C.↗

Flight Experiment to Study Double-Diffusive Instabilities in Silver-Doped Lead Bromide Crystals

A detailed study on the effect of convection on crystal quality was carried out by growing lead bromide crystals in transparent Bridgman furnace. Direct observations were made on the solid-liquid interface and a new kind of instability was observed. This could be explained on the basis of toroidal flow in the AgBr-doped lead bromide sample. With the increasing translation velocity, the interface changed from flat to depressed, and then formed a cavity in the center of the growth tube. The crystal grown at the lowest thermal Rayleigh number showed the highest quality and crystal grown at the largest thermal Rayleigh number showed the worst quality. Numerical studies were carried out to provide a framework for interpreting the observed convective and morphological instabilities, and to determine the critical (limiting) concentration of dopant for a particular growth velocity and gravity level. Theoretical instability diagrams were compared with data obtained from the experimental studies. These studies provided basic data on convective behavior in doped lead bromide crystals grown by the commercially important Bridgman process.

Singh, N. B.↗

Time-dependent solutions of multimode convection equations

Truncated modal equations are used to study the time evolution of thermal convection. In the Boussinesq approximation these nonlinear equations are obtained by expanding the fluctuating velocity and temperature fields in a finite set of planforms of the horizontal coordinates. Numerical studies dealing with two or three modes with triad interactions are discussed. Rich time dependence was found in these cases: periodic and aperiodic solutions can be obtained, along with various steady solutions. Three-mode solutions reproduce the qualitative appearance of spoke-pattern convection as observed in experiments at high Prandtl numbers. Though the values of the periods of the time-dependent solutions do not agree with those of the experiments, their variation with Rayleigh number compares favorably. Except at the highest Rayleigh number considered (10,000,000), the theoretical Nusselt numbers agree well with experiment.

Toomre, J.↗

Steady bimodal convection in a cylinder at large Prandtl numbers

Steady bimodal convection of an infinite Prandtl-number Boussinesq fluid in a cylinder is considered. An asymptotic analysis similar to the one used by Buell and Catton (1986) for axisymmetric convection yields a solvability condition that determines the radial wavenumber. The analysis is valid for convection far away from the origin, the lateral boundary, and any pattern dislocations. The azimuthal wave number is treated as a parameter, although in real systems it is dependent on the initial and boundary conditions. Results are presented for Rayleigh numbers between 14,000 and 60,000, and for azimuthal wave numbers between 5 and 7. It is shown that for increasing Rayleigh numbers, the selected radial wave number and the heat transfer tend to become independent of the azimuthal wave number. No quantitative experimental data are available, but one qualitative comparison is good.

Buell, Jeffrey C.↗

Effect of enclosure shape on natural convection velocities

A numerical analysis was performed to compare natural convection velocities in two dimensional enclosures of various shape. The following shapes were investigated: circle, square, horizontal and upright 2 x 1 aspect ratio rectangles, horizontal and upright half circles, diamond. In all cases, the length scale in the various dimensionless parameters, such as Rayleigh number, is defined as the diameter of the equal area circle. Natural convection velocities were calculated for Rayleigh numbers of 1000 and 5000 with the temperature difference taken to be across (1) the maximum horizontal dimension, (2) the median horizontal line (line through centroid) and (3) the horizontal distance such that the temperature gradient is the same for shapes of equal area. For the class of shapes including the square, upright half circle and upright rectangle, the computed velocities were found to agree very closely with that of the equal area circle when the temperature difference is taken to be across the maximum horizontal dimension (condition (a)). The velocities for the horizontal rectangle and half circle were found to be approximately one half that of the equal area circle for the same condition. Better overall agreement among all shapes was obtained by setting the temperature difference across a distance such that the temperature gradients were equal for shapes of equal area.

Robertson, S. J.↗

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.↗