Search NASA⌕ Search

SEARCH · Search NASA

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 199 records · Page 11

Convection during Thermally Unstable Solidification of Pb-Sn in a Magnetic Field

Convection and macrosegregation in directionally solidified hypoeutectic Pb-38 wt pct Sn and hypereutectic Pb-64.5 wt pct Sn have been examined during upward and downward growth. Temperature fluctuations are observed along the length of the melt column during downward growth. With increasing Rayleigh number, these fluctuations change from none, to cyclic, to time periodic having multiple harmonics, and finally to random. At the higher convective driving force of 350 K temperature inversion, the transverse magnetic field decreased convective levels, strong random temperature fluctuations (flows) becoming smaller and periodic. The maximum field of 0.45 T was unable to completely eliminate convection. For the lower convective driving force of 150 K temperature inversion, the 0.05 T magnetic field decreased flows, and at 0.15 T, the field caused a dramatic decrease in the characteristic frequency of the temperature fluctuations, indicating a change in the nature of the flow, the waveform of the temperature fluctuations changing from sinusoidal to a pulsed wave. Temperature fluctuations and time delays between thermocouples were used to estimate flow velocities. Irrespective of the convection in the bulk melt (ahead of the mushy zone), longitudinal macrosegregation occurs only if the interdendritic melt mixes with the bulk melt.

Song, H.↗

Microgravity Experiments and Numerical Modeling of Rotating Buoyant Convection in a Spherical Shell with Latitudinal Thermal Gradients

The results of numerical model calculations are compared with space laboratory experiments for rotating, thermally driven flow in a hemispherical shell, and various flow regimes are described. A spherically symmetric body force, analogous to gravity, is imposed in the radial direction through the use of a dielectric fluid and an electrostatic potential difference across the gap. The spherical boundaries are maintained at constant temperature profiles (with the pole and the inner sphere being the warmer surfaces), and the equatorial wall is an insulator. Typical parameter combinations result in highly nonlinear, but laminar, flow. For weak enough buoyant forcing, the flow is axisymmetric. It consists of a single meridional cell, rising in warm latitudes (in this case, near the pole) and sinking in cool latitudes, with prograde flow in the equatorial region and near the inner hemisphere for other latitudes, and retrograde motion elsewhere. For fast rotation, the first transition due to the instability of this simple flow is also axisymmetric, consisting of rings of convection in the polar region, which propagate poleward. The first nonaxisymmetric convection occurs at Rayleigh numbers which increase with rotation rate. The form of the convection near the transition also depends upon the rate of rotation. Selected flow patterns near the transition as well as those beyond it are studied numerically. For those cases where there exist laboratory experiments with which to compare, the numerical and experimental results agree very well.

Leslie, Fred W.↗

Stability of a Horizontal Fluid Layer with Unsteady Heating from Below and Time-Dependent Body Force

The stability of a horizontal layer of fluid in an accelerated container heated unsteadily from below was investigated theoretically, assuming an incompressible fluid with small density changes resulting from heating. The critical Rayleigh numbers based on the over-all density differential are much higher than for the static case, are dependent only on the density distribution, and instantaneous value of the acceleration, and are independent of Prandtl number and rate of change of temperature and body force field. The initial motion corresponds to approximately the same cell shape as for the static case. Rate of transition of temperature perturbations from a stable to an unstable condition is proportional to the rate of increase of the temperature gradient in the region well removed from the walls; rate of transition of the slow motion is proportional to the Prandtl number and rate of increase of the body force field.

Goldstein, Arthur W.↗

Experimental Observations on a Low Strain Counter-Flow Diffusion Flame: Flow and Bouyancy Effects

