[论文解读] 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，采用通量限制扩散方案进行辐射传热，包含查表得到的Rosseland平均消光截面和基于边缘检测的黑体辐射模型。
• AMR代码FLASH用于对比，采用自适应网格加密，并设置高初始网格分辨率以解析主导螺旋模态的波长。
• 通过基于立体角的边缘检测算法识别边缘粒子，辐射损失通过包含表面积和边缘度因子的修正能量方程进行建模。
• 能量方程中引入通量限制项以防止出现非物理解的热流，冷却 timescale 与轨道 timescale 对比以评估碎片化可能性。
• 模拟采用等温与辐射方程状态，结果在不同代码和分辨率水平间进行对比。
• 关键准则为：初始网格或光滑长度必须能解析最强螺旋模态在刚进入非线性状态时的波长。

实验结果

研究问题

• RQ1 为何不同数值代码在R < 50 AU处报告的盘碎片化结果存在矛盾？
• RQ2 初始网格或光滑长度的分辨率在多大程度上决定了盘碎片化模拟的准确性？
• RQ3 不同的辐射传热与热力学处理方式如何影响原行星盘碎片化的可能性？
• RQ4 若初始分辨率不足以解析主导螺旋模态，SPH结合自适应软化能否产生可靠结果？
• RQ5 AMR与SPH代码在存在自引力与辐射冷却的情况下，正确模拟盘碎片化的最低分辨率要求是什么？

主要发现

• 当分辨率相近且初始网格或光滑长度能解析最强螺旋模态波长时，SPH与AMR代码在盘碎片化模拟中结果极为相似。
• AMR中大量网格加密是必要的但不充分；初始分辨率必须足够高，以从一开始就捕捉到最不稳定的模态的波长。
• 若初始分辨率无法解析主导螺旋模态波长，即使采用自适应软化，SPH模拟的精度也不优于固定软化。
• 必须解析主导螺旋模态波长的要求，比基于局部Jeans长度的传统Bate & Burkert（1997）准则更为严格。
• GASOLINE与FLASH在相近分辨率下表现出极佳的一致性，表明当分辨率足够时，SPH在自引力与热力学复杂系统中具有鲁棒性。
• 本研究为未来代码比较建立了基准，强调在盘碎片化模拟中必须采用相同的初始条件与分辨率标准。