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.     

Heat transport by variable viscosity convection and implications for the Earth's thermal evolution

The case is presented that the efficiency of variable viscosity convection in the Earth's mantle to remove heat may depend only very weakly on the internal viscosity or temperature. An extensive numerical study of the heat transport by 2-D steady state convection with free boundaries and temperature dependent viscosity was carried out. The range of Rayleigh numbers (Ra) is 10⁴?10⁷ and the viscosity contrast goes up to 250000. Although an absolute or relative maximum of the Nusselt number (Nu) is obtained at long wavelength in a certain parameter range, at sufficiently high Rayleigh number optimal heat transport is achieved by an aspect ratio close to or below one. The results for convection in a square box are presented in several ways. With the viscosity ratio fixed and the Rayleigh number defined with the viscosity at the mean of top and bottom temperature the increase of Nu with Ra is characterized by a logarithmic gradient β = ?ln(Nu)/? ln(Ra) in the range of 0.23–0.36, similar to constant viscosity convection. More appropriate for a cooling planetary body is a parameterization where the Rayleigh number is defined with the viscosity at the actual average temperature and the surface viscosity is fixed rather than the viscosity ratio. Now the logarithmic gradient β falls below 0.10 when the viscosity ratio exceeds 250, and the velocity of the surface layer becomes almost independent of Ra. In an end-member model for the Earth's thermal evolution it is assumed that the Nusselt number becomes virtually constant at high Rayleigh number. In the context of whole mantle convection this would imply that the present thermal state is still affected by the initial temperature, that only 25–50% of the present-day heat loss is balanced by radiogenic heat production, and the plate velocities were about the same during most of the Earth's history.     

Convection with heat flux prescribed on the boundaries of the system. I. The effect of temperature dependence of material properties

Abstract

We study the bifurcation to steady two-dimensional convection with the heat flux prescribed on the fluid boundaries. The fluid is weakly non-Boussinesq on account of a slight temperature dependence of its material properties. Using expansions in the spirit of shallow water theory based on the preference for large horizontal scales in fixed flux convection, we derive an evolution equation for the horizontal structure of convective cells. In the steady state, this reduces to a simple nonlinear ordinary differential equation. When the horizontal scales of the cells exceed a certain critical size, the bifurcation to steady convection is subcritical and the degree of subcriticality increases with increasing cell size.     

The effect of stratification and compressibility on anelastic convection in a rotating plane layer

We study the effect of stratification and compressibility on the threshold of convection and the heat transfer by developed convection in the nonlinear regime in the presence of strong background rotation. We consider fluids both with constant thermal conductivity and constant thermal diffusivity. The fluid is confined between two horizontal planes with both boundaries being impermeable and stress-free. An asymptotic analysis is performed in the limits of weak compressibility of the medium and rapid rotation (τ^?1/12???|θ|???1, where τ is the Taylor number and θ is the dimensionless temperature jump across the fluid layer). We find that the properties of compressible convection differ significantly in the two cases considered. Analytically, the correction to the characteristic Rayleigh number resulting from small compressibility of the medium is positive in the case of constant thermal conductivity of the fluid and negative for constant thermal diffusivity. These results are compared with numerical solutions for arbitrary stratification. Furthermore, by generalizing the nonlinear theory of Julien and Knobloch [Fully nonlinear three-dimensional convection in a rapidly rotating layer. Phys. Fluids 1999, 11, 1469–1483] to include the effects of compressibility, we study the Nusselt number in both cases. In the weakly nonlinear regime we report an increase of efficiency of the heat transfer with the compressibility for fluids with constant thermal diffusivity, whereas if the conductivity is constant, the heat transfer by a compressible medium is more efficient than in the Boussinesq case only if the specific heat ratio γ is larger than two.     

Modeling of Rapidly Rotating Thermal Convection Using Vorticity and Vector Potential

