The Influence of the Prandtl Number on the Style of Vigorous Thermal Convection

The effect of the Prandtl number on convection in a planar three-dimensional geometry is investigated in this study. We have employed a numerical scheme to integrate the governing equations. Differently from previous studies we have chosen stress-free boundaries. Experiments have been performed at a Rayleigh number of Ra = 10 6 for Prandtl numbers (Pr) ranging from 0.025 to 100. We have further conducted one experiment in the limiting case of infinite Prandtl number. Despite the differences in the geometry and the boundary conditions, as compared to other studies, we find a similar transition in the dynamics of the flow when the Prandtl number is increased. While the velocity and the temperature structure show diffusive character at low Pr, sharp thermal boundary layers form at high Pr. The heat transport efficiency increases with Pr until a transition value is reached, from there on Nu behaves almost asymptotically. The transition can not be caused by a change in hierarchies between velocity and thermal boundary layers, as suggested in other studies. Due to the stress-free boundaries, a velocity boundary layer does not exist. We observe that the toroidal part of the flow is strong at low Pr and looses its strength with increasing Pr, thus it is likely to be responsible for the transition. In a further chapter we demonstrate that due to the neglect of the toroidal part in two-dimensional calculations at low Pr results are obtained which are misleading, even in a qualitative sense. Infinite Pr results from 2D calculations closely resemble the dynamics of fully 3D flows.     

Plume Mode of Thermal Convection in the Earth’s Mantle

In our previous works, based on numerical models, it was shown that under certain conditions a hot material can rise in portions in the tails of thermal mantle plumes. The spectrum of these pulsations can correspond to the observed spectra of catastrophic hotspot eruptions. Since most of the existing numerical models of thermal convection for the mantle of the present Earth do not reveal these pulsations, in this work, we analyze the physical cause and initiation conditions of pulsations of thermal plumes. The results of a numerical solution of the thermal convection equations for a material with varying parameters in the extended Boussinesq approximation are presented. It is shown how the structure of the convection is transformed with the increase of convection intensity. At the Rayleigh numbers Ra > 10⁶, convection becomes unsteady, and the configuration of the ascending and descending flows changes. The new flow emerging at the mantle bottom acquires a mushroom shape with a head and a tail. After the rise of the plume’s head to the surface, the tail remains in the mantle in the form of a quasi-stationary hot steam. It turns out that at Ra ~ 5 × 10⁷, the thermal mantle plume becomes pulsating and its tail is in fact a heated channel through which the hot material rises in successive portions. At the Rayleigh numbers Ra > 5 × 10⁸, the tail of the thermal plume breaks and the plume becomes a regular conveyor of separate ascending portions of the hot material, which are referred to as thermals. Thus, thermal convection with pulsating plumes takes place at the transitional stage from the regime of quasi-stationary plumes to the regime of thermals.     

On the upper bounding approach to thermal convection at moderate rayleigh numbers

Abstract

Two upper bounding problems for thermal convection in a layer of fluid contained between perfectly conducting stress-free boundaries are treated numerically. Since the Euler equations resulting from this variational approach are simpler than the Navier-Stokes equations, they allow numerical calculations to be carried out economically to fairly large values of the Rayleigh number. The upper bounding problem formulated by Howard (1963), which yields a Nusselt number independent of Prandtl number, diverges from the correct behavior as the Rayleigh number increases. In hopes of coming closer to results of previous investigations of the Boussinesq equations of motion, a more restrictive upper bounding problem is formulated. For large Prandtl numbers the momentum equation is linearized and is used as an explicit side constraint on the variational problem, thereby forcing the solutions to more closely resemble the solutions of the Boussinesq equations. Numerical calculations at values of the Rayleigh number up to 1.5 × 10⁵ indicate that the additional constraint decreases the upper bound on the Nusselt number; it appears that this upper bound differs by only a multiplicative factor from that calculated from solutions of the full equations of motion and may be a reasonable approximation for large Rayleigh numbers.     

Two-dimensional density currents in a confined basin

