《含激波的黑洞吸积流的稳定性分析课件.pptx》由会员分享，可在线阅读，更多相关《含激波的黑洞吸积流的稳定性分析课件.pptx（19页珍藏版）》请在三一办公上搜索。
1、1,主 要 内 容,研究背景数值模拟 Molteni et al. (1999)线性化数值计算 Gu & Lu, MNRAS (2006),1主 要 内 容研究背景,2,Mach,1,0,黑洞,激波,含激波的黑洞吸积,声速点,rs,rsh,2Mach10黑洞激波含激波的黑洞吸积声速点rsrsh,3,？,激波对轴对称扰动的稳定性,两种方法得出相同的结论：内激波不稳定，外激波稳定。1. 考虑轴对称扰动下 的计算和模拟2. 简单地分析激波前后 的压强关系,3？激波对轴对称扰动的稳定性两种方法得出相同的结论：,4,激波对非轴对称扰动的稳定性 ？,数值模拟绝热Molteni et al. (1999),
2、线性化数值计算等温Gu & Foglizzo (2003),线性化数值计算绝热Gu & Lu (2006),4激波对非轴对称扰动的稳定性 ？数值模拟线性化数值计算线性化,5,Molteni et al. (1999),5Molteni et al. (1999),6,Molteni et al. (1999),6Molteni et al. (1999),7,数值模拟的结论,轴对称的无粘滞绝热吸积流的外激波对于非轴对称的扰动一般是不稳定的，但是激波不会被扰动彻底破坏，最终形成一个扭曲的，非轴对称的，绕中心黑洞一直转动的激波面。 应用： 解释准周期振荡(QPO) 数值模拟发现激波面的转动周期P在
3、一个范围内，60P300 (rg/c 10-5M/M) 取P=100，则对于10和106个太阳质量的黑洞，振荡周期分别为0.01s和1000s,7数值模拟的结论 轴对称的无粘滞绝热吸积流的外激波对于,8,数值模拟现象的物理解释？,为什么一般情况下都是不稳定的，而对极少数情况是稳定的 ？为什么吸积流参数变化不大的情况下，数值模拟的最终结果有可能相差很大 ?该不稳定性的物理机制 ?,8数值模拟现象的物理解释？为什么一般情况下都是不稳定的，而对,9,我们的工作 对非轴对称扰动的线性计算,吸积流：无粘滞，绝热，有角动量激波：外激波扰动：非轴对称稳定性：整体稳定性分析求解过程： 1. 取扰动量正比于e-
4、it+im，对原含时的偏微分方程组进行线性化分析，得到关于扰动量的常微分方程组 2. 求解声速点和激波面上的边界条件 3. 利用龙格-库塔法数值求解常微分方程组，得出特征频率的值,9我们的工作 对非轴对称扰动的线性计算吸积流：无粘滞，绝,10,10,11,一. 与数值模拟的比较,关于非轴对称扰动下的稳定性，我们的计算结果与数值模拟的结果符合得很好：激波对非轴对称扰动一般都是不稳定的（只有非常弱的激波才可能稳定）。,11一. 与数值模拟的比较关于非轴对称扰动下的稳定性，我们的,12,外激波对轴对称扰动（m=0）的确是稳定的，而对非轴对称扰动是不稳定的。在这个例子中，m=1的扰动增长最快，与该例子
5、在数值模拟中形成m=1的扭曲激波面相符合。,Case 5,12外激波对轴对称扰动（m=0）的确是稳定的，而对非轴对称扰,13,数值模拟的最终状态与线性计算出的具有最快的扰动增长率所对应的m有关。,返回,13数值模拟的最终状态与线性计算出的具有最快的扰动增长率所对,14,二. 不稳定性的物理机制 ？,该不稳定性是基于激波面与共转半径之间的纯声波循环，而与径移声波循环关系不大。,Case 12,共转半径,14二. 不稳定性的物理机制 ？该不稳定性是基于激波面与共转,15,三. 经验公式,15三. 经验公式,16,四. 应用：解释黑洞X射线双星的高频QPO现象,McClintock & Remill
6、ard, astro-ph/0306213,16 四. 应用：解释黑洞X射线双星的高频QPO现象McC,17,结 论,外激波对非轴对称扰动一般是不稳定的我们的线性化数值计算结果能很好地解释 之前的数值模拟现象这种不稳定性的物理机制是基于激波面与共转半径之间的声波循环，其实质是Papaloizou-Pringle 不稳定性的一种具体表现这种不稳定性可能与双星X射线源的准周期振荡 有关,17结 论外激波对非轴对称扰动一般是不稳定的,18,讨 论,激波的存在并不是发生这种现象的必要条件。当吸积流中存在一个相对比较陡峭的转变面时，如标准薄盘到径移主导吸积流（ADAF）的转变面，也可能出现类似的不稳定性。因此我们认为由这种不稳定性引起的非轴对称的转动面应该是有角动量吸积系统的一个常见现象。,18讨 论激波的存在并不是发生这种现象的必要条件。当吸积流,19,Thanks !,19Thanks !,
