Volume 439, Number 1, August III 2005
|Page(s)||153 - 158|
|Section||Interstellar and circumstellar matter|
|Published online||22 July 2005|
Monte Carlo radiative transfer in SPH density fields
School of Physics & Astronomy, Cardiff University, 5 The Parade, Cardiff CF24 3YB, Wales, UK e-mail: D.Stamatellos@astro.cf.ac.uk
Accepted: 25 April 2005
We combine a Monte Carlo radiative transfer code with an SPH code, so that – assuming thermal equilibrium – we can calculate dust-temperature fields, spectral energy distributions, and isophotal maps, for the individual time-frames generated by an SPH simulation. On large scales, the radiative transfer cells (RT cells) are borrowed from the tree structure built by the SPH code, and are chosen so that their size – and hence the resolution of the calculated temperature field – is comparable with the resolution of the density field. We refer collectively to these cubic RT cells as the global grid. The code is tested and found to treat externally illuminated dust configurations very well. However, when there are embedded discrete sources, i.e. stars, these produce very steep local temperature gradients which can only be modelled properly if – in the immediate vicinity of, and centred on, each embedded star – we supplement the global grid with a star grid of closely spaced concentric RT cells.
Key words: ISM: clouds / ISM: structure / ISM: dust, extinction / methods: numerical / radiative transfer / hydrodynamics
© ESO, 2005
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.