A review of: “Problems of Solar and Stellar Oscillations. Iau Colloquium No. 66. Edited by D. O. Gough. D. Reidel Publishing Company, 493 pp. Df1195.00/U.S. $84.50 (ISBN 90-277-1554-8). 1983.”

Abstract

The model equations describing two-dimensional thermohaline convection of a Boussinesq fluid in a rotating horizontal layer are known to support multiple instabilities, depending on the values of certain control parameters (Arneodo et al., 1985). Most of these multiple instabilities have already been studied for double or triple diffusive convection, where behaviours ranging from simple steady to irregular motions have been found. Here we consider the one remaining bifurcation mentioned by Arneodo et al. (1985): the interaction between a steady and an oscillatory convection roll when the linear spectrum for a single wavenumber comprises one zero and one pair of purely imaginary eigenvalues. The method of centre manifolds and normal forms is used to derive evolution equations for the amplitudes of the convection rolls close to bifurcation and the behaviours associated with the equations is discussed.     

Nonlinear convection in an imposed horizontal magnetic field

Abstract

Nonlinear two-dimensional magnetoconvection, with a Boussinesq fluid driven across the field-lines, is taken as a model for giant-cell convection in the sun and late-type stars. A series of numerical experiments shows the sensitivity of the horizontal scale of convection to the applied field and to the Rayleigh number R. Overstable oscillations occur in cells as broad as they are deep, but increasing R leads to steady motions of much greater wavelength. Purely geometrical effects can cause oscillation: this work implies that strong horizontal field will in general lead to time-dependent convection.     

Magnetoconvection

Abstract

Cowling investigated the effect of an imposed magnetic field on convection in order to explain the origin of sunspots. After summarizing the classical linear theory of Boussinesq magnetoconvection, this review proceeds to more recent nonlinear results. Weakly nonlinear theory is used to establish the relevant bifurcation structure, which involves steady, oscillatory and chaotic solutions. Behaviour found in numerical experiments can then be related to these analytical results. Thereafter, attention is focused on the astrophysically relevant problem of fully compressible magnetoconvection. Steady two-dimensional nonlinear solutions show two important effects: stratification introduces an asymmetry between rising and falling fluid, while compressibility leads to evacuated magnetic flux sheets. Time-dependent behaviour includes transitions between standing waves and travelling waves, as well as changes in horizontal scale, leading to the development of more complicated spatial structures. Work on three-dimensional models, which is now in progress, will lead to a better understanding of the structure of a sunspot.     

The kinematics of hexagonal magnetoconvection

Abstract

Calculations are presented for the evolution of a magnetic field which is subject to the effect of three-dimensional motions in a convecting layer of highly conducting fluid with hexagonal symmetry. The back reaction of the field on the motions via the Lorentz force is neglected. We consider cases where the imposed field is either vertical or horizontal. In the former case, flux accumulates at cell centres, with subsidiary concentrations at the vertices of the pattern. In the latter, topological asymmetries between up- and down-moving fluid regions generate positive flux at the base of the layer and negative flux at the top, though the system is actually an amplifier rather than a self-excited dynamo. Spiral field lines form in the interiors of the cells, and the phenomenon of “flux expulsion” found in two-dimensional solutions is somewhat altered when the imposed field is horizontal. Applications for stellar magnetic fields include a possible mechanism for burying flux at the base of a convection zone.     

Long wavelength thermal convection between non-conducting boundaries

Non-linear Rayleigh-Bénard convection in a fluid layer is considered as a model of convection in the Earth's upper mantle. Previous studies have shown that when the temperature is held fixed at one of the boundaries of the layer, convection takes place in cells of width of the order of the layer depth or less. We investigate the effects of a different thermal boundary condition, in which the flux of heat is held fixed on both layer boundaries; then if this flux is just greater than that required for the onset of convection, motion takes place on horizontal scales much greater than the layer depth. An analytical treatment of the equations, based on an expansion in the depth-to-width ratio of the cells, shows that cells of a definite horizontal scale are the fastest growing according to linearised theory, but that these cells are unstable to ones of larger wavelength than themselves. Thus the dominant wavelength lengthens with time. The results hold whether the heat flux is generated internally of comes from beneath the layer. These results produce flow patterns similar to those found when the heat flux is much greater than the critical value. The results have important consequences for the understanding of mantle convection.     

