原恒星吸积盘的数值模拟.pdfVIP

原恒星吸积盘的数值模拟 张 辉 南京大学天文系 摘 要 这篇论文的主要目的是用流体力学方法模拟在一个原恒星的吸积盘中加入一颗行星后吸 积盘的演化过程。由于背景知识的限制，我们没有对吸积盘的具体物理过程或者具体物理量 进行讨论，我们把这个程序作为一种训练以及日后做深入研究的一个工具。为了保证其可靠 性，我们将得到的结果与几个已经比较成熟的结果(例如.Kley 1999 )进行了比较。最后我们得 出结论，我们的程序是可靠的。 关键词：原行星 吸积盘 数值模拟 一．简 介 自从1995 年，Mayor 和 Queloz 发现第一个地外行星系统以来，外太阳系的研究已经成为 天体物理与天体力学研究的热点。最近的观测显示，地外行星系统与太阳系相比，有很多奇 怪之处。例如：有些行星质量与木星相当(0.4MjMp11Mj)，然而它们的公转轨道却十分接近 中心天体，最近的只有大约 0.04AU ，小于水星到太阳的距离。有的行星具有很大的偏心率 0.0e0.67 。它们是如何演化成这种与太阳系大相径庭的状态？到底它们是“怪胎”还是太阳 系自己是特例？要解释这些问题,研究行星形成后的动力学演化已经不能满足要求了，而需要 上溯到行星形成的早期――原恒星吸积盘的研究。通过研究吸积盘的演化，我们才能更深入 的了解行星形成的过程和机制，对外太阳系和太阳系的研究都有重大意义。而研究原恒星吸 积盘的一种有效工具，就是计算机数值模拟。 现在比较流行的方法有两大类，流体力学方法和多粒子方法（N 体问题），我们选用了前 者作为物理基础，将在第二部分中介绍。基本算法我们参考了 ZEUS2D （Stone Norman 1992） 和其改进算法 FarGo （Masset 2000 ），其基本思想是用有限差分的方法来解欧拉流体力学偏微 分方程组，我们对其中一些地方做了修改（例如：粘性的处理），在第三部分会详细介绍。最 后我们将得到的结果与一些比较成熟的结果（Kley 1999 ）做了比较，证明我们的程序是可靠 的和有效的。对具体演化过程和物理量的具体讨论已经超出本论文的范围。 1 二． 物 理 模 型 1 ．基本方程组 首先，我们研究的模型是一个原恒星的吸积盘，在盘中嵌入一颗原行星。取柱坐 标系（r ，φ，z ），坐标原点在恒星中心，计算平面取在z ＝0 处。坐标系固定，不随盘转动。 不考虑 z 轴方向的变化，因此是一个二维的模拟。考虑压力梯度；考虑中心天体和行星对盘 上物质的引力势；考虑中心天体以及行星对盘的吸积，为简化模型不考虑盘对行星的影响； 同时考虑由于曲线坐标而产生的惯性力（如：离心力）和盘上物质的粘滞力。 在吸积盘中，物质近似做 Keplerian 运动。根据薄盘的假设，即 H/r1 。因此我们可以沿 z 轴方向积分流体力学方程组，我们讨论的变量都是垂直积分后的平均值，这也就是不考虑 z 方向变化的原因。H/r 的取值一般为0.04—0.05，本程序取了 0.05 。 基本的流体力学方程组为： 连续性方程 ∂ρ +∇=⋅ ( ) 0ρU （1） ∂t 动量方程，r 方向 2