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

Abstract

Radiation is an important contributor to the energetics of the interstellar medium, yet its transport is difficult to solve numerically. We present a novel approach towards solving radiative transfer of diffuse sources via backwards ray tracing. Here we focus on the radiative transfer of infrared radiation and the radiation pressure on dust. The new module, \textsc{TreeRay/RadPressure}, is an extension to the novel radiative transfer method \textsc{TreeRay} implemented in the grid-based MHD code {\sc Flash}. In \textsc{TreeRay/RadPressure}, every cell and every star particle is a source of infrared radiation. We also describe how gas, dust and radiation are coupled via a chemical network. This allows us to compute the local dust temperature in thermal equilibrium, leading to a significantly improvement over the classical grey approximation. In several tests, we demonstrate that the scheme produces the correct radiative intensities as well as the correct momentum input by radiation pressure. Subsequently, we apply our new scheme to model massive star formation from a collapsing, turbulent core of 150 ${\mathrm{M}}_{\odot }$. We include the effects of both, ionizing and infrared radiation on the dynamics of the core. We find that the newborn massive star prevents fragmentation in its proximity due to radiative heating. Over time, dust and radiation temperature equalize, while the gas temperature can be either warmer due to shock heating or colder due to insufficient dust-gas coupling. Compared to gravity, the effects of radiation pressure are insignificant for the stellar mass on the simulated time scale in this work.

25 pages, 19 figures, accepted by MNRAS

References in corpus (11)

• Modeling Collapse and Accretion in Turbulent Gas Clouds: Implementation and Comparison of Sink Particles in AMR and SPH
• Simulating the formation of molecular clouds. I. Slow formation by gravitational collapse from static initial conditions
• A scheme for radiation pressure and photon diffusion with the M1 closure in RAMSES-RT
• Equations and Algorithms for Mixed Frame Flux-Limited Diffusion Radiation Hydrodynamics
• Tree-based solvers for adaptive mesh refinement code FLASH - I: gravity and optical depths
• How fast do young star clusters expel their natal gas?: Estimating the upper limit of the gas expulsion time-scale
• Hybrid Adaptive Ray-Moment Method (HARM${}^2$): A Highly Parallel Method for Radiation Hydrodynamics on Adaptive Grids
• Makemake + Sedna: A Continuum Radiation Transport and Photoionization Framework for Astrophysical Newtonian Fluid Dynamics
• Tree-based solvers for adaptive mesh refinement code FLASH -- II: radiation transport module TreeRay
• VETTAM: A scheme for radiation hydrodynamics with adaptive mesh refinement using the variable Eddington tensor method
• TREVR: A general $Nlo{g}^{2}N$ radiative transfer algorithm

Cited by in corpus (3)

• Radiation transport methods in star formation simulations
• The origin and evolution of the [CII] deficit in HII regions and star-forming molecular clouds
• 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