Markov Chain Monte Carlo solutions for radiative transfer problems
Institut f. Angewandte Mathematik,
2 Interdisziplinäres Zentrum für Wissenschaftliches Rechnen, INF 368, 19120 Heidelberg, Germany
3 ZAH, Institut f. Theoret. Astrophysik, Albert-Ueberle-Str. 2, 69120 Heidelberg, Germany
Accepted: 19 September 2010
Algorithms for the solution of the radiative transfer equation that are simultaneously accurate, fast, and easy to implement are still highly desirable, in particular for multidimensional media. We present a stochastic algorithm that solves the transfer equation by assuming that the transfer of radiation can be modelled as a Markov process. It is a generalization of the classical Monte Carlo method and can be applied to the solution of the Milne equation.
Key words: radiative transfer / methods: numerical / methods: statistical / scattering / line: formation
© ESO, 2010