You are here: Home / Groups / Dynamo / Benchmarks


Code/Group Name Authors Description Accuracy Results
ASH Mob Nick Featherstone, Mark Miesch, Sacha Brun    
Calypso Hiroaki Matsui Spherical harmonics expansion on sphere and 2nd-order FDM in the radial direction. Crank-Nicolson Scheme is used for the diffusion terms, and 2nd-order Adams-Bashforth scheme is used for the other terms. Vorticity equation and Poisson equation for the toroidal vorticity is used for the time integration of the fluid motion. Received
ETH Code Andrey Sheyko Spectral simulation using spherical harmonics for the angular component and finite differences in radius. The incompressibility condition is guaranteed by the use of a toroidal/poloidal decomposition of the vector fields. A second order predictor-corrector scheme is used for the time integration Received
GeoFEM-MHD Hiroaki Matsui Spatial discretization: Finite-element methods with tri-linear hexahedral elements. Time integration: Fractional step scheme is applied. Crank-Nicolson scheme is used for the diffusion terms, and the other terms are solved by Adams-Bashforth scheme. Pressure and electric potential are solved to satisfy mass conservation and Coulomb gauge for magnetic vector potential. Received
GFD-dennou group Shin-ichi Takehiro, Youhei Sasaki, Yoshi-Yuki Hayashi Fully spectral method. Spherical harmonics and Chebyshev polinomials are used to expand the variables in the horizontal and radial directions, respectively. The nonlinear terms and the Coriolis terms are evaluated in the physical space and convert back to the spectral space (so called transform method). Diffusion terms are integrated with Crank-Nicolson scheme while other terms are done with Adams-Bashforth scheme. Received
Goddard Code Weijia Kuang, Weiyuan Jiang For the momentum equation, we solve the radial component of the velocity, the radial component of the vorticity, and the modified pressure simultaneously. For the induction equation, we solve the radial components of the magnetic field and current simultaneously. In the model, the velocity and the magnetic field are described by the poloidal and toroidal scalars. On the spherical surfaces, spherical harmonic expansions are used. In the radial direction, both compact finite difference algorithms and the Chebyshev collocation method are employed. For this benchmark effort, the Chebyshev collocation method is used. Received
H2000 Rainer Hollerbach Poloidal-toroidal decomposition, then expanded in spherical harmonics in angle (fully de-aliased). In r quantities are expanded in Chebyshev polynomials. The time-stepping is 2nd order Runge-Kutta, with the diffusive terms treated implicitly. Received
Helsinki Pencil Code Petri Käpylä, Joern Warnecke  
IPGP / JHU team Maylis Landeau, Julien Aubert The code uses a spherical harmonic expansion in lateral directions, and finite differences in the radial direction (Dormy et al., 1998). The radial mesh interval decreases in geometrical progression towards the boundaries. Three-point stencil are used for second-order derivatives and five-point stencil for biharmonic operators (second-order accurate). Time integration involves a Crank-Nicolson scheme for diffusion terms and a second-order Adams-Bashforth scheme for other terms. To ensure numerical stability the time step is chosen as the minimum between the characteristic times of advection and Alfven wave propagation in one grid. Received
Leeds Spherical Dynamo (LSD) code David Gubbins, Ashley Willis, Chris Davies, Maggie Avery, Chris Jones The Leeds Spherical Dynamo (LSD) code solves the Boussinesq dynamo equations by representing velocity, v, and magnetic field, B, as poloidal and toroidal scalars. It is pseudospectral;the theta and phi variations are expanded in spherical harmonics, and r variations by variable order, variable mesh, finite differences with non-equidistant grid using Chebyshev zeros as grid points. The nonlinear terms are evaluated by the transform method. Time stepping is by a predictor-corrector method and the time step controlled first by the CFL (Δt < Δx/v) condition, and then modified using information from the predictor step. Received
MagIC Johannes Wicht    
SFEMaNS Jean-Luc Guermond, Francky Luddens, and Caroline Nore Uses quadratic finite elements in the meridian section and Fourier expansions in the azimuthal direction. The time stepping is done with second-order backward Euler. The diffusive terms are treated implicitly. The divergence of the magnetic induction is controlled through a magnetic pressure. Parallelism is done in the meridian section using Petsc and with respect to the Fourier modes using FFTW. Received
Simitev-Busse Dynamo Code Radostin Simitev, Friedrich Busse This is a pseudospectral numerical code for the solution of the governing equations using a poloidal-toroidal representation for the velocity and magnetic fields. Unknown scalar fields are expanded in spherical harmonics in the angular variables and Chebyshev polynomials in radius. Linear and diffusion terms are computed in spectral space, while Coriolis and non-linear terms are computed in physical space and transformed to spectral space at every time step. Time stepping is implemented by a combination of an implicit Crank–Nicolson scheme for the diffusion terms and an explicit Adams–Bashforth scheme for the Coriolis and the non-linear terms; both schemes are second order accurate. Adaptive time-stepping is possible through a modification of this scheme. Received
UCSC Code Gary Glatzmaier Spectral: spherical harmonics and Chebyshev polynomials. Semi-implicit time integration: Adams-Bashforth and Crank-Nicolson. Received
xshells Nathanael Schaeffer xshells uses finite differences in radius and spherical harmonic expansion. The diffusive terms are treated using Cranck-Nicholson scheme, while the other terms are treated explicitely with a second order Adams-Bashford scheme. Non-linear terms are computed in physical space using the SHTns spherical harmonic transform library. Received
Created by Lorraine Hwang Last Modified Thu February 24, 2022 11:09 pm by Lorraine Hwang