A note on the heat transfer by the axisymmetric thermal convection in a rotating fluid annulus

Abstract

Some new measurements are presented of the axisymmetric heat transport in a differentially heated rotating fluid annulus. Both rigid and free upper surface cases are studied, for Prandtl numbers of 7 and 45, from low to high rotation rates. The rigid lid case is extended to high rotation rates by suppressing the baroclinic waves, that would normally develop at some intermediate rotation rate, with the use of sloping endwalls.

A parameter P is defined as the square of the ratio of the (non-rotating) thermal sidewall layer thickness to the Ekman layer thickness. For small P the heat transport remains unaffected by the rotation, but as P increases to order unity the Ekman layer becomes thin enough to inhibit the radial mass transport, and hence the heat flux. No explicit Prandtl number dependence is observed. Also this scaling allows the identification of the region in which the azimuthal velocity reaches its maximum. Direct comparisons are drawn with previous experimental and numerical results, which show what can be interpreted as an inhibiting effect of increasing curvature on the heat transport.     

Numerical simulations of thermal convection in a rapidly rotating spherical shell cooled inhomogeneously from above

Abstract

Numerical simulations of thermal convection in a rapidly rotating spherical fluid shell with and without inhomogeneous temperature anomalies on the top boundary have been carried out using a three-dimensional, time-dependent, spectral-transform code. The spherical shell of Boussinesq fluid has inner and outer radii the same as those of the Earth's liquid outer core. The Taylor number is 10⁷, the Prandtl number is 1, and the Rayleigh number R is 5R_c (R_c is the critical value of R for the onset of convection when the top boundary is isothermal and R is based on the spherically averaged temperature difference across the shell). The shell is heated from below and cooled from above; there is no internal heating. The lower boundary of the shell is isothermal and both boundaries are rigid and impermeable. Three cases are considered. In one, the upper boundary is isothermal while in the others, temperature anomalies with (l,m) = (3,2) and (6,4) are imposed on the top boundary. The spherically averaged temperature difference across the shell is the same in all three cases. The amplitudes of the imposed temperature anomalies are equal to one-half of the spherically averaged temperature difference across the shell. Convective structures are strongly controlled by both rotation and the imposed temperature anomalies suggesting that thermal inhomogeneities imposed by the mantle on the core have a significant influence on the motions inside the core. The imposed temperature anomaly locks the thermal perturbation structure in the outer part of the spherical shell onto the upper boundary and significantly modifies the velocity structure in the same region. However, the radial velocity structure in the outer part of the shell is different from the temperature perturbation structure. The influence of the imposed temperature anomaly decreases with depth in the shell. Thermal structure and velocity structure are similar and convective rolls are more columnar in the inner part of the shell where the effects of rotation are most dominant.     

Magnetic instabilities in rapidly rotating spherical geometries I. from cylinders to spheres

Abstract

The solution of the full nonlinear hydromagnetic dynamo problem is a major numerical undertaking. While efforts continue, supplementary studies into various aspects of the dynamo process can greatly improve our understanding of the mechanisms involved. In the present study, the linear stability of an electrically conducting fluid in a rigid, electrically insulating spherical container in the presence of a toroidal magnetic field B_o(r,θ)lø and toroidal velocity field U_o(r,θ)lø, [where (r,θ,ø) are polar coordinates] is investigated. The system, a model for the Earth's fluid core, is rapidly rotating, the magnetostrophic approximation is used and thermal effects are excluded. Earlier studies have adopted a cylindrical geometry in order to simplify the numerical analysis. Although the cylindrical geometry retains the fundamental physics, a spherical geometry is a more appropriate model for the Earth. Here, we use the results which have been found for cylindrical systems as guidelines for the more realistic spherical case. This is achieved by restricting attention to basic states depending only on the distance from the rotation axis and by concentrating on the field gradient instability. We then find that our calculations for the sphere are in very good qualitative agreement both with a local analysis and with the predictions from the results of the cylindrical geometry. We have thus established the existence of field gradient modes in a realistic (spherical) model and found a sound basis for the study of various other, more complicated, classes of magnetically driven instabilities which will be comprehensively investigated in future work.     

The deepening of the wind-Mixed layer

Abstract

A simple model is given that describes the response of the upper ocean to an imposed wind stress. The stress drives both mean and turbulent flow near the surface, which is taken to mix thoroughly a layer of depth h, and to erode the stably stratified fluid below. A marginal stability criterion based on a Froude number is used to close the problem, and it is suggested that the mean momentum has a strong role in the mixing process. The initial deepening is predicted to obey

where u. is the friction velocity of the imposed stress, N the ambient buoyancy frequency, and t the time.

