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 (59)

• Observational Evidence of AGN Feedback
• X-ray Properties of Black-Hole Binaries
• Advection-Dominated Accretion: A Self-Similar Solution
• Advection-Dominated Accretion: Underfed Black Holes and Neutron Stars
• A luminous quasar at a redshift of z = 7.085
• Advection-Dominated Accretion and the Spectral States of Black Hole X-Ray Binaries: Application to Nova Muscae 1991
• Modelling the behaviour of accretion flows in X-ray binaries
• General Relativistic Magnetohydrodynamic Simulations of Magnetically Choked Accretion Flows around Black Holes
• HARM: A Numerical Scheme for General Relativistic Magnetohydrodynamics
• Nonthermal Electrons in Radiatively Inefficient Accretion Flow Models of Sagittarius A*
• Thermal Equilibria of Accretion Disks
• The Canada-France High-z Quasar Survey: nine new quasars and the luminosity function at redshift 6
• GRMHD Simulations of Magnetized Advection Dominated Accretion on a Non-Spinning Black Hole: Role of Outflows
• A Survey of z>5.7 Quasars in the Sloan Digital Sky Survey IV: Discovery of Seven Additional Quasars
• Intermediate-Mass Black Holes
• GRS 1915+105 and the Disc-Jet Coupling in Accreting Black Hole Systems
• Global Structure of Three Distinct Accretion Flows and Outflows around Black Holes through Two-Dimensional Radiation-Magnetohydrodynamic Simulations
• Magnetically Driven Accretion Flows in the Kerr Metric I: Models and Overall Structure
• Primitive Variable Solvers for Conservative General Relativistic Magnetohydrodynamics
• ECHO: an Eulerian Conservative High Order scheme for general relativistic magnetohydrodynamics and magnetodynamics
• Simulations of Magnetized Disks Around Black Holes: Effects of Black Hole Spin, Disk Thickness, and Magnetic Field Geometry
• Spectral Transitions in Cyg X-1 and Other Black Hole X-Ray Binaries
• Global Radiation-Magnetohydrodynamic Simulations of Black Hole Accretion Flow and Outflow: Unified Model of Three States
• Measuring the Spins of Accreting Black Holes
• Radiation-Dominated Disks Are Thermally Stable
• Broad iron K-alpha emission lines as a diagnostic of black hole spin
• Prograde and Retrograde Black Holes: Whose Jet is More Powerful?
• GYOTO: a new general relativistic ray-tracing code
• Three-Dimensional Simulations of Magnetized Thin Accretion Disks around Black Holes: Stress in the Plunging Region
• Slim Disk Model for Ultra-Luminous X-Ray Sources
• Sagittarius A* Accretion Flow and Black Hole Parameters from General Relativistic Dynamical and Polarized Radiative Modeling
• Slim disks around Kerr black holes revisited
• Equations and Algorithms for Mixed Frame Flux-Limited Diffusion Radiation Hydrodynamics
• Iron Emission in the z=6.4 Quasar SDSS J114816.64+525150.3
• Cosmos++: Relativistic Magnetohydrodynamics on Unstructured Grids with Local Adaptive Refinement
• The Exact Solution of the Riemann Problem in Relativistic MHD
• A Godunov Method for Multidimensional Radiation Magnetohydrodynamics based on a variable Eddington tensor
• A Radiation Transfer Solver for Athena using Short Characteristics
• General relativistic radiation hydrodynamics of accretion flows. I: Bondi-Hoyle accretion
• Beyond Flux-Limited Diffusion: Parallel Algorithms for Multidimensional Radiation Hydrodynamics
• Are quasars accreting at super-Eddington rates?
• Local Three-Dimensional Simulations of Magneto-Rotational Instability in Radiation-Dominated Accretion Disks
• Relativistic Radiation Magnetohydrodynamics in Dynamical Spacetimes: Numerical Methods and Tests
• Light Curves from an MHD Simulation of a Black Hole Accretion Disk
• General Relativistic Magnetohydrodynamic Simulations of the Hard State as a Magnetically-Dominated Accretion Flow
• GRMHD simulations of accretion onto Sgr A*: How important are radiative losses?
• Imaging of SDSS z>6 quasar fields: gravitational lensing, companion galaxies and the host dark matter halos
• Turbulent Stresses in Local Simulations of Radiation-Dominated Accretion Disks, and the Possibility of the LIghtman-Eardley Instability
• Radiative efficiency and thermal spectrum of accretion onto Schwarzschild black holes
• Numerical Simulations of Optically Thick Accretion onto a Black Hole - I. Spherical Case
• Dissipation and Vertical Energy Transport in Radiation-Dominated Accretion Disks
• Testing slim-disk models on the thermal spectra of LMC X-3
• Two-dimensional radiation-hydrodynamic model for limit-cycle oscillations of luminous accretion disks
• General relativistic radiation hydrodynamics of accretion flows: II. Treating stiff source terms and exploring physical limitations
• Thermodynamics of an Accretion Disk Annulus with Comparable Radiation and Gas Pressure
• Spectral Energy Distribution in Supercritical Disk Accretion Flows through Photon-trapping Effects
• Surface Structure in an Accretion Disk Annulus with Comparable Radiation and Gas Pressure
• Radiation Magnetohydrodynamics for Black Hole-Torus System in Full General Relativity: A Step toward Physical Simulation
• The formation of the coronal flow/ADAF

