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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

Radial segregation induced by natural convection and melt/solid interface shape in vertical Bridgman growth

The roles of natural convection in the melt and the shape of the melt/solid interface on radial dopant segregation are analyzed for a prototype of vertical Bridgman crystal growth system by finite element methods that solve simultaneously for the velocity field in the melt, the shape of the solidification isotherm, and the temperature distribution in both phases. Results are presented for crystal and melt with thermophysical properties similar to those of gallium-doped germanium in Bridgman configurations with melt below (thermally destabilizing) and above (stabilizing) the crystal. Steady axisymmetric flow are classified according to Rayleigh number as either being nearly the growth velocity, having a weak cellular structure or having large amplitude cellular convention. The flows in the two Bridgman configurations are driven by different temperature gradients and are in opposite directions. Finite element calculations for the transport of a dilute dopant by these flow fields reveal radial segregation levels as large as sixty percent of the mean concentration. Segregation is found most severe at an intermediate value of Rayleigh number above which the dopant distribution along the interface levels as the intensity of the flow increases.

Chang, C. J.↗
NASA NTRS1975

Studies of heat source driven natural convection

Natural convection energy transport in a horizontal layer of internally heated fluid was measured for Rayleigh numbers from 1890 to 2.17 x 10 to the 12th power. The fluid layer is bounded below by a rigid zero-heat-flux surface and above by a rigid constant-temperature surface. Joule heating by an alternating current passing horizontally through the layer provides the uniform volumetric energy source. The overall steady-state heat transfer coefficient at the upper surface was determined by measuring the temperature difference across the layer and power input to the fluid. The correlation between the Nusselt and Rayleigh numbers for the data of the present study and the data of the Kulacki study is given.

Kulacki, F. A.↗
NASA NTRS2014

Transient Convection Due to Imposed Heat Flux: Application to Liquid-Acquisition Devices

A model problem is considered that addresses the effect of heat load from an ambient laboratory environment on the temperature rise of liquid nitrogen inside an enclosure. This model has applications to liquid acquisition devices inside the cryogenic storage tanks used to transport vapor-free propellant to the main engine. We show that heat loads from Q = 0.001 to 10 W, with corresponding Rayleigh numbers from Ra = 109 to 1013, yield a range of unsteady convective states and temperature rise in the liquid. The results show that Q = 1 to 10 W (Ra = 1012 to 1013) yield temperature distributions along the enclosure height that are similar in trend to experimental measurements. Unsteady convection, which shows selfsimilarity in its planforms, is predicted for the range of heat-load conditions. The onset of convection occurs from a free-convection-dominated base flow that becomes unstable against convective instability generated at the bottom of the enclosure while the top of the enclosure is convectively stable. A number of modes are generated with small-scale thermals at the bottom of the enclosure in which the flow selforganizes into two symmetric modes prior to the onset of the propagation of the instability. These symmetric vertical modes transition to asymmetric modes that propagate as a traveling-wave-type motion of convective modes and are representative of the asymptotic convective state of the flow field. Intense vorticity production is created in the core of the flow field due to the fact that there is shear instability between the vertical and horizontal modes. For the higher Rayleigh numbers, 1012 to 1013, there is a transition from a stationary to a nonstationary response time signal of the flow and temperature fields with a mean value that increases with time over various time bands and regions of the enclosure.

Duval, Walter M. B.↗
NASA NTRS2003

Flow Visualization of Low Prandtl Number Fluids using Electrochemical Measurements

