Tensor calculus in spherical coordinates using Jacobi polynomials. II: Implementation and examples
From MaRDI portal
Publication:6145374
Abstract: We present a simulation code which can solve broad ranges of partial differential equations in a full sphere. The code expands tensorial variables in a spectral series of spin-weighted spherical harmonics in the angular directions and a scaled Jacobi polynomial basis in the radial direction, as described in Part-I. Nonlinear terms are calculated by transforming from the coefficients in the spectral series to the value of each quantity on the physical grid, where it is easy to calculate products and perform other local operations. The expansion makes it straightforward to solve equations in tensor form (i.e., without decomposition into scalars). We propose and study several unit tests which demonstrate the code can accurately solve linear problems, implement boundary conditions, and transform between spectral and physical space. We then run a series of benchmark problems proposed in Marti et al (2014), implementing the hydrodynamic and magnetohydrodynamic equations. We are able to calculate more accurate solutions than reported in Marti et al 2014 by running at higher spatial resolution and using a higher-order timestepping scheme. We find the rotating convection and convective dynamo benchmark problems depend sensitively on details of timestepping and data analysis. We also demonstrate that in low resolution simulations of the dynamo problem, small changes in a numerical scheme can lead to large changes in the solution. To aid future comparison to these benchmarks, we include the source code used to generate the data, as well as the data and analysis scripts used to generate the figures.
Cites work
- A fully spectral methodology for magnetohydrodynamic calculations in a whole sphere
- A spectral method for polar coordinates
- Chebyshev and Fourier spectral methods.
- Comparing seven spectral methods for interpolation and for solving the Poisson equation in a disk: Zernike polynomials, Logan-Shepp ridge polynomials, Chebyshev-Fourier series, cylindrical Robert functions, Bessel-Fourier expansions, square-to-disk confor
- Cross-code comparisons of mixing during the implosion of dense cylindrical and spherical shells
- Fourier transform summation of Legendre series and \(D\)-functions
- scientific article; zbMATH DE number 1455322 (Why is no real title available?)
- Implicit-Explicit Methods for Time-Dependent Partial Differential Equations
- New tensor spherical harmonics, for application to the partial differential equations of mathematical physics
- Representation of the Elastic - Gravitational Excitation of a Spherical Earth Model by Generalized Spherical Harmonics
- Spectral radial basis functions for full sphere computations
- Tensor calculus in polar coordinates using Jacobi polynomials
- Trigonometric Interpolation of Empirical and Analytical Functions
- Variable step-size implicit-explicit linear multistep methods for time-dependent partial differential equations
Cited in
(6)- Building hierarchies of semiclassical Jacobi polynomials for spectral methods in annuli
- On nonlinear transitions, minimal seeds and exact solutions for the geodynamo
- Rational function approximation with normalized positive denominators
- Waves in the Earth's core. I: Mildly diffusive torsional oscillations
- Thermal convection in the internally heated sphere
- Weakly nonlinear analysis of the onset of convection in rotating spherical shells
This page was built for publication: Tensor calculus in spherical coordinates using Jacobi polynomials. II: Implementation and examples
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6145374)
