吸积盘的研究进展

1、王建成云南天文台l球对称吸积qBondi（1952）研究定常态球对称吸积，天体对周围气体的影响。q模型给出吸积率、吸积半径、吸积流跨声速等物理量的关系。qParker（1969）等人在Bondi解的基础上研究了球对称星风和吸积的过程，进一步发展了理论。l薄吸积盘q流体具有角动量，吸积过程需要角动量的转移。qShakura&Sunyaev（1973）提出粘性的模型，发展了薄盘理论。q粘滞作用导致角动量沿径向向外转移。q粘滞耗散产生的能量以辐射方式转移出去。q吸积流形成一个几何薄、光学厚的吸积盘。q辐射谱是不同温度黑体谱的叠加。q广泛用于解释高能天体的红外、光学、紫外、X射线的辐射谱。l离

2、子主导吸积盘（Shapiro et al. 1976）q吸积流形成双温等离子体，离子1011K，电子108-109K。q流体是光学薄的、能产生X和波段的非热辐射。q热不稳定l超爱丁顿吸积盘（Kato et al. 1977）q流体是光学厚的，大部分辐射被俘获。q耗散产生的能量被物质内流拖曳进黑洞。q吸积率大、光度小。l薄盘模型的局限性q低光度的天体：低态的X射线双星、低光度活动星系、Sgr A等。q不能解释宽波段的辐射谱lADAF的一维模型（Narayan& Yi 1994）q定常态轴对称吸积流 q动力学性质由四个高度积分的微分方程描述：吸积流的质量、径向动量、角动量和能量的守恒 。2

3、221KSdvdvRRcdRdR 0dR H vdR231dRddvR HdRR H dRdR2222dsddvTqqRqfRdRdRdRl参数f是平流能量与粘性产生的能量的比，f=1为平流主导，则f=0为辐射冷却主导。l参数是运动学粘性系数l假设为一个与半径R无关的量2SSKcc Hl质量守恒方程导致吸积率为常数l以开普勒角速度 和自由落体速度 做自相似假设： l密度轮廓为：22MRHv常数3/KG MR/ffvG MR1/2ffvvR1/2SffcvR3/2KR 3/2/ 4/KSMR H vR v cR l利用自相似和f与R无关的假设，可得盘的解： 252,3ffv Rgv 1/2225

4、2,9KRg 2222 52,9SffcRgv 1 5/31ff 1/22218,1152g lADAF有趣的特征q对于高粘性ADAF （0.2-0.3），径向速度与自由落体速度相当（V 0.1Vff）q流体以低于开普勒的角速度旋转，离心力只起部分支撑作用，剩下的支撑来自于径向压强梯度。当 流体几乎没有旋转( )。 5/30q辐射率低，耗散能量转变为热能，流体温度达维里温度。q盘被热压力顶起，标高H Cs/K R，ADAF在几何形状上更像球吸积。qADAF中的流体具有正的Bernoulli参数 （正的比能），有可能产生喷流或是某种形式的物质外流。q流体的熵随着半径的减小而增加，ADAF是对流不

5、稳定的。 l二维ADAF模型 （Narayan&Yi 1995）qADAF具有准球吸积性质，使用高度积分方程可能是一种过度的简化。q为证明高度积分的有效性，Narayan和Yi（1995）考察了ADAF在球坐标极角方向上的结构。q他们考虑了球坐标下的一个轴对称无子午流( )的ADAF二维结构。 0vq自相似假设 3/2r rKG Mvvrvr Kvr SKScrc q代入流体力学方程，得到四个关于函数 的六阶常微分方程组。q方程组需要六个边界条件才能求解：在赤道面和旋转轴解没有奇异性，满足对称性要求。q利用解常微分方程双边界问题的驰豫方法，求解盘的结构。 ,Svc 和:02Sdcdvd

