[论文解读] Selfgravitating disks in binary systems: an SPH approach -- I. Implementation of the code and reliability tests

Luis Diego Pinto, R. Capuzzo‐Dolcetta+1 more|arXiv (Cornell University)|Jun 30, 2019

Astrophysics and Star Formation Studies参考文献 59被引用 3

一句话总结

本文提出 GaSPH，一种基于树结构引力求解器的新型光滑粒子流体动力学（SPH）代码，旨在模拟双星系统中自引力原行星盘的动力学行为。该代码通过了严格的稳定性与性能测试，包括 Sedov-Taylor 爆炸波验证，并在双星扰动下展现出盘面偏心率与近日点角演化过程中的收敛性，实现了对非孤立环境中盘动力学的可靠建模。

ABSTRACT

The study of the stability of massive gaseous disks around a star in a non-isolated context is not a trivial issue and becomes a more complicated task for disks hosted by binary systems. The role of self-gravity is thought to be significant, whenever the ratio of the disk to the star mass is non-negligible. To tackle these issues we implemented, tested and applied our own Smoothed Particle Hydrodynamics (SPH) algorithm. The code (named GaSPH) passed various quality tests and shows good performances, so to be reliably applied to the study of disks around stars accounting for self-gravity. This work aims to introduce and describe the algorithm, making some performance and stability tests. It constitutes the first part of a series of studies in which self-gravitating disks in binary systems will be let evolve in larger environments such as Open Clusters.

研究动机与目标

• 开发一种稳健的自引力SPH代码，能够模拟非孤立环境中的原行星盘，特别是双星系统中的情况。
• 实现并测试一种基于树结构的引力求解器，以精确处理近距离粒子相互作用与远距离力。
• 利用标准流体动力学与引力基准测试，验证代码的数值性能、稳定性和准确性。
• 通过建模气体-恒星-双星相互作用并包含自引力，为未来研究开放星团中盘体演化奠定基础。
• 为在高密度恒星环境中研究恒星对盘体的反馈作用提供理论基础。

提出的方法

• 代码采用拉格朗日型SPH格式，使用三次样条核函数进行流体插值，确保平滑长度下具有二阶精度。
• 引力场通过基于树结构的多极展开方案计算，远距离粒子簇近似为单极与四极项。
• 近距离粒子引力相互作用直接计算，避免软化带来的误差，从而在高密度区域保持高精度。
• 点质量（恒星或行星）的运动采用高阶显式时间积分器进行积分。
• 代码实现了气体自引力与气体-恒星间相互引力的自洽处理，不包含辐射传热过程。
• 通过偏心率矢量与单位质量的角动量，计算盘体特性如偏心率与近日点角。

实验结果

研究问题

• RQ1 基于树结构引力求解器的SPH代码在多大程度上能准确模拟双星系统中自引力盘体的演化？
• RQ2 在改变粒子分辨率时，盘体偏心率与近日点角的数值收敛性达到何种水平？
• RQ3 该代码能否以高保真度重现已知的流体动力学基准，如Sedov-Taylor爆炸波？
• RQ4 当双星伴星偏心率 e = 0.4 时，盘体轨道参数如何演化？
• RQ5 与二维模型相比，垂直引力与三维结构在盘体动力学中的作用是什么？

主要发现

• GaSPH代码成功以高精度重现Sedov-Taylor爆炸波，证实其具备处理强激波动力学的能力。
• 在初始混沌阶段后，盘体偏心率稳定在约 0.075 的平均值，与先前研究结果一致。
• 盘体偏心率表现出由双星轨道周期（Pbin ≈ 134 年）调制的小幅度振荡，与伴星近日点通过相关联。
• 盘体近日点角 ωdisk 在粒子数从 N = 20,000 增加至 100,000 的模拟中表现出收敛性，表明数值稳定性良好。
• 该代码在缓慢与剧烈动力学区域均展现出良好的数值性能、稳定性和准确性。
• 三维引力处理可实现垂直方向的角动量输运，引入了二维模型中所不具备的新自由度。