Linear Stability of Thermal Convection in Rotating Systems with Fixed Heat Flux Boundaries

Linear stability of rotating thermal convection in a horizontal layer of Boussinesq fluid under the fixed heat flux boundary condition is examined by the use of a vertically truncated system up to wavenumber one. When the rotation axis is in the vertical direction, the asymptotic behavior of the critical convection for large rotation rates is almost the same as that under the fixed temperature boundary condition. However, when the rotation axis is horizontal and the lateral boundaries are inclined, the mode with zero horizontal wavenumber remains as the critical mode regardless of the rotation rate. The neutral curve has another local minimum at a nonzero horizontal wavenumber, whose asymptotic behavior coincides with the critical mode under the fixed temperature condition. The difference of the critical horizontal wavenumber between those two geometries is qualitatively understood by the difference of wave characteristics; inertial waves and Rossby waves, respectively.     

Large aspect ratio cells in two-dimensional thermal convection

Numerical experiments have been carried out on two-dimensional thermal convection, in a Boussinesq fluid with infinite Prandtl number, at high Rayleigh numbers. With stress free boundary conditions and fixed heat flux on upper and lower boundaries, convection cells develop with aspect ratios (width/depth) λ? 5, if heat is supplied either entirely from within or entirely from below the fluid layer. The preferred aspect ratio is affected by the lateral boundary conditions. If the temperature, rather than the heat flux, is fixed on the upper boundary the cells haveλ ≈ 1. At Rayleigh numbers of 2.4 × 10⁵ and greater, small sinking sheets are superimposed on the large aspect ratio cells, though they do not disrupt the circulation. Similar two-scale flows have been proposed for convection in the earth's mantle. The existence of two scales of flow in two-dimensional numerical experiments when the viscosity is constant will allow a variety of geophysically important effects to be investigated.     

Intermittent convection

Abstract

A theoretical analysis of pseudo two-dimensional, finite-amplitude, thermal convection is made for an infinite Prandtl number fluid which is subjected to a constant heat flux out of the top boundary and insulated at the bottom. For large Rayleigh numbers the convective flow becomes intermittent and the system is characterized by the following cyclic process: the formation of a thermal boundary layer by diffusion, the instability of this layer when it becomes sufficiently thick, the destruction of the layer by the convective flow, the dying down of the convection, and the reforming of the thermal boundary layer by diffusion. The periodicity and the horizontal wave number of the intermittent convective flow are found to be independent of the depth of the fluid layer but depend on the rate of cooling and the properties of the fluid.     

The effect of radiative heat loss on steady cellular convection

Summary In the atmosphere there may be layers undergoing cellular convection with a much larger heat flux through the base of the layer than through the top. This may be either because there is a steady loss of heat by radiation from the body of the fluid or because the temperature is everywhere rising. In this latter case the temperature gradients could remain constant so that the mechanics would be the same as if the heat were being lost and the temperature kept steady. The fluid is considered incompressible as in the classical theory of cellular convection, and we determine the critical Rayleigh number for the onset of convection and the width to height ratio of the cells as functions of the heat loss. The problem, is in some respects analogous to that of the motion of a viscous fluid between rotating cylinders but in this case there are two non-dimensional-numbers-the Rayleigh number (g h ⁴/K v) and a number representing the ratio of the heat loss by radiation to the heat flux. It is found that the critical Rayleigh number is decreased and the cells widened as had already been found for the case of a fluid with transfer coefficients having a spatial variation, with free boundaries, but the cells are made more narrow if the boundaries are rigid.     

Large wavenumber convection in the rotating annulus

Abstract

