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.     

Topology of Rayleigh–Bénard convection and magnetoconvection in plane layer

ABSTRACT

The present study aims to link the dynamics of geophysical fluid flows with their vortical structures in physical space and to study the transition of these structures due to the control parameters. The simulations are carried in a rectangular box filled with liquid gallium for three different cases, namely, Rayleigh–Bénard convection (RBC), magnetoconvection (MC) and rotating magnetoconvection (RMC). The physical setup and material properties are similar to those considered by Aurnou and Olson in their experimental work. The simulated results are validated with theoretical results of Chandrasekhar and experimental results of Aurnou and Olson. The results are also topologically verified with the help of Euler number given by Ma and Wang. For RBC, the onset is obtained at Ra greater than 1708 and at this Ra, the symmetric rolls are orientated in/along a horizontal axis. As the value of Ra increases further, the width of the horizontal rolls starts to amplify. It is observed that these two-dimensional rolls are nothing but the cross-sections of three-dimensional (3D) cylindrical rolls with wave structures. When the vertically imposed magnetic field is added to RBC, the onset of convection is delayed due to the effect of Lorentz force on the thermal buoyancy force. The presence of 3D rectangular structures is highlighted and analysed. When the magnetically influenced rectangular box rotates about vertical axis at low rotation rates in magnetoconvection model, the onset of convection gets further delayed by magnetic field, which is in general agreement with the theoretical predictions. The critical Ra increases linearly with magnetic field intensity. Coherent thermal oscillations are detected near the onset of convection, at moderate rotation rates.     

Small-scale magnetic buoyancy and magnetic pumping effects in a turbulent convection

We determine the nonlinear drift velocities of the mean magnetic field and nonlinear turbulent magnetic diffusion in a turbulent convection. We show that the nonlinear drift velocities are caused by three kinds of the inhomogeneities; i.e., inhomogeneous turbulence, the nonuniform fluid density and the nonuniform turbulent heat flux. The inhomogeneous turbulence results in the well-known turbulent diamagnetic and paramagnetic velocities. The nonlinear drift velocities of the mean magnetic field cause the small-scale magnetic buoyancy and magnetic pumping effects in the turbulent convection. These phenomena are different from the large-scale magnetic buoyancy and magnetic pumping effects which are due to the effect of the mean magnetic field on the large-scale density stratified fluid flow. The small-scale magnetic buoyancy and magnetic pumping can be stronger than these large-scale effects when the mean magnetic field is smaller than the equipartition field. We discuss the small-scale magnetic buoyancy and magnetic pumping effects in the context of the solar and stellar turbulent convection. We demonstrate also that the nonlinear turbulent magnetic diffusion in the turbulent convection is anisotropic even for a weak mean magnetic field. In particular, it is enhanced in the radial direction. The magnetic fluctuations due to the small-scale dynamo increase the turbulent magnetic diffusion of the toroidal component of the mean magnetic field, while they do not affect the turbulent magnetic diffusion of the poloidal field.     

Incorporating velocity shear into the magneto-Boussinesq approximation

Motivated by consideration of the solar tachocline, we derive, via an asymptotic procedure, a new set of equations incorporating velocity shear and magnetic buoyancy into the Boussinesq approximation. We demonstrate, by increasing the magnetic field scale height, how these equations are linked to the magneto-Boussinesq equations of Spiegel and Weiss (Magnetic buoyancy and the Boussinesq approximation. Geophys. Astrophys. Fluid Dyn. 1982, 22, 219–234).     

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.     

Generation of magnetic fields by convection in a rotating sphere,I

The magnetohydrodynamic dynamo problem is solved for an electrically conducting spherical fluid shell with spherically symmetric distributions of gravity and heat sources. The dynamics of motions generated by thermal buoyancy are dominated by the effects of rotation of the fluid shell. Dynamos are found for low and intermediate values of the Taylor number, T ? 10⁵, if the scale of the nonaxisymmetric component of the velocity field is sufficiently small. The generation of magnetic fields of quadrupolar symmetry is preferred at Rayleigh numbers close to the critical value R_c for onset of convection. As the Rayleigh number increases, the generation of dipolar magnetic fields becomes preferred.     

A simple model of convection in the terrestrial planets

Abstract

In this paper we study analytically the simplest fluid mechanical model which can mimic the convective behavior which is thought to occur in the solid mantles of the terrestrial planets. The convecting materials are polycrystalline rocks, whose creep behavior depends very strongly on temperature and probably also on pressure. As a simple model of this situation, we consider the flow of a Newtonian viscous fluid, whose viscosity depends strongly on temperature (only), and in fact has an infinite viscosity below a certain temperature, and a constant viscosity above this temperature. This model would also be directly relevant to the convection of a melt beneath its own solid phase (e.g. water below ice, though in that case there are other physical complications).

