3D continuum radiative transfer

Free access article

Issue		A&A Volume 456, Number 1, September II 2006


Page(s)		1 - 12
Section		Astrophysical processes
DOI		http://dx.doi.org/10.1051/0004-6361:20054410

A&A 456, 1-12 (2006)
DOI: 10.1051/0004-6361:20054410

A new adaptive grid construction algorithm based on the Monte Carlo method

G. Niccolini and J. Alcolea

Observatorio Astronómico Nacional (IGN), Alfonso XII 3, 28014 Madrid, Spain
e-mail: g.niccolini@oan.es

(Received 25 October 2005 / Accepted 12 May 2006 )

Abstract
Aims. The continuum radiative transfer is a common problem in Astrophysics, that must be solved to describe the physics and interpret the observations of objects embedded (for instance) in a dusty medium (e.g. evolved stars, young stellar objects, $\ldots$ ). Extreme high angular resolution observations reveal in general complex geometries (e.g. mediums with clumps of matter, warps $\ldots$ ). It is thus necessary to solve the radiative transfer problem making no simplifying assumptions on the geometry, i.e. resolve the full spatial three dimensional problem.
Methods. We use the flexible Monte Carlo (MC) method to solve the continuum radiative transfer problem. We focus on the computation of the temperature structure of a dusty medium and the construction of the grid that samples this structure. We propose to use a natural by-product of the MC method, the locations of the photon package absorptions in the medium, to build an adaptive grid to reach the finest resolution in the regions of interest.
Results. We present the details of the new algorithm used for the construction of the grid. We also propose modifications to a widely used MC propagation framework. The reliability of our method is numerically tested against well established numerical results in one and two dimensions. In order to illustrate the capabilities of our three dimensional methods, we examine the observational evidence for the presence of warps in circumstellar discs.

Key words: radiative transfer -- methods: numerical -- stars: AGB and post-AGB -- stars: planetary systems: protoplanetary disks

What is OpenURL?

The OpenURL standard is a protocol for transmission of metadata describing the resource that you wish to access. An OpenURL link contains article metadata and directs it to the OpenURL server of your choice. The OpenURL server can provide access to the resource and also offer complementary services (specific search engine, export of references...). The OpenURL link can be generated by different means.

If your librarian has set up your subscription with an OpenURL resolver, OpenURL links appear automatically on the abstract pages.
You can define your own OpenURL resolver with your EDPS Account. In this case your choice will be given priority over that of your library.
You can use an add-on for your browser (Firefox or I.E.) to display OpenURL links on a page (see http://www.openly.com/openurlref/). You should disable this module if you wish to use the OpenURL server that you or your library have defined.