Tree-based solvers for adaptive mesh refinement code FLASH -- II: radiation transport module TreeRay

Abstract

The treatment of radiative transfer with multiple radiation sources is a critical challenge in simulations of star formation and the interstellar medium. In this paper we present the novel TreeRay method for solving general radiative transfer problems, based on reverse ray tracing combined with tree-based accelerated integration. We implement TreeRay in the adaptive mesh refinement code FLASH, as a module of the tree solver developed by Wünsch et al. However, the method itself is independent of the host code and can be implemented in any grid based or particle based hydrodynamics code. A key advantage of TreeRay is that its computational cost is independent of the number of sources, making it suitable for simulations with many point sources (e.g. massive star clusters) as well as simulations where diffuse emission is important. A very efficient communication and tree-walk strategy enables TreeRay to achieve almost ideal parallel scalings. TreeRay can easily be extended with sub-modules to treat radiative transfer at different wavelengths and to implement related physical processes. Here, we focus on ionising (EUV) radiation and use the On-the-Spot approximation to test the method and its parameters. The ability to set the tree solver time step independently enables the speedy calculation of radiative transfer in a multi-phase interstellar medium, where the hydrodynamic time step is typically limited by the sound speed of the hot gas produced in stellar wind bubbles or supernova remnants. We show that complicated simulations of star clusters with feedback from multiple massive stars become feasible with TreeRay.

26 pages, 21 figures, accepted by MNRAS

References in corpus (13)

• Modeling Collapse and Accretion in Turbulent Gas Clouds: Implementation and Comparison of Sink Particles in AMR and SPH
• Radiation-Hydrodynamic Simulations of Collapse and Fragmentation in Massive Protostellar Cores
• Cosmological Radiative Transfer Codes Comparison Project I: The Static Density Field Tests
• The SILCC project --- IV. Impact of dissociating and ionising radiation on the interstellar medium and Halpha emission as a tracer of the star formation rate
• iVINE - Ionization in the parallel tree/SPH code VINE: First results on the observed age-spread around O-stars
• Smoothed Particle Hydrodynamics simulations of expanding HII regions. I. Numerical methods and tests
• Tree-based solvers for adaptive mesh refinement code FLASH - I: gravity and optical depths
• Effects of Magnetic Fields on Photoionised Pillars and Globules
• Modeling UV Radiation Feedback from Massive Stars: I. Implementation of Adaptive Ray Tracing Method and Tests
• Hybrid Adaptive Ray-Moment Method (HARM${}^2$): A Highly Parallel Method for Radiation Hydrodynamics on Adaptive Grids
• Numerical Simulations of Turbulent Molecular Clouds Regulated by Radiation Feedback Forces II: Radiation-Gas Interactions and Outflows
• A general hybrid radiation transport scheme for star formation simulations on an adaptive grid
• TREVR: A general $Nlo{g}^{2}N$ radiative transfer algorithm