We present new experimental results on the mechanisms through which steady two-dimensional density currents lead to the formation of a stratification in a closed basin. A motivation for this work is to test the underlying assumptions in a diffusive “filling box” model that describes the oceanic thermohaline circulation (Hughes, G.O. and Griffiths, R.W., A simple convective model of the global overturning circulation, including effects of entrainment into sinking regions, Ocean Modeling, 2005, submitted.). In particular, they hypothesized that a non-uniform upwelling velocity is due to weak along-slope entrainment in density currents associated with a large horizontal entrainment ratio of E _eq ?～?0.1. We experimentally measure the relationship between the along-slope entrainment ratio, E, of a density current to the horizontal entrainment ratio, E _eq, of an equivalent vertical plume. The along-slope entrainment ratios show the same quantitative decrease with slope as observed by Ellison and Turner (, 6, 423–448.), whereas the horizontal entrainment ratio E _eq appears to asymptote to a value of E eq?=?0.08 at low slopes. Using the measured values of E _eq we show that two-dimensional density currents drive circulations that are in good agreement with the two-dimensional filling box model of Baines and Turner (Baines, W.D. and Turner, J.S., Turbulent buoyant convection from a source in a confined region, J. Fluid. Mech., 1969, 37, 51–80.). We find that the vertical velocities of density fronts collapse onto their theoretical prediction that U =-2^?2/3 B ^1/3 E _eq ^2/3 (H/R) ζ, where U is the velocity, H the depth, B the buoyancy flux, R the basin width, E _eq the horizontal entrainment ratio and?ζ?= z/H the dimensionless depth. The density profiles are well fitted with?Δ?= 2^?1/3 B ^2/3 E _eq ^?2/3 H _-1 [ln(ζ )?+?τ ], where?τ?is the dimensionless time. Finally, we provide a simple example of a diffusive filling box model, where we show how the density stratification of the deep Caribbean waters (below 1850?m depth) can be described by a balance between a steady two-dimensional entraining density current and vertical diffusion in a triangular basin.     

Resistive instabilities in rapidly rotating fluids: Linear theory of the convective modes

Abstract

This paper analyzes the linear stability of a rapidly-rotating, stratified sheet pinch in a gravitational field, g, perpendicular to the sheet. The sheet pinch is a layer (O ? z ? d) of inviscid, Boussinesq fluid of electrical conductivity σ, magnetic permeability μ, and almost uniform density ρ _o; z is height. The prevailing magnetic field. B _o(z), is horizontal at each z level, but varies in direction with z. The angular velocity, Ω, is vertical and large (Ω ? V_A/d, where V_A = B₀√(μρ₀) is the Alfvén velocity). The Elsasser number, Λ = σB² ₀/2Ωρ₀, measures σ. A (modified) Rayleigh number, R = gβd²/ρ₀V² _A, measures the buoyancy force, where β is the imposed density gradient, antiparallel to g. A Prandtl number, P_K = μσK, measures the diffusivity, k, of density differences.     

Symmetric baroclinic instability for small ekman numbers

Abstract

The stability of a zonal shear flow to symmetric baroclinic perturbations is examined when the Ekman number, E, is asymptotically small. It is assumed, following Antar and Fowlis (1982), that the zonal Row is generated by imposing a constant horizontal temperature gradient γ* at the horizontal boundaries, and by maintaining a constant temperature difference δT* between them. The boundaries are at rest relative to a rotating frame.

Features of the neutral stability curve are determined for several ranges of values of δT/E ^1/3, where δT = δT*/Hγ* and H is the depth of the fluid layer, and all values of the Prandtl number, [sgrave]. In some cases it is possible to determine the whole curve analytically. The most important feature of the results is that the neutral stability curve is closed.

The results are compared to the numerical integrations of Antar and Fowlis (1982). The qualitative features of the solutions are in accord and the quantitative results are, in most cases, as good as can be expected for E only as small as ～ 10^?4. The implications of the results for experimental observations of symmetric baroclinic instability are explored.     

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.     

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.     

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.     