The annulus model considers convection between concentric cylinders with sloping endwalls. It is used as a simplified model of convection in a rapidly rotating sphere. Large azimuthal wavenumbers are preferred in this problem, and this has been exploited to develop an asymptotic approach to nonlinear convection in the annulus. The problem is further reduced because the Taylor-Proudman constraint simplifies the dependence in the direction of the rotation vector, so that a nonlinear system dependent only on the radial variable and time results. As Rayleigh number is increased a sequence of bifurcations is found, from steady solutions to periodic solutions and 2-tori, typically ending in chaotic behaviour. Both the magnetic (MHD convection) and non-magnetic problem has been considered, and in the non-magnetic case our bifurcation sequence can be compared with those found by previous two-dimensional numerical simulations.     

Transitions to broad cells in a nonlinear thermal convection system

Abstract

This paper explores the properties of a two-dimensional, Boussinesq convection model with an ad hoc term in the buoyancy tendency equation that represents a positive external feedback process acting on the buoyancy fluctuations. Linear stability analyses and nonlinear integrations are presented for the case of constant heat flux boundary conditions. Although the large wavenumber modes grow the fastest from a state of rest, the nonlinear solutions progressively evolve to cells of small wavenumber. Applications to mesoscale cellular convection in the atmosphere are discussed.     

Distribution of magnetic energy in αΩ-dynamos,I: The method

Abstract

In this paper a method for solving the equation for the mean magnetic energy <BB> of a solar type dynamo with an axisymmetric convection zone geometry is developed and the main features of the method are described. This method is referred to as the finite magnetic energy method since it is based on the idea that the real magnetic field B of the dynamo remains finite only if <BB> remains finite. Ensemble averaging is used, which implies that fields of all spatial scales are included, small-scale as well as large-scale fields. The method yields an energy balance for the mean energy density ε ≡ B ²/8π of the dynamo, from which the relative energy production rates by the different dynamo processes can be inferred. An estimate for the r.m.s. field strength at the surface and at the base of the convection zone can be found by comparing the magnetic energy density and the outgoing flux at the surface with the observed values. We neglect resistive effects and present arguments indicating that this is a fair assumption for the solar convection zone. The model considerations and examples presented indicate that (1) the energy loss at the solar surface is almost instantaneous; (2) the convection in the convection zone takes place in the form of giant cells; (3) the r.m.s. field strength at the base of the solar convection zone is no more than a few hundred gauss; (4) the turbulent diffusion coefficient within the bulk of the convection zone is about 10¹⁴cm²s^?1, which is an order of magnitude larger than usually adopted in solar mean field models.     

Convection in a rotating cylinder with its sidewall having a finite thermal conductance

Abstract

An investigation is made of steady thermal convection of a Boussinesq fluid confined in a vertically-mounted rotating cylinder. The top and bottom endwall disks are thermal conductors at temperatures T_t and T_b with δT = T_t ? T_b >0. The vertical sidewall has a finite thermal conductance. A Newtonian heat flux condition is adopted at the sidewall. The Rayleigh number of the fluid system is large to render a boundary layer-type flow. Finite-difference numerical solutions to the full Navier-Stokes equations are obtained. The vertical motions within the buoyancy layer along the sidewall induce weak meridional flows in the interior. Because of the Coriolis acceleration, the meridional flows give rise to azimuthal flows relative to the rotating container. Strong vertical gradients of azimuthal flows exist in the regions near the endwalls. As the stratification effect increases, concentration of flow gradients in thin endwall boundary layers becomes more pronounced. The azimuthal flow field exhibits considerable horizontal gradients. The temperature field develops horizontal variations superposed on the dominant vertical distribution. As either the sidewall thermal conductance or the stratification effect decreases, the temperature distribution tends to the profile varying linearly with height. Comparisons of the sizes of the dynamic effects demonstrate that, in the bulk of flow field, the vertical shear of azimuthal velocity is supported by the horizontal temperature gradient, resulting in a thermal-wind relation.     

Convection at the melting point: A thermal history of the earth's core

Abstract

