Search NASASearch

Engineering topics

Fowlis, W. W.

Publications and source records attributed to Fowlis, W. W..

At least 19 records

A laboratory model of planetary and stellar convection

Experiments on thermal convection in a rotating, differentially-heated spherical shell with a radial buoyancy force were conducted in an orbiting microgravity laboratory. A variety of convective structures, or planforms, were observed depending on the magnitude of the rotation and the nature of the imposed heating distribution. The results are in agreement with numerical simulations that can be conducted at modest parameter values, and suggest possible regimes of motion in rotating planets and stars.

Hart, J. E.

Flow regimes in a shallow rotating cylindrical annulus with temperature gradients imposed on the horizontal boundaries

Experimental flow regime diagrams are determined for a new rotating cylindrical annulus configuration which permits a measure of control over the internal vertical temperature gradient. The new annulus has radial temperature gradients imposed on plane horizontal thermally conducting endwalls (with the cylindrical sidewalls as insulators) and is considered to be more relevant to atmospheric dynamics studies than the classical cylindrical annulus. Observations have revealed that, in addition to the axisymmetric flow and nonaxisymmetric baroclinic wave flow which occur in the classical annulus, two additional nonaxisymmetric flow types occur in the new annulus: boundary-layer thermal convection and deep thermal convection. Flow regime diagrams for three different values of the imposed vertical temperature difference are presented, and explanations for the flow transitions are offered. The new annulus provides scientific backup for the proposed Atmospheric General Circulation Experiment for Spacelab. The apparatus diagram is included.

Hathaway, D. H.

Laboratory experiments on planetary and stellar convection performed on Spacelab 3

Experiments on thermal convection in a rotating, differentially heated hemispherical shell with a radial buoyancy force were conducted in an orbiting microgravity laboratory. A variety of convective structures, or planforms, were observed, depending on the magnitude of the rotation and the nature of the imposed heating distribution. The results are compared with numerical simulations that can be conducted at the more modest heating rates, and suggest possible regimes of motion in rotating planets and stars.

Hart, J. E.

Confinement of thermocapillary floating zone flow by uniform rotation

Smith and Greenspan (1984) examined theoretically the idea of using a uniform rotation of the floating zone system to confine the thermocapillary flow in crystal growth experiments to the melt sidewall, leaving the interior of the melt passive. Here, that model is extended to a full zone with a more realistic temperature distribution imposed on the sidewall, and both linear and nonlinear thermocapillary flows are theoretically studied. Linearized, analytical solutions are found using singular perturbation theory and the various sidewall boundary layers described by Greenspan (1969) for rotating fluids. The analytical and linearized numerical results are compared, and the linear and nonlinear flows are discussed. The results demonstrate that the thermocapillary flow is strong and that rotation cannot confine the flow. Temperature advection by strongly nonlinear flow is significant even for the small Prandtl number of silicon.

Fowlis, W. W.

Analytical and numerical studies of the thermocapillary flow in a uniformly floating zone

The microgravity environment of an orbiting vehicle permits crystal growth experiments in the presence of greatly reduced buoyant convection in the liquid melt. Crystals grown in ground-based laboratories do not achieve their potential properties because of dopant variations caused by flow in the melt. The floating zone crystal growing system is widely used to produce crystals of silicon and other materials. However, in this system the temperature gradient on the free sidewall surface of the melt is the source of a thermocapillary flow which does not disappear in the low-gravity environment. The idea of using a uniform rotation of the floating zone system to confine the thermocapillary flow to the melt sidewall leaving the interior of the melt passive is examined. A cylinder of fluid with an axial temperature gradient imposed on the cylindrical sidewall is considered. A half zone and the linearized, axisymmetric flow in the absence of crystal growth is examined. Rotation is found to confine the linear thermocapillary flow. A simplified model is extended to a full zone and both linear and nonlinear thermocapillary flows are studied theoretically. Analytical and numerical methods are used for the linear flows and numerical methods for the nonlinear flows. It was found that the linear flows in the full zone have more complicated and thicker boundary layer structures than in the half zone, and that these flows are also confined by the rotation. However, for the simplified model considered and for realistic values for silicon, the thermocapillary flow is not linear. The fully nonlinear flow is strong and unsteady (a weak oscillation is present) and it penetrates the interior. Some non-rotating flow results are also presented. Since silicon as a large value of thermal conductivity, one would expect the temperature fields to be determined by conduction alone. This is true for the linear and weakly nonlinear flows, but for the stronger nonlinear flow the results show that temperature advection is also important. Uniform rotation may still be a means of confining the flow and the results obtained define the procedure to be used to examine this hypothesis.

Fowlis, W. W.

Laboratory experiments in a baroclinic annulus with heating and cooling on the horizontal boundaries