Diffusion flames are of great interest in fire safety and many industrial processes. The counter-flow configuration provides a constant strain flow, and therefore is ideal to study the structure of diffusion flames. Most studies have concentrated on the high velocity, high strain limit, since buoyantly induced instabilities will disintegrate the planar flame as the velocity decreases. Only recently, experimental studies in microgravity conditions have begun to explore the low strain regimes. Numerical work has shown the coupling between gas phase reaction rates, soot reaction rates, and radiation. For these programs, size, geometry and experimental conditions have been chosen to keep the flame unaffected by the physical boundaries. When the physical boundaries can not be considered infinitely far from the reaction zone discrepancies arise. A computational study that includes boundary effects and accounts for the deviations occurring when the major potential flow assumptions are relaxed was presented by Borlik et al. This development properly incorporates all heat loss terms and shows the possibility of extinction in the low strain regime. A major constraint of studying the low strain regime is buoyancy. Buoyant instabilities have been shown to have a significant effect on the nature of reactants and heat transport, and can introduce instabilities on the flow that result in phenomena such as flickering or fingering. The counter-flow configuration has been shown to provide a flame with no symmetry disrupting instabilities for inlet velocities greater than 50 mm/s. As the velocity approaches this limit, the characteristic length of the experiment has to be reduced to a few millimetres so as to keep the Rayleigh number (Ra L = (β g 0 L 3 ∆T)/(𝑎v)) below 2000. In this work, a rectangular counter-flow burner was used to study a two-dimensional counter-flow diffusion flame. Flow visualisation and Particle Image Velocimetry served to describe the nature of the stagnation plane for strain rates smaller than 100 (1/s). These experiments were conducted with a non-reacting flow. Video images of a propane air diffusion flame were used to describe the behaviour of a diffusion flame in this regime. Flame geometry and pulsation frequency are described.

J. A. Sutula↗

Geoid Anomalies and Dynamic Topography from Time Dependent, Spherical Axisymmetric Mantle Convection

Geoid anomalies and dynamic topography are two important diagnostics of mantle convection. We present geoid and topography results for several time-dependent convection models in spherical axisymmetric geometry for Rayleigh numbers between 10(exp 6) and 10(exp 7) with depth-dependent viscosity and mixtures of bottom and internal heating. The models are strongly chaotic, with boundary layer instabilities erupting out of both thermal boundary layers. In some instances, instabilities from one boundary layer influence the development of instabilities in the other boundary layer. Such coupling between events at the top and bottom of the mantle has been suggested to play a role in a mid-Cretaceous episode of enhanced volcanism in the Pacific. These boundary layer instabilities produce large temporal variations in the geoid anomalies and dynamic nd to the topography associated with the convection. The amplitudes of these fluctuations depend on the detailed model parameter,.% it of this but fluctuations of 30-50% relative to the time-averaged geoid and topography are common. The convective planform is strongly sensitive to the specific initial conditions. Convection cells with larger aspect ratio tend to have larger fractional fluctuations in their geoid and topography amplitudes, because boundary layer instabilities have more time to develop in long cells. In some instances, we observe low-amplitude topographic highs adjacent to the topographic lows produced by cold downwellings. We discuss applications of these results to several situations, including the temporal variability of m basis. hotspots such as Hawaii, the topography of subduction zone outer rises, and the topography of coronae on Venus.

Kiefer, Walter S.↗

Crystal Growth and Fluid Mechanics Problems in Directional Solidification

Our work in directional solidification has been in the following areas: (1) Dynamics of dendrites including rigorous mathematical analysis of the resulting equations; (2) Examination of the near-structurally unstable features of the mathematically related Hele-Shaw dynamics; (3) Numerical studies of steady temperature distribution in a vertical Bridgman device; (4) Numerical study of transient effects in a vertical Bridgman device; (5) Asymptotic treatment of quasi-steady operation of a vertical Bridgman furnace for large Rayleigh numbers and small Biot number in 3D; and (6) Understanding of Mullins-Sererka transition in a Bridgman device with fluid dynamics is accounted for.

Tanveer, Saleh A.↗

Microgravity Science Glovebox Investigations SUBSA

