Browse by author
Lookup NU author(s): Dr Dr DG MMillan McMillan, Dr Graeme Sarson
Full text for this publication is not currently held within this repository. Alternative links are provided below where available.
Most numerical simulations of the geodynamo are cast in spherical geometry, using a spherical harmonic representation for lateral variations, and an expansion in Chebyshev polynomials, or discretisation, in radius. A number of research groups have produced three-dimensional, time-dependent, self-consistent simulations of the geodynamo using this pseudospectral methodology. Computational limitations place a practical bound on the parameter regime that can be explored in this context, however, and Earth-like values currently remain several orders of magnitude out of reach. For the spherically pseudospectral codes, the absence of a fast Legendre transform is a strong factor contributing to this limitation. As a first step toward alternative computational methods for geodynamo modelling, we adapt a pre-existing, efficiently parallelised magnetohydrodynamics (MHD) code, to model spherical shell problems. The Pencil-Code is a Cartesian code that uses sixth-order finite differences, applied to 'pencils' (i.e. array sections) in the x direction in a cache-efficient way. The domain is tiled in the y and z directions, with the communication of boundary elements handled by Message Passing Interface (MPI). Time stepping is via a third order Runge-Kutta method. The code's modular structure allows a flexible selection of various physical processes and variables, making it easily adaptable for many types of compressible MHD problems. We demonstrate dynamo action driven by thermal convection in a spherical shell of compressible ideal gas, for comparison with the ideal gas dynamos of Kageyama and Sato [Kageyama, A., Sato, T., the Complexity Simulation Group, 1995. Computer simulation of a magnetohydrodynamic dynamo. II. Phys. Fluids, B 2, 2793-2805]. As with the previous results, our solutions are characterised by steady state thermal convection with axial columnar cells, and exponential growth and subsequent saturation of magnetic energy. Minor differences in the details of the benchmark and current simulations are attributed to differences in the computational methods and geometry of the domains involved. We discuss the further modifications still needed for a fuller adaptation to geodynamo problems. © 2005 Elsevier B.V. All rights reserved.
Author(s): McMillan DG, Sarson GR
Publication type: Article
Publication status: Published
Journal: Physics of the Earth and Planetary Interiors
Year: 2005
Volume: 153
Issue: 1-3
Pages: 124-135
ISSN (print): 0031-9201
ISSN (electronic): 1872-7395
Publisher: Elsevier BV
URL: http://dx.doi.org/10.1016/j.pepi.2005.03.012
DOI: 10.1016/j.pepi.2005.03.012
Altmetrics provided by Altmetric
