Semi-implicit scheme for treating radiation under M1 closure in general relativistic conservative fluid dynamics codes

Abstract

A numerical scheme is described for including radiation in multi-dimensional general-relativistic conservative fluid dynamics codes. In this method, a covariant form of the M1 closure scheme is used to close the radiation moments, and the radiative source terms are treated semi-implicitly in order to handle both optically thin and optically thick regimes. The scheme has been implemented in a conservative general relativistic radiation hydrodynamics code KORAL. The robustness of the code is demonstrated on a number of test problems, including radiative relativistic shock tubes, static radiation pressure supported atmosphere, shadows, beams of light in curved spacetime, and radiative Bondi accretion. The advantages of M1 closure relative to other approaches such as Eddington closure and flux-limited diffusion are discussed, and its limitations are also highlighted.

19 pages, 15 figures, as published in MNRAS

References in corpus (12)

• X-ray Properties of Black-Hole Binaries
• A luminous quasar at a redshift of z = 7.085
• Modelling the behaviour of accretion flows in X-ray binaries
• ECHO: an Eulerian Conservative High Order scheme for general relativistic magnetohydrodynamics and magnetodynamics
• Broad iron K-alpha emission lines as a diagnostic of black hole spin
• Three-Dimensional Simulations of Magnetized Thin Accretion Disks around Black Holes: Stress in the Plunging Region
• Slim disks around Kerr black holes revisited
• Equations and Algorithms for Mixed Frame Flux-Limited Diffusion Radiation Hydrodynamics
• Relativistic Radiation Magnetohydrodynamics in Dynamical Spacetimes: Numerical Methods and Tests
• Light Curves from an MHD Simulation of a Black Hole Accretion Disk
• Thermodynamics of an Accretion Disk Annulus with Comparable Radiation and Gas Pressure
• Surface Structure in an Accretion Disk Annulus with Comparable Radiation and Gas Pressure