[論文レビュー] Rossby Wave Instability in three dimensional discs

H. Méheut, Cong Yu+1 more|arXiv (Cornell University)|Mar 2, 2012

Astrophysics and Star Formation Studies被引用数 38

ひとこと要約

本稿は、3次元的層状化・等温な原始惑星系円盤におけるロスビー波不安定性（RWI）の線形解析を提示し、3次元渦が顕著な垂直速度構造を示すことを示している。3Dの増幅率は2Dと比較してわずかに低下するにとどまり、数値シミュレーションは線形固有関数を正確に再現しており、円盤のデッドゾーンにおける非線形的・長期的渦の進化を研究する際にその有用性が裏付けられている。

ABSTRACT

The Rossby wave instability (RWI) is a promising mechanism for producing large-scale vortices in protoplanetary discs. The instability operates around a density bump in the disc, and the resulting vortices may facilitate planetesimal formation and angular momentum transfer in the disc dead zone. Most previous works on the RWI deal with two-dimensional (height-integrated) discs. However, vortices may have different dynamical behaviours in 3D than in 2D. Recent numerical simulations of the RWI in 3D global discs by Meheut et al. have revealed intriguing vertical structure of the vortices, including appreciable vertical velocities. In this paper we present a linear analysis of the RWI in 3D global models of isothermal discs. We calculate the growth rates of the Rossby modes (of various azimuthal wave numbers m = 2 - 6) trapped around the fiducial density bump and carry out 3D numerical simulations to compare with our linear results. We show that the 3D RWI growth rates are only slightly smaller than the 2D growth rates, and the velocity structures seen in the numerical simulations during the linear phase are in agreement with the velocity eigenfunctions obtained in our linear calculations. This numerical benchmark shows that numerical simulations can accurately describe the instability. The angular momentum transfer rate associated with Rossby vortices is also studied.

研究の動機と目的

• 原始惑星系円盤におけるロスビー波不安定性（RWI）によって生成されるロスビー渦の3次元的速度構造を理解すること。
• 垂直方向の層状化が2Dモデルと比較して、RWIの増幅率および固有モードに与える影響を評価すること。
• 線形摂動理論と比較することで、3次元的グローバル数値シミュレーションのRWIの妥当性を検証すること。
• 角運動量の輸送と垂直速度の役割が、渦のダイナミクスおよび地球型惑星形成に与える影響を検討すること。

提案手法

• 摂動法を用いて、3次元的層状化・等温円盤におけるRWIのグローバル線形安定性解析を実施する。
• 基準密度バンプを持つ円盤において、方位角波番号 $m = 2$ から $6$ のロスビーモードの固有値問題を解く。
• 垂直方向にヘルミート多項式展開を用いて、3次元的速度および密度構造をモデル化する。
• アングル方向フーリエ変換を用いて3次元非線形シミュレーションから抽出した固有関数（振幅、位相、垂直構造）と線形固有関数を比較する。
• ロスビー波が運ぶ角運動量フラックスを計算し、不安定性のメカニズムおよび密度バンプへの影響を理解する。

実験結果

リサーチクエスチョン

• RQ1 3次元的垂直速度構造は、原始惑星系円盤におけるロスビー渦の増幅および安定性にどのように影響を与えるか？
• RQ2 垂直方向の層状化の影響により、3次元的RWI増幅率は2次元的増幅率と比べてどの程度低下するか？
• RQ3 3次元的数値シミュレーションは、特にモードの垂直構造を的確に再現できるか？
• RQ4 ロスビー波による角運動量輸送が、密度バンプの進化および渦の発達に果たす役割は何か？
• RQ5 垂直速度の存在が、渦内でのダスト集積および地球型惑星形成にどのように寄与するか？

主な発見

• ロスビーモードの3次元的線形増幅率は、2次元的対応物と比較してわずかに小さいが、その減少は垂直構造に起因するコリレーション共鳴領域での波の吸収に起因する。
• 数値的シミュレーションは、複雑な径方向および垂直速度構造を含めた線形固有関数を的確に再現しており、非線形的解析に於いてその信頼性が裏付けられている。
• 垂直速度成分は、特にコリレーション領域において顕著であり、非無視可能な大きさを示しており、線形理論の予測と整合的である。
• 正規化されたアングル方向速度固有関数について、シミュレーションと線形計算の結果は良好に一致しており、シミュレーションフレームワークの妥当性が検証された。
• 角運動量フラックスは、渦度最大値の両側で相互作用する2つのロスビー波によって運ばれ、不安定性を駆動し、初期の密度バンプの侵食を促進する。
• $m \geq 4$ モードでは、非線形シミュレーションにおける増幅率が線形理論の予測よりわずかに低く、これは線形モデルに含まれない数値的粘性または高次垂直モード寄与の影響による可能性がある。