We have used a numerical scheme based on higher-order finite differences to investigate effects of adiabatic heating and viscous dissipation on 3-D rapidly rotating thermal convection in a Cartesian box with an aspect-ratio of 221. Although we omitted coupling with the magnetic field, which can play a key role in the dynamics of the Earth's core, the understanding of non-linear rotating convection including realistic thermodynamic effects is a necessary prerequisite for understanding the full complexity of the Earth's core dynamics. The system of coupled partial differential equations has been solved in terms of the principal variables vorticity , vector potential A and temperature T. The use of the vector potential A allows the velocity field to be calculated with one spatial differentiation in contrast to the spheroidal and toroidal function approach. The temporal evolution is governed by a coupled time-dependent system consisting of and T. The equations are discretized in all directions by using an eighth-order, variable spaced scheme. Rayleigh number Ra of 10⁶, Taylor number Ta of 10⁸ and a Prandtl number Pr of 1 have been employed. The dissipation number of the outer core was taken to be 0.2. A stretched grid has been employed in the vertical direction for resolving the thin shear boundary layers at the top and bottom. This vertical resolution corresponds to around 240 regularly spaced points with an eighth-order accuracy. For the regime appropriate to the Earth's outer core, the dimensionless surface temperature T ₀ takes a large value, around 4. This large value in the adiabatic heating/cooling term is found to cause stabilization of both the temperature and velocity fields.     

Numerical simulations of mantle convection: Time-dependent,three-dimensional,compressible, spherical shell

Abstract

We describe nonlinear time-dependent numerical simulations of whole mantle convection for a Newtonian, infinite Prandtl number, anelastic fluid in a three-dimensional spherical shell for conditions that approximate the Earth's mantle. Each dependent variable is expanded in a series of 4,096 spherical harmonics to resolve its horizontal structure and in 61 Chebyshev polynomials to resolve its radial structure. A semiimplicit time-integration scheme is used with a spectral transform method. In grid space there are 61 unequally-spaced Chebyshev radial levels, 96 Legendre colatitudinal levels, and 192 Fourier longitudinal levels. For this preliminary study we consider four scenarios, all having the same radially-dependent reference state and no internal heating. They differ by their radially-dependent linear viscous and thermal diffusivities and by the specified temperatures on their isothermal, impermeable, stress-free boundaries. We have found that the structure of convection changes dramatically as the Rayleigh number increases from 10⁵ to 10⁶ to 10⁷. The differences also depend on how the Rayleigh number is increased. That is, increasing the superadiabatic temperature drop, δT, across the mantle produces a greater effect than decreasing the diffusivities. The simulation with a Rayleigh number of 10⁷ is approximately 10,000 times critical, close to estimates of that for the Earth's mantle. However, although the velocity structure for this highest Rayleigh number scenario may be adequately resolved, its thermodynamic structure requires greater horizontal resolution. The velocity and thermodynamic structures of the scenarios at Rayleigh numbers of 10⁵ and 10⁶ appear to be adequately resolved. The 10⁵ Rayleigh number solution has a small number of broad regions of warm upflow embedded in a network of narrow cold downflow regions; whereas, the higher Rayleigh number solutions (with large δT) have a large number of small hot upflow plumes embedded in a broad weak background of downflow. In addition, as would be expected, these higher Rayleigh number solutions have thinner thermal boundary layers and larger convective velocities, temperatures perturbations, and heat fluxes. These differences emphasize the importance of developing even more realistic models at realistic Rayleigh numbers if one wishes to investigate by numerical simulation the type of convection that occurs in the Earth's mantle.     

Convective Instability of a Mushy Layer - I: Uniform Permeability

Mushy layers arise and are significant in a number of geophysical contexts, including freezing of sea ice, solidification of magma chambers and inner-core solidification. A mushy layer is a region of solid and liquid in phase equilibrium which commonly forms between the liquid and solid regions of a solidifying system composed of two or more constituents. We consider the convective instability of a plane mushy layer which advances steadily upwards as heat is withdrawn at a uniform rate from the bottom of a eutectic binary alloy. The solid which forms is assumed to be composed entirely of the denser constituent, making the residual liquid within the mush compositionally buoyant and thus prone to convective motion. In this article we focus on the large-scale mush mode of instability, arguing that the 'boundary-layer' mode is not amenable to the standard stability analysis, because convective motions occur on that scale for any non-zero value of the Rayleigh number. We quantify the minimum critical Rayleigh number and determine the structure of the convective modes of motion within the mush and the associated deflections of the mush-melt and mush-solid boundaries. This study of convective perturbations differs from previous analyses in two ways; the inhibition of motion and deformation of the mush-melt interface by the stable stratification of the overlying melt is properly quantified and deformation of the mush-solid interface is permitted and quantified. We find that the mush-melt interface is almost unaffected by convection while significant deformation of the mush-solid interface occurs. We show that each of these effects causes significant (unit-order) changes in the predicted critical Rayleigh number. The marginal modes depend on three dimensionless parameters: a scaled eutectic temperature, τ _e (which characterizes the eutectic temperature relative to the depression of the liquidus), a scaled superheat, τ_∞ (which measures the amount by which the temperature of the incoming melt exceeds the liquidus temperature) and the Stefan number, S (which measures the latent heat of crystallization). To survey parameter space, we focus on seven cases, a standard case having S = τ_∞ = τ _e = 1, and six others in which one of the parameters is either large or small compared with unity: a nearly pure case (τ _e = 100; having little of the light constituent), the large superheat limit (τ_∞→ ∞), a case of large latent heat (S = 100), the near eutectic limit (τ _e → 0), a case of small superheat (τ_∞ = 0.01) and the case of zero latent heat (S = 0). The critical Rayleigh number and the associated wavelength of the convection pattern are determined in each case. The eigenvector for each case is presented in terms of the streamlines and the isolines of the perturbation temperature and solid fraction.     

