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

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

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

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

Asymptotic analysis and analytical solutions of a model of cardiac excitation.

The original publication is available at www.springerlink.com - http://link.springer.com/article/10.1007/s11538-007-9267-0Journal ArticleCopyright © SpringerWe describe an asymptotic approach to gated ionic models of single-cell cardiac excitability. It has a form essentially different from the Tikhonov fast-slow form assumed in standard asymptotic reductions of excitable systems. This is of interest since the standard approaches have been previously found inadequate to describe phenomena such as the dissipation of cardiac wave fronts and the shape of action potential at repolarization. The proposed asymptotic description overcomes these deficiencies by allowing, among other non-Tikhonov features, that a dynamical variable may change its character from fast to slow within a single solution. The general asymptotic approach is best demonstrated on an example which should be both simple and generic. The classical model of Purkinje fibers (Noble in J. Physiol. 160:317-352, 1962) has the simplest functional form of all cardiac models but according to the current understanding it assigns a physiologically incorrect role to the Na current. This leads us to suggest an "Archetypal Model" with the simplicity of the Noble model but with a structure more typical to contemporary cardiac models. We demonstrate that the Archetypal Model admits a complete asymptotic solution in quadratures. To validate our asymptotic approach, we proceed to consider an exactly solvable "caricature" of the Archetypal Model and demonstrate that the asymptotic of its exact solution coincides with the solutions obtained by substituting the "caricature" right-hand sides into the asymptotic solution of the generic Archetypal Model. This is necessary, because, unlike in standard asymptotic descriptions, no general results exist which can guarantee the proximity of the non-Tikhonov asymptotic solutions to the solutions of the corresponding detailed ionic model

How far can minimal models explain the solar cycle?

A physically consistent model of magnetic field generation by convection in a rotating spherical shell with a minimum of parameters is applied to the Sun. Despite its unrealistic features the model exhibits a number of properties resembling those observed on the Sun. The model suggests that the large scale solar dynamo is dominated by a non-axisymmetric $m=1$ component of the magnetic field.Comment: Accepted for publication in the Astrophysical Journal on 2012/01/3

A spherical shell numerical dynamo benchmark with pseudo vacuum magnetic boundary conditions

It is frequently considered that many planetary magnetic fields originate as a result of convection within planetary cores. Buoyancy forces responsible for driving the convection generate a fluid flow that is able to induce magnetic fields; numerous sophisticated computer codes are able to simulate the dynamic behaviour of such systems. This paper reports the results of a community activity aimed at comparing numerical results of several different types of computer codes that are capable of solving the equations of momentum transfer, magnetic field generation and heat transfer in the setting of a spherical shell, namely a sphere containing an inner core. The electrically conducting fluid is incompressible and rapidly rotating and the forcing of the flow is thermal convection under the Boussinesq approximation. We follow the original specifications and results reported in Harder & Hansen to construct a specific benchmark in which the boundaries of the fluid are taken to be impenetrable, non-slip and isothermal, with the added boundary condition for the magnetic field <b>B</b> that the field must be entirely radial there; this type of boundary condition for <b>B</b> is frequently referred to as ‘pseudo-vacuum’. This latter condition should be compared with the more frequently used insulating boundary condition. This benchmark is so-defined in order that computer codes based on local methods, such as finite element, finite volume or finite differences, can handle the boundary condition with ease. The defined benchmark, governed by specific choices of the Roberts, magnetic Rossby, Rayleigh and Ekman numbers, possesses a simple solution that is steady in an azimuthally drifting frame of reference, thus allowing easy comparison among results. 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 among codes