Double-diffusive convection in a rotating cylindrical annulus with conical caps

Double-diffusive convection driven by both thermal and compositional buoyancy in a rotating cylindrical annulus with conical caps is considered with the aim to establish whether a small fraction of compositional buoyancy added to the thermal buoyancy (or vice versa) can significantly reduce the critical Rayleigh number and amplify convection in planetary cores. It is shown that the neutral surface describing the onset of convection in the double-buoyancy case is essentially different from that of the well-studied purely thermal case, and does indeed allow the possibility of low-Rayleigh number convection. In particular, isolated islands of instability are formed by an additional "double-diffusive" eigenmode in certain regions of the parameter space. However, the amplitude of such low-Rayleigh number convection is relatively weak. At similar flow amplitudes purely compositional and double-diffusive cases are characterized by a stronger time dependence compared to purely thermal cases, and by a prograde mean zonal flow near the inner cylindrical surface. Implications of the results for planetary core convection are briefly discussed.Comment: Accepted for publication in Physics of the Earth and Planetary Interiors on 20 April 201

Prandtl-number dependence of convection-driven dynamos in rotating spherical fluid shells

The value of the Prandtl number P exerts a strong influence on convection-driven dynamos in rotating spherical shells filled with electrically conducting fluids. Low Prandtl numbers promote dynamo action through the shear provided by differential rotation, while the generation of magnetic fields is more difficult to sustain in high-Prandtl-number fluids where higher values of the magnetic Prandtl number Pm are required. The magnetostrophic approximation often used in dynamo theory appears to be valid only for relatively high values of P and Pm. Dynamos with a minimum value of Pm seem to be most readily realizable in the presence of convection columns at moderately low values of P. The structure of the magnetic field varies strongly with P in that dynamos with a strong axial dipole field are found for high values of P while the energy of this component is exceeded by that of the axisymmetric toroidal field and by that of the non-axisymmetric components at low values of P. Some conclusions are discussed in relation to the problem of the generation of planetary magnetic fields by motions in their electrically conducting liquid cores

Inertial convection in rotating fluid spheres

The onset of convection in the form of inertial waves in a rotating fluid sphere is studied through a perturbation analysis in an extension of earlier work by Zhang (1994). Explicit expressions for the dependence of the Rayleigh number on the azimuthal wavenumber are derived and new results for the case of a nearly thermally insulating boundary are obtained

Magneto-inertial convection in rotating fluid spheres

The onset of convection in the form of magneto-inertial waves in a rotating fluid sphere permeated by a constant axial electric current is studied through a perturbation analysis. Explicit expressions for the dependence of the Rayleigh number on the azimuthal wavenumber are derived in the limit of high thermal diffusivity. Results for the cases of thermally infinitely conducting and of nearly thermally insulating boundaries are obtained.Comment: 10 pages, 5 figures, to be submitted for publicatio

Asymptotic properties of mathematical models of excitability

We analyse small parameters in selected models of biological excitability, including Hodgkin-Huxley (1952) model of nerve axon, Noble (1962) model of heart Purkinje fibres, and Courtemanche et al. (1998) model of human atrial cells. Some of the small parameters are responsible for differences in the characteristic timescales of dynamic variables, as in the traditional singular perturbation approaches. Others appear in a way which makes the standard approaches inapplicable. We apply this analysis to study the behaviour of fronts of excitation waves in spatially-extended cardiac models. Suppressing the excitability of the tissue leads to a decrease in the propagation speed, but only to a certain limit; further suppression blocks active propagation and leads to a passive diffusive spread of voltage. Such a dissipation may happen if a front propagates into a tissue recovering after a previous wave, e.g. re-entry. A dissipated front does not recover even when the excitability restores. This has no analogy in FitzHugh-Nagumo model and its variants, where fronts can stop and then start again. In two spatial dimensions, dissipation accounts for break-ups and self-termination of re-entrant waves in excitable media with Courtemanche et al. (1998) kinetics.Comment: 15 pages, 8 figures, to appear in Phil Trans Roy Soc London

Kinetic energy cascades in quasi-geostrophic convection in a spherical shell

We consider triadic nonlinear interaction in the Navier-Stokes equation for quasi-geostrophic convection in a spherical shell. This approach helps understanding the origin of kinetic energy transport in the system and the particular scheme of mode interaction, as well as the locality of the energy transfer. The peculiarity of convection in the sphere, concerned with excitation of Rossby waves, is considered. The obtained results are compared with our previous study in Cartesian geometry