It is well established that residual flows exist in contained liquid metal processes. In 1-g processing, buoyancy forces often drive these flows and their magnitudes can be substantial. It is also known that residual flows can exist during microgravity processing, and although greatly reduced in magnitude, they can influence the properties of the processed materials. Unfortunately, there are very few techniques to visualize flows in opaque, high temperature liquid metals, and those available are not easily adapted to flight investigation. In this study, a novel technique is developed that uses liquid tin as the model fluid and solid-state electrochemical cells constructed from Yttria-Stabilized Zirconia (YSZ) to establish and measure dissolved oxygen boundary conditions. The melt serves as a common electrode for each of the electrochemical cells in this design, while independent reference electrodes are maintained at the outside surfaces of the electrolyte. By constructing isolated electrochemical cells at various locations along the container walls, oxygen is introduced or extracted by imposing a known electrical potential or passing a given current between the melt and the reference electrode. This programmed titration then establishes a known oxygen concentration boundary condition at the selected electrolyte-melt interface. Using the other cells, the concentration of oxygen at the electrolyte-melt interface is also monitored by measuring the open-circuit potentials developed between the melt and reference electrodes. Thus the electrochemical cells serve to both establish boundary conditions for the passive tracer and sense its path. Rayleigh-Benard convection was used to validate the electrochemical approach to flow visualization. Thus, a numerical characterization of the second critical Rayleigh numbers in liquid tin was conducted for a variety of Cartesian aspect ratios. The extremely low Prandtl number of tin represents the lowest value studied numerically. Additionally, flow field oscillations are visualized and the effect of tilt on convecting systems is quantified. Experimental studies of the effect of convection in liquid tin are presented. Three geometries are studied: (1) double electrochemical cell with vertical concentration gradients; (2) double cell with horizontal concentration gradients; and (3) multiple cells with vertical temperature gradients. The first critical Rayleigh number transition is detected with geometry (1) and it is concluded that current measurements are not as affected by convection as EMF measurements. The system is compared with numerical simulations in geometry (2), and oscillating convection is detected with geometry (3).

Crunkleton, D.↗
NASA NTRS1996

Flow Transitions in a Rotating Magnetic Field

Critical Rayleigh numbers have been measured in a liquid metal cylinder of finite height in the presence of a rotating magnetic field. Several different stability regimes were observed, which were determined by the values of the Rayleigh and Hartmann numbers. For weak rotating magnetic fields and small Rayleigh numbers, the experimental observations can be explained by the existence of a single non-axisymmetric meridional roll rotating around the cylinder, driven by the azimuthal component of the magnetic field. The measured dependence of rotational velocity on magnetic field strength is consistent with the existence of laminar flow in this regime.

Volz, M. P.↗
NASA NTRS1990

Boundary layer roll circulations during FIRE

The probable mechanism underlying the development of boundary layer roll circulations are studied using wind and temperature profiles measured by the National Center for Atmospheric Research (NCAR) Electra during the stratocumulus phase of the First ISCCP Regional Experiment (FIRE). The expected, or preferred, roll orientations, horizontal wavelengths, and propagation periods are determined by finding the minimum values of the dynamic and thermodynamic forcing parameters, which here are the eddy Reynolds number (Re) and moist Rayleigh number (Ra sub m). These minimum values depend on the height z sub T of the capping temperature inversion and on the values of the Fourier coefficients of the background height-dependent vector wind profile. As input to our nonlinear spectral model, descent and ascent runs by the Electra provide for initial estimates of the inversion height and the wind profiles. In the first phase of the investigation presented here, a mechanism is said to be a probable contributor to the development of roll circulations within the stratocumulus-topped boundary layer if the modeled roll orientation and wavelengths agree with their observed values. Preliminary results using the 14-coefficient model of Haack-Hirschberg (1988) are discussed for the 7 July 1987 Electra Mission 188-A (Flight 5). This mission was flown across a sharp cloud boundary that was within a LANDSAT/SPOT scene. The stratocumulus deck was relatively solid in the eastern part of the scene, while there was a rapid decrease in cloud cover to scattered cumulus clouds aligned in streets to the west. These cloud streets were oriented nearly parallel to the mean wind direction in the layer, which was approximately 340 degrees. The hypothesis that roll circulations occurred in both the relatively clear and the cloudy regions is investigated using as model input a descent profile obtained in the relatively clear air and an ascent profile obtained in the cloudy air. Initial results for the clear air case are that the pure inflection point mode is not possible and the pure thermal mode was oriented 35 degrees to the right of the mean wind direction. The origin of this unacceptably large discrepancy between the observed and modeled results will be investigated further and the conclusions reported at the next FIRE workshop.

