We introduce arepo-rt, a novel radiation hydrodynamic (RHD) solver for the unstructured moving-mesh code arepo. Our method solves the moment-based radiative transfer equations using the M1 closure relation. We achieve second-order convergence by using a slope-limited linear spatial extrapolation and a first-order time prediction step to obtain the values of the primitive variables on both sides of the cell interface. A Harten–Lax–van Leer flux function, suitably modified for moving meshes, is then used to solve the Riemann problem at the interface. The implementation is fully conservative and compatible with the individual time-stepping scheme of arepo. It incorporates atomic hydrogen (H) and helium (He) thermochemistry, which is used to couple the ultraviolet radiation field to the gas. Additionally, infrared (IR) radiation is coupled to the gas under the assumption of local thermodynamic equilibrium between the gas and the dust. We successfully apply our code to a large number of test problems, including applications such as the expansion of HII regions, radiation pressure-driven outflows, and the levitation of optically thick layer of gas by trapped IR radiation. The new implementation is suitable for studying various important astrophysical phenomena, such as the effect of radiative feedback in driving galactic scale outflows, radiation-driven dusty winds in high-redshift quasars, or simulating the reionization history of the Universe in a self-consistent manner.

arepo-rt: radiation hydrodynamics on a moving mesh

Marinacci, Federico;
2019

Abstract

We introduce arepo-rt, a novel radiation hydrodynamic (RHD) solver for the unstructured moving-mesh code arepo. Our method solves the moment-based radiative transfer equations using the M1 closure relation. We achieve second-order convergence by using a slope-limited linear spatial extrapolation and a first-order time prediction step to obtain the values of the primitive variables on both sides of the cell interface. A Harten–Lax–van Leer flux function, suitably modified for moving meshes, is then used to solve the Riemann problem at the interface. The implementation is fully conservative and compatible with the individual time-stepping scheme of arepo. It incorporates atomic hydrogen (H) and helium (He) thermochemistry, which is used to couple the ultraviolet radiation field to the gas. Additionally, infrared (IR) radiation is coupled to the gas under the assumption of local thermodynamic equilibrium between the gas and the dust. We successfully apply our code to a large number of test problems, including applications such as the expansion of HII regions, radiation pressure-driven outflows, and the levitation of optically thick layer of gas by trapped IR radiation. The new implementation is suitable for studying various important astrophysical phenomena, such as the effect of radiative feedback in driving galactic scale outflows, radiation-driven dusty winds in high-redshift quasars, or simulating the reionization history of the Universe in a self-consistent manner.
2019
Kannan, Rahul; Vogelsberger, Mark; Marinacci, Federico; McKinnon, Ryan; Pakmor, Rüdiger; Springel, Volker
File in questo prodotto:
File Dimensione Formato  
stz287_compressed.pdf

accesso aperto

Tipo: Versione (PDF) editoriale
Licenza: Licenza per accesso libero gratuito
Dimensione 2.01 MB
Formato Adobe PDF
2.01 MB Adobe PDF Visualizza/Apri

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11585/683207
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 67
  • ???jsp.display-item.citation.isi??? 60
social impact