Solidification Using a Baffle in Sealed Ampoules (SUBSA) is a Microgravity Science Glovebox Investigation manifested for the UF2 flight, on the U.S. Orbiter 111, to the International Space Station (ISS). SUBSA complements the "parent" flight investigation CG13 (Space-and Groundbased Crystal Growth Using a Baffle). During directional solidification, the disk-shaped baffle acts as a partition, creating a small melt zone at the solid-liquid interface. As a result, the level of buoyancy-driven convection at the interface is significantly reduced. In space, the baffle will reduce convection driven by residual micro acceleration. The baffle reduces the Rayleigh number (Ra) of the melt by a factor of 103. The combined effect of the baffle and microgravity will yield a reduction in Ra by a factor of 107 to 109 approaching effectively the acceleration conditions in "nanogravity". The results of ground based tests and numerical modeling will be presented. The furnace for directional solidification (flight hardware and the ground unit) was developed by Tec-Masters Inc. The flight ampoules were produced jointly at Rensselaer, Tec-Masters Inc. and Crystallod Inc.

Ostrogorsky, A.↗

Differential Flatness and Cooperative Tracking in the Lorenz System

In this paper the control of the Lorenz system for both stabilization and tracking problems is studied via feedback linearization and differential flatness. By using the Rayleigh number as the control, only variable physically tunable, a barrier in the controllability of the system is incidentally imposed. This is reflected in the appearance of a singularity in the state transformation. Composite controllers that overcome this difficulty are designed and evaluated. The transition through the manifold defined by such a singularity is achieved by inducing a chaotic response within a boundary layer that contains it. Outside this region, a conventional feedback nonlinear control is applied. In this fashion, the authority of the control is enlarged to the whole. state space and the need for high control efforts is mitigated. In addition, the differential parametrization of the problem is used to track nonlinear functions of one state variable (single tracking) as well as several state variables (cooperative tracking). Control tasks that lead to integrable and non-integrable differential equations for the nominal flat output in steady-state are considered. In particular, a novel numerical strategy to deal with the non-integrable case is proposed. Numerical results validate very well the control design.

Crespo, Luis G.↗

Suppression of Buoyancy in Gaseous Media at High Temperatures

Consider a rectangular box filled with a fluid having a heated bottom and a cold top surface, and insulated side-walls (Benard problem). As the temperature difference between the horizontal top and bottom surfaces increases, a critical condition, defined quantitatively by the Rayleigh number, is reached beyond which density stratification can no longer be sustained by conduction and the fluid disrupts from its stable, quiescent state into an unstable, convective mode in which lighter and heavier gas mix. This paper suggests that such a statement is not necessarily true for gaseous media under normalized temperature differences that are much larger than justifiable for the Boussinesq approximation! In fact, there may be situations where a system cannot ever be made unstable with respect to the onset on buoyant convection no matter how large the temperature (density) difference becomes at a given pressure even under normal gravity! This unexpected behavior is primarily attributed to highly temperature-sensitive kinematic viscosity which counteracts the tendency toward instability and dampens convection by making the gas more viscous at higher temperatures. This compensation of the buoyant force by the viscous force exhibits itself by the formation of a peak hot-surface temperature beyond which a system will tend to be more stable as the hot-surface temperature increases.

Gokoglu, Suleyman A.↗

Mass Redistribution in the Core and Time-varying Gravity at the Earth's Surface

The Earth's liquid outer core is in convection, as suggested by the existence of the geomagnetic field in much of the Earth's history. One consequence of the convection is the redistribution of mass resulting from relative motion among fluid parcels with slightly different densities. This time dependent mass redistribution inside the core produces a small perturbation on the gravity field of the Earth. With our numerical dynamo solutions, we find that the mass redistribution (and the resultant gravity field) symmetric about the equator is much stronger than that anti-symmetric about the equator. In particular, J(sub 2) component is the strongest. In addition, the gravity field variation increases with the Rayleigh number that measures the driving force for the geodynamo in the core. With reasonable scaling from the current dynamo solutions, we could expect that at the surface of the Earth, the J(sub 2) variation from the core is on the order of l0(exp -16)/year relative to the mean (i.e. spherically symmetric) gravity field of the Earth. The possible shielding effect due to core-mantle boundary pressure variation loading is likely much smaller and is therefore negligible. Our results suggest that time-varying gravity field perturbation due to core mass redistribution may be measured with modem space geodetic observations, which will result a new means of detecting dynamical processes in the Earth's deep interior.

Kuang, Wei-Jia↗

