宇宙里有风么？

三年前的问题，尽管我现在才看到，但还是想答一下。

首先，宇宙中有风吗？

如果你是指"Wind"这个词的话，那答案是有，而且还有很多种，比如星风（Stellar Wind），以及问题描述中提到的盘风（disk wind）。既然问题描述是关于盘风的，那我就介绍一下盘风好了。

什么是盘风？它与地球上的风又有什么不同？

从吸积流的温度的角度进行划分的话，黑洞吸积盘被分为"冷"和"热"两种。其中热吸积盘模型被认为可以很好的诠释低光度活动星系核的机制以及黑洞X射线双星的低/硬态。

研究表明黑洞热吸积流中存在着很强的外流，外流产生的物理原因是由于吸积流的熵随着半径的减小而增大，因此吸积流是对流不稳定的，会产生对流的外流。因此，与地球上风的定义不同，热吸积流中的风可以被理解为是一种非准直的外流。

也就是说，气体在被黑洞吸积、落入黑洞的同时，会产生风。没错，和黑洞热吸积流中产生的喷流的很像，但区别在于以下两点：

  1. 喷流是准直性很强的结构，而风是非准直性的，风的立体角很大；
  2. 喷流的速度可以达到相对论速度或接近光速，而相比来说，风的速度很慢。

而按照问题描述中提到的袁峰老师本人的解释，风"指的是非常发散（即张开的立体角很大，如几十度），运动速度远低于光速的物质外流"。

那么天文学家是怎么发现风的存在的？

这要从1999年说起了，那一年萨比星人入侵...咳咳不好意思拿错剧本了。

1999年，普林斯顿大学的James.Stone利用他们自己开发的Zeus-2D程序完成了第一个黑洞吸积流的整体数值模拟工作。这一开创性的工作最重要的发现是，吸积率并不是大家之前认为的是半径的常数，而是随着半径的减小而减小。

在这项工作中，Stone首先在Zeus-2D中采用球坐标系 r,\theta,\phi 设置黑洞的初始条件。设定一个中心位于1r_{g} 处的具有一定角动量和密度分布的厚的环面，其所处环境气体压强 P=A\rho^{\gamma} ， A 和 \gamma 为给定的系数， \rho 为密度，且有如下关系式：

  • \frac{P}{\rho}=\frac{GM}{(n+1)R_{0}}[\frac{R_{0}}{r}-\frac{1}{2}(\frac{R_{0}}{rsin\theta})^{2}-\frac{1}{2d}]
  • \rho_{max}=[\frac{1}{(n+1)A}(\frac{d-1}{2d})]

其中 n=(\gamma-1)^{-1} 是多方指数， R_{0} 是环面的中心也是密度最大的位置， d=1.125 是给定的畸变参数，取 \rho_{max}=1 、 G=M=R_{0}=1 以及 \gamma=\frac{5}{3} ，并通过计算的出环面中星际物质所具有的角动量是一个常数，该常数为 v_{\theta}=\frac{1}{rsin\theta} 。

设置好参数之后利用超级计算机计算便得出了环面密度在吸积过程随时间演化的图像：

随后Stone通过计算得出了质量内流率、外流率和净吸积率的表达式：

  • \dot{M}_{acc}(r)=2\pi r^{2}\int^{\pi}_{0}\rho v_{r}sin\theta d\theta
  • \dot{M}_{in}(r)=-2\pi r^{2}\int^{\pi}_{0}\rho min(v_{r},0)sin\theta d\theta
  • \dot{M}_{out}(r)=2\pi r^{2}\int^{\pi}_{0}\rho max(v_{r},0)sin\theta d\theta
  • \dot{M}_{net}=\dot{M}_{in}-\dot{M}_{out}

然后根据之前的数值模拟所得出的数据计算质量内流率、外流率和净吸积率，得到吸积率的图像：

