Accurate and efficient Jones-Worland spectral transforms for planetary applications

Spectral transforms between physical space and spectral space are needed for fluid dynamical calculations in the whole sphere, representative of a planetary core. In order to construct a representation that is everywhere smooth, regular and differentiable, special polynomials called Jones-Worland polynomials, based on a type of Jacobi polynomial, are used for the radial expansion, coupled to spherical harmonics in angular variables. We present an exact, efficient transform that is partly based on the FFT and which remains accurate in finite precision. Application is to high-resolution solutions of the Navier-Stokes equation, possibly coupled to the heat transfer and induction equations. Expected implementations would be in simulations with P3 degrees of freedom, where P may be greater than 103. Memory use remains modest at high spatial resolution, indeed typically P times lower than competing algorithms based on quadrature.

[1]  Taylor state dynamos found by optimal control: axisymmetric examples , 2018, Journal of Fluid Mechanics.

[2]  Daniel Potts,et al.  Fast algorithms for discrete polynomial transforms on arbitrary grids , 2003 .

[3]  Jie Shen,et al.  Fast structured Jacobi-Jacobi transforms , 2018, Math. Comput..

[4]  A. Jackson,et al.  A fully spectral methodology for magnetohydrodynamic calculations in a whole sphere , 2016, J. Comput. Phys..

[5]  Keaton J. Burns,et al.  Tensor calculus in spherical coordinates using Jacobi polynomials. Part-II: Implementation and examples , 2018, J. Comput. Phys. X.

[6]  Nicholas Hale,et al.  A Fast, Simple, and Stable Chebyshev-Legendre Transform Using an Asymptotic Formula , 2014, SIAM J. Sci. Comput..

[7]  C. Barenghi,et al.  Magnetoconvection in a rapidly rotating sphere , 1997 .

[8]  Andrew Jackson,et al.  An optimal Galerkin scheme to solve the kinematic dynamo eigenvalue problem in a full sphere , 2010, J. Comput. Phys..

[9]  J. Aurnou,et al.  The influence of magnetic fields in planetary dynamo models , 2012 .

[10]  Jan S. Hesthaven,et al.  Computation of connection coefficients and measure modifications for orthogonal polynomials , 2012 .

[11]  Ulrich R. Christensen,et al.  Numerical modelling of the geodynamo: a systematic parameter study , 1999 .

[12]  Daniel Lecoanet,et al.  Tensor calculus in spherical coordinates using Jacobi polynomials. Part-I: Mathematical analysis and derivations , 2018, J. Comput. Phys. X.

[13]  A. Jackson,et al.  Precession-driven dynamos in a full sphere and the role of large scale cyclonic vortices , 2016, 1606.03230.

[14]  W. Verkley,et al.  A Spectral Model for Two-Dimensional Incompressible Fluid Flow in a Circular Basin , 1997 .

[15]  J. Guermond,et al.  Full sphere hydrodynamic and dynamo benchmarks , 2014, Geophysical Journal International.

[16]  A. Hofmann,et al.  Geodynamo, Solar Wind, and Magnetopause 3.4 to 3.45 Billion Years Ago , 2010, Science.

[17]  K. Li,et al.  The optimal kinematic dynamo driven by steady flows in a sphere , 2018, Journal of Fluid Mechanics.

[18]  Richard Mikael Slevinsky,et al.  On the use of Hahn's asymptotic formula and stabilized recurrence for a fast, simple, and stable Chebyshev--Jacobi transform , 2016, 1602.02618.

[19]  Philip W. Livermore,et al.  Spectral radial basis functions for full sphere computations , 2007, J. Comput. Phys..

[20]  G. Glatzmaier,et al.  A three-dimensional convective dynamo solution with rotating and finitely conducting inner core and mantle , 1995 .