As a consequence of this assumption, there is a vigorous convection zone overlain by a stagnant lid, as also observed in analogous laboratory experiments (Nataf and Richter, 1982). The analysis is then very similar to that of Roberts (1979), but the extension to variable viscosity introduces important differences, most notably that the boundary between the lid and the convecting zone is unknown, and not horizontal. The resulting buoyancy induced stresses near this boundary are much larger than the stresses produced by buoyancy in the side-wall plumes, and mean that the dynamics of this region, and hence also the heat flux, are independent of the rest of the cell. We give a first order approximation for the Nusselt number-Rayleigh number relationship.     

Axially asymmetric velocities in the boundary layer of the nearly symmetric hydromagnetic dynamo

Abstract

An attempt has been made to include the axially asymmetric velocities into the calculation of Braginsky's Z-model of the nearly symmetric hydromagnetic dynamo. In this axisymmetric non-linear model dominated by Lorentz and Coriolis forces and maintained by a specified convection, the α-effect is prescribed. An example is shown of the axially asymmetric Archimedean buoyancy, which can imply an arbitrary alpha effect in the model with viscous core-mantle coupling. The formalisms of Tough and Roberts (1968) is also discussed and a modified α-effect in the Z-model is suggested.     

Nonlinear Convection in a Rotating Annulus with a Finite Gap

Thermal convection in a fluid-filled gap between the two corotating, concentric cylindrical sidewalls with sloping curved ends driven by radial buoyancy was first studied by Busse (Busse, F.H., "Thermal instabilities in rapidly rotating systems", J. Fluid Mech . 44 , 441-460 (1970)). The annulus model captures the key features of rotating convection in full spherical geometry and has been widely employed to study convection, magnetoconvection and dynamos in planetary systems, usually in connection with the small-gap approximation neglecting the effect of azimuthal curvature of the annulus. This article investigates nonlinear thermal convection in a rotating annulus with a finite gap through numerical simulations of the full set of nonlinear convection equations. Three representative cases are investigated in detail: a large-gap annulus with the ratio of the radii ( s i and s o ) of the sidewalls ξ = s i / o s = 0.1, a medium-gap annulus with ξ = 0.35 and a small-gap annulus with ξ = 0.8. Near the onset of convection, the effect of rapid rotation through the sloping ends forces the first (Hopf) bifurcation in the form of small-scale, steadily drifting rolls (thermal Rossby waves). At moderately large Rayleigh numbers, a variety of different convection patterns are found, including mixed-mode steadily drifting, quasi-periodic (vacillating) and temporally chaotic convection in association with various temporal and spatial symmetry-breaking bifurcations. Our extensive simulations suggest that competition between nonlinear and rotational effects with increasing Rayleigh number leads to an unusual sequence of bifurcation characterized by enlarging the spatial scale of convection.     

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.     

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.     

Convective dynamos with intermediate and strong fields

Abstract

The weak-field Benard-type dynamo treated by Soward is considered here at higher levels of the induced magnetic field. Two sources of instability are found to occur in the intermediate field regime M ～ T ^1/12, where M and T are the Hartmann and Taylor numbers. On the time scale of magnetic diffusion, solutions may blow up in finite time owing to destabilization of the convection by the magnetic field. On a faster time scale a dynamic instability related to MAC-wave instability can also occur. It is therefore concluded that the asymptotic structure of this dynamo is unstable to virtual increases in the magnetic field energy.

In an attempt to model stabilization of the dynamo in a strong-field regime we consider two approximations. In the first, a truncated expansion in three-dimensional plane waves is studied numerically. A second approach utilizes an ad hoc set of ordinary differential equations which contains many of the features of convection dynamos at all field energies. Both of these models exhibit temporal intermittency of the dynamo effect.     

Numerical Simulation of Deep-Water Convection in the Gulf of Lion

The unsteady-state process of deep-water convection in the Gulf of Lion has been observed and investigated in recent decades. However, the mechanisms of the uncertainty and irregularity of the deep-water convection in this region have not yet been fully understood. In this study, the effects of time variation of the surface buoyancy flux on the formation of the deep-water convection are examined. Numerical simulations using the Coupled Ocean/Atmosphere Mesoscale Prediction System (COAMPS^®) with the NRL Coastal Ocean Model (NCOM) as an oceanic component were conducted for the period from October 1998 to September 2000 to cover two winters, 1998–1999 and 1999–2000, over the Gulf of Lion region. The results show large differences in the deep-water convection between the two winters, even though the total surface heat fluxes during the two winter seasons are similar. The differences are related to the time variation of the surface buoyancy flux that causes large differences in the preconditioning and mixing stages of the convection.     

