Large eddy simulations of two-dimensional turbulent convection in a density-stratified fluid

Large eddy simulations (LES) of two-dimensional turbulent convection within the anelastic approximation are presented for Rayleigh number Ra?=?10⁹, Prandtl number Pr?=?1 with free-slip boundary conditions. Various subgrid-scale (SGS) models are investigated such as a similarity model, a dynamic similarity model, a dynamic eddy-viscosity model, a hyperdiffusion model and a hybrid model (dynamic similarity hyperdiffusion model). To study the effects of density stratification on the models, we have also carried out simulations for a Boussinesq flow. The SGS models are compared to direct numerical simulation (DNS) data on the basis of kinetic energy and entropy variance spectra, mean entropy profiles, r.m.s. entropy profiles and r.m.s. kinetic energy density profiles. The results show that for the Boussinesq flow, all the SGS models agree fairly well with the high resolution DNS data. However, for the strongly density-stratified flow, only the hyperdiffusion and the hybrid model show good performance.     

Transitions to chaotic thermal convection in a rapidly rotating spherical fluid shell

Abstract

Numerical simulations of thermal convection in a rapidly rotating spherical fluid shell heated from below and within have been carried out with a nonlinear, three-dimensional, time-dependent pseudospectral code. The investigated phenomena include the sequence of transitions to chaos and the differential mean zonal rotation. At the fixed Taylor number T _a =10⁶ and Prandtl number Pr=1 and with increasing Rayleigh number R, convection undergoes a series of bifurcations from onset of steadily propagating motions SP at R=R _c = 13050, to a periodic state P, and thence to a quasi-periodic state QP and a non-periodic or chaotic state NP. Examples of SP, P, QP, and NP solutions are obtained at R = 1.3R _c , R = 1.7 R _c , R = 2R _c , and R = 5 R _c , respectively. In the SP state, convection rolls propagate at a constant longitudinal phase velocity that is slower than that obtained from the linear calculation at the onset of instability. The P state, characterized by a single frequency and its harmonics, has a two-layer cellular structure in radius. Convection rolls near the upper and lower surfaces of the spherical shell both propagate in a prograde sense with respect to the rotation of the reference frame. The outer convection rolls propagate faster than those near the inner shell. The physical mechanism responsible for the time-periodic oscillations is the differential shear of the convection cells due to the mean zonal flow. Meridional transport of zonal momentum by the convection cells in turn supports the mean zonal differential rotation. In the QP state, the longitudinal wave number m of the convection pattern oscillates among m = 3,4,5, and 6; the convection pattern near the outer shell has larger m than that near the inner shell. Radial motions are very weak in the polar regions. The convection pattern also shifts in m for the NP state at R = 5R _c , whose power spectrum is characterized by broadened peaks and broadband background noise. The convection pattern near the outer shell propagates prograde, while the pattern near the inner shell propagates retrograde with respect to the basic rotation. Convection cells exist in polar regions. There is a large variation in the vigor of individual convection cells. An example of a more vigorously convecting chaotic state is obtained at R = 50R _c . At this Rayleigh number some of the convection rolls have axes perpendicular to the axis of the basic rotation, indicating a partial relaxation of the rotational constraint. There are strong convective motions in the polar regions. The longitudinally averaged mean zonal flow has an equatorial superrotation and a high latitude subrotation for all cases except R = 50R _c , at this highest Rayleigh number, the mean zonal flow pattern is completely reversed, opposite to the solar differential rotation pattern.     

Anisotropic turbulent thermal diffusion and thermal convection in a rapidly rotating fluid sphere

Estimates of the molecular values of magnetic, viscous and thermal diffusion suggest that the state of the Earth’s core is turbulent and that complete numerical simulation of the geodynamo is not realizable at present. Large eddy simulation of the geodynamo with modelling of the sub-grid scale turbulence must be used. Current geodynamo models effectively model the sub-grid scale turbulence with isotropic diffusivities larger than the molecular values appropriate for the core. In the Braginsky and Meytlis (1990) picture of core turbulence the thermal and viscous diffusivities are enhanced up to the molecular magnetic diffusivity in the directions of the rotation axis and mean magnetic field. We neglect the mean magnetic field herein to isolate the effects of anisotropic thermal diffusion, enhanced or diminished along the rotation axis, and explore the instability of a steady conductive basic state with zero mean flow in the Boussinesq approximation. This state is found to be more stable (less stable) as the thermal diffusion parallel to the rotation axis is increased (decreased), if the transverse thermal diffusion is fixed. To examine the effect of simultaneously varying the diffusion along and transverse to the rotation axis, the Frobenius norm is used to control for the total thermal diffusion. When the Frobenius norm of the thermal diffusion tensor is fixed, it is found that increasing the thermal diffusion parallel to the rotation axis is destabilising. This result suggests that, for a fixed total thermal diffusion, geodynamo codes with anisotropic thermal diffusion may operate at lower modified Rayleigh numbers.     

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.     