Magnetic field generation by convective flows in a plane layer: the dependence on the Prandtl numbers

Investigation of magnetic field generation by convective flows is carried out for three values of kinematic Prandtl number: P = 0.3, 1 and 6.8. We consider Rayleigh–Bénard convection in Boussinesq approximation assuming stress-free boundary conditions on horizontal boundaries and periodicity with the same period in the x and y directions. Convective attractors are modelled for increasing Rayleigh numbers for each value of the kinematic Prandtl number. Linear and non-linear dynamo action of these attractors is studied for magnetic Prandtl numbers P _m ≤ 100. Flows, which can act as magnetic dynamos, have been found for all the three considered values of P, if the Rayleigh number R is large enough. The minimal R, for which of magnetic field generation occurs, increases with P. The minimum (over R) of critical P_m for magnetic field generation in the kinematic regime is admitted for P = 0.3. Thus, our study indicates that smaller values of P are beneficial for magnetic field generation.     

How to anchor hotspots in a convecting mantle?

Laboratory experiments were performed to study the influence of density and viscosity layering on the formation and stability of plumes. Viscosity ratios ranged from 0.1 to 6400 for buoyancy ratios between 0.3 and 20, and Rayleigh numbers between 10⁵ and 2.10⁸. The presence of a chemically stratified boundary layer generates long-lived thermochemical plumes. These plumes first develop from the interface as classical thermal boundary layer instabilities. As they rise, they entrain by viscous coupling a thin film of the other layer and locally deform the interface into cusps. The interfacial topography and the entrainment act to further anchor the plumes, which persist until the chemical stratification disappears through entrainment, even for Rayleigh numbers around 10⁸. The pattern of thermochemical plumes remains the same during an experiment, drifting only slowly through the tank. Scaled to an Earth’s mantle without plate tectonics, our results show that: (1) thermochemical plumes are expected to exist in the mantle, (2) they could easily survive hundreds of millions of years, depending on the size and magnitude of the chemical heterogeneity on which they are anchored, and (3) their drift velocity would be at most 1-2 mm/yr. They would therefore produce long-lived and relatively fixed hotspots on the lithosphere. However, the thermochemical plumes would follow any large scale motion imposed on the chemical layer. Therefore, the chemical heterogeneity acts more as a ‘floating anchor’ than as an absolute one.     

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.     

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.     

Transitions in the style of mantle convection at high Rayleigh numbers

The pattern and style of mantle convection govern the thermal evolution, internal dynamics, and large-scale surface deformation of the terrestrial planets. In order to characterize the nature of heat transport and convective behaviour at Rayleigh numbers, Ra, appropriate for planetary mantles (between 10⁴ and 10⁸), we perform a set of laboratory experiments. Convection is driven by a temperature gradient imposed between two rigid surfaces, and there is no internal heating. As the Rayleigh number is increased, two transitions in convective behaviour occur. First we observe a change from steady to time-dependent convection at Ra≈10⁵. A second transition occurs at higher Rayleigh numbers, Ra≈5×10⁶, with large-scale time-dependent flow being replaced by isolated rising and sinking plumes. Corresponding to the latter transition, the exponent β in the power law relating the Nusselt number Nu to the Rayleigh number (NuRa^β) is reduced. Both rising and sinking plumes always consist of plume heads followed by tails. There is no characteristic frequency for the formation of plumes.     

Models of solar differential rotation

Abstract

Models of a differentially rotating compressible convection zone are calculated, considering the inertial forces in the poloidal components of the equations of motion. Two driving mechanisms have been considered: latitude dependent heat transport and anisotropic viscosity. In the former case a meridional circulation is induced initially which in turn generates differential rotation, whereas in the latter case differential rotation is directly driven by the anisotropic viscosity, and the meridional circulation is a secondary effect.

In the case of anisotropic viscosity the choice of boundary conditions has a big influence on the results: depending on whether or not the conditions of vanishing pressure perturbation are imposed at the bottom of the convection zone, one obtains differential rotation with a fast (≥ 10 ms^?1) or a slow (～ 1 ms^?1) circulation. In the latter case the rotation law is mainly a function of radius and the rotation rate increases inwards if the viscosity is larger in radial direction than in the horizontal directions.

