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

A. Jackson, A. Sheyko, P. Marti, A. Tilgner, D. Cébron, S. Vantieghem, R. Simitev, F. Busse, X. Zhan, G. Schubert, S. Takehiro, Y. Sasaki, Y. Y. Hayashi, A. Ribeiro, C. Nore, J. L. Guermond

Research output: Contribution to journalArticle

13 Citations (Scopus)

Abstract

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 that the field must be entirely radial there; this type of boundary condition for B is frequently referred to as 'pseudovacuum'. 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.

Original languageEnglish
Pages (from-to)712-723
Number of pages12
JournalGeophysical Journal International
Volume196
Issue number2
DOIs
Publication statusPublished - 21 Nov 2013
Externally publishedYes

Fingerprint

spherical shells
boundary condition
Boundary conditions
Vacuum
shell
boundary conditions
Magnetic fields
vacuum
computer programs
convection
magnetic field
planetary cores
planetary magnetic fields
magnetic fields
Boussinesq approximation
conducting fluids
expansion
Fluids
Momentum transfer
spherical harmonics

Keywords

  • Dynamo
  • Electromagnetic theory
  • Non-linear differential equations
  • theories and simulations

ASJC Scopus subject areas

  • Geochemistry and Petrology
  • Geophysics

Cite this

Jackson, A., Sheyko, A., Marti, P., Tilgner, A., Cébron, D., Vantieghem, S., ... Guermond, J. L. (2013). A spherical shell numerical dynamo benchmark with pseudo-vacuum magnetic boundary conditions. Geophysical Journal International, 196(2), 712-723. https://doi.org/10.1093/gji/ggt425

A spherical shell numerical dynamo benchmark with pseudo-vacuum magnetic boundary conditions. / Jackson, A.; Sheyko, A.; Marti, P.; Tilgner, A.; Cébron, D.; Vantieghem, S.; Simitev, R.; Busse, F.; Zhan, X.; Schubert, G.; Takehiro, S.; Sasaki, Y.; Hayashi, Y. Y.; Ribeiro, A.; Nore, C.; Guermond, J. L.

In: Geophysical Journal International, Vol. 196, No. 2, 21.11.2013, p. 712-723.

Research output: Contribution to journalArticle

Jackson, A, Sheyko, A, Marti, P, Tilgner, A, Cébron, D, Vantieghem, S, Simitev, R, Busse, F, Zhan, X, Schubert, G, Takehiro, S, Sasaki, Y, Hayashi, YY, Ribeiro, A, Nore, C & Guermond, JL 2013, 'A spherical shell numerical dynamo benchmark with pseudo-vacuum magnetic boundary conditions' Geophysical Journal International, vol. 196, no. 2, pp. 712-723. https://doi.org/10.1093/gji/ggt425
Jackson, A. ; Sheyko, A. ; Marti, P. ; Tilgner, A. ; Cébron, D. ; Vantieghem, S. ; Simitev, R. ; Busse, F. ; Zhan, X. ; Schubert, G. ; Takehiro, S. ; Sasaki, Y. ; Hayashi, Y. Y. ; Ribeiro, A. ; Nore, C. ; Guermond, J. L. / A spherical shell numerical dynamo benchmark with pseudo-vacuum magnetic boundary conditions. In: Geophysical Journal International. 2013 ; Vol. 196, No. 2. pp. 712-723.
@article{beff1d3bdd4f4668a3d86b1eb2bd141d,
title = "A spherical shell numerical dynamo benchmark with pseudo-vacuum magnetic boundary conditions",
abstract = "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 that the field must be entirely radial there; this type of boundary condition for B is frequently referred to as 'pseudovacuum'. 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.",
keywords = "Dynamo, Electromagnetic theory, Non-linear differential equations, theories and simulations",
author = "A. Jackson and A. Sheyko and P. Marti and A. Tilgner and D. C{\'e}bron and S. Vantieghem and R. Simitev and F. Busse and X. Zhan and G. Schubert and S. Takehiro and Y. Sasaki and Hayashi, {Y. Y.} and A. Ribeiro and C. Nore and Guermond, {J. L.}",
year = "2013",
month = "11",
day = "21",
doi = "10.1093/gji/ggt425",
language = "English",
volume = "196",
pages = "712--723",
journal = "Geophysical Journal International",
issn = "0956-540X",
publisher = "Oxford University Press",
number = "2",

}

TY - JOUR

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

AU - Jackson, A.

AU - Sheyko, A.

AU - Marti, P.

AU - Tilgner, A.

AU - Cébron, D.

AU - Vantieghem, S.

AU - Simitev, R.

AU - Busse, F.

AU - Zhan, X.

AU - Schubert, G.

AU - Takehiro, S.

AU - Sasaki, Y.

AU - Hayashi, Y. Y.

AU - Ribeiro, A.

AU - Nore, C.

AU - Guermond, J. L.

PY - 2013/11/21

Y1 - 2013/11/21

N2 - 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 that the field must be entirely radial there; this type of boundary condition for B is frequently referred to as 'pseudovacuum'. 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.

AB - 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 that the field must be entirely radial there; this type of boundary condition for B is frequently referred to as 'pseudovacuum'. 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.

KW - Dynamo

KW - Electromagnetic theory

KW - Non-linear differential equations

KW - theories and simulations

U2 - 10.1093/gji/ggt425

DO - 10.1093/gji/ggt425

M3 - Article

VL - 196

SP - 712

EP - 723

JO - Geophysical Journal International

JF - Geophysical Journal International

SN - 0956-540X

IS - 2

ER -