A time-stepping dynamically-consistent spherical-shell dynamo code

D. J. Ivers

Abstract


A pseudo-spectral dynamo code, developed as a computational laboratory, is described. The magnetic, heat and Boussinesq Navier-Stokes equations, with inertia, non-linear advection, buoyancy with asymmetric gravity, Coriolis, viscous and Lorentz forces, are solved numerically in a rotating conducting fluid shell. The convection is thermally driven by prescribed boundary temperatures. The equations are discretised using toroidal-poloidal fields, Chebychev collocation in radius and spherical harmonic expansion in angles. Derivatives are performed spectrally. Products are evaluated in physical space for efficiency. Fields are transformed between physical and spectral spaces by fast Fourier and Gauss-Legendre methods. Linear terms are time-stepped implicitly and product terms explicitly using an Adams predictor/corrector. Results are presented for two benchmark models.

Full Text:

PDF BibTeX


DOI: https://doi.org/10.21914/anziamj.v44i0.688



Remember, for most actions you have to record/upload into this online system
and then inform the editor/author via clicking on an email icon or Completion button.
ANZIAM Journal, ISSN 1446-8735, copyright Australian Mathematical Society.