After one-half inertial period the deepening is arrested by rotadeon at a depth h = 2^2/4 u.{(Nf)⁺

where f is the Coriolis frequency. The flow is then a “mixed Ekman” layer, with strong inertial oscillations superimposed on it. Three quarters of the mean energy of the deepening layer is found to be kinetic, and only one-quarter potential.

Heating and cooling are included in the model, but stress dominates for time-scales of a day or less. Non-uniform stratification and currents existing prior to the onset of the wind are easily included.

Agreement between the first formula above and laboratory experiments of Kato and Phillips is very satisfactory; the second formula is consistent with observations of Francis and Stommel, though a more thorough test is needed. Oceanic observations in general support the assumption of slab-like mean profiles and direct response of the fluid to local winds.     

The structure of separated flow past a circular cylinder in a rotating frame

Abstract

This paper examines the detailed E ^1/4-layer structure of separated flow past a circular cylinder in a low-Rossby-number rotating fluid as the Ekman number E tends to zero. This structure is based on an initial proposal by Page (1987) but with some modifications in response to further evidence, outlined both in this paper and elsewhere, on the behaviour of E ^1/4-layer flows in this context. Numerical calculations for flow in an E ^1/4 shear layer along the separated free streamline are described and the mass flux from this layer is then used to calculate the higher-order flow within the separation bubble. The flow structure is found to have two forms, depending on the value of the O(1) parameter λ, and these are compared with results from published “Navier-Stokes” type calculations for the flow at small but finite values of E.     

The vortex created by mass transfer between layers of a rotating fluid

Abstract

A simple way to model stratification of the ocean or atmosphere is in terms of two superposed homogeneous layers of different density. Effects of cooling of the upper layer, such as that which occurs during bottom-water formation in the ocean, can be simulated by mass transfer from the upper layer to the lower layer. A model is constructed to see What effect such a mass transfer has on the flow when the mass transfer is confined to a limited region. The main effects are (i) doming of the interface, which maintains pressure gradients in balance with the velocity field, (ii) cyclonic rotation in the upper layer due to conservation of angular momentum of particles king drawn toward the sink, yet anticyclonic vorticity for those particles outside the mass transfer region due to shrinking of vortex lines drawn up over the dome. (iii) generally anticyclonic rotation in the lower layer due to particles tending to maintain their angular momentum while being pushed outwards, but some cyclonic rotation near the centre of mass transfer, due to momentum transfer from the upper layer. Similar effects to these are seen in the Greenland Sea where bottom water formation occurs. Results of the same sort are also found in a laboratory model of the process.     

Magnetic Rossby waves in a stably stratified layer near the surface of the Earth's outer core

Abstract

The stratification profile of the Earth's magnetofluid outer core is unknown, but there have been suggestions that its upper part may be stably stratified. Braginsky (1984) suggested that the magnetic analog of Rossby (planetary) waves in this stable layer (the ‘H’ layer) may be responsible for a portion of the short-period secular variation. In this study, we adopt a thin shell model to examine the dynamics of the H layer. The stable stratification justifies the thin-layer approximations, which greatly simplify the analysis. The governing equations are then the Laplace's tidal equations modified by the Lorentz force terms, and the magnetic induction equation. We linearize the Lorentz force in the Laplace's tidal equations and the advection term in the magnetic induction equation, assuming a zeroth order dipole field as representative of the magnetic field near the insulating core-mantle boundary. An analytical β-plane solution shows that a magnetic field can release the equatorial trapping that non-magnetic Rossby waves exhibit. A numerical solution to the full spherical equations confirms that a sufficiently strong magnetic field can break the equatorial waveguide. Both solutions are highly dissipative, which is a consequence of our necessary neglect of the induction term in comparison with the advection and diffusion terms in the magnetic induction equation in the thin-layer limit. However, were one to relax the thin-layer approximations and allow a radial dependence of the solutions, one would find magnetic Rossby waves less damped (through the inclusion of the induction term). For the magnetic field strength appropriate for the H layer, the real parts of the eigenfrequencies do not change appreciably from their non-magnetic values. We estimate a phase velocity of the lowest modes that is rather rapid compared with the core fluid speed typically presumed from the secular variation.     

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.     

Penetrative convection in rotating fluids: A model for the base of stellar convection zones

Abstract

To model penetrative convection at the base of a stellar convection zone we consider two plane parallel, co-rotating Boussinesq layers coupled at their fluid interface. The system is such that the upper layer is unstable to convection while the lower is stable. Following the method of Kondo and Unno (1982, 1983) we calculate critical Rayleigh numbers R_c for a wide class of parameters. Here, R_c is typically much less than in the case of a single layer, although the scaling R_c～T^2/3 as T → ∞ still holds, where T is the usual Taylor number. With parameters relevant to the Sun the helicity profile is discontinuous at the interface, and dominated by a large peak in a thin boundary layer beneath the convecting region. In reality the distribution is continuous, but the sharp transition associated with a rapid decline in the effective viscosity in the overshoot region is approximated by a discontinuity here. This source of helicity and its relation to an alpha effect in a mean-field dynamo is especially relevant since it is a generally held view that the overshoot region is the location of magnetic field generation in the Sun.     

On the onset of thermal convection in slowly rotating fluid shells

Abstract

The problem of the removal of the degeneracy of the patterns of convective motion in a spherically symmetric fluid shell by the effects of rotation is considered. It is shown that the axisymmetric solution is preferred in sufficiently thick shells where the minimum Rayleigh number corresponds to degree l = 1 of the spherical harmonics. In all cases with l > 1 the solution described by sectional spherical harmonics Y^l _l (θ,φ) is preferred.     

Linear theory of rotating fluids using spherical harmonics part I: Steady flows

Abstract

It is shown that a systematic development of physical quantities using spherical harmonics provides analytical solutions to a whole class of linear problems of rotating fluids.

These solutions are regular throughout the whole domain of the fluid and are not much affected by the equatorial singularity of steady boundary layers in spherical geometries.

A comparison between this method and the one based on boundary layer theory is carried out in the case of the steady spin-up of a fluid inside a sphere.     

The destabilising nature of differential rotation

Abstract

In a rapidly rotating, electrically conducting fluid we investigate the thermal stability of the fluid in the presence of an imposed toroidal magnetic field and an imposed toroidal differential rotation. We choose a magnetic field profile that is stable. The familiar role of differential rotation is a stabilising one. We wish to examine the less well known destabilising effect that it can have. In a plane layer model (for which we are restricted to Roberts number q = 0) with differential rotation, U = sΩ(z)1 _?, no choice of Ω(z) led to a destabilising effect. However, in a cylindrical geometry (for which our model permits all values of q) we found that differential rotations U = sΩ(s)1 _? which include a substantial proportion of negative gradient (dΩ/ds ≤ 0) give a destabilising effect which is largest when the magnetic Reynolds number R _m = O(10); the critical Rayleigh number, Ra _c, is about 7% smaller at minimum than at R_m = 0 for q = 10⁶. We also find that as q is reduced, the destabilising effect is diminished and at q = 10^?6, which may be more appropriate to the Earth's core, the effect causes a dip in the critical Rayleigh number of only about 0.001%. This suggests that we see no dip in the plane layer results because of the q = 0 condition. In the above results, the Elsasser number A = 1 but the effect of differential rotation is also dependent on A. Earlier work has shown a smooth transition from thermal to differential rotation driven instability at high A [A = O(100)]. We find, at intermediate A [A = O(10)], a dip in the Ra_c vs. R_m curve similar to the A = 1 case. However, it has Ra_c ≤ 0 at its minimum and unlike the results for high A, larger values of R_m result in a restabilisation.     

Scattering of internal waves in a two-layer fluid flowing through a channel with small undulations

Using two-dimensional linear water wave theory, we consider the problem of normal water wave (internal wave) propagation over small undulations in a channel flow consisting of a two-layer fluid in which the upper layer is bounded by a fixed wall, an approximation to the free surface, and the lower one is bounded by a bottom surface that has small undulations. The effects of surface tension at the surface of separation is neglected. Assuming irrotational motion, a perturbation analysis is employed to calculate the first-order corrections to the velocity potentials in the two-layer fluid by using Green’s integral theorem in a suitable manner and the reflection and transmission coefficients in terms of integrals involving the shape function c(x) representing the bottom undulation. Two special forms of the shape function are considered for which explicit expressions for reflection and transmission coefficients are evaluated. For the specific case of a patch of sinusoidal ripples having the same wave number throughout, the reflection coefficient up to the first order is an oscillatory function in the quotient of twice the interface wave number and the ripple wave number. When this quotient approaches one, the theory predicts a resonant interaction between the bed and the interface, and the reflection coefficient becomes a multiple of the number of ripples. High reflection of the incident wave energy occurs if this number is large. Again, when a patch of sinusoidal ripples having two different wave numbers for two consecutive stretches is considered, the interaction between the bed and the interface near resonance attains in the neighborhood of two (singular) points along the x-axis (when the ripple wave number of the bottom undulation become approximately twice as large as the interface wave number). The theoretical observations are presented in graphical form.     

An analytical model for ice-edge upwelling

Abstract

This paper presents an analytical, two-dimensional model of the wind-induced homogeneous circulation near the edge of an ice pack floating on the ocean surface. It is shown that a vertical shear layer arises under the ice edge, by which the wind-driven geostrophic motion in the open ocean is matched to the flow region underneath the ice. As in coastal upwelling models, this shear layer consists of a thin E ^1/2-layer inside a thicker E ^1/4-layer (E being the Ekman number). Under certain conditions the shear layer produces a vertical mass flux from the bottom to the surface Ekman layer. Near the surface this upwelling flux is concentrated in the narrow E ^1/2-layer. Comparison with observations of upwelling at the edge of a polar ice pack shows good agreement.     

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.     

Spin-up of a two-layer fluid in a rotating cylinder

Abstract

we report the results of experiments on the spin-up of two layers of immiscible fluid with a free upper surface in a rotating cylinder over a wide range of internal Froude numbers. Observations of the evolution of the velocity field by particle tracking indicates that spin-up of the azimuthal velocity in the upper layer take much longer than in a homogeneous fluid. Initially, spin-up occurs at a rate comparable to that of homogeneous fluid but, at high internal Froude number, a second phase follows in which the remaining lative motion decays much more slowly. Quantitative comparison of these measurements to the theory of Pedlosky (1967) shows good agreement.

Visualization of the interface displacement during spin-up detected the presence of transient azimuthal variations in the interface elevation over a wide range of Froude (F), Ekman (E), and Rossby (ε) number. nalysis of the occurrence of the asymmetric variations using the parameter space (Q, F), where Q = E ^1/2/ε, suggested by the baroclinic instability theory and experiments of Hart (1972), showed that the flow was stable for Q > 0.06 with no discernable dependence on F. This result, together with the prediction of Pedlosky's theory that radial gradient of potential vorticity in the two layers have opposite signs, suggests at the baroclinic instability mechanism was responsible for the asymmetries. The location and timing of these instabilities may account for the discrepancies between the observations and the Pedlosky (1967) theory.     

Jet instability and the stabilizing effect of topography on jets in two-layer rotating systems

Abstract

Laboratory experiments concerning azimuthal jets in two-layer rotating systems in the absence and presence of bottom topography aligned along the jets have been conducted. The jets were forced by the selective withdrawal of fluid from the upper layer of a two-fluid system contained in a circular dishpan geometry. The principal parameters measured in the experiments were the jet Rossby number, Ro, and a stratification parameter F = r ₁/(λ₁λ₂)^1/2 where r ₁ is the radius of the circular disc used for the selective withdrawal (i.e., r ₁ is the approximate radius of curvature of the jet) and λ₁,λ₂ are the internal Rossby radii of deformation in the upper and lower fluids, respectively.

The no-topography experiments show that for a sufficiently small F, the particular value depending on Ro, the jet is stable for the duration of the experiment. For sufficiently large F, again as a function of Ro, the jet becomes unstable, exhibiting horizontal wave disturbances from modes three to seven. An Ro against F flow regime diagram is presented.

Experiments are then conducted in the presence of a bottom topography having constant cross-section and extending around a mid-radius of the dishpan. The axis of the topography is in the vicinity of the jet axis forced in the no-topography experiments and the crest of the topography is in the vicinity of the interface between the two fluids (i.e., the front associated with the jet). The experiments show that in all cases investigated the jet tends to be stabilized by the bottom topography. Experiments with the topography in place, but with the interface between the fluids being above the topography crest, are shown to be unstable but more irregular than their no-topography counterparts.

Various quantitative measurements of the jet are presented. It is shown, for example, that the jet Rossby number defined in terms of the fluid withdrawal rate from the tank. Q, can be well correlated with a dimensionless vorticity gradient, V_G , across the upper layer jet. This allows for an assessment of the stability characteristics of a jet based on a knowledge of V_G (which can be estimated given a jet profile) and F.     

Turbulent mixing in a rotating,stratified fluid

Abstract

An experimental study was carried out to investigate the effect of rotation on turbulent mixing in a stratified fluid when the turbulence in the mixed layer is generated by an oscillating grid. Two types of experiments were carried out: one of them is concerned with the deepening of the upper mixed layer in a stable, two-fluid system, and the other deals with the interaction between a stabilizing buoyancy flux and turbulence.

In the first type of experiments, it was found that rotation suppresses entrainment at larger Rossby numbers. As the Rossby number becomes smaller (Ro 0.1), the entrainment rate increases with rotation—the onset of this phenomenon, however, was found to coincide with the appearance of coherent vortices within the mixed layer. The radiation of energy from the mixed layer to the lower non-turbulent layer was found to occur and the magnitude of the energy flux was found to be increased with the rotational frequency. It was also observed that vortices are generated, rather abruptly, in the lower layer as the mixed layer deepens.

In the second set of experiments a quasi-steady mixed layer was found to develop of which the thickness varies with rotation in a fashion that is consistent with the result of the first experiment. Also the rotation was found to delay the formation of a pycnocline.     

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.     

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.