A pseudospectral matrix method for time-dependent tensor fields on a spherical shell
arXiv:1104.3408 · doi:10.1016/j.jcp.2012.11.007
Abstract
We construct a pseudospectral method for the solution of time-dependent, non-linear partial differential equations on a three-dimensional spherical shell. The problem we address is the treatment of tensor fields on the sphere. As a test case we consider the evolution of a single black hole in numerical general relativity. A natural strategy would be the expansion in tensor spherical harmonics in spherical coordinates. Instead, we consider the simpler and potentially more efficient possibility of a double Fourier expansion on the sphere for tensors in Cartesian coordinates. As usual for the double Fourier method, we employ a filter to address time-step limitations and certain stability issues. We find that a tensor filter based on spin-weighted spherical harmonics is successful, while two simplified, non-spin-weighted filters do not lead to stable evolutions. The derivatives and the filter are implemented by matrix multiplication for efficiency. A key technical point is the construction of a matrix multiplication method for the spin-weighted spherical harmonic filter. As example for the efficient parallelization of the double Fourier, spin-weighted filter method we discuss an implementation on a GPU, which achieves a speed-up of up to a factor of 20 compared to a single core CPU implementation.
33 pages, 9 figures
References in corpus (7)
- Calibration of Moving Puncture Simulations
- Simulation of Binary Black Hole Spacetimes with a Harmonic Evolution Scheme
- Testing outer boundary treatments for the Einstein equations
- An explicit harmonic code for black-hole evolution using excision
- Implementation of higher-order absorbing boundary conditions for the Einstein equations
- Black hole evolution with the BSSN system by pseudo-spectral methods
- Outer boundary conditions for Einstein's field equations in harmonic coordinates
Cited by in corpus (21)
- Continuum and Discrete Initial-Boundary-Value Problems and Einstein's Field Equations
- Numerical solution of the 2+1 Teukolsky equation on a hyperboloidal and horizon penetrating foliation of Kerr and application to late-time decays
- A Pseudospectral Method for Gravitational Wave Collapse
- Asymptotic gravitational wave fluxes from a spinning particle in circular equatorial orbits around a rotating black hole
- Hyperbolic Relaxation Method for Elliptic Equations
- Critical phenomena in the collapse of gravitational waves
- The new discontinuous Galerkin methods based numerical relativity program Nmesh
- Numerical evolutions of fields on the 2-sphere using a spectral method based on spin-weighted spherical harmonics
- The Hyperboloidal Numerical Evolution of a Good-Bad-Ugly Wave Equation
- Evolution of Brill waves with an adaptive pseudospectral method
- Numerical solutions of Einstein's equations for cosmological spacetimes with spatial topology S3 and symmetry group U(1)
- Adaptive hp-Refinement for Spectral Elements in Numerical Relativity
- Boson star head-on collisions with constraint-violating and constraint-satisfying initial data
- Twist-free axisymmetric critical collapse of a complex scalar field
- An Implementation of DF-GHG with Application to Spherical Black Hole Excision
- Solving the Hamiltonian constraint for 1+log trumpets
- Analytical and numerical treatment of perturbed black holes in horizon-penetrating coordinates
- GPU-Accelerated Simulations of Isolated Black Holes
- Well-balanced high order finite difference WENO schemes for a first-order Z4 formulation of the Einstein field equations
- Dual Foliation Formulations of General Relativity
- Twist and higher modes of a complex scalar field at the threshold of collapse