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.

arXiv comment25 pages, 19 figures, accepted by MNRAS

References in corpus (34)

• HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
• Dust Grain Size Distributions and Extinction in the Milky Way, LMC, and SMC
• Radiation Pressure Supported Starburst Disks and AGN Fueling
• Rosseland and Planck Mean Opacities for Protoplanetary Discs
• Self-Regulated Star Formation in Galaxies via Momentum Input from Massive Stars
• Modeling Collapse and Accretion in Turbulent Gas Clouds: Implementation and Comparison of Sink Particles in AMR and SPH
• The SILCC (SImulating the LifeCycle of molecular Clouds) project: I. Chemical evolution of the supernova-driven ISM
• On the Role of Massive Stars in the Support and Destruction of Giant Molecular Clouds
• Simulating the formation of molecular clouds. I. Slow formation by gravitational collapse from static initial conditions
• Modeling UV Radiation Feedback from Massive Stars: II. Dispersal of Star-Forming Giant Molecular Clouds by Photoionization and Radiation Pressure
• Mocassin: A fully three-dimensional Monte Carlo photoionization code
• On Radiation Pressure in Static, Dusty HII Regions
• A scheme for radiation pressure and photon diffusion with the M1 closure in RAMSES-RT
• Direct Numerical Simulation of Radiation Pressure-Driven Turbulence and Winds in Star Clusters and Galactic Disks
• Equations and Algorithms for Mixed Frame Flux-Limited Diffusion Radiation Hydrodynamics
• Adaptive Ray Tracing for Radiative Transfer around Point Sources
• First hydrodynamics simulations of radiation forces and photoionization feedback in massive star formation
• AREPO-RT: Radiation hydrodynamics on a moving mesh
• Filamentary fragmentation in a turbulent medium
• The relative impact of photoionizing radiation and stellar winds on different environments
• StarBench: The D-type expansion of an HII region
• Tree-based solvers for adaptive mesh refinement code FLASH - I: gravity and optical depths
• Ionizing Radiation in Smoothed Particle Hydrodynamics
• Fervent: Chemistry-coupled, ionising and non-ionising radiative feedback in magnetohydrodynamical simulations
• SILCC-ZOOM: The early impact of ionizing radiation on forming molecular clouds
• Massive Star Formation via the Collapse of Subvirial and Virialized Turbulent Massive Cores
• 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
• Simulating protostellar evolution and radiative feedback in the cluster environment
• Makemake + Sedna: A Continuum Radiation Transport and Photoionization Framework for Astrophysical Newtonian Fluid Dynamics
• Urchin: A Reverse Ray Tracer for Astrophysical Applications
• 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 (4)

• 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
• Linking simulations and observations of star-forming regions I. Masses and temperatures of ALMAGAL-like cores