1、1主 要 内 容 研究背景 数值模拟 Molteni et al.(1999)线性化数值计算 Gu&Lu,MNRAS(2006)2Mach10黑洞激波含激波的黑洞吸积含激波的黑洞吸积声速点rsrsh3?激波对轴对称扰动的稳定性激波对轴对称扰动的稳定性两种方法得出相同的结论:内激波不稳定,外激波稳定。内激波不稳定,外激波稳定。1.考虑轴对称扰动下 的计算和模拟2.简单地分析激波前后 的压强关系4激波对激波对非非轴对称扰动的稳定性轴对称扰动的稳定性?数值模拟绝热Molteni et al.(1999)线性化数值计算等温Gu&Foglizzo(2003)线性化数值计算绝热Gu&Lu(2006)5Mo
2、lteni et al.(1999)6Molteni et al.(1999)7数值模拟的结论 轴对称的无粘滞绝热吸积流的外激波对于非轴对称的扰动一般是不稳定的,但是激波不会被扰动彻底破坏,最终形成一个扭曲的,非轴对称的,绕中心黑洞一直转动的激波面。应用:应用:解释准周期振荡(QPO)数值模拟发现激波面的转动周期P在一个范围内,60P300 (rg/c 10-5M/M)取P=100,则对于10和106个太阳质量的黑洞,振荡周期分别为0.01s和1000s8数值模拟现象的物理解释?为什么一般情况下都是不稳定的,而对极少数情况是稳定的?为什么吸积流参数变化不大的情况下,数值模拟的最终结果有可能相差
3、很大?该不稳定性的物理机制?91011关于非轴对称扰动下的稳定性,我们的计算结果与数值模拟的结果符合得很好:激波对非轴对称扰动激波对非轴对称扰动一般都是不稳定的一般都是不稳定的(只有非常弱的激波(只有非常弱的激波才可能稳定)。才可能稳定)。12外激波对轴对称扰动(m=0)的确是稳定的,而对非轴对称扰动是不稳定的。在这个例子中,m=1的扰动增长最快,与该例子在数值模拟中形成m=1的扭曲激波面相符合。Case 513数值模拟的最终状态与线性计算出的具有最快的扰动增长率所对应的m有关。返回返回14该不稳定性是基于激波面与共转半径之间的纯声波循环,而与径移声波循环关系不大。Case 12共转半径151
4、6 四.应用:解释黑洞X射线双星的高频QPO现象McClintock&Remillard,astro-ph/030621317结 论 外激波对非轴对称扰动一般是不稳定的 我们的线性化数值计算结果能很好地解释 之前的数值模拟现象 这种不稳定性的物理机制是基于激波面与共转半径之间的声波循环,其实质是Papaloizou-Pringle 不稳定性的一种具体表现 这种不稳定性可能与双星X射线源的准周期振荡 有关18讨 论 激波的存在并不是发生这种现象的必要条件。当吸积流中存在一个相对比较陡峭的转变面时,如标准薄盘到径移主导吸积流(ADAF)的转变面,也可能出现类似的不稳定性。因此我们认为由这种不稳定性引起的非轴对称的转动面应该是有角动量吸积系统的一个常见现象。19