Physical and numerical experiments on layered convection in a density-Stratified fluid

Abstract

The formation and growth of horizontal layered convection cells in a density stratified solution of salt water subject to an impulsively applied lateral temperature gradient is investigated with physical and numerical experiments. Results indicate that lyers are induced by two mechanisms. One is the successive formation of layers due to the presence of the top and bottom boundaries. The other is the spontaneous occurrence of layers when a suitably defined Rayleigh number exceeds a critical value. It is found that well established layers are homogeneous in temperature and salinity and are separated by sharp gradients in density. Lateral heat transfer is of a periodic nature. Numerical experiments were carried out for finite and infinite geometry cases. For the finite geometry case, convection cells are generated successively inward from the horizontal boundaries. For the infinite geometry case, periodic conditions in the vertical direction are assumed. With continuous input of small perturbations, simultaneous occurrence of the convection cells is obtained at supercritical Rayleigh numbers. Criteria for determining the onset of spontaneous cells numerically are explored.     

Dynamos with weakly convecting outer layers: implications for core-mantle boundary interaction

Convection in the Earth's core is driven much harder at the bottom than the top. This is partly because the adiabatic gradient steepens towards the top, partly because the spherical geometry means the area involved increases towards the top, and partly because compositional convection is driven by light material released at the lower boundary and remixed uniformly throughout the outer core, providing a volumetric sink of buoyancy. We have therefore investigated dynamo action of thermal convection in a Boussinesq fluid contained within a rotating spherical shell driven by a combination of bottom and internal heating or cooling. We first apply a homogeneous temperature on the outer boundary in order to explore the effects of heat sinks on dynamo action; we then impose an inhomogeneous temperature proportional to a single spherical harmonic Y ₂² in order to explore core-mantle interactions. With homogeneous boundary conditions and moderate Rayleigh numbers, a heat sink reduces the generated magnetic field appreciably; the magnetic Reynolds number remains high because the dominant toroidal component of flow is not reduced significantly. The dipolar structure of the field becomes more pronounced as found by other authors. Increasing the Rayleigh number yields a regime in which convection inside the tangent cylinder is strongly affected by the magnetic field. With inhomogeneous boundary conditions, a heat sink promotes boundary effects and locking of the magnetic field to boundary anomalies. We show that boundary locking is inhibited by advection of heat in the outer regions. With uniform heating, the boundary effects are only significant at low Rayleigh numbers, when dynamo action is only possible for artificially low magnetic diffusivity. With heat sinks, the boundary effects remain significant at higher Rayleigh numbers provided the convection remains weak or the fluid is stably stratified at the top. Dynamo action is driven by vigorous convection at depth while boundary thermal anomalies dominate in the upper regions. This is a likely regime for the Earth's core.     

Finite-Amplitude thermal convection within a self-gravitating fluid sphere

Abstract

Finite-difference calculations have been carried out to determine the structure of finite-amplitude thermal convection within a self-gravitating fluid sphere with uniform heat release. For a fixed-surface boundary condition single-cell convection breaks up into double-cell convection at a Rayleigh number of 3 × 10⁴, at a Rayleigh number of 5 × 10⁵ four-cell convection is observed. With a free-surface boundary condition only single cell convection is obtained up to a Rayleigh number of 5 × 10⁶.     

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.     

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.     