Shirer, Hampton N.↗
NASA NTRS1993

Plume structures in the hard-turbulent regime of three-dimensional infinite Prandtl number convection

Numerical simulations of three-dimensional infinite Prandtl number thermal convection with Rayleigh number (Ra) up to 10 exp 8 are reported. Convection with Ra higher than 10 exp 7 is characterized by the appearance of disconnected thermal plumes. The smaller plumes are detached by the currents produced by the larger plumes. The low wavenumber portion of a thermal power spectrum near the boundary layer becomes flat at high Ra, while the spectrum measured in the interior shows a positive slope for low wavenumbers. Differences are found in the thermal spectra between 2D and 3D models.

Malevsky, A. V.↗
NASA NTRS1993

Steady, three-dimensional, internally heated convection

Numerical calculations have been carried out of steady, symmetric, 3D modes of convection in internally heated, infinite Prandtl number, Boussinesq fluids at a Rayleigh number of 1.4 x 10 exp 4 in a spherical shell with inner/outer radius of 0.55 and in a 3 x 3 x 1 rectangular box. Multiple patterns of convection occur in both geometries. In the Cartesian geometry, the patterns are dominated by cylindrical cold downflows and a broad hot upwelling. In the spherical geometry the patterns consist of cylindrical cold downwellings centered either at the vertices of a tetrahedron or the centers of the faces of a cube. The cold downflow cylinders are immersed in a background of upwelling within which there are cylindrical hot concentrations (plumes) and hot halos around the downflows. The forced hot upflow return plumes of internally heated spherical convection are fundamentally different from the buoyancy-driven plumes of heated from below convection.

Schubert, G.↗
DOE OSTI2026

Flow dynamics and heat transfer in simplified battery energy storage systems with heated battery modules

Large-scale energy storage systems (ESSs) composed of batteries show promise in addressing current energy challenges, but dissipation of generated heat is important. Here, this paper focuses on buoyant convective flows in simplified ESS battery racks. Natural convection is not generally the primary cooling strategy but can be important in abnormal scenarios where there is module overheat or potentially thermal runaway. We use computational fluid dynamics to investigate the flow dynamics and heat transfer mechanisms in a simplified parameterized rack design. Despite its simplicity, this configuration produces many of the relevant features expected in real ESSs without details of module geometry or hardware, allowing broad conclusions independent of manufacture-specific designs. We start by providing visualizations of the flowfield and measurements of entrainment, heat flux, and pressure. To characterize the dependence on the system parameters, we develop an integral-scale analysis of the average temperature equation to highlight the dominant source terms. We use results from this analysis to derive a steady network model composed of simple algebraic expressions to provide first-order predictions of entrainment through the rack. The network model leads to a linear scaling of the Reynolds number based on convective mass flux with respect to the Grashof number based on the heat source. We deduce empirical relationships that relate the heat exchanged between modules using a surface-averaged Nusselt number as a function of the local Reynolds and Rayleigh numbers. Lastly, we investigate how space between the modules and rack in the spanwise direction creates flow bypass, resulting in different flow pathways.

Battery thermal management↗
NASA NTRS1990

Three-dimensional compressible convection at low Prandtl numbers

Numerical simulations are used to study fully compressible thermal convection at large Rayleigh numbers. Results are presented from a sequence of three-dimensional simulations that reveal a transition from gradually-evolving laminar convection to nearly turbulent convection as the Prandtl number is reduced from a value of unity to one-tenth. The convective flows form irregular cellular patterns near the upper surface, possesing a network of fast downflow at cell peripheries and gentler upflow at cell centers. At greater depths the curving sheets of downflow collapse into plumes which may twist and possess substantial vertical vorticity. For the lowest Prandtl number, the convection near the bottom of the layer appears to be turbulent, yet the rapidly varying small-scale flow structure there is accompanied by more ordered sites of wavering upflow, with the latter able to penetrate all the way to the upper boundary. Thus a significant component of the flow is able to extend over multiple density scale heights, in contrast to what is argued in formulating mixing-length models for stellar convection. Results are also shown from two-dimensional simulations carried out with very high spatial resolution, which reveal that supersonic convection with fluttering shock systems can be realized.

