uu.seUppsala University Publications
Change search

Magnetohydrodynamics of Turbulent Accretion Disks and Planet Formation
Uppsala University, Teknisk-naturvetenskapliga vetenskapsområdet, Faculty of Science and Technology, Astronomical observatory.
2007 (English)Licentiate thesis, monograph (Other scientific)
Abstract [en]

Planets have long been believed to form in disks of gas and dust around young

stars, interacting with their surroundings via a set

of complex and highly nonlinear processes. In the core accretion scenario for

giant planet formation dust coagulates first into km-sized icy

and rocky planetesimals that further collide, forming

progressively larger solid bodies that eventually give rise to cores of several

Earth masses. If a critical mass is attained, these cores become gas giant

planets by undergoing runaway accretion of gas. Otherwise, just a small amount of

nebular gas is retained by the core, which ends up as an ice giant.

Although explaining the overall architecture of the solar system, this picture was shaken by the discovery of the extra-solar planets, which present a diversity far greater than the standard scenario above could predict. Planet-disk interaction seems to be one of the obvious candidates to account

for this diversity, since planets exchange angular momentum with the disk,

leading to either inward or outward migration. An understanding of the

physical state of accretion disks is therefore essential to provide a detailed picture of

the effect of migration on planetary orbits. Such understanding can be achieved through hydrodynamical modelling of the environment where planets are formed.

On the same token, the study of accretion disks is in itself an interesting field of research, since our

current understanding of the physics of accretion is only marginal: the accretion rates observed in circumstellar environments are far too great to be explained by molecular viscosity alone. A milestone was achieved 15 years ago with the realization that the magneto-rotational instability (MRI) generates a vigorous turbulence that appears to be present under a wide range of conditions in accretion disks. Much work was done since in an effort to understand the properties of turbulent magnetized accretion disks and its implications to stellar accretion and planet formation.

The equations of hydrodynamics are nonlinear and we must therefore resort to numerical

simulations. For all simulations we employ the \pencilcodec, a 3D high-order finite-difference MHD code

using Cartesian coordinates. A cylindrical formulation of {\sc Pencil} is presented, but was not yet used in a publication. Planets and stars are treated as particles

evolved with an $N$-body scheme, which is also introduced. In paper I we test the reliability of the numerical setup by contributing

to a code comparison project that consisted of performing numerical simulations of a disk-planet system

using various grid-based and

smoothed particle hydrodynamics codes. Surface density contours, potential vorticity

and smoothed radial profiles were compared at several times, as well as

the disk mass and gravitational torque time evolution. Our results were consistent between the many

codes and outperformed the only other Cartesian code used in the comparison.

In paper II we develop our models into three-dimensional magnetohydrodynamical simulations

of disks of gas and solids, studying the turbulence generated by the magneto-rotational instability

that is believed to be the source of the aforementioned anomalous viscosity. We find that Cartesian grids

are well-suited for accretion disk problems. The

disk-in-a-box models based on Cartesian grids presented in paper II develop and

sustain MHD turbulence, in good agreement with published results achieved with

cylindrical codes. We investigate the dependence of the magneto-rotational

instability on disk scale height, finding evidence that the turbulence

generated by the magneto-rotational instability grows with thermal pressure. The

turbulent stresses depend on the thermal pressure obeying a power law with index $0.24\pm0.03$, compatible with the value of $0.25$ found in shearing box

calculations. The ratio of stresses decreased with increasing temperature. We

also study the dynamics of boulders in the hydromagnetic turbulence. The

vertical turbulent diffusion of the embedded boulders is comparable to the

turbulent viscosity of the flow. Significant overdensities arise in the solid

component as boulders concentrate in high pressure regions, with interesting implications for planet formation.

We went further, on paper III, to study whether the predictions based on local shearing box simulations

are also valid for global reference frames. By doing so, we help on the disentangling of physics and

numerics of simulated hydromagnetic turbulence, since it has been shown that the saturated state of the MRI depends not only on physical parameters, but also on numerical choices such as box length and resolution. We find that the saturation predictor derived from shearing box data also applies for the case of global disks in the minimum saturation regime where the fastest growing wavelength of the MRI is not resolved. A simulation where this wavelength is resolved shows that the scaling law is also accurate for global disks in this regime. Analytical considerations of the ratio of stresses, however, do not agree with the values found in the non-linear saturated state of the turbulence in the global simulations reported here. A brief summary of ongoing projects is presented in the last chapter.

2007.
National Category
Astronomy, Astrophysics and Cosmology
Identifiers
OAI: oai:DiVA.org:uu-85843DiVA: diva2:116652
Available from: 2008-02-23 Created: 2008-02-23

No full text

Search in DiVA

By organisation
Astronomical observatory
On the subject
Astronomy, Astrophysics and Cosmology