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

Numerical study of baroclinic instability

The effect of a power law gravity field on baroclinic instability is examined with emphasis on the case of inverse fifth power gravity, since this is the power law produced when terrestrial gravity is simulated in spherical geometry by electrostatic means. Growth rates of unstable normal modes were obtained as a function of parameters of the problem by solving a second order differential equation numerically. Results are compared with those from an earlier study where gravity was a constant. The conclusion is that, over the range of parameter space explored here, there is no significant change in the character of theoretical regime diagrams if the vertically averaged gravity is used as a parameter

A study of the expected effects of latitude-dependent rotation rate on laboratory geophysical flow experiments

Results of a theoretical model study of some of the expected effects of spherical geometry on laboratory simulations of the type of geophysical flow that dominates the general circulation of the earth's troposphere are reported

The wave structures of the Eady model of baroclinic instability

By solving the linear quasi-geostrophic set of equations pertinent to the Eady model, the complex eigenvalues and the eigenfunctions are obtained. The propagation speed and the growth rate are computed. Quantitative information is provided about the wave structures for several unstable models, a marginally stable mode, and a stable mode. The peculiarities concerning the amplitude and the phase variations of the waves are noted as the wavenumber varies from the unstable region to the stable region. Physical interpretations of the interrelationships among the dynamical variables are given, with a view toward revealing important aspects of the energy transfer from the basic state to the growing waves

The Numerical Studies Program for the Atmospheric General Circulation Experiment (AGCE) for Spacelab Flights

The atmospheric general circulation experiment (AGCE) numerical design for Spacelab flights was studied. A spherical baroclinic flow experiment which models the large scale circulations of the Earth's atmosphere was proposed. Gravity is simulated by a radial dielectric body force. The major objective of the AGCE is to study nonlinear baroclinic wave flows in spherical geometry. Numerical models must be developed which accurately predict the basic axisymmetric states and the stability of nonlinear baroclinic wave flows. A three dimensional, fully nonlinear, numerical model and the AGCE based on the complete set of equations is required. Progress in the AGCE numerical design studies program is reported

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

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

Theoretical regime diagrams for thermally driven flows in a beta-plane channel in the presence of variable gravity

The effect of a power law gravity field on baroclinic instability is examined, with a focus on the case of inverse fifth power gravity, since this is the power law produced when terrestrial gravity is simulated in spherical geometry by a dielectric force. Growth rates are obtained of unstable normal modes as a function of parameters of the problem by solving a second order differential equation numerically. It is concluded that over the range of parameter space explored, there is no significant change in the character of theoretical regime diagrams if the vertically averaged gravity is used as parameter

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

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