FORMING PLANETESIMALS BY GRAVITATIONAL INSTABILITY. I. THE ROLE OF THE RICHARDSON NUMBER IN TRIGGERING THE KELVIN-HELMHOLTZ INSTABILITY

被引:39
作者
Lee, Aaron T. [1 ]
Chiang, Eugene [1 ,2 ]
Asay-Davis, Xylar [3 ]
Barranco, Joseph [4 ]
机构
[1] Univ Calif Berkeley, Dept Astron, Berkeley, CA 94720 USA
[2] Univ Calif Berkeley, Dept Earth & Planetary Sci, Berkeley, CA 94720 USA
[3] Los Alamos Natl Lab, Ctr Nonlinear Studies, Los Alamos, NM 87545 USA
[4] San Francisco State Univ, Dept Phys & Astron, San Francisco, CA 94132 USA
基金
美国国家科学基金会;
关键词
hydrodynamics; instabilities; planets and satellites: formation; protoplanetary disks; PROTOPLANETARY DISKS; SOLAR NEBULA; PROTOSTELLAR DISKS; DUST LAYER; TURBULENCE; CONVECTION; PARTICLES; GROWTH; SHEAR; COAGULATION;
D O I
10.1088/0004-637X/718/2/1367
中图分类号
P1 [天文学];
学科分类号
0704 ;
摘要
Gravitational instability (GI) of a dust-rich layer at the midplane of a gaseous circumstellar disk is one proposed mechanism to form planetesimals, the building blocks of rocky planets and gas giant cores. Self-gravity competes against the Kelvin-Helmholtz instability (KHI): gradients in dust content drive a vertical shear which risks over-turning the dusty subdisk and forestalling GI. To understand the conditions under which the disk can resist the KHI, we perform three-dimensional simulations of stratified subdisks in the limit that dust particles are small and aerodynamically well coupled to gas, thereby screening out the streaming instability and isolating the KHI. Each subdisk is assumed to have a vertical density profile given by a spatially constant Richardson number Ri. We vary Ri and the midplane dust-to-gas ratio mu(0) and find that the critical Richardson number dividing KH-unstable from KH-stable flows is not unique; rather, Ri(crit) grows nearly linearly with mu(0) for mu(0) = 0.3-10. Plausibly, a linear dependence arises for mu(0) << 1 because in this regime the radial Kepler shear replaces vertical buoyancy as the dominant stabilizing influence. Why this dependence should persist at mu(0) > 1 is a new puzzle. The bulk (height-integrated) metallicity is uniquely determined by Ri and mu(0). Only for disks of bulk solar metallicity is Ri(crit) approximate to 0.2, which is close to the classical value. Our empirical stability boundary is such that a dusty sublayer can gravitationally fragment and presumably spawn planetesimals if embedded within a solar metallicity gas disk similar to 4 x more massive than the minimum-mass solar nebula; or a minimum-mass disk having similar to 3 x solar metallicity; or some intermediate combination of these two possibilities. Gravitational instability seems possible without resorting to the streaming instability or to turbulent concentration of particles.
引用
收藏
页码:1367 / 1377
页数:11
相关论文
共 54 条
[1]  
[Anonymous], 2005, ARXIVASTROPH0508659
[2]  
[Anonymous], 2004, Hydrodynamic Stability
[3]  
[Anonymous], 2009, ARXIV09060854
[4]   HEAT AND DUST IN ACTIVE LAYERS OF PROTOSTELLAR DISKS [J].
Bai, Xue-Ning ;
Goodman, Jeremy .
ASTROPHYSICAL JOURNAL, 2009, 701 (01) :737-755
[5]  
Barranco J.A., 2000, Center for Turbulence Research - Proceedings of the 2000 Summer Program, P97
[6]   Three-dimensional vortices in stratified protoplanetary disks [J].
Barranco, JA ;
Marcus, PS .
ASTROPHYSICAL JOURNAL, 2005, 623 (02) :1157-1170
[7]   A 3D spectral anelastic hydrodynamic code for shearing, stratified flows [J].
Barranco, Joseph A. ;
Marcus, Philip S. .
JOURNAL OF COMPUTATIONAL PHYSICS, 2006, 219 (01) :21-46
[8]   THREE-DIMENSIONAL SIMULATIONS OF KELVIN-HELMHOLTZ INSTABILITY IN SETTLED DUST LAYERS IN PROTOPLANETARY DISKS [J].
Barranco, Joseph A. .
ASTROPHYSICAL JOURNAL, 2009, 691 (02) :907-921
[9]   The growth mechanisms of macroscopic bodies in protoplanetary disks [J].
Blum, Juergen ;
Wurm, Gerhard .
ANNUAL REVIEW OF ASTRONOMY AND ASTROPHYSICS, 2008, 46 :21-56
[10]  
Chandrasekhar S, 1981, HYDRODYNAMIC HYDROMA