EDP Sciences
Free access
Issue A&A
Volume 452, Number 3, June IV 2006
Page(s) 751 - 762
Section Astrophysical processes
DOI http://dx.doi.org/10.1051/0004-6361:20054612

A&A 452, 751-762 (2006)
DOI: 10.1051/0004-6361:20054612

Dust settling in local simulations of turbulent protoplanetary disks

S. Fromang1, 2 and J. Papaloizou2

1  Astronomy Unit, Queen Mary, University of London, Mile End Road, London E1 4NS, UK
2  Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Centre for Mathematical Sciences, Wilberforce Road, Cambridge, CB3 0WA, UK
    e-mail: S.Fromang@damtp.cam.ac.uk

(Received 30 November 2005 / Accepted 25 February 2006)

Aims.In this paper, we study the effect of MHD turbulence on the dynamics of dust particles in protoplanetary disks. We vary the size of the particles and relate the dust evolution to the turbulent velocity fluctuations.
Methods.We performed numerical simulations using two Eulerian MHD codes, both based on finite difference techniques: ZEUS-3D and NIRVANA. These were local shearing box simulations incorporating vertical stratification. Both ideal and non ideal MHD simulations with midplane dead zones were carried out. The codes were extended to incorporate different models for the dust as an additional fluid component. Good agreement between results obtained using the different approaches was obtained.
Results.The simulations show that a thin layer of very small dust particles is diffusively spread over the full vertical extent of the disk. We show that a simple description obtained using the diffusion equation with a diffusion coefficient simply expressed in terms of the velocity correlations accurately matches the results. Dust settling starts to become apparent for particle sizes of the order of 1 to 10 centimeters for which the gas begins to decouple in a standard solar nebula model at 5.2 AU. However, for particles which are 10 centimeters in size, complete settling toward a very thin midplane layer is prevented by turbulent motions within the disk, even in the presence of a midplane dead zone of significant size.
Conclusions.These results indicate that, when present, MHD turbulence affects dust dynamics in protoplanetary disks. We find that the evolution and settling of the dust can be accurately modelled using an advection diffusion equation that incorporates vertical settling. The value of the diffusion coefficient can be calculated from the turbulent velocity field when that is known for a time of several local orbits.

Key words: accretion, accretion disks -- magnetohydrodynamics (MHD) -- methods: numerical -- planets and satellites: formation

© ESO 2006

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.

Editor-in-Chief: T. Forveille
Letters Editor-in-Chief: J. Alves
Managing Editor: C. Bertout

ISSN: 0004-6361 ; e-ISSN: 1432-0746
Frequency: 12 volumes per year
Published by: EDP Sciences

Mirror sites: CDS | EDP Sciences
  RSS feeds
© The European Southern Observatory (ESO)