Reynolds stresses and mean fields generated by pure waves: applications to shear flows and convection in a rotating shell

A general reformulation of the Reynolds stresses created by two-dimensional waves breaking a translational or a rotational invariance is described. This reformulation emphasizes the importance of a geometrical factor: the slope of the separatrices of the wave flow. Its physical relevance is illustrated by two model systems: waves destabilizing open shear flows; and thermal Rossby waves in spherical shell convection with rotation. In the case of shear-flow waves, a new expression of the Reynolds–Orr amplification mechanism is obtained, and a good understanding of the form of the mean pressure and velocity fields created by weakly nonlinear waves is gained. In the case of thermal Rossby waves, results of a three-dimensional code using no-slip boundary conditions are presented in the nonlinear regime, and compared with those of a two-dimensional quasi-geostrophic model. A semi-quantitative agreement is obtained on the flow amplitudes, but discrepancies are observed concerning the nonlinear frequency shifts. With the quasi-geostrophic model we also revisit a geometrical formula proposed by Zhang to interpret the form of the zonal flow created by the waves, and explore the very low Ekman-number regime. A change in the nature of the wave bifurcation, from supercritical to subcritical, is found

Full sphere hydrodynamic and dynamo benchmarks

Convection in planetary cores can generate fluid flow and magnetic fields, and a number of sophisticated codes exist to simulate the dynamic behaviour of such systems. We report on the first community activity to compare numerical results of computer codes designed to calculate fluid flow within a whole sphere. The flows are incompressible and rapidly rotating and the forcing of the flow is either due to thermal convection or due to moving boundaries. All problems defined have solutions that allow easy comparison, since they are either steady, slowly drifting or perfectly periodic. The first two benchmarks are defined based on uniform internal heating within the sphere under the Boussinesq approximation with boundary conditions that are uniform in temperature and stress-free for the flow. Benchmark 1 is purely hydrodynamic, and has a drifting solution. Benchmark 2 is a magnetohydrodynamic benchmark that can generate oscillatory, purely periodic, flows and magnetic fields. In contrast, Benchmark 3 is a hydrodynamic rotating bubble benchmark using no slip boundary conditions that has a stationary solution. Results from a variety of types of code are reported, including codes that are fully spectral (based on spherical harmonic expansions in angular coordinates and polynomial expansions in radius), mixed spectral and finite difference, finite volume, finite element and also a mixed Fourier–finite element code. There is good agreement between codes. It is found that in Benchmarks 1 and 2, the approximation of a whole sphere problem by a domain that is a spherical shell (a sphere possessing an inner core) does not represent an adequate approximation to the system, since the results differ from whole sphere results

Asymptotics of conduction velocity restitution in models of electrical excitation in the heart.

Copyright © Springer 2011Journal ArticleThe original publication is available at www.springerlink.com - http://link.springer.com/article/10.1007/s11538-010-9523-6We extend a non-Tikhonov asymptotic embedding, proposed earlier, for calculation of conduction velocity restitution curves in ionic models of cardiac excitability. Conduction velocity restitution is the simplest non-trivial spatially extended problem in excitable media, and in the case of cardiac tissue it is an important tool for prediction of cardiac arrhythmias and fibrillation. An idealized conduction velocity restitution curve requires solving a non-linear eigenvalue problem with periodic boundary conditions, which in the cardiac case is very stiff and calls for the use of asymptotic methods. We compare asymptotics of restitution curves in four examples, two generic excitable media models, and two ionic cardiac models. The generic models include the classical FitzHugh-Nagumo model and its variation by Barkley. They are treated with standard singular perturbation techniques. The ionic models include a simplified "caricature" of Noble (J. Physiol. Lond. 160:317-352, 1962) model and Beeler and Reuter (J. Physiol. Lond. 268:177-210, 1977) model, which lead to non-Tikhonov problems where known asymptotic results do not apply. The Caricature Noble model is considered with particular care to demonstrate the well-posedness of the corresponding boundary-value problem. The developed method for calculation of conduction velocity restitution is then applied to the Beeler-Reuter model. We discuss new mathematical features appearing in cardiac ionic models and possible applications of the developed method