[論文レビュー] Protoplanetary disk fragmentation with varying radiative physics, initial conditions and numerical techniques

Lucio Mayer, A. Gawryszczak|ArXiv.org|Oct 19, 2007

Astro and Planetary Science参考文献 2被引用数 6

ひとこと要約

本論文は、放射物理、初期条件、数値的手法を変えて、原始惑星系円盤における重力不安定性および破壊を、SPHおよびAMRシミュレーションを用いて調査している。高解像度のシミュレーションで初期グリッド解像度を適切に設定した場合、コードに依存せず一貫して円盤の破壊を捉えることができ、SPHで適応的ソフトニングを用いることで、数十AUにおける破壊に関する長年の文献上の矛盾を解消する結果が得られた。

ABSTRACT

We review recent results of SPH simulations of gravitational instability in gaseous protoplanetary disks,emphasizing the role of thermodynamics in both isolated and binary systems. Contradictory results appeared in the literature regarding disk fragmentation at tens of AU from the central star are likely due to the different treatment of radiation physics as well as reflecting different initial conditions. Further progress on the subject requires extensive comparisons between different codes with the requirement that the same initial conditions are adopted. It is discussed how the local conditions of the disks undergoing fragmentation at $R < 25$ AU in recent SPH simulations are in rough agreement with the prediction of analytical models, with small differences being likely related to the inability of analytical models to account for the dynamics and thermodynamics of three-dimensional spiral shocks. We report that radically different adaptive hydrodynamical codes, SPH and adaptive mesh refinement (AMR), yield very similar results on disk fragmentation at comparable resolution in the simple case of an isothermal equation of state. A high number of refinements in AMR codes is necessary but not sufficient to correctly follow fragmentation, rather an initial resolution of the grid high enough to capture the wavelength of the strongest spiral modes when they are still barely nonlinear is essential. These tests represent a useful benchmark and a starting point for a forthcoming code comparison with realistic radiation physics.

研究の動機と目的

• 異なる数値処理および初期条件による、数十AUにおける円盤破壊に関する文献上の矛盾を解消すること。
• 同一の初期条件および放射物理を用いた場合に、SPHおよびAMRコードの円盤破壊シミュレーションに対する耐久性を検証すること。
• 特に非軸対称モードの最も強い波長を正確に捉えるために必要な最小解像度要件を特定すること。
• SPHにおける適応的ソフトニングが正確性を向上させるのか、それとも初期解像度の方がより重要なのかを評価すること。
• 原始惑星系円盤のシミュレーションにおいて、現実的な放射物理を含む今後のコード比較のベンチマークを確立すること。

提案手法

• シミュレーションは、フラックス制限付き拡散法を用いたSPHコードGASOLINEを用い、表形式のローゼラン平均透過率と、黒体放射のためのエッジ検出法を含む放射移動を扱う。
• 比較のためAMRコードFLASHを用い、メッシュの適応的細分化と、支配的スパイラルモードの波長を解像できる高い初期グリッド解像度を確保する。
• エッジ粒子は、立体角に基づくエッジ検出アルゴリズムを用いて特定され、放射損失は表面積およびエッジ度要因を含む修正されたエネルギー方程式でモデル化される。
• 物理的でない熱フラックスを防ぐために、エネルギー方程式にフラックスリミッター項が含まれており、冷却時定数と軌道時定数を比較して破壊の可能性を評価する。
• 等温および放射的状態方程式を用いたシミュレーションを実行し、コード間および解像度レベル間での結果を比較する。
• 重要な基準として、初期グリッドまたはスムージング長が、非線形にわずかに達した段階での最も強いスパイラルモードの波長を解像している必要がある。

実験結果

リサーチクエスチョン

• RQ1 なぜ異なる数値コードは、R < 50 AUにおける円盤破壊に関して、矛盾する結果を報告するのか？
• RQ2 初期グリッド解像度またはスムージング長が、円盤破壊シミュレーションの正確性をどの程度決定づけるのか？
• RQ3 放射移動および熱力学の異なる取り扱い方は、原始惑星系円盤における破壊の可能性にどのように影響するのか？
• RQ4 初期解像度が支配的スパイラルモードの波長を解像できない場合、SPHで適応的ソフトニングを用いても信頼できる結果が得られるのか？
• RQ5 AMRおよびSPHコードが、自己重力および放射冷却を伴う状況で円盤破壊を正しくシミュレートするための最小解像度要件は何か？

主な発見

• SPHおよびAMRコードは、解像度が同等で、初期グリッドまたはスムージング長が最も強いスパイラルモードの波長を解像している場合、円盤破壊に関して非常に類似した結果を示す。
• AMRでは多数の細分化が必要だが、それだけでは不十分であり、初期解像度が最も不安定なモードの波長を最初から捉えられるほど高くていなければならない。
• 初期解像度が支配的スパイラルモードの波長を解像できない場合、適応的ソフトニングを用いたSPHシミュレーションも、固定ソフトニングと比べて本質的に正確性に優れない。
• 支配的スパイラルモードの波長を解像するという要件は、局所的ジェイン長に基づく従来のBate & Burkert (1997)の基準よりも厳密である。
• GASOLINEとFLASHが同等の解像度で優れた一致を示したことは、解像度が十分にあれば、SPHが自己重力および熱力学的に複雑な系に対しても頑健であることを示唆している。
• 本研究は、今後のコード比較のためのベンチマークを確立した。初期条件と解像度基準が同一であることが、円盤破壊シミュレーションにおいて不可欠であることを強調している。