A&A 496, 597-608 (2009)
DOI: 10.1051/0004-6361/200811220

Global MHD simulations of stratified and turbulent protoplanetary discs

II. Dust settling

S. Fromang^{1, 2} and R. P. Nelson<sup>3

1</sup>  CEA, Irfu, SAp, Centre de Saclay, 91191 Gif-sur-Yvette, France
    e-mail: sebastien.fromang@cea.fr
²  UMR AIM, CEA-CNRS-Univ. Paris VII, Centre de Saclay, 91191 Gif-sur-Yvette, France
³  Astronomy Unit, Queen Mary, University of London, Mile End Road, London E1 4NS, UK 


Received 24 October 2008 / Accepted 15 January 2009

Abstract
Aims. The aim of this paper is to study the vertical profile of small dust particles in protoplanetary discs in which angular momentum transport is due to MHD turbulence driven by the magnetorotational instability. We consider particle sizes that range from approximately 1 micron up to a few millimeters.
Methods. We use a grid-based MHD code to perform global two-fluid simulations of turbulent protoplanetary discs which contain dust grains of various sizes.
Results. In quasi-steady state, the gravitational settling of dust particles is balanced by turbulent diffusion. Simple and standard models of this process fail to describe accurately the vertical profile of the dust density. The disagreement is larger for small dust particles (of a few microns in size), especially in the disc upper layers (Z>3H, where H is the scale-height). Here there can be orders of magnitude in the disagreement between the simple model predictions and the simulation results. This is because MHD turbulence is not homogeneous in accretion discs, since velocity fluctuations increase significantly in the disc upper layer where a strongly magnetized corona develops. We provide an alternative model that gives a better fit to the simulations. In this model, dust particles are diffused away from the midplane by MHD turbulence, but the diffusion coefficient varies vertically and is everywhere proportional to the square of the local turbulent vertical velocity fluctuations.
Conclusions. The spatial distribution of dust particles can be used to trace the properties of MHD turbulence in protoplanetary discs, such as the amplitude of the velocity fluctuations. In the future, detailed and direct comparison between numerical simulations and observations should prove a useful tool for constraining the properties of turbulence in protoplanetary discs.

Key words: accretion, accretion disks -- Magnetohydrodynamics (MHD) -- methods: numerical -- turbulence

© ESO 2009

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.