Convection driven by centrifugal bouyancy in a rotating annulus

Abstract

Drift rates and amplitudes of convection columns driven by centrifugal bouyancy in a cylindrical fluid annulus rotating about a vertical axis have been measured by thermistor probes. Conical top and bottom boundaries of the annular fluid region are responsible for the prograde Rossby wave like dynamics of the convection columns. A constant positive temperature difference between the outer and the inner cylindrical boundaries is generated by the circulation of thermostatically controled water. Mercury and water have been used as converting fluids. The measurements extend the earlier visual observations of Busse and Carrigan (1974) and provide quantitative data for an eventual comparison with nonlinear theories of thermal Rossby waves. The measured drift frequencies are in general agreement with linear theory. Of particular interest is the decline of the amplitude of convection with increasing Rayleigh number in a region beyond the onset of convection.     

Transmission of Rayleigh waves over the surface of a heterogeneous medium

Summary The frequency equation of Rayleigh waves propagating over the free surface of an isotropic, perfectly elastic, heterogeneous semi-infinite medium with material properties varying as = ₀ e ^az , = ₀ e ^az , = ₀ e ^az (a>0) has been obtained. Solution of the frequency equation in closed form is obtained in two cases (i) =0, (ii) =, and the Rayleigh wave dispersion curves for phase and group velocities drawn. In both the cases the medium yields single Rayleigh modes which cannot propagate below certain cut-off frequencies. It is found that in case (i), <c<c ₀ and 0.87500 <c _g <c ₀, and in case (ii), 1.03082 <c<c ₁ and 0.90850 <c _g <c ₁, wherec andc _g denote phase nad group velocities respectively, is the constant shear wave velocity of the mediumc ₀ andc ₁ are the corresponding Rayleigh wave velocities of the homogeneous medium of the same Poisson's ratio. The motion of the surface particles is found to be retrograde elliptical as in the homogeneous case, but the ratic of the major and minor axes now becomes frequency dependent and is plotted against frequency. In both the cases (i) and (ii), the ratio starts at a lower value at the cut-off frequency and approaches the corresponding value of the homogeneous medium at high frequencies.     

Numerical model of convection in sills

The numerical model of convection in magma sills is developed. The model is based on a full system of equations of fluid dynamics and includes heat transfer, buoyancy effects and diffusion of some minor component (marker). Solidification is treated as a phase transition. The results indicate that there are some qualitative differences between very thin sills with Rayleigh number Ra = 10⁵ and thin sills with Ra = 10⁶. For a basaltic magma the first case corresponds to the thickness of the sills of approximately 30 cm and the second case corresponds to the thickness of 60 cm. In the first case mixing is inefficient and conduction is the dominant form of heat transfer. In the second case mixing is efficient and convection is the dominant form of heat transfer. Some of the results can be scaled for the more viscous magmas in thicker sills.     

Flow of a stratified fluid bounded by walls of finite thermal conductance

Abstract

A study is made of the behavior of a thermally stratified fluid in a container when the non-horizontal boundaries have finite thermal conductance. The theory of Rahm and Walin is briefly recounted. Numerical solutions to the Navier-Stokes equations for a Boussinesq fluid in a cylinder, adopting a Newtonian heat flux condition at the vertical sidewall, are presented. Results on the details of flow and temperature fields are given over ranges of the Rayleigh number Ra, the container aspect ratio H, and the sidewall conductance S. As S increases, the isotherms in the meridional plane are horizontal at small radii but they diverge at large radii. This creates temperature nonuniformilies in the horizontal direction, and convective motions result. The salient features of the interior temperature profiles are captured by the theoretical model. The velocity field is characterized by two oppositely-directed circulations. As Ra or S varies, the qualitative circulation patterns remain substantially unchanged, but the magnitudes of the convective flows differ by large amounts. The effects of the externally-imposed parameters on the flow and temperature structures are examined.     

Thermal convection in a twisted horizontal magnetic field

The onset of convection in a layer of an electrically conducting fluid heated from below is considered in the case when the layer is permeated by a horizontal magnetic field of strength B ₀ the orientation of which varies sinusoidally with height. The critical value of the Rayleigh number for the onset of convection is derived as a function of the Chandrasekhar number Q. With increasing Q the height of the convection rolls decreases, while their horizontal wavelength slowly increases. Potential applications to the penumbral filaments of sunspots are briefly discussed.     

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.