Numerical Methods for General Relativistic Radiative Transfer

Kyle Felker, Princeton University

The numerical solution of radiation hydrodynamics in an arbitrarily curved spacetime is a computationally demanding challenge with applications such as black hole accretion disk formation and compact object mergers. We formulate the time-dependent radiation transfer equation using the 3+1 ADM decomposition of spacetime. Flux-conservative monotonic upwinding is used to advect the specific intensity along spatial and angular meshes. Results for cylindrical, spherical, and toy curved coordinates on Minkowski space are presented and compared to standard flat-space Cartesian solvers.

Abstract Author(s): Kyle Gerard Felker, James Stone, Frans Pretorius