Experiments have been performed in a cylindrical annulus with horizontal temperature gradients imposed upon the horizontal boundaries and in which the vertical depth was smaller than the width of the annulus. Qualitative observations were made by the use of small, suspended, reflective flakes in the liquid (water). Four basic regimes of flow were observed: (1) axisymmetric flow, (2) deep cellular convection, (3) boundary layer convective rolls, and (4) baroclinic waves. In some cases there was a mix of baroclinic and convective instabilities present. As a 'mean' interior Richardson number was decreased from a value greater than unity to one less than zero, axisymmetric baroclinic instability of the Solberg type was never observed. Rather, the transition was from non-axisymmetric baroclinic waves, to a mix of baroclinic and convective instability, to irregular cellular convection.

Miller, T. L.

Theoretical and experimental design studies for the Atmospheric General Circulation Experiment

The major criterion for the Atmospheric General Circulation Experiment (AGCE) design is that it be possible to realize strong baroclinic instability in the spherical configuration chosen. A configuration was selected in which a hemispherical shell of fluid is subjected to latitudinal temperature gradients on its spherical boundaries and the latitudinal boundaries are insulators. Work in the laboratory with a cylindrical version of this configuration revealed more instabilities than baroclinic instability. Since researchers fully expect these additional instabilities to appear in the spherical configuration also, they decided to continue the laboratory cylindrical annulus studies. Four flow regimes were identified: an axisymmetric Hadley circulation, boundary layer convection, baroclinic waves and deep thermal convection. Regime diagrams were prepared.

Fowlis, W. W.

Three-dimensional baroclinic instability of a Hadley cell for small Richardson number

A three-dimensional, linear stability analysis of a baroclinic flow for Richardson number, Ri, of order unity is presented. The model considered is a thin horizontal, rotating fluid layer which is subjected to horizontal and vertical temperature gradients. The basic state is a Hadley cell which is a solution of the complete set of governing, nonlinear equations and contains both Ekman and thermal boundary layers adjacent to the rigid boundaries; it is given in a closed form. The stability analysis is also based on the complete set of equations; and perturbation possessing zonal, meridional, and vertical structures were considered. Numerical methods were developed for the stability problem which results in a stiff, eighth-order, ordinary differential eigenvalue problem. The previous work on three-dimensional baroclinic instability for small Ri was extended to a more realistic model involving the Prandtl number, sigma, and the Ekman number, E, and to finite growth rates and a wider range of the zonal wavenumber.

Antar, B. N.

Theoretical analyses of baroclinic flows

Completed and ongoing research activities are discussed briefly, including a three-dimensional, linear stability analysis of the baroclinic Hadley cell and a numerical model of the baroclinic flow between two rotating concentric spheres. This model simulates axisymmetric flow in the Atmospheric General Circulation Experiment configuration. A computer code designed to solve the strongly nonlinear stability problem for the Eady basic state is mentioned.

Antar, B. N.

Theoretical flow regime diagrams for the AGCE

The major criterion for the design of the Atmospheric General Circulation Experiment is that it be possible to realize strong baroclinic instability in the apparatus. A spherical annulus configuration which allows only steady basic state flows was chosen for the first set of stability analyses. Baroclinic instability was found for this configuration and few results suggest a regime diagram very different from the cylindrical annulus regime diagram.

Fowlis, W. W.

Laboratory studies of baroclinic instability at small Richardson number

As part of the support program for the Atmospheric General Circulation Experiment, laboratory studies of baroclinic and other convective instabilities were performed for a thin layer of fluid between thermally conducting horizontal discs. There were three types of modes identified. The first has a spiral-arm appearance, and exists for large enough horizontal thermal forcing, weak enough static stability, and large enough rotation. The source of this wave is shown to be the Eady mode of instability. The second mode is due to convective instability in the thermal boundary layers which exist due to the thermally conducting horizontal boundaries. Finally, for strong enough negative static stability, thermal convection of the Benard type appears. The most significant result is that the symmetric (Solberg) mode was not found, even though the infinite-plane theory predicts this mode under certain experimental conditions.

Miller, T. L.

Finite-difference fluid dynamics computer mathematical models for the design and interpretation of experiments for space flight

Numerical methods are used to design a spherical baroclinic flow model experiment of the large scale atmosphere flow for Spacelab. The dielectric simulation of radial gravity is only dominant in a low gravity environment. Computer codes are developed to study the processes at work in crystal growing systems which are also candidates for space flight. Crystalline materials rarely achieve their potential properties because of imperfections and component concentration variations. Thermosolutal convection in the liquid melt can be the cause of these imperfections. Such convection is suppressed in a low gravity environment. Two and three dimensional finite difference codes are being used for this work. Nonuniform meshes and implicit iterative methods are used. The iterative method for steady solutions is based on time stepping but has the options of different time steps for velocity and temperature and of a time step varying smoothly with position according to specified powers of the mesh spacings. This allows for more rapid convergence. The code being developed for the crystal growth studies allows for growth of the crystal as the solid-liquid interface. The moving interface is followed using finite differences; shape variations are permitted. For convenience in applying finite differences in the solid and liquid, a time dependent coordinate transformation is used to make this interface a coordinate surface.

