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.

arXiv comment26 pages, 21 figures, accepted by MNRAS

References in corpus (49)

• HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
• Towards a complete accounting of energy and momentum from stellar feedback in galaxy formation simulations
• 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
• The SILCC (SImulating the LifeCycle of molecular Clouds) project: I. Chemical evolution of the supernova-driven ISM
• RAMSES-RT: Radiation hydrodynamics in the cosmological context
• Massive stars in massive clusters II: Disruption of bound clusters by photoionization
• Three-phase Interstellar medium in Galaxies Resolving Evolution with Star formation and Supernova feedback (TIGRESS): Algorithms, Fiducial model, and Convergence
• The SILCC (SImulating the LifeCycle of molecular Clouds) project - II. Dynamical evolution of the supernova-driven ISM and the launching of outflows
• Dispersal of molecular clouds by ionising radiation
• Cosmological Radiative Transfer Codes Comparison Project I: The Static Density Field Tests
• The SILCC project: III. Regulation of star formation and outflows by stellar winds and supernovae
• Enzo+Moray: Radiation Hydrodynamics Adaptive Mesh Refinement Simulations with Adaptive Ray Tracing
• Modelling CO formation in the turbulent interstellar medium
• TreeCol: a novel approach to estimating column densities in astrophysical simulations
• Radiation driven implosion and triggered star formation
• Dynamical HII Region Evolution in Turbulent Molecular Clouds
• Cosmological Radiative Transfer Comparison Project II: The Radiation-Hydrodynamic Tests
• Detailed Numerical Simulations on the Formation of Pillars around HII-regions
• AREPO-RT: Radiation hydrodynamics on a moving mesh
• Fast and accurate frequency-dependent radiation transport for hydrodynamics simulations in massive star formation
• Radiation Driven Implosion of Molecular Cloud Cores
• A Two-moment Radiation Hydrodynamics Module in Athena Using a Time-explicit Godunov Method
• 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
• The relative impact of photoionizing radiation and stellar winds on different environments
• Hybrid Characteristics: 3D radiative transfer for parallel adaptive mesh refinement hydrodynamics
• iVINE - Ionization in the parallel tree/SPH code VINE: First results on the observed age-spread around O-stars
• StarBench: The D-type expansion of an HII region
• Clumps and triggered star formation in ionised molecular clouds
• 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
• Dynamical Models for the Formation of Elephant Trunks in H II Regions
• The Pillars of Creation revisited with MUSE: gas kinematics and high-mass stellar feedback traced by optical spectroscopy
• An Investigation on the Morphological Evolution of Bright-Rimmed Clouds
• Effects of Magnetic Fields on Photoionised Pillars and Globules
• The Elephant Trunk Nebula and the Trumpler 37 cluster: Contribution of triggered star formation to the total population of an HII region
• Fervent: Chemistry-coupled, ionising and non-ionising radiative feedback in magnetohydrodynamical simulations
• Modeling UV Radiation Feedback from Massive Stars: I. Implementation of Adaptive Ray Tracing Method and Tests
• SILCC-ZOOM: The early impact of ionizing radiation on forming molecular clouds
• Globules and pillars seen in the [CII] 158 micron line with SOFIA
• Hybrid Adaptive Ray-Moment Method (HARM${}^2$): A Highly Parallel Method for Radiation Hydrodynamics on Adaptive Grids
• Evolution of Blister-Type HII Regions in a Magnetized Medium
• SimpleX2: radiative transfer on an unstructured, dynamic grid
• 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
• Non-Equilibrium Chemistry and Destruction of CO by X-ray Flares
• Urchin: A Reverse Ray Tracer for Astrophysical Applications
• Radiative feedback and cosmic molecular gas: numerical method
• TREVR: A general $Nlo{g}^{2}N$ radiative transfer algorithm

Cited by in corpus (16)

• Introducing TIGRESS-NCR: I. Co-Regulation of the Multiphase Interstellar Medium and Star Formation Rates
• SILCC VII -- Gas kinematics and multiphase outflows of the simulated ISM at high gas surface densities
• Photochemistry and Heating/Cooling of the Multiphase Interstellar Medium with UV Radiative Transfer for Magnetohydrodynamic Simulations
• Why are thermally- and cosmic ray-driven galactic winds fundamentally different?
• Radiation transport methods in star formation simulations
• Catastrophic Cooling in Superwinds. III. Non-equilibrium Photoionization
• Tree-based solvers for adaptive mesh refinement code FLASH -- III: a novel scheme for radiation pressure on dust and gas and radiative transfer from diffuse sources
• The impact of cosmic ray heating on the cooling of the low-metallicity interstellar medium
• Multidimensional half-moment multigroup radiative transfer. Improving moment-based thermal models of circumstellar disks
• A Discontinuous Galerkin Solver in the FLASH Multi-Physics Framework
• Tree-based solvers for adaptive mesh refinement code FLASH -- IV: An X-ray radiation scheme to couple discrete and diffuse X-ray emission sources to the thermochemistry of the interstellar medium
• The origin and evolution of the [CII] deficit in HII regions and star-forming molecular clouds
• A multi-ion non-equilibrium solver for ionised astrophysical plasmas with arbitrary elemental abundances
• Predicting interstellar radiation fields from chemical evolution models
• SILCC -- IX. The multi-phase interstellar medium at low metallicity
• Linking simulations and observations of star-forming regions I. Masses and temperatures of ALMAGAL-like cores