Multigroup Radiation Diffusion on a Moving Mesh: Implementation in RICH and Application to Tidal Disruption Events
arXiv:2601.05120 · doi:10.3847/1538-4365/ae4874
Abstract
Radiation-hydrodynamics (RHD) determines the bulk evolution and observable emission in a wide variety of high-energy astrophysical phenomena. Due to their complexity, RHD problems must usually be studied through numerical simulation. We have extended the publicly available RICH code, which previously solved the equations of RHD in the limit of grey flux-limited diffusion (FLD), to operate with a multigroup FLD solver. RICH is a semi-Lagrangian code that solves the equations of RHD on an unstructured moving mesh, and is the first multigroup RHD moving mesh code, making it uniquely applicable to problems with extreme dynamic range and dynamically important radiation forces. We validate our multigroup module against multiple analytic benchmarks, including a novel test of the RHD Doppler term. The computational efficiency of the code is aided by a novel scheme to accelerate convergence in optically thick cells by limiting the absorption coefficients. Finally, we apply multigroup \textsc{rich} in a pilot three dimensional study of a stellar tidal disruption event (TDE), using a intermediate-mass black hole. Our simulations self-consistently produce a bright early-time X-ray flash prior to peak optical/UV light, in qualitative agreement with post-processing of (grey) RICH simulations of supermassive black hole TDEs, as well as X-ray observations of the TDE AT 2022dsb.
References in corpus (43)
- HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
- Hot Accretion Flows Around Black Holes
- Radiation-Hydrodynamic Simulations of Collapse and Fragmentation in Massive Protostellar Cores
- A unified model for tidal disruption events
- Moving mesh cosmology: tracing cosmological gas accretion
- General Relativistic Hydrodynamic Simulation of Accretion Flow from a Stellar Tidal Disruption
- PS1-10jh: The Disruption of a Main-Sequence Star of Near-Solar Composition
- Super-Eddington Accretion Disks around Supermassive black Holes
- Optical-Ultraviolet Tidal Disruption Events
- Semi-implicit scheme for treating radiation under M1 closure in general relativistic conservative fluid dynamics codes
- Finite, Intense Accretion Bursts from Tidal Disruption of Stars on Bound Orbits
- STARFORGE: Toward a comprehensive numerical model of star cluster formation and feedback
- Dynamics of Dusty Radiation Pressure Driven Shells and Clouds: Fast Outflows from Galaxies, Star Clusters, Massive Stars, and AGN
- Equations and Algorithms for Mixed Frame Flux-Limited Diffusion Radiation Hydrodynamics
- A radiative transfer scheme for cosmological reionization based on a local Eddington tensor
- A Godunov Method for Multidimensional Radiation Magnetohydrodynamics based on a variable Eddington tensor
- AREPO-RT: Radiation hydrodynamics on a moving mesh
- Beyond Flux-Limited Diffusion: Parallel Algorithms for Multidimensional Radiation Hydrodynamics
- The Prospects of Observing Tidal Disruption Events with the LSST
- Monte Carlo Radiative Transfer
- Continuum-Fitting the X-ray Spectra of Tidal Disruption Events
- Tidal Disruption Disks Formed and Fed by Stream-Stream and Stream-Disk Interactions in Global GRHD Simulations
- First light from tidal disruption events
- On the Missing Energy Puzzle of Tidal Disruption Events
- The Multi-Dimensional Structure of Radiative Shocks: Suppressed Thermal X-rays and Relativistic Ion Acceleration
- Rich: Open Source Hydrodynamic Simulation on a Moving Voronoi Mesh
- Time Dependent Radiation Hydrodynamics on a Moving Mesh
- Solutions of the imploding shock problem in a medium with varying density
- Analytic solutions of the nonlinear radiation diffusion equation with an instantaneous point source in non-homogeneous media
- Piston driven shock waves in non-homogeneous planar media
- A New Discrete Implicit Monte Carlo Scheme for Simulating Radiative Transfer Problems
- Self-similar solutions for the non-equilibrium nonlinear supersonic Marshak wave problem
- Simulating cold shear flows on a moving mesh
- Solutions of the converging and diverging shock problem in a medium with varying density
- Frequency-dependent Discrete Implicit Monte-Carlo Scheme for the Radiative Transfer Equation
- Grid Noise in Moving Mesh Codes: Fixing The Volume Inconsistency Problem
- The Primary Flare Following a Stellar Collision in a Galactic Nucleus
- An accurate reaction-diffusion limit to the spherical-symmetric Boltzmann equation
- Converged simulations of the nozzle shock in tidal disruption events
- Future Simulations of Tidal Disruption Events
- Nonrelativistic grey Sn-transport radiative-shock solutions
- Balancing the load: A Voronoi based scheme for parallel computations
- A novel numerical method for mixed-frame multigroup radiation-hydrodynamics with GPU acceleration implemented in the QUOKKA code