An efficient numerical scheme for a 3D spherical dynamo equation

We develop an efficient numerical scheme for the 3D mean-field spherical dynamo equation. The scheme is based on a semi-implicit discretization in time and a spectral method in space based on the divergence-free spherical harmonic functions. A special semi-implicit approach is proposed such that at each time step one only needs to solve a linear system with constant coefficients. Then, using expansion in divergence-free spherical harmonic functions in the transverse directions allows us to reduce the linear system at each time step to a sequence of one-dimensional equations in the radial direction, which can then be efficiently solved by using a spectral-element method. We show that the solution of fully discretized scheme remains bounded independent of the number of unknowns, and present numerical results to validate our scheme.

[1]  M D Nornberg,et al.  Numerical simulations of current generation and dynamo excitation in a mechanically forced turbulent flow. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[2]  Mohammad M. Rahman,et al.  A spectral solution of nonlinear mean field dynamo equations: With inertia , 2009, Comput. Math. Appl..

[3]  J. Nédélec.  Acoustic and electromagnetic equations , 2001 .

[4]  Ohannes A. Karakashian,et al.  A Nonconforming Finite Element Method for the Stationary Navier--Stokes Equations , 1998 .

[5]  T. Clune,et al.  Numerical simulation of a spherical dynamo excited by a flow of von Kármán type , 2010 .

[6]  J. Zou,et al.  A nonlinear vacillating dynamo induced by an electrically heterogeneous mantle , 2001 .

[7]  Rainer Hollerbach.  ON THE THEORY OF THE GEODYNAMO , 1996 .

[8]  Direct numerical simulation of pipe flow using a solenoidal spectral method , 2012 .

[9]  Ohannes A. Karakashian,et al.  Piecewise solenoidal vector fields and the Stokes problem , 1990 .

[10]  Claus-Dieter Munz,et al.  Divergence Correction Techniques for Maxwell Solvers Based on a Hyperbolic Model , 2000 .

[11]  K. Yee.  Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media , 1966 .

[12]  D. Moss.  Numerical simulation of the Gailitis dynamo , 2006 .

[13]  F. Busse,et al.  Convection driven magnetohydrodynamic dynamos in rotating spherical shells , 1989 .

[14]  Jeremy Bloxham,et al.  Numerical Modeling of Magnetohydrodynamic Convection in a Rapidly Rotating Spherical Shell , 1999 .

[15]  Jun Zou,et al.  A Three-dimensional Spherical Nonlinear Interface Dynamo , 2003 .

[16]  C. Guervilly,et al.  Numerical simulations of dynamos generated in spherical Couette flows , 2010, 1010.3859.

[17]  E. L. Hill.  The Theory of Vector Spherical Harmonics , 1954 .

[18]  Roger Temam,et al.  Some mathematical questions related to the MHD equations , 1983 .

[19]  C. Jones,et al.  Planetary Magnetic Fields and Fluid Dynamos , 2011 .

[20]  Rainer Hollerbach,et al.  A spectral solution of the magneto-convection equations in spherical geometry , 2000 .

[21]  Jun Zou,et al.  Spherical Interface Dynamos: Mathematical Theory, Finite Element Approximation, and Application , 2006, SIAM J. Numer. Anal..

[22]  Edward Crisp Bullard,et al.  Homogeneous dynamos and terrestrial magnetism , 1954, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.