A spectral method for the wave equation of divergence-free vectors and symmetric tensors inside a sphere
arXiv:0905.2048 · doi:10.1016/j.jcp.2009.09.033
Abstract
The wave equation for vectors and symmetric tensors in spherical coordinates is studied under the divergence-free constraint. We describe a numerical method, based on the spectral decomposition of vector/tensor components onto spherical harmonics, that allows for the evolution of only those scalar fields which correspond to the divergence-free degrees of freedom of the vector/tensor. The full vector/tensor field is recovered at each time-step from these two (in the vector case), or three (symmetric tensor case) scalar fields, through the solution of a first-order system of ordinary differential equations (ODE) for each spherical harmonic. The correspondence with the poloidal-toroidal decomposition is shown for the vector case. Numerical tests are presented using an explicit Chebyshev-tau method for the radial coordinate.
25 pages, 4 figures, to appear in Journal of Computational Physics
References in corpus (7)
- Spectral Methods for Numerical Relativity
- A constrained scheme for Einstein equations based on Dirac gauge and spherical coordinates
- Implementation of higher-order absorbing boundary conditions for the Einstein equations
- Excised black hole spacetimes: quasi-local horizon formalism applied to the Kerr example
- Poloidal-toroidal decomposition in a finite cylinder. I. Influence matrices for the magnetohydrodynamic equations
- Poloidal-toroidal decomposition in a finite cylinder. II. Discretization, regularization and validation
- MHD of rotating compact stars with spectral methods: description of the algorithm and tests