6、ddddd0:0,0Sdcdvddvdddd q图展示了三个典型解（ ）的角速度、径向速度、密度和声速平方的角向分布轮廓。q对于ADAF盘( )， 在同一半径的球壳层上几乎都是常数，径向速度在旋转轴上为零，在赤道面达到极大值。q解( )的轮廓显示出薄盘解的特征，密度在赤道面上达到最大，随着纬度的增加密度迅速下降，说明物质主要集中在赤道面上，而角速度接近开普勒值。 0.1,0.1,1,10 0.1,1 2,S和c10lADAF的应用 q辐射以ADAF中的高能电子非热辐射为主，主要有同步辐射、轫致辐射和逆康普顿辐射。q辐射谱从射电延展到硬X射线波段。q质子质子碰撞产生的中性介子衰变能产生 辐射。l

7、有物质外流的ADAF（Xu & Chen 1997）q经典ADAF没有子午流动，吸积内流的物质不可能改变方向变成外流物质。q具有子午流动（ ）的情况下，求解自相似的二维ADAF结构，研究物质外流的影响。0vq自相似假设：密度的自相似幂律指数由3/2换成任意参数n，其他物理量假设不变。q穿越整个球面的净吸积率是： 若质量守恒，净吸积率为常数，有两种情况：幂律指数净吸积率为零： 23/22sin,2sinrnrMrrvrdG M rvd 。3/2n sin0rvd q第一种情况对应经典的ADAF，只有物质吸积，没有物质外流。q第二种情况对应修正的ADAF，内流的物质和外流的物质相等，使净吸

8、积率为零。 两种情况都是特例，一种是没有物质外流；一种是吸积物质全部逃逸出来，没有物质落入中心天体。 q模型对 进行傅立叶展开，代入流体力学方程组，将一个以一定的边界条件求解微分方程组的问题转变为一个非线性代数方程组的求根问题。q经过相应的数值计算，他们发现两个类型的解吸积外流和抛射外流。 ,rpvvv 和 左上图是密度等值线和子午面内的速度矢量场；左下图是温度等值线；右图是相应的物理量的角分布图。 吸积外流解处处都是正能量，有逃逸到无穷远的潜力。抛射外流解处处都是负能量，外流的物质最终都能返回来，不产生外流。 l绝热内流外流模型（Adiabatic Inflow-Outflow Soluti

9、ons, ADIOS） （Blandford& Belegman 1999）q假设辐射冷却无效q假设径向速度远小于旋转速度 （只适合小粘滞情况，即0.01）q离子主导的物态方程：q引入三个自相似参量p、5/3物质吸积率满足：角动量内流满足：能量外流满足：G是半径r处内层物质对外层物质的力矩。物质外流带走的角动量和能量： ;01pmrp21/2;0lFm rGm r222115022EamFGmrrr1/21/21;lEprpdFdFdmpdmpr。径向运动方程： Bernoulli常数（物质的比能量）： 利用以上方程组，可求解具有外流的盘结构22215/20rp ar 。2222221

10、51,222raB eparrl三个参数 决定吸积盘的性质q 情况，对应无外流无旋转的Bondi球对称吸积。q 情况，对应无外流但有辐射损失的吸积流，因此只有能量外流而无角动量外流。q 情况，对应磁主导的风，盘上的物质流动是守恒的，所有的角动量和能量被风所带走，在盘中不存在耗散过程，而且盘是冷和薄的。p， 和0p0,1,1/2p0,1,1/2pq 情况，对应纯气体动力学风，即风只带走它自己在出发点的角动量，而不对吸积盘的其它流体产生力矩作用。q 情况，对应Bernoulli常数（比能量）为零的临界束缚吸积流。q 情况，对应 的中间解。 1/2221044/ 21 4815pppppp 232/

11、21/5/2rarpp0.75,0.75,0.5p0.35/Ber l粗线所围四边形为满足(1) ，(2)G0，(3)lwl，(4)Be0四个限制条件的允许区域。 ,0a H q质量守恒、动量守恒和能量守恒方程q自相似假设q物态和粘性假设2221sssKpccecq可得到五个关于函数 的常微分方程。q它们分别对应着质量守恒、三个方向的动量守恒和能量守恒方程。q这五个方程组成的常微分方程组是一个八阶的方程组，必须给定八个边界条件才能求解。 ,rsvvvc 和l边界条件q赤道面边界条件赤道面边界条件（ ） 物理量在赤道面光滑过渡，并且具有镜物理量在赤道面光滑过渡，并且具有镜像对称的性质。像对称的性