On the dynamics of 3-D single thermal plumes at various Prandtl numbers and Rayleigh numbers

Three-dimensional (3-D) numerical simulations of single turbulent thermal plumes in the Boussinesq approximation are used to understand more deeply the interaction of a plume with itself and its environment. In order to do so, we varied the Rayleigh and Prandtl numbers from Ra?～?10⁵ to Ra?～?10⁸ and from Pr?～?0.025 to Pr?～?70. We found that thermal dissipation takes place mostly on the border of the plume. Moreover, the rate of energy dissipation per unit mass ε _T has a critical point around Pr?～?0.7. The reason is that at Pr greater than ～0.7, buoyancy dominates inertia and thermal advection dominates wave formation whereas this trend is reversed at Pr less than ～0.7. We also found that for large enough Prandtl number (Pr?～?70), the velocity field is mostly poloidal although this result was known for Rayleigh–Bénard convection (see Schmalzl et al. [On the validity of two-dimensional numerical approaches to time-dependent thermal convection. Europhys. Lett. 2004, 67, 390--396]). On the other hand, at small Prandtl numbers, the plume has a large helicity at large scale and a non-negligible toroidal part. Finally, as observed recently in details in weakly compressible turbulent thermal plume at Pr?=?0.7 (see Plourde et al. [Direct numerical simulations of a rapidly expanding thermal plume: structure and entrainment interaction. J. Fluid Mech. 2008, 604, 99--123]), we also noticed a two-time cycle in which there is entrainment of some of the external fluid to the plume, this process being most pronounced at the base of the plume. We explain this as a consequence of calculated Richardson number being unity at Pr?=?0.7 when buoyancy balance inertia.     

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.     

Finite amplitude convection and magnetic field generation in a rotating spherical shell

Abstract

Finite amplitude solutions for convection in a rotating spherical fluid shell with a radius ratio of η=0.4 are obtained numerically by the Galerkin method. The case of the azimuthal wavenumber m=2 is emphasized, but solutions with m=4 are also considered. The pronounced distinction between different modes at low Prandtl numbers found in a preceding linear analysis (Zhang and Busse, 1987) is also found with respect to nonlinear properties. Only the positive-ω-mode exhibits subcritical finite amplitude convection. The stability of the stationary drifting solutions with respect to hydrodynamic disturbances is analyzed and regions of stability are presented. A major part of the paper is concerned with the growth of magnetic disturbances. The critical magnetic Prandtl number for the onset of dynamo action has been determined as function of the Rayleigh and Taylor numbers for the Prandtl numbers P=0.1 and P=1.0. Stationary and oscillatory dynamos with both, dipolar and quadrupolar, symmetries are close competitors in the parameter space of the problem.     

A nonlinear dynamo driven by rapidly rotating convection

In Kim et al. (Kim, E., Hughes, D.W. and Soward, A.M., “An investigation into high conductivity dynamo action driven by rotating convection”, Geophys. Astrophys. Fluid Dynam. 91, 303–332 ().) we investigated kinematic dynamo action driven by rapidly rotating convection in a cylindrical annulus. Here we extend this work to consider self-consistent nonlinear dynamo action in which the back-reaction of the Lorentz force on the flow is taken into account. In particular, we investigate, as a function of magnetic Prandtl number, the evolution of an initially weak magnetic field in two different types of convective flow – one chaotic and the other integrable. On saturation, the latter shows a systematic dependence on the magnetic Prandtl number whereas the former appears not to. In addition, we show how, in keeping with the findings of Cattaneo et al. (Cattaneo, F., Hughes, D.W. and Kim, E., “Suppression of chaos in a simplified nonlinear dynamo model”, Phys. Rev. Lett. 76, 2057–2060 ().), saturation of the growth of the magnetic field is brought about, for the originally chaotic flow, by a strong suppression of chaos.     

Local turbulence in the Earth's core

Abstract

It is shown that in the Earth's core, where the geodynamo is at work (and is supplied with energy by the prevailing unstable density stratification), a buoyancy instability of a local character exists which is highly supercritical. This instability results in fully developed turbulence dominated by small scale vortices. The influence of the Earth's rotation and of the magnetic field produced by the geodynamo makes this small scale turbulence highly anisotropic. A qualitative picture of this local anisotropic turbulence is devised and the main parameters characterizing it are estimated. Expressions for the turbulent diffusivity are developed and discussed.