Patterns of thermal convection in slowly rotating,weakly magnetic fluid shells

The effects of rotation and a toroidal magnetic field on the preferred pattern of small amplitude convection in spherical fluid shells are considered. The convective motions are described in terms of associated Legendre functions P_l^|m| (cos θ). For a given pair of Prandtl number P and magnetic Prandtl number $\text{P}$_m the physically realized solution is represented either by m = 0 or |m| = l depending on the ratio of the rotation rate Λ to the magnetic field amplitude H. The case of m = 0 is preferred if this ratio ranges below a critical value, which is a function of the shell thickness, and |m| = l otherwise.     

Solitary waves in a compressible fluid

Evolution equations for long nonlinear internal waves in a compressible fluid are derived, with the aim of comparing these equations with their counterparts in an incompressible fluid. Both the Korteweg-de Vries equation, and the deep fluid equation are discussed, for both dry and moist atmospheres. It is shown that the effects of compressibility, or non-Boussinesq terms, are generally small, but measurable, and are manifested mainly in the nonlinear term of the evolution equation. For the case of a moist atmosphere the effect of a gain in energy by latent heat release is compared with the energy lost by radiation damping.     

Steady-state convection in a rotating long cavity

Abstract

Convection experiments in a rotating long cavity were conducted to investigate the changes introduced by Coriolis-force to the well understood convection in the (nonrotating) long cavity. In the far convective regime, the results indicated that large-scale features of the convection were dominated by geostrophically balanced confined primary currents. Close to the transition regime, where the intruding currents gain a thickness close to the vertical dimension of the tank, secondary features became predominant. In the transition regime, these secondary features were not present and the simple picture of two cavity filling counter currents with laterally tilted isotherms could be confirmed.     

Travelling wave convection in a rotating layer

Abstract

Small amplitude two-dimensional Boussinesq convection in a plane layer with stress-free boundaries rotating uniformly about the vertical is studied. A horizontally unbounded layer is modelled by periodic boundary conditions. When the centrifugal force is balanced by an appropriate pressure gradient the resulting equations are translation invariant, and overstable convection can take the form of travelling waves. In the Prandtl number regime 0.53 < [sgrave] < 0.68 such solutions are preferred over the more usual standing waves. For [sgrave] < 0.53, travelling waves are stable provided the Taylor number is sufficiently large.     

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.     

Non-linear marginal convection in a rotating magnetic system

Abstract

An inviscid, electrically conducting fluid is contained between two rigid horizontal planes and bounded laterally by two vertical walls. The fluid is permeated by a strong uniform horizontal magnetic field aligned with the side wall boundaries and the entire system rotates rapidly about a vertical axis. The ratio of the magnitudes of the Lorentz and Coriolis forces is characterized by the Elsasser number, A, and the ratio of the thermal and magnetic diffusivities, q. By heating the fluid from below and cooling from above the system becomes unstable to small perturbations when the adverse density gradient as measured by the Rayleigh number, R, is sufficiently large.

With the viscosity ignored the geostrophic velocity, U, which is aligned with the applied magnetic field, is independent of the coordinate parallel to the rotation axis but is an arbitrary function of the horizontal cross-stream coordinate. At the onset of instability the value of U taken ensures that Taylor's condition is met. Specifically the Lorentz force, which results from marginal convection must not cause any acceleration of the geostrophic flow. It is found that the critical Rayleigh number characterising the onset of instability is generally close to the corresponding value for the usual linear problem, in which Taylor's condition is ignored and U is chosen to vanish. Significant differences can occur when q is small owing to a complicated flow structure. There is a central interior region in which the local magnetic Reynolds number, R_m , based on U is small of order q and on exterior region in which R_m is of order unity.     

Two-and three-dimensional linear convection in a rotating annulus