从图中我们可以发现，质量内流率和外流率都随着半径的减小而降低。质量外流率和内流率可以近似看做半径的幂律函数。这是很令人吃惊的结果，因为以前的解析研究几乎无一例外全都假设黑洞吸积的质量吸积率是半径的常数。

为解释这一数值模拟结果，国际上两个著名的研究组提出了两个竞争性的模型。这两个模型分别为斯坦福大学的Blandford与他的合作者提出的ADIOS模型 (adiabaticinflow-outflow solution)，以及哈佛大学的Narayan与他的合作者们提出的CDAF模型(convection-dominated accretion flow)。

在 ADIOS 模型中，作者假设由于风的存在，使得吸积气体逐渐损失，导致吸积率随着半径减小而减小。至于风是怎样产生的，ADIOS模型没有给出解释。

CDAF模型则是基于基于热吸积流是对流不稳定的假设而建立的。在CDAF模型中，对流向内转移角动量的效率和粘滞向外转移角动量的效率相同，互相抵消。因此，在黑洞附近会形成一个对流元包。对流元包中，气体始终在向内和向外做对流运动，气体逐渐被锁在对流元包中而没有被吸积进入黑洞，这导致了吸积率随半径的变化。

这两个模型相互竞争了十几年的时间，期间邵逸夫天文奖获得者Hawley教授和原作者Stone教授等著名吸积盘专家也都加入了对该问题的讨论，但一直没有得出明确的结论。

直到2012年，问题描述中提到的中国科学院上海天文台袁峰老师在与导师Narayan教授合写一篇ARA&A（Annual Review of Astronomy and Astrophysics）的约稿时发现，自己老师的模型也许是错误的！

随后袁峰老师等人运用磁流体力学数值模拟方法，对这一问题进行了研究。问题的关键在于，数值模拟中，真正的外流 (风) 到底有多少比例。在外吸积率公式中，只要速度 v_{r}>0，我们就将这部分气体看做外流。在吸积流中，真正的外流和对流元包中向外运动的成分同时存在，因此，外吸积流公式计算的外流率既包含真正的外流也包含对流成分。他随后计算和比较了外流和内流的各种物理性质，如温度、径向速度、转动速度、熵等。他发现，外流和内流的性质差别巨大。如果内吸积流公式和外吸积流公式给出的内外流率是对流元包对流运动所主导，可以预期，内流和外流的性质差别不大。随后他还研究了含磁场的吸积流的对流稳定性，发现吸积流是对流稳定的，也就是说CDAF模型的基础假设存在问题。基于这两方面的原因，袁峰老师得出结论，Stone给出的外流率中，真正的外流，即风，应该是主导的。吸积率随半径的变化原因在于外流而不是对流，即ADIOS模型是正确的。这一研究解决了黑洞吸积的一个基本问题，结束了在该问题上的争论。

同时，幸运的是，2013年，国际上一个由60余人组成的国际合作研究团队一起申请到了针对银河系中心黑洞吸积流的钱德拉X射线太空望远镜的长达300万秒的观测时间。这次国际合作成果发表在了Science上，主要成果就是证明热吸积流的确存在外流，证实了袁峰老师上一年的理论研究结果。

吾爱吾师，吾更爱真理。

参考文献：

[1]. Note：关于黑洞吸积率的数值模拟工作

[2]. Stone, J., Pringle, J. & Begelman, M. 1999, MNRAS, 310,1002.

[3]. 袁峰：黑洞周围发生了什么？

[4]. Yuan, F. & Narayan, R. 2014, ARA&A, 52, 529.



来源：知乎 www.zhihu.com
作者：不拿诺奖不改名

【知乎日报】千万用户的选择，做朋友圈里的新鲜事分享大牛。 点击下载

此问题还有 9 个回答，查看全部。
延伸阅读：
目前宇宙中热门的研究对象有哪些？
如何测定太阳系外恒星的质量？
时间：

没有评论:

发表评论

订阅： 博文评论 (Atom)