Toomre, Juri↗
NASA NTRS1989

Numerical modeling of enclosure convection

A numerical study on the steady and unsteady natural convection in two-dimensional rectangular enclosures has been performed by a time-accurate ADI finite difference scheme. The study covered a range of Rayleigh numbers between 1000 and 10 to the 7th, aspect ratios between 0.2 and 10.0, and tilt angles between -90 (heating from bottom) and +90 deg (heating from top). Various Prandtl numbers have been studied, but only the results of water (Pr = 7.0) are reported here due to space limitations. The physics revealed, however, includes the convection phenomena and the Rayleigh-Benard stability, as well as the combined mechanism of these two. The onset of secondary cells is determined by using a velocity map, which is simpler and cleaner, instead of a streamline plot. The critical Ra number for the occurrence of these secondary cells is shown to be lower than can be detected by experimental studies. On the Rayleigh-Benard stability part, a second transition from stable single-cell convection to periodic multicellular convection is disclosed.

Duh, J. C.↗
NASA NTRS1991

Characteristics-based methods applied to infinite Prandtl number thermal convection in the hard turbulent regime

Characteristics-based methods for the advection-diffusion equation are presented and directly applied to study thermal convection with extremely large Rayleigh number (Ra). It is shown that the operator-splitting method for advection-diffusion problems is very accurate for determining the advected field at extremely high Peclet number (Pe). The technique presented is considered to have great potential for solving advection-dominated problems, while the Langrangian method is more accurate for lower Pe. It is noted that the accuracy of these characteristics-based methods strongly depends on the quality of interpolation. The computational time for the operator-splitting method grows with the number of time steps employed. The Langrangian method was used for simulations of convection at very high Ra, up to 3 x 10 to the 9th, and time-dependent, thermal convection solutions were obtained for infinite Prandtl number.

Malevsky, A. V.↗
NASA NTRS1999

Thermocapillary Flows with Low Frequency g-Jitter

A thermocapillary parallel flow is established in a fluid filled slot with an applied temperature gradient. Low frequency jitter is imposed in arbitrary directions. Vertical jitter proves to be relatively uninteresting, merely augmenting or opposing the basic thermocapillary flow. Streamwise jitter still produces parallel flows, but these now exhibit boundary layers or layered cellular structures for large Rayleigh number as the applied stratification alternates between stable and unstable. Runaways are possible for unstable stratification and these correspond to resonant excitation of stationary long wave Rayleigh-Benard modes. Spanwise jitter produces fully three dimensional motion. A spanwise-streamwise circulation results for weak spanwise jitter, which advects the interfacial temperature establishing a subsidiary spanwise thermocapillary flow. This flow is strong at small Biot number when advected temperature is trapped in the slot, and has a counter-intuitive dependence on the spanwise-streamwise aspect ratio.

Grassia, P.↗
NASA NTRS1981

Effects of truncation in modal representations of thermal convection

The Galerkin equations for convection in a sphere are examined to determine which physical processes are neglected by the severe truncation of the equations of motion. It is demonstrated that the gross features of the flow are affected by truncation in the horizontal direction, with all of the models considered being well resolved in the vertical direction. One of the effects of truncation is to enhance the high-wave number end of the kinetic energy and thermal variance spectra. The examples cited indicate that as long as the kinetic energy spectrum decreases with wave number, a truncation gives a qualitatively correct solution. Conclusions are tested by calculating solutions to the equations of motion for several values of the Rayleigh number and the limit of horizontal spatial resolution.

Marcus, P. S.↗