Higgins and Kennedy (1971) concluded that the Earth's fluid core has a stable stratification if it is at its melting point. Busse (1972) and Elsasser suggested as an alternative that a hydrostatic-isentropic distribution of particulate solid can produce neutral stability in a partially molten core. Here this suggestion is quantified and a determination is made of the efficiency of the production of fluid motion from the heat flux. This is used to establish that macroscopic convection can exist only if the particulate solid is of sufficiently small size. A thermal history of the core compatible with upper mantle heat flux is advanced in which it is suggested that the inner core is a fairly recent feature. The implication of these results for convection-driven and precession-driven dynamos is that both can function for a small enough suspended particulate, that the convection-dynamo will fail for particles greater than one micron in diameter, and that the precession-driven dynamo probably cannot survive particles greater than ten microns in diameter.     

Testing density-dependent groundwater models: two-dimensional steady state unstable convection in infinite,finite and inclined porous layers

This study proposes the use of several problems of unstable steady state convection with variable fluid density in a porous layer of infinite horizontal extent as two-dimensional (2-D) test cases for density-dependent groundwater flow and solute transport simulators. Unlike existing density-dependent model benchmarks, these problems have well-defined stability criteria that are determined analytically. These analytical stability indicators can be compared with numerical model results to test the ability of a code to accurately simulate buoyancy driven flow and diffusion. The basic analytical solution is for a horizontally infinite fluid-filled porous layer in which fluid density decreases with depth. The proposed test problems include unstable convection in an infinite horizontal box, in a finite horizontal box, and in an infinite inclined box. A dimensionless Rayleigh number incorporating properties of the fluid and the porous media determines the stability of the layer in each case. Testing the ability of numerical codes to match both the critical Rayleigh number at which convection occurs and the wavelength of convection cells is an addition to the benchmark problems currently in use. The proposed test problems are modelled in 2-D using the SUTRA [SUTRA––A model for saturated–unsaturated variable-density ground-water flow with solute or energy transport. US Geological Survey Water-Resources Investigations Report, 02-4231, 2002. 250 p] density-dependent groundwater flow and solute transport code. For the case of an infinite horizontal box, SUTRA results show a distinct change from stable to unstable behaviour around the theoretical critical Rayleigh number of 4π² and the simulated wavelength of unstable convection agrees with that predicted by the analytical solution. The effects of finite layer aspect ratio and inclination on stability indicators are also tested and numerical results are in excellent agreement with theoretical stability criteria and with numerical results previously reported in traditional fluid mechanics literature.     

Nonlinear internal waves in the upper atmosphere

To investigate the mechanism of mixing in oscillatory doubly diffusive (ODD) convection, we truncate the horizontal modal expansion of the Boussinesq equations to obtain a simplified model of the process. In the astrophysically interesting case with low Prandtl number (traditionally called semiconvection), large-scale shears are generated as in ordinary thermal convection. The interplay between the shear and the oscillatory convection produces intermittent overturning of the fluid with significant mixing. By contrast, in the parameter regime appropriate to sea water, large-scale flows are not generated by the convection. However, if such flows are imposed externally, intermittent overturning with enhanced mixing is observed.     

The onset of magnetoconvection at large Prandtl number in a rotating layer I. Finite magnetic diffusion

Abstract

This paper develops further a convection model that has been studied several times previously as a very crude idealization of planetary core dynamics. A plane layer of electrically-conducting fluid rotates about the vertical in the presence of a magnetic field. Such a field can be created spontaneously, as in the Childress—Soward dynamo, but here it is uniform, horizontal and externally-applied. The Prandtl number of the fluid is large, but the Ekman, Elsasser and Rayleigh numbers are of order unity, as is the ratio of thermal to magnetic diffusivity. Attention is focused on the onset of convection as the temperature difference applied across the layer is increased, and on the preferred mode, i.e., the planform and time-dependence of small amplitude convection. The case of main interest is the layer confined between electrically-insulating no-slip walls, but the analysis is guided by a parallel study based on illustrative boundary conditions that are mathematically simpler.