[论文解读] 3D-radiation hydro simulations of disk-planet interactions: I. Numerical algorithm and test cases

Hubert Klahr, W. Kley|ArXiv.org|Oct 13, 2005

Astrophysics and Star Formation Studies参考文献 29被引用 64

一句话总结

本研究提出了一种三维辐射流体动力学代码（TRAMP），用于模拟原行星-盘相互作用并实现完整的热力学处理，表明正确的能量传输会导致形成厚实的、由压力支撑的原行星包层，其H/R ≈ 0.5，而非薄的开普勒盘。关键结果是，这种结构抑制了原行星物质中卫星的形成，原因在于固态物质的快速内向漂移，同时与等温模型相比，扭矩波动显著增加。

ABSTRACT

We study the evolution of an embedded protoplanet in a circumstellar disk using the 3D-Radiation Hydro code TRAMP, and treat the thermodynamics of the gas properly in three dimensions. The primary interest of this work lies in the demonstration and testing of the numerical method. We show how far numerical parameters can influence the simulations of gap opening. We study a standard reference model under various numerical approximations. Then we compare the commonly used locally isothermal approximation to the radiation hydro simulation using an equation for the internal energy. Models with different treatments of the mass accretion process are compared. Often mass accumulates in the Roche lobe of the planet creating a hydrostatic atmosphere around the planet. The gravitational torques induced by the spiral pattern of the disk onto the planet are not strongly affected in the average magnitude, but the short time scale fluctuations are stronger in the radiation hydro models. An interesting result of this work lies in the analysis of the temperature structure around the planet. The most striking effect of treating the thermodynamics properly is the formation of a hot pressure--supported bubble around the planet with a pressure scale height of H/R ~ 0.5 rather than a thin Keplerian circumplanetary accretion disk. We also observe an outflow of gas above and below the planet during the gap opening phase.

研究动机与目标

• 开发并测试一种用于模拟原行星-盘相互作用的三维辐射流体动力学代码（TRAMP），实现完整的热力学处理。
• 评估辐射传输、数值分辨率和质量吸积建模对间隙形成与行星迁移的影响。
• 研究在原恒星盘中质量为木星质量的行星周围的三维温度与密度结构。
• 评估真实热力学对原行星吸积与卫星形成的影响。
• 将辐射流体动力学与常用的局部等温近似进行比较。

提出的方法

• TRAMP代码求解三维流体动力学方程，并采用通量限制扩散近似处理辐射传输。
• 状态方程包含热能与内能演化，可实现对吸积光度加热的正确处理。
• 模拟采用固定质量为木星质量的行星位于5 AU处，在三维盘中使用嵌套网格加密围绕行星。
• 模型包含可变的质量吸积率，并追踪行星吸积过程中的能量沉积。
• 边界和初始条件设定为模拟一个安静的原行星盘，包含中心恒星和大质量行星。
• 系统性地改变数值参数（如分辨率和吸积处理方式）以检验结果的稳健性。

实验结果

研究问题

• RQ1 辐射传输如何影响三维行星-盘模拟中的间隙形成过程与扭矩演化？
• RQ2 当真实处理热力学时，大质量原行星周围的三维温度与密度结构是怎样的？
• RQ3 原行星气体的特性（如标高与旋转速度分布）与薄盘近似相比有何不同？
• RQ4 数值分辨率与质量吸积建模在多大程度上影响模拟结果？
• RQ5 是否能形成围绕原行星的压力支撑包层？其对卫星形成有何影响？

主要发现

• 原行星气体形成一个厚实的、由压力支撑的包层，其压力标高H/R ≈ 0.5，远厚于薄盘近似。
• 原行星区域的旋转速度仅为开普勒速度的约50%，表明存在强烈的热压力支撑。
• 厚包层结构使薄盘近似失效，表明应为原行星云而非开普勒吸积盘。
• 若在0.1罗氏半径处释放大质量固态物质，其在一个月内即可落入行星，表明原行星物质中难以形成卫星。
• 在具有真实吸积的辐射流体动力学模型中，扭矩波动比等温模型大一个数量级。
• 在早期间隙形成阶段观测到气体的垂直喷流，持续数万年，可能在散射光中可探测。