The models with latitude dependent heat transport exhibit a strong dependence on the Prandtl number. For values of the Prandtl number less than 0.2 the pole-equator temperature difference and the surface velocity of the meridional circulation are compatible with observations. For sufficiently small values of the Prandtl number the convection zone becomes globally unstable like a layer of fluid for which the critical Rayleigh number is exceeded.     

Effect of a Viscosity Interface at 1000 km Depth on Mantle Circulation

We present the preliminary results of axisymmetric numerical simulations of thermal convection in the mantle with a phase transition boundary at 660 km depth and a viscosity interface at 1000 km depth. The results, obtained for Ra = 2 × 10 ⁶ , are compared with the case when both the phase and the viscosity boundaries are located at the same depth of 660 km.     

Enhanced scavenging ofPb andPo by processes associated with the East Pacific Rise near 8°45′N

In 1977 and 1981 a hydrothermal plume was detected at the East Pacific Rise (EPR) near 8°45′N and at MANOP Site M, 25 km east of the EPR, by anomalous²²²Rn and Mn concentrations. In 1981, samples were also taken for²¹⁰Pb,²¹⁰Po and²²⁶Ra analyses to determine if enhanced scavenging of these elements occurred in the plume. At both the ridge crest and at Site M, the ²¹⁰Pb/²²⁶Ra ratios range from 0.09 to 0.35, which are among the lowest values ever measured. It appears that removal of the²¹⁰Pb is occurring by processes operating at or near the seafloor. There is also significant ²¹⁰Po/²¹⁰Pb disequilibrium at both locations, which appears to increase away from the seafloor towards the plume (as indicated by elevated Mn concentrations). An in-situ water column scavenging process is suggested by correlation of [²¹⁰Pb-²¹⁰Po] and Mn concentrations in the plume. The residence time of the²¹⁰Po is only about 1 year, which is close to its residence time in surface waters but quite short compared to typical deep sea values. Thus rapid scavenging of the polonium in the plume and relatively rapid settling of particles from the plume is suggested.     

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.     

Numerical simulation of Martian historical dynamo: Impact of the Rayleigh number on the dynamo state

The observed Mars remnant magnetism suggests that there was an active dynamo in the Martian core. We use the MoSST core dynamics model to simulate the Martian historical dynamo, focusing on the variation of the dynamo states with the Rayleigh number Ra (a non-dimensional parameter describing the buoyancy force in the core). Our numerical results show that the mean field length scale does not vary monotonically with the Rayleigh number, and the field morphology at the core mantle boundary changes with Rayleigh number. In particular, it drifts westward with a speed decreasing with Rayleigh number. Supported by National Natural Science Foundation of China (Grant No. 40328006)     

Chemical reaction rates and entrainment within the Endeavour Ridge hydrothermal plume

The aging of the hydrothermal plume over the Endeavour segment of the Juan de Fuca Ridge was estimated by measuring the²²²Rn³He ratio in the plume as it dispersed. Despite uncertainties in the source function of hydrothermal input, it wa determined that the relative sequence of removal from the plume isH₂ > Δc >²²²Rn>CH₄ Mn, whereΔc is a measure of particle concentration and the mean life of²²²Rn is 5.5 days. H₂ is removed from the plume within hours of input while Mn is not removed within the two-week timescale of the radon-helium clock.Entrainment of bottom water within the buoyant plume may introduce additional chemical signatures into the spreading effluent layer over that which would be introduced by hydrothermal discharge alone. This is particularly significant for those chemical species which are not greatly enriched in the vent fluids relative to bottom water concentration and which display a nutrient-like profile in the deep ocean. Thus we found that significant fractions of the Si and²²⁶Ra anomalies in the plume were not of hydrothermal origin but were derived from entrained bottom water which has a higher concentration of these elements than ambient water at plume height.