Issue |
A&A
Volume 689, September 2024
|
|
---|---|---|
Article Number | A180 | |
Number of page(s) | 9 | |
Section | Numerical methods and codes | |
DOI | https://doi.org/10.1051/0004-6361/202451423 | |
Published online | 12 September 2024 |
DIES: Parallel dust radiative transfer program with the immediate re-emission method
Department of Physics,
PO Box 64, University of Helsinki,
00014
Helsinki,
Finland
Received:
8
July
2024
Accepted:
10
August
2024
Context. Radiative transfer (RT) modelling is a necessary tool in the interpretation of observations of the thermal emission of interstellar dust. It is also often part of multi-physics modelling. In this context, the efficiency of radiative transfer calculations is important, even for one-dimensional models.
Aims. We investigate the use of the so-called immediate re-emission (IRE) method for fast calculation of one-dimensional spherical cloud models. We wish to determine whether weighting methods similar to those used in traditional Monte Carlo simulations can speed up the estimation of dust temperature.
Methods. We present the program DIES, a parallel implementation of the IRE method, which makes it possible to do the calculations also on graphics processing units (GPUs). We tested the program with externally and internally heated cloud models, and examined the potential improvements from the use of different weighted sampling schemes.
Results. The execution times of the program compare favourably with previous programs, especially when run on GPUs. On the other hand, weighting schemes produce only limited improvements. In the case of an internal radiation source, the basic IRE method samples the re-emission well, while traditional Monte Carlo requires the use of spatial importance sampling. Some noise reduction could be achieved for externally heated models by weighting the initial photon directions. Only in optically very thin models does weighting – such as the proposed method of forced first interaction – result in noise reduction by a factor of several.
Conclusions. The IRE method performs well for both internally and externally heated models, typically without the need for any additional weighting schemes. With run times of the order of one second for our test models, the DIES program is suitable even for larger parameter studies.
Key words: radiative transfer / methods: numerical / ISM: clouds / dust, extinction
© The Authors 2024
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.