An optimal Galerkin scheme to solve the kinematic dynamo eigenvalue problem in a full sphere (Q613383)

From the authors' abstract: ``The kinematic dynamo approximation describes the generation of the magnetic field in a prescribed flow of electrically-conducting liquid. In this paper, we demonstrate the numerical superiority of a Galerkin scheme in solving the kinematic dynamo eigenvalue problem in a full sphere. After adopting a poloidal-toroidal decomposition and expanding in spherical harmonics, we express the radial dependence in terms of a basis of exponentially convergent orthogonal polynomials. Each basis function is constructed from a terse sum of one-sided Jacobi polynomials that not only satisfies the boundary conditions of matching to an electrically insulating exterior, but is everywhere infinitely differentiable, including at the origin. This Galerkin method exhibits more rapid convergence, for a given problem size, than any other scheme hitherto reported, as demonstrated by a benchmark of the magnetic diffusion problem and by comparison to numerous kinematic dynamos from the literature.''