Abstract

An exceptional case to the model-independent theory of Knobloch (1995) is presented, by investigating a rotating cylindrical annulus of height H and side wall radii r _o and r _i, with non-slip, perfectly thermally conducting side walls and thermally insulating stress-free ends. Radial heating permits the possibility of either two- or three-dimensional convective solutions being the preferred mode. An analytical solution is obtained for the two-dimensional case and a numerical solution for the three-dimensional solution, which is also applied to the two-dimensional solution. It is shown that both two- and three-dimensional solutions can be realized depending on the aspect ratio, γ = H/d, where d = r _o-r _i is the thickness of the annulus, the radii ratio λ = r _i/r _o and the rotation rate of the model. For γ = O(1) and λ = 0.4, the preferred convective solution is three-dimensional when the Taylor number, T < 10² and two-dimensional for T > 10². For small aspect ratios, γ ? 1, the preferred mode is two-dimensional for all rotation rates.     

The catastrophe structure of thermohaline convection in a two-dimensional fluid model and a comparison with low-order box models

Abstract

We impose a surface forcing on the 2D, Boussinesq, thermohaline equations in a rectangular domain, in the form of equatorially symmetric cosine distributions of salinity flux and temperature. This system may be seen as an idealization of the ocean thermohaline circulation on the global scale over intervals of centuries or millenia. Multiple steady states are found numerically. They reflect the competition between the opposite signs of the temperature and salinity-driven equatorially symmetric circulations. There are also pole-to-pole, equatorially asymmetric circulations. In the control space of the temperature and salinity-flux forcing amplitudes, these equilibria form two cusp catastrophes, and transitions between stable equilibria occur through several distinct bifurcations. These catastrophes can be reproduced in simple box models connecting stirred reservoirs through capillary pipes. This steady-state analysis may provide a framework for a better understanding of climatic transitions between different stable regimes of the ocean-atmosphere system.     

Thermal convection in differentially rotating systems

Abstract

The onset of convection in a cylindrical fluid annulus is analyzed in the case when the cylindrical walls are rotating differentially, a temperature gradient in the radial direction is applied, and the centrifugal force dominates over gravity. The small gap approximation is used and no-slip conditions on the cylindrical walls are assumed. It is found that over a considerable range of the parameter space either convection rolls aligned with the axis of rotation or rolls in the perpendicular (azimuthal) direction are preferred. It is shown that by a suitable redefinition of parameters, results for finite amplitude Taylor vortices and for convection rolls in the presence of shear can be applied to the present problem. Weakly nonlinear results for transverse rolls in a Couette flow indicate the possibility of subcritical bifurcation for Prandtl numbers P less than 0.82. Heat and momentum transports are derived as functions of P and the problem of interaction between transverse and longitudinal rolls is considered. The relevance of the analysis for problems of convection in planetary and stellar atmospheres is briefly discussed.     

Transition to two-dimensional turbulent convection in a rapidly-rotating annulus

Abstract

We study a semi-analytical model of convection in a rapidly-rotating, differentially-heated annulus with sloping top and bottom lids. Rapid rotation leads to a preservation of relatively simple, two-dimensional (2-D) structure in the experimentally-observed flow, while temporal complexity increases with the Rayleigh number. The model is, therefore, two-dimensional; it exhibits a sequence of bifurcations from steadily-drifting, azimuthally-periodic convection columns, also called thermal Rossby waves, through vacillation and a period-doubling cascade, to aperiodic, weakly-turbulent solutions.

Our semi-analytical results match to within a few percent previous numerical results with a limited-resolution 2-D model, and extend these results, due to the greater flexibility of the model presented here. Two types of vacillation are obtained, which we call, by analogy with classical nomenclature of the baroclinic annulus with moderate rotation rates, amplitude vacillation and tilted-trough vacillation. Their properties and dependence on the problem's nondimensional parameters are investigated. The period-doubling cascade for each type of vacillation is studied in some detail.     

Wave propagation in a uniformly rotating fluid

