Numerical evolutions of fields on the 2-sphere using a spectral method based on spin-weighted spherical harmonics

Abstract

Many applications in science call for the numerical simulation of systems on manifolds with spherical topology. Through use of integer spin weighted spherical harmonics we present a method which allows for the implementation of arbitrary tensorial evolution equations. Our method combines two numerical techniques that were originally developed with different applications in mind. The first is Huffenberger and Wandelt's spectral decomposition algorithm to perform the mapping from physical to spectral space. The second is the application of Luscombe and Luban's method, to convert numerically divergent linear recursions into stable nonlinear recursions, to the calculation of reduced Wigner d-functions. We give a detailed discussion of the theory and numerical implementation of our algorithm. The properties of our method are investigated by solving the scalar and vectorial advection equation on the sphere, as well as the 2+1 Maxwell equations on a deformed sphere.

New version matches the published version

References in corpus (3)

• Fast and Exact Spin-s Spherical Harmonic Transforms
• On the Use of Multipole Expansion in Time Evolution of Non-linear Dynamical Systems and Some Surprises Related to Superradiance
• A spectral solver for evolution problems with spatial S3-topology

Cited by in corpus (13)

• Spin memory effect for compact binaries in the post-Newtonian approximation
• The Poincaré and BMS flux-balance laws with application to binary systems
• Coadjoint representation of the BMS group on celestial Riemann surfaces
• Gravitational-wave memory effects in Brans-Dicke theory: Waveforms and effects in the post-Newtonian approximation
• Numerical solutions of Einstein's equations for cosmological spacetimes with spatial topology S3 and symmetry group U(1)
• The numerical initial boundary value problem for the generalized conformal field equations
• COFFEE -- An MPI-parallelized Python package for the numerical evolution of differential equations
• Differentiable and accelerated spherical harmonic and Wigner transforms
• The Teukolsky--Starobinsky constants: facts and fictions
• Asymptotically flat vacuum initial data sets from a modified parabolic-hyperbolic formulation of the Einstein vacuum constraint equations
• Asymptotically hyperboloidal initial data sets from a parabolic-hyperbolic formulation of the Einstein vacuum constraints
• Criticality of inhomogeneous Nariai-like cosmological models
• Waveform models for the gravitational-wave memory effect: Extreme mass-ratio limit and final memory offset