12、质。 /2/2/2/2/20rsdvddvdcdddd /20v/21 l外流边界条件q= 0作为外流的边界q在边界上物质流从内流开始转向外流，它们的径向速度是零（ ）q盘的厚度满足条件cs/K(/2- 0)r或 cs /2- 0q模型参量：f=1，=1.43（能均分），另外三个参数 做为自由参数。00rv0, ,n qn表征物质内流的强弱，越大则物质的内流越强，当n=3/2时只有内流没有外流，对应典型的ADAF；当=1/2时对应对流主导吸积流，1/2n3/2。q表征粘性大小的参数，0 1。l三个典型解：1.25,1,0.75n l 与BB99模型的对比 BB99模型是一维垂直方向高度积分模型

13、，新模型是轴对称的二维模型，对模型做方向的积分也可以给出如BB99模型的质量吸积率、角动量内流率和能量外流率。 q质量吸积率（内流为正）：q角动量内流率（内流为正）： 0/23/204sin3/2nrpMrG Mvdrrpn 。 002/22130/22221/21/2002sinsin34sinsin2lrnrsnpF rrv v r T rd dG Mvvvcd drrrM r 。q能量外流率（向外为正）： 002/2211121302/23/222222222sin141213322222cotsin33ErrsrrrsrrrFrrv B eT vT vT vd dcG Mvvvvdvc

14、vvvvddvvvvdd 1/21/210nnrrM r。qBB99模型的参数p对应于我们的3/2-n，而他们的 对应于我们的 q定义单位质量外流物质所带走的角动量和能量:, 1/2/2/3/2lldFdFdrnrdMdMdrn1/1/2/3/2EEdFdFdrnrdMdMdrnl我们的解是不是ADIOS的一个二维扩展呢？ q平均比能对于上述三个典型解我们有： 随着的n减小，平均比能在减小。这说明物质外流的减小将导致能量在吸积气体中的迅速堆积，使这些气体更加接近非束缚的状态。 BB edd 0.0627(1.25)Bn 0.124(1)Bn 0.128(0.75)Bn q物质的外流有将吸积盘中

15、过多的能量带走减小能量堆积的作用，这也是ADIOS（BB99模型）引入物质外流以解决ADAF中能量的过度堆积问题的基本思想。 q我们的解与ADIOS有很多相似之处，但是我们的解却不是ADIOS的一个扩展：ADIOS对应低粘性1(0.01)，VrV我们的解对应高粘性0.3， VrVq我们的解和ADIOS的允许区域不同l热前（Thermal Front） 讨论模型中流体的耗散性质q定义四个物理量在不同纬度上的分布情况： (压强梯度 ) (粘性耗散率) (比动能) (流体的焓) 21npArr 5/2:nBru T 222112rvvvCr 211scDr 子午压强梯度轮廓（实线），粘性耗散率轮廓（

16、虚线），比动能轮廓（点划线），气体比焓轮廓（点线）。 q子午压强梯度的方向在低纬度的地方（大值）是负方向，与 的方向一致，对子午流做正功。因此，推动物质外流的源动力是气体的压强梯度。它在高纬度（ ）情况相反，对流体做负功，起稳定吸积盘的作用。这个区域在耗散性质上十分类似于Blandford & Begeleman 2004（BB04）的BB99模型二维扩展解中假设的“热前（Thermal Front）”区域：运动能量快速耗散为热能。v5765l外流对吸积盘质量、角动量和能量的影响nB/ / 0.751.551.330.37-0.13 0.860.24 1.44-0.081.001.591.320.37-0.13 0.830.23 1.66-0.231.251.421.150.2