Summary An unsteady flow generated by a harmonically oscillating pressure distribution of frequency acting on the paraboloidal free surface of an inviscid, incompressible fluid rotating with uniform angular velocity has been investigated. It is shown that case (i), >2 , corresponds to the usual surface waves, and case (ii), <2 , in contrast to the surface waves, corresponds to the inertial waves which are originated entirely due to rotation and have no counterpart in a non-rotating fluid motion. An explicit solution of the problem related to the above cases are obtained by the joint Laplace and Hankel transforms treatment in conjunction with asymptotic methods. The significant effects of the Coriolis force and the curvature of the free surface on the wave motions have been investigated. A comparison is made between the solutions of the problems with the horizontal and the paraboloidal free surface curvature. The analysis is concluded by exihibiting the characteristic features of the wave motions.     

Frontal upwelling in a rotating two-layer fluid

Abstract

This paper describes the source-sink driven flow in a two-layer fluid confined in a rotating annulus. Light fluid is injected at the inner wall, while denser fluid is withdrawn at the outer wall. The interface between the immiscible fluids intersects the bottom and thus produces a front. The net transport from the source to the sink is carried by Ekman layers at the bottom and at the interface, and by Stewartson layers at the side walls. A detached Stewartson layer arises at the front, leading to a pronounced upwelling circulation.     

Nonlinear dynamics of boussinesq convection in a deep rotating spherical shell-i

Abstract

Convection in a rotating spherical shell has wide application for understanding the dynamics of the atmospheres and interiors of many celestial bodies. In this paper we review linear results for convection in a shell of finite depth at substantial but not asymptotically large Taylor numbers, present nonlinear multimode calculations for similar conditions, and discuss the model and results in the context of the problem of solar convection and differential rotation. Detailed nonlinear calculations are presented for Taylor number T = 10⁵, Prandtl number P = 1, and Rayleigh number R between 1 |MX 10⁴ and 4 |MX 10⁴ (which is between about 4 and 16 times critical) for a shell of depth 20% of the outer radius. Sixteen longitudinal wave numbers are usually included (all even wave numbers m between 0 and 30) the amplitudes of which are computed on a staggered grid in the meridian plane.

The kinetic energy spectrum shows a peak in the wave number range m = 12–18 at R = 10⁴, which straddles the critical wave number m = 14 predicted by linear theory. These are modes which peak near the equator. The spectrum shows a second strong peak at m = 0, which represents the differential rotation driven by the peak convective modes. As R is increased, the amplitude of low wave numbers increases relative to high wave numbers as convection fills in in high and middle latitudes, and as the longitudinal scale of equatorial convection grows. By R = 3 |MX 10⁴, m = 8 is the peak convective mode. There is a clear minimum in the total kinetic energy at middle latitudes relative to low and high, well into the nonlinear regime, representing the continued dominance of equatorial and polar modes found in the linear case. The kinetic energy spectrum for m > 0 is maintained primarily by buoyancy work in each mode, but with substantial nonlinear transfer of kinetic energy from the peak modes to both lower and higher wave numbers.

For R = 1 to 2 |MX 10⁴, the differential rotation takes the form of an equatorial acceleration, with angular velocity generally decreasing with latitude away from the equator (as on the sun) and decreasing inwards. By R = 4 |MX 10⁴, this equatorial profile has completely reversed, with angular velocity increasing with depth and latitude. Also, a polar vortex which has positive rotation relative to the reference frame (no evidence of which has been seen on the sun) builds up as soon as polar modes become important. Meridional circulation is quite weak relative to differential rotation at R = 10⁴, but grows relative to it as R is increased. This circulation takes the farm of a single cell of large latitudinal extent in equatorial regions, with upward flow near the equator, together with a series of narrower cells in high latitudes. It is maintained primarily by axisymmetric buoyancy forces. The differential rotation is maintained at all R primarily by Reynolds stresses, rather than meridional circulation. Angular momentum transport toward the equator for R = 1–2 |MX 10⁴ maintains the equatorial acceleration while radially inward transport maintains the opposite profile at R = 4 |MX 10⁴.

The total heat flux out the top of the convective shell always shows two peaks for the range of R studied, one at the equator and the other near the poles (no significant variation with latitude is seen on the sun), while heat flux in at the bottom shows only a polar peak at large R. The meridional circulation and convective cells transport heat toward the equator to maintain this difference.

The helicity of the convection plus the differential rotation produced by it suggest the system may be capable of driving a field reversing dynamo, but the toroidal field may migrate with lime in each cycle toward the poles and equator, rather than just toward the equator as apparently occurs on the sun.

We finally outline additions to the physics of the model to make it more realistic for solar application.     

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.     

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.