Global numerical simulations of differentially rotating disks with free eccentricity

Free access article

Issue		A&A Volume 432, Number 3, March IV 2005


Page(s)		757 - 769
Section		Astrophysical processes
DOI		http://dx.doi.org/10.1051/0004-6361:20041948

A&A 432, 757-769 (2005)
DOI: 10.1051/0004-6361:20041948

Global numerical simulations of differentially rotating disks with free eccentricity

J. C. B. Papaloizou

Astronomy Unit, Queen Mary, University of London, Mile End Rd, London E1 4NS, UK
e-mail: jcbp@maths.qmw.ac.uk

(Received 3 September 2004 / Accepted 10 December 2004)

Abstract
We study the nonlinear evolution of global m=1 modes with low pattern speeds in a differentially rotating non magnetic disk, by means of two and three dimensional numerical simulations. The modes make disk streamlines eccentric with maximum eccentricity in the range 0.13 - 0.24.

We found that long lived patterns corresponding to eccentricities ~0.1 lasted for the duration of the simulations of ~64 orbits at the disk outer boundary. They had slow retrograde precession with period ~10 orbits at the outer boundary, in good agreement with that found from linear normal mode analysis.

As expected from linear stability analysis, which leads one to expect a parametric instability associated with the non circular streamlines, we also found that three dimensional simulations showed the growth of a local instability to producing vertical motions on a local scale that eventually results in small amplitude turbulence with root mean square vertical velocity typically $\sim$ $0.03 c_{\rm s},$ $c_{\rm s}$ being the sound speed. This turbulence together with much more effective shocks produced by the interaction of the eccentric disk with the inner boundary were responsible for damping the disk eccentricity. Our first estimate of the damping time associated with turbulence for eccentricities ~0.1, corresponded to that associated with a turbulent viscosity acting to circularize the eccentric streamlines with " $\alpha$ " parameter ~10^-3. For parameters appropriate to protostellar disks at 5 AU this corresponds to a decay time $\sim$ $10^{5{-}6} ~\rm y.$

Thus although the indication is that disks with eccentric streamlines are long lived, there is associated turbulence even in the non magnetic case leading to estimated decay times that may be significant when the possibility of the growth of orbital eccentricity for extrasolar planets is considered through disk planet interaction on sufficiently long timescales.

Key words: accretion, accretion disks -- instabilities -- hydrodynamics -- celestial mechanics -- planets and satellites: formation

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.