Roberts, G. O.

Three-dimensional baroclinic instability of a Hadley cell for small Richardson number

For the case of a baroclinic flow whose Richardson number, Ri, is of order unity, a three-dimensional linear stability analysis is conducted on the basis of a model for a thin, horizontal, rotating fluid layer which is subjected to horizontal and vertical temperature gradients. The Hadley cell basic state and stability analysis are both based on the Navier-Stokes and energy equations, and perturbations possessing zonal, meridional, and vertical structures are considered. An attempt is made to extend the previous theoretical work on three-dimensional baroclinic instability for small Ri to a more realistic model involving the Prandtl and Ekman numbers, as well as to finite growth rates and a wider range of the zonal wavenumber. In general, it is found that the symmetric modes of maximum growth are not purely symmetric, but have a weak zonal structure.

Antar, B. N.

Numerical solutions for spin-up from rest in a cylinder

A set of three-dimensional flow-field data for the region around a cylinder impulsively spun-up from rest was derived with a numerical model based on the Navier-Stokes equations. Laser-Doppler anemometer data in the azimuthal direction was employed to test the model predictions, and data was developed for a flowfield with Ekman numbers from 9.18/1,000,000 to 9.18/10,000. The contributions of inviscid and viscous terms were determined as functions of radius and time. It was found that immediately after start-up viscous diffusion is the dominant factor, which is replaced by nonlinear radial advection. The Coriolis force dominates in the later stages of spin-up. The inward radial flow is a maximum near the front, where the vertical velocity is small, but features strong radial gradients, as it does at the edge of the Ekman layer.

Hyun, J. M.

Symmetric baroclinic instability of a Hadley cell

A symmetric baroclinic instability is examined in terms of a Boussinesq fluid contained between two horizontal plates to determine the effects of the Ekman and thermal layers. Governing equations are written for a rotating reference frame, taking into account the Rossby, Ekman, and Prandtl numbers. Equations are defined for the perturbation functions, treated as an eigenvalue problem, and a numerical integration of the full eighth order differential system is performed by a shooting technique. An instability is found to occur in the Hadley cell containing both Ekman and thermal boundary layers when the Richardson number is close to unity. If the Prandtl number is fixed the critical Richardson number decreases with an increasing Ekman number until the Ekman number reaches a certain value, at which time the fluid is stable.

Antar, B. N.

Numerical solutions for the spin-up of a stratified fluid

The model of Warn-Varnas et al. (1978) is used to numerically examine the spin-up flow of a thermally stratified fluid in a cylinder with an insulating side wall, and comparison of the results with the laser-Doppler measurements of Lee (1975) shows excellent agreement. It is shown that flow gradients are created in the interior of the fluid during the meridional circulation spin-up phase, and that the azimuthal flow decayed faster than has been predicted by Wallin (1969). It is established that viscous diffusion in the interior, arising from the interior-flow gradients, is the cause of the discrepancy with Wallin's theory.

Hyun, J. M.

The numerical design of a spherical baroclinic experiment for Spacelab flights

The near-zero G environment of Spacelab is the basis of a true spherical experimental model of synoptic scale baroclinic atmospheric processes, using a radial dielectric body force analogous to gravity over a volume of liquid within two concentric spheres. The baroclinic motions are generated by corotating the spheres and imposing thermal boundary conditions, such that the liquid is subjected to a stable radial gradient and a latitudinal gradient. Owing to mathematical difficulties associated with the spherical geometry, quantitative design criteria can be acquired only by means of numerical models. The procedure adopted required the development of two computer codes based on the Navier-Stokes equations. The codes, of which the first calculates axisymmetric steady flow solutions and the second determines the growth or decay rates of linear wave perturbations with different wave numbers, are combined to generate marginal stability curves.

Fowlis, W. W.

Geophysical fluid dynamics

Systematic scaling or dimensional analysis reveals that certain scales of geophysical fluid flows (such as stellar, ocean, and planetary atmosphere circulations) can be accurately modeled in the laboratory using a procedure which differs from conventional engineering modeling. Rather than building a model to obtain numbers for a specific design problem, the relative effects of the significant forces are systematically varied in an attempt to deepen understanding of the effects of these forces. Topics covered include: (1) modeling a large-scale planetary atmospheric flow in a rotating cylindrical annulus; (2) achieving a radial dielectric body force; (3) spherical geophysical fluid dynamics experiments for Spacelab flights; (4) measuring flow and temperature; and (5) the possible effect of rotational or precessional disturbances on the flow in the rotating spherical containers.

Fowlis, W. W.