Numerical Treatment of Anisotropic Radiation Field Coupling with the Relativistic Resistive Magnetofluids

Abstract

We develop a numerical scheme for solving a fully special relativistic resistive radiation magnetohydrodynamics. Our code guarantees conservations of total mass, momentum and energy. Radiation energy density and radiation flux are consistently updated using the M-1 closure method, which can resolve an anisotropic radiation fields in contrast to the Eddington approximation as well as the flux-limited diffusion approximation. For the resistive part, we adopt a simple form of the Ohm's law. The advection terms are explicitly solved with an approximate Riemann solver, mainly HLL scheme, and HLLC and HLLD schemes for some tests. The source terms, which describe the gas-radiation interaction and the magnetic energy dissipation, are implicitly integrated, relaxing the Courant-Friedrichs-Lewy condition even in optically thick regime or a large magnetic Reynolds number regime. Although we need to invert $4\times 4$ (for gas-radiation interaction) and $3\times 3$ (for magnetic energy dissipation) matrices at each grid point for implicit integration, they are obtained analytically without preventing massive parallel computing. We show that our code gives reasonable outcomes in numerical tests for ideal magnetohydrodynamics, propagating radiation, and radiation hydrodynamics. We also applied our resistive code to the relativistic Petschek type magnetic reconnection, revealing the reduction of the reconnection rate via the radiation drag.

16 pages, 1 table, 13 Figures, accepted for publication in ApJ

References in corpus (12)

• Stability of Relativistic Jets from Rotating, Accreting Black Holes via Fully Three-Dimensional Magnetohydrodynamic Simulations
• MHD simulations of the magnetorotational instability in a shearing box with zero net flux. II. The effect of transport coefficients
• Impact of dimensionless numbers on the efficiency of MRI-induced turbulent transport
• Equation of State in Relativistic Magnetohydrodynamics: variable versus constant adiabatic index
• Semi-implicit scheme for treating radiation under M1 closure in general relativistic conservative fluid dynamics codes
• An HLLC Solver for Relativistic Flows -- II. Magnetohydrodynamics
• Multi-dimensional Numerical Scheme for Resistive Relativistic MHD
• Tearing instability in relativistic magnetically dominated plasmas
• Two-dimensional MHD simulations of relativistic magnetic reconnection
• Resistive Magnetohydrodynamic Simulations of Relativistic Magnetic Reconnection
• Saturation of the MRI in Strongly Radiation Dominated Accretion Disks
• Explicit-Implicit Scheme for Relativistic Radiation Hydrodynamics

Cited by in corpus (26)

• Numerical simulations of super-critical black hole accretion flows in general relativity
• Three-Dimensional General Relativistic Radiation Magnetohydrodynamical Simulation of Super-Eddington Accretion, using a new code HARMRAD with M1 Closure
• Global simulations of axisymmetric radiative black hole accretion disks in general relativity with a sub-grid magnetic dynamo
• Formation of Overheated Regions and Truncated Disks around Black Holes; Three-dimensional General Relativistic Radiation-magnetohydrodynamics Simulations
• Radiation Hydrodynamic Simulations of Line-Driven Disk Winds for Ultra Fast Outflows
• Radiation hydrodynamics simulations of massive star cluster formation in giant molecular clouds
• Three-dimensional structure of clumpy outflow from supercritical accretion flow onto black holes
• Hyper-Eddington accretion flows onto black holes accompanied by powerful outflows
• Rapid growth of black holes accompanied with hot or warm outflows exposed to anisotropic super-Eddington radiation
• Magnetohydrodynamic Simulation Code CANS+: Assessments and Applications
• Radiative Magnetic Reconnection in Astrophysics
• A radiative transfer module for relativistic magnetohydrodynamics in the PLUTO code
• The thermal-radiative wind in low mass X-ray binary H 1743-322; Radiation hydrodynamic simulations
• Super-Eddington growth of black holes in the early Universe: effects of disk radiation spectra
• A two-moment radiation hydrodynamics scheme applicable to simulations of planet formation in circumstellar disks
• 2D radiaition-hydrodynamic simulations of supernova shock breakout in bipolar explosions of a blue supergiant progenitor
• General-relativistic radiation transport scheme in $\text{Gmunu}$ I: Implementation of two-moment based multi-frequency radiative transfer and code tests
• A Monte-Carlo based relativistic radiation hydrodynamics code with a higher-order scheme
• CAFE-R a code that solves the special relativistic radiation hydrodynamics equations
• Radiation Magnetohydrodynamic Simulations of Sub-Eddington accretion Flows in AGN: Origin of Soft X-ray Excess and Rapid Time Variabilities
• Oscillations of radiation pressure supported tori near black holes
• Time-dependent Multi-group Multidimensional Relativistic Radiative Transfer Code Based On Spherical Harmonic Discrete Ordinate Method
• On the application of Jacobian-free Riemann solvers for relativistic radiation magnetohydrodynamics under M1 closure
• Approaching the black hole by numerical simulations
• Model for an optically thick torus in local thermodynamic equilibrium around a black hole
• Self-similar structure of resistive ADAFs with outflow and large-scale magnetic field