Cool Flames in Propane-Oxygen Premixtures at Low and Intermediate Temperatures at Reduced-Gravity

The Cool Flame Experiment aims to address the role of diffusive transport on the structure and the stability of gas-phase, non-isothermal, hydrocarbon oxidation reactions, cool flames and auto-ignition fronts in an unstirred, static reactor. These reactions cannot be studied on Earth where 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. On Earth, reactions with associated Rayleigh numbers (Ra) less than the critical Ra for onset of convection (Ra(sub cr) approx. 600) cannot be achieved in laboratory-scale vessels for conditions representative of nearly all low-temperature reactions. In fact, the Ra at 1g ranges from 10(exp 4) - 10(exp 5) (or larger), while at reduced-gravity, these values can be reduced two to six orders of magnitude (below Ra(sub cr)), depending on the reduced-gravity test facility. Currently, laboratory (1g) and NASA s KC-135 reduced-gravity (g) aircraft studies are being conducted in parallel with the development of a detailed chemical kinetic model that includes thermal and species diffusion. Select experiments have also been conducted at partial gravity (Martian, 0.3gearth) aboard the KC-135 aircraft. This paper discusses these preliminary results for propane-oxygen premixtures in the low to intermediate temperature range (310- 350 C) at reduced-gravity.

Pearlman, Howard↗

Convection in a Very Compressible Fluid: Comparison of Simulations With Experiments

The time profile (Delta)T(t) of the temperature difference, measured across a very compressible fluid layer of supercritical He-3 after the start of a heat flow, shows a damped oscillatory behavior before steady state convection is reached. The results for (Delta)T(t) obtained from numerical simulations and from laboratory experiments are compared over a temperature range where the compressibility varies by a factor of approx. = 40. First the steady-state convective heat current j(sup conv) as a function of the Rayleigh number R(alpha) is presented, and the agreement is found to be good. Second, the shape of the time profile and two characteristic times in the transient part of (Delta)T(t) from simulations and experiments are compared, namely: 1) t(sub osc), the oscillatory period and 2) t(sub p), the time of the first peak after starting the heat flow. These times, scaled by the diffusive time tau(sub D) versus R(alpha), are presented. The agreement is good for t(sup osc)/tau(sub D), where the results collapse on a single curve showing a powerlaw behavior. The simulation hence confirms the universal scaling behavior found experimentally. However for t(sub p)/tau(sub D), where the experimental data also collapse on a single curve, the simulation results show systematic departures from such a behavior. A possible reason for some of the disagreements, both in the time profile and in t(sub p) is discussed. In the Appendix a third characteristic time, t(sub m), between the first peak and the first oscillation minimum is plotted and a comparison between the results of experiments and simulations is made.

Meyer, H.↗

New Species of Fire Discovered: Fingering Flamelets Form a Dynamic Population

Poets and artists have long used fire as a metaphor for life. At the NASA Glenn Research Center, recent experiments in a subcritical Rayleigh number flow channel demonstrated that this analogy holds up surprisingly well when tools developed to characterize a biological population are applied to a class of fire that occurs in near-extinction, weakly convective environments (such as microgravity) or in vertically confined spaces (such as our apparatus). Under these conditions, the flame breaks into numerous 'flamelets" that form a Turing-type reaction-diffusion fingering pattern as they spread across the fuel. It is standard practice on U.S. spacecraft for the astronaut crew to turn off the ventilation to help extinguish a fire, both to eliminate the fresh oxygen supply and to reduce the distribution of the smoke. When crew members think that the fire is fully extinguished, they reactivate the ventilation system to clear the smoke. However, some flamelets can survive, and our experiments have demonstrated that flamelets quickly grow into a large fire when ventilation increases.

Olson, Sandra L.↗

Parallel Computation for Natural Convection in Cavities

Parallel computation for thermal convective flows in cavities with adiabatic horizontal boundaries and driven by differential heating of the two vertical end walls, is investigated using supercomputers. A parallel computation code has been implemented using a finite-difference method with a multigrid elliptic solver and a Dufort-Frankel scheme. The domain decomposition techniques are discussed in detail. The parallel code is numerically stable, computationally efficient, and portable to various parallel architectures which support either PVM or NX libraries for communications. Finally, numerical results for various Rayleigh numbers and Prandtl numbers are presented.

