paper

Multi-Resolution Radiative Transfer for Line Emission

arXiv:astro-ph/0408321 · doi:10.1086/426112

Abstract

We present describe a new computer code that solves the radiative transfer problem on multi-resolution grids. If the cloud model is from an MHD simulation on a regular cartesian grid, criteria based for example on local density or velocity gradients are used to refine the grid by dividing selected cells into sub-cells. Division can be repeated hierarchically. Alternatively, if the cloud model is from an MHD simulations with adaptive mesh refinement, the same multi-resolution grid used for the MHD simulation is adopted in the radiative transfer calculations. High discretization is often needed only in a small fraction of the total volume. This makes it possible to simulate spectral line maps with good accuracy, also minimizing the total number of cells and the computational cost (time and memory). Multi-resolution models are compared with models on regular grids. In the case of moderate optical depths (tau ~ a few) an accuracy of 10% can be reached with multi-resolution models where only 10% of the cells of the full grid are used. For optically thick species (tau~100), the same accuracy is achieved using 15% of the cells. The relation between accuracy and number of cells is not found to be significantly different in the two MHD models we have studied. The new code is used to study differences between LTE and non-LTE spectra and between isothermal and non-isothermal cloud models. We find significant differences in line ratios and individual spectral line profiles of the isothermal and LTE models relative to the more realistic non-isothermal case. The slope of the power spectrum of integrated intensity is instead very similar in all models.

Submitted to ApJ, 34 pages, 14 figures

References in corpus (6)

Cited by in corpus (8)