Cited by in corpus (142)

• First M87 Event Horizon Telescope Results. V. Physical Origin of the Asymmetric Ring
• First M87 Event Horizon Telescope Results. VIII. Magnetic Field Structure near The Event Horizon
• Comprehensive nucleosynthesis analysis for ejecta of compact binary mergers
• First Sagittarius A* Event Horizon Telescope Results. V. Testing Astrophysical Models of the Galactic Center Black Hole
• The Event Horizon General Relativistic Magnetohydrodynamic Code Comparison Project
• 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
• The Black Hole Accretion Code
• Magnetic Reconnection and Hot Spot Formation in Black Hole Accretion Disks
• Three-dimensional simulations of super-critical black hole accretion disks --- luminosities, photon trapping and variability
• Observing the Inner Shadow of a Black Hole: A Direct View of the Event Horizon
• Energy, momentum and mass outflows and feedback from thick accretion discs around rotating black holes
• Global simulations of axisymmetric radiative black hole accretion disks in general relativity with a sub-grid magnetic dynamo
• Jets in Magnetically Arrested Hot Accretion Flows: Geometry, Power and Black Hole Spindown
• Two-temperature, Magnetically Arrested Disc simulations of the jet from the supermassive black hole in M87
• Powerful radiative jets in super-critical accretion disks around non-spinning black holes
• A new multidimensional, energy-dependent two-moment transport code for neutrino-hydrodynamics
• Magnetohydrodynamic Simulations of Active Galactic Nucleus Disks and Jets
• The role of electron heating physics in images and variability of the Galactic Centre black hole Sagittarius A*
• Radiative, two-temperature simulations of low luminosity black hole accretion flows in general relativity
• The Size, Shape, and Scattering of Sagittarius A* at 86 GHz: First VLBI with ALMA
• Hydrodynamic Simulations of the Inner Accretion Flow of Sagittarius A* Fueled By Stellar Winds
• Magnetohydrodynamical simulations of a tidal disruption in general relativity
• An Algorithm for Radiation Magnetohydrodynamics Based on Solving the Time-dependent Transfer Equation
• Thin accretion disks are stabilized by a strong magnetic field
• H-AMR: A New GPU-accelerated GRMHD Code for Exascale Computing With 3D Adaptive Mesh Refinement and Local Adaptive Time-stepping
• Formation of Magnetically Truncated Accretion Disks in 3D Radiation-Transport Two-Temperature GRMHD Simulations
• Mildly Super-Eddington Accretion Onto Slowly-Spinning Black Holes Explains the X-Ray Weakness of the Little Red Dots
• Spectra of Black Hole Accretion Models of Ultra-Luminous X-ray Sources
• Formation of Overheated Regions and Truncated Disks around Black Holes; Three-dimensional General Relativistic Radiation-magnetohydrodynamics Simulations
• General relativistic resistive magnetohydrodynamics with robust primitive variable recovery for accretion disk simulations
• GRRMHD Simulations of Tidal Disruption Event Accretion Disks around Supermassive Black Holes: Jet Formation, Spectra, and Detectability
• Can we make sense of dissipation without causality?
• Neutrino transport in general relativistic neutron star merger simulations
• Energy flows in thick accretion disks and their consequences for black hole feedback
• Numerical Simulations of Optically Thick Accretion onto a Black Hole - II. Rotating Flow
• Flux eruption events drive angular momentum transport in magnetically arrested accretion flows
• An Implicit Finite Volume Scheme to Solve the Time Dependent Radiation Transport Equation Based on Discrete Ordinates
• Analytic Closures for M1 Neutrino Transport
• Puffy accretion disks: sub-Eddington, optically thick, and stable
• Kinetic and radiative power from optically thin accretion flows
• Evolving Nonthermal Electrons in Simulations of Black Hole Accretion
• X-Ray Polarimetry with the Polarization Spectroscopic Telescope Array (PolSTAR)
• General Overview of Black Hole Accretion Theory
• Relativistic Viscous Radiation Hydrodynamic Simulations of Geometrically Thin Disks: I. Thermal and Other Instabilities
• Photon-conserving Comptonization in simulations of accretion disks around black holes
• Astrophysics of Super-massive Black Hole Mergers
• Iharm3D: Vectorized General Relativistic Magnetohydrodynamics
• HEROIC: 3D General Relativistic Radiative Postprocessor with Comptonization for Black Hole Accretion Discs
• The Physics of Accretion Discs, Winds And Jets in Tidal Disruption Events
• Three-dimensional, global, radiative GRMHD simulations of a thermally unstable disc
• Black Hole Polarimetry I: A Signature of Electromagnetic Energy Extraction
• Radiation Transport Two-Temperature GRMHD Simulations of Warped Accretion Disks
• Numerical Treatment of Anisotropic Radiation Field Coupling with the Relativistic Resistive Magnetofluids
• GRRMHD Simulations of MAD Accretion Disks Declining from Super-Eddington to Sub-Eddington Accretion Rates
• Applying the Gibbs stability criterion to relativistic hydrodynamics
• Multi-frequency Black Hole Imaging for the Next-Generation Event Horizon Telescope
• Supercritical Accretion onto a Non-Magnetized Neutron Star: Why is it Feasible?
• Beyond moments: relativistic Lattice-Boltzmann methods for radiative transport in computational astrophysics
• A radiative transfer module for relativistic magnetohydrodynamics in the PLUTO code
• Two-moment scheme for general-relativistic radiation hydrodynamics: a systematic description and new applications
• Observational Signatures of Frame Dragging in Strong Gravity
• GRMHD Simulations and Modeling for Jet Formation and Acceleration Region in AGNs
• Observational properties of puffy disks: radiative GRMHD spectra of mildly sub-Eddington accretion
• An Extension of the Athena++ Code Framework for Radiation-Magnetohydrodynamics in General Relativity Using a Finite-Solid-Angle Discretization
• Theoretical modeling of Comptonized X-ray spectra of super-Eddington accretion flow: origin of hard excess in Ultraluminous X-Ray Sources
• Testing General Relativity with the Black Hole Shadow Size and Asymmetry of Sagittarius A*: Limitations from Interstellar Scattering
• Multi-Frequency General Relativistic Radiation-Hydrodynamics with ${\mathbf{M}}_1$ Closure
• Magnetized discs and photon rings around Yukawa-like black holes
• Rotation and toroidal magnetic field effects on the stability of two-component jets
• Multifluid Modelling of Relativistic Radiation Hydrodynamics
• Location of the bow shock ahead of cloud G2 at the Galactic Center
• A two-moment radiation hydrodynamics scheme applicable to simulations of planet formation in circumstellar disks
• Shock wave structure in astrophysical flows with an account of photon transfer
• Radiative GRMHD simulations of accretion and outflow in non-magnetized neutron stars and ultraluminous X-ray sources
• Beamed emission from a neutron-star ULX in a GRRMHD simulation
• GRay2: A General Purpose Geodesic Integrator for Kerr Spacetimes
• Global Simulations of Tidal Disruption Event Disk Formation via Stream Injection in GRRMHD
• 2D radiaition-hydrodynamic simulations of supernova shock breakout in bipolar explosions of a blue supergiant progenitor
• Rank-3 Moment Closures in General Relativistic Neutrino Transport
• Accretion Flow Morphology in Numerical Simulations of Black Holes from the ngEHT Model Library: The Impact of Radiation Physics
• General-relativistic radiation transport scheme in $\text{Gmunu}$ I: Implementation of two-moment based multi-frequency radiative transfer and code tests
• Structured, relativistic jets driven by radiation
• MOCMC: Method of Characteristics Moment Closure, a Numerical Method for Covariant Radiation Magnetohydrodynamics
• Levitating atmospheres of Eddington-luminosity neutron stars
• Numerical simulations of the Cosmic Battery in accretion flows around astrophysical black holes
• Blazar Variability From Turbulence in Jets Launched by Magnetically Arrested Accretion Flows
• Covariant Radiative Transfer for Black Hole Spacetimes
• Simulating the effect of the Sgr A* accretion flow on the appearance of G2 after pericenter
• The Guided Moments formalism: a new efficient full-neutrino treatment for astrophysical simulations
• Energy-dependent and energy-integrated two-moment general-relativistic neutrino transport simulations of hypermassive neutron star
• Component of energy flow from supercritical accretion disks around rotating stellar mass black holes
• Radiation transport methods in star formation simulations
• Survey of Radiative, Two-Temperature Magnetically Arrested Simulations of the Black Hole M87* I: Turbulent Electron Heating
• General Relativistic Implicit Monte Carlo Radiation-Hydrodynamics
• Breathing Oscillations in a Global Simulation of a Thin Accretion Disk
• Optically thin outbursts of rotating neutron stars can not be spherical
• Relativistic, Viscous, Radiation Hydrodynamic Simulations of Geometrically Thin Disks. II. Disk Variability
• Performance-Portable Binary Neutron Star Mergers with AthenaK
• Modeling of Thermal Emission from ULX Pulsar Swift J0243.6+6124 with General Relativistic Radiation MHD simulations
• The general relativistic equations of radiation hydrodynamics in the viscous limit
• Hybrid GRMHD and Force-Free Simulations of Black Hole Accretion
• Jet Signatures in the Spectra of Accreting Black Holes
• Optical and X-ray luminosity of expanding nebulae around ultraluminous X-ray sources
• Choked accretion onto a Kerr black hole
• CAFE-R a code that solves the special relativistic radiation hydrodynamics equations
• Bulk Comptonization by Turbulence in Accretion Disks
• Simulations of recoiling black holes: adaptive mesh refinement and radiative transfer
• Relativistic radiation hydrodynamics in a reference-metric formulation
• Jets from SANE Super-Eddington Accretion Disks: Morphology, Spectra, and Their Potential as Targets for ngEHT
• cuHARM : a new GPU accelerated GR-MHD code and its application to ADAF disks
• Global high-order numerical schemes for the time evolution of the general relativistic radiation magneto-hydrodynamics equations
• Time-dependent Multi-group Multidimensional Relativistic Radiative Transfer Code Based On Spherical Harmonic Discrete Ordinate Method
• Oscillations of radiation pressure supported tori near black holes
• A High-Frequency Doppler Feature in the Power Spectra of Simulated GRMHD Black Hole Accretion Disks
• Modeling long GRBs using a single shock with relativistic radiation hydrodynamics
• On the application of Jacobian-free Riemann solvers for relativistic radiation magnetohydrodynamics under M1 closure
• Infinite Order Hydrodynamics: An Analytical Example
• GR-RMHD Simulations of Super-Eddington Accretion Flows onto a Neutron Star with Dipole and Quadrupole Magnetic Fields
• ULX collimation by outflows in moderately magnetized neutron stars
• Spectra of Puffy Accretion Discs: the kynbb Fit
• Deep learning inference with the Event Horizon Telescope II. The Zingularity framework for Bayesian artificial neural networks
• Evolution of jets driven by relativistic radiation hydrodynamics as Long and Low Luminosity GRBs
• Multidimensional half-moment multigroup radiative transfer. Improving moment-based thermal models of circumstellar disks
• Approaching the black hole by numerical simulations
• Dispersion relations of relativistic radiation hydrodynamics
• 3D Photon Conserving Code for Time-dependent General Relativistic Radiative Transfer : CARTOON
• Future Simulations of Tidal Disruption Events
• Relativistic heat conduction in the large-flux regime
• Model for an optically thick torus in local thermodynamic equilibrium around a black hole
• Powerful Jets from Radiatively Efficient Disks, a Decades-Old Unresolved Problem in High Energy Astrophysics
• Energy distribution and substructure formation in astrophysical MHD simulations
• Current Status of Simulations
• Multigroup Radiation Diffusion on a Moving Mesh: Implementation in RICH and Application to Tidal Disruption Events
• Radiative GRMHD simulations of accretion disks
• The radiation hydrodynamics of relativistic shear flows
• Non-Newtonian corrections to radiative viscosity: Israel-Stewart theory as a viscosity limiter
• PATOKA: Simulating Electromagnetic Observables of Black Hole Accretion
• A multifluid perspective on multimessenger modelling
• Global Magnetohydrodynamic Simulations of Monster Shocks in Neutron Star Magnetospheres
• Interactions of Type I X-ray Bursts with Thin Accretion Discs
• GRRMHD Simulations of State Transitions in Non-Jetted Tidal Disruption Events