convection simulation parallel computation↗

Thermal Vibrational Convection in a Two-phase Stratified Liquid

The response of a two-phase stratified liquid system subject to a vibration parallel to an imposed temperature gradient is analyzed using a hybrid thermal lattice Boltzmann method (HTLB). The vibrations considered correspond to sinusoidal translations of a rigid cavity at a fixed frequency. The layers are thermally and mechanically coupled. Interaction between gravity-induced and vibration-induced thermal convection is studied. The ability of applied vibration to enhance the flow, heat transfer and interface distortion is investigated. For the range of conditions investigated, the results reveal that the effect of vibrational Rayleigh number and vibrational frequency on a two-phase stratified fluid system is much different than that for a single-phase fluid system. Comparisons of the response of a two-phase stratified fluid system with a single-phase fluid system are discussed.

Chang, Qingming↗

Gravity Scaling of a Power Reactor Water Shield

A similarity analysis on a water-based reactor shield examined the effect of gravity on free convection between a reactor shield inner and outer vessel boundaries. Two approaches established similarity between operation on the Earth and the Moon: 1) direct scaling of Rayleigh number equating gravity-surface heat flux products, 2) temperature difference between the wall and thermal boundary layer held constant. Nusselt number for natural convection (laminar and turbulent) is assumed of form Nu = CRa(sup n).

Reid, Robert S.↗

Gravity Scaling of a Power Reactor Water Shield

Water based reactor shielding is being considered as an affordable option for use on initial lunar surface power systems. Heat dissipation in the shield from nuclear sources must be rejected by an auxiliary thermal hydraulic cooling system. The mechanism for transferring heat through the shield is natural convection between the core surface and an array of thermosyphon radiator elements. Natural convection in a 100 kWt lunar surface reactor shield design has been previously evaluated at lower power levels (Pearson, 2007). The current baseline assumes that 5.5 kW are dissipated in the water shield, the preponderance on the core surface, but with some volumetric heating in the naturally circulating water as well. This power is rejected by a radiator located above the shield with a surface temperature of 370 K. A similarity analysis on a water-based reactor shield is presented examining the effect of gravity on free convection between a radiation shield inner vessel and a radiation shield outer vessel boundaries. Two approaches established similarity: 1) direct scaling of Rayleigh number equates gravity-surface heat flux products, 2) temperature difference between the wall and thermal boundary layer held constant on Earth and the Moon. Nussult number for natural convection (laminar and turbulent) is assumed of form Nu = CRa(sup n). These combined results estimate similarity conditions under Earth and Lunar gravities. The influence of reduced gravity on the performance of thermosyphon heat pipes is also examined.

Reid, Robert S.↗

Decadal Polar Motion of the Earth Excited by the Convective Outer Core From Geodynamo Simulations

Long time geodetic observation records show that the orientation of the Earth's rotation axis with respect to the terrestrial reference frame, or polar motion, changes on a broad range of timescales. Apart from external torques from the luni-solar tides, these changes are excited by interactions among different components of the Earth system. The convective fluid outer core has long been conjectured a likely contributor to the observed polar motion on timescales upward of decades, such as the 30 year Markowitz wobble. We investigated the electromagnetic coupling scenario across the core-mantle boundary via numerical geodynamo simulation for different geodynamo parameters (Rayleigh numbers and magnetic Rossby numbers). Our simulated polar motion varies strongly with the dynamo parameters, while its excitation on decadal timescales appear to converge asymptotically within the adopted range of numerical Rossby numbers. Three strongest asymptotic modes emerge from numerical results, with periods around 30, 40, and 60 years for the prograde excitation and around 24, 30, and 60 years for the retrograde excitation. Their amplitudes are all larger than 5 times 10 (sup minus 8), or approximately 10 milliseconds of arc.The results suggest that the electromagnetic core-mantle coupling could explain a substantial portion, if not all, of the observed decadal polar motion. In particular, the predicted 60 year polar motion deserves special attention for future observations and studies.

Kuang, W.↗