降着円盤磁気流体シミュレータの開発(2) 輻射磁気流体ソルバー J02b 小田寛(千葉大)町田真美(国立天文台)、中村賢仁(松江高専) 小川崇之、川島朋尚、松元亮治(千葉大) 降着円盤磁気流体シミュレータ 磁気流体 HLLD + CT法 (小川;J01a) 輻射磁気流体 FLD近似 (本講演) (相対論的磁気流体) 並列化:スカラー計算機用に最適化(松元;J03c)
磁場も輻射輸送も両方重要 輻射輸送を考慮せず 磁場を考慮せず J02b 輻射圧、光子捕獲 など輻射輸送が重要 光学的に薄い円盤 光学的に厚い円盤 角運動量輸送や 円盤加熱には 磁場が重要 磁気圧優勢な 熱平衡状態も存在 磁場も輻射輸送も両方重要 輻射輸送を考慮せず 磁場を考慮せず
計算方法 基礎方程式 J02b 残りの項は陽解法(修正Lax-Wendroff法) FLD近似 & Variable Eddington法 質量保存 運動量保存 ガスの エネルギー式 輻射の エネルギー式 誘導方程式 FLD近似 & Variable Eddington法 計算方法 ガスと輻射の相互作用の項 のみ陰解法 (e.g., Turner & Stone ‘01) 残りの項は陽解法(修正Lax-Wendroff法)
一次元テスト計算結果 光学的に厚い物質中の 輻射の拡散 光学的に薄い物質中の Radiation frontの伝播 J02b 輻射エネルギー x x 初期条件: 物質は一様 輻射エネルギーのみ原点付近で大 初期条件: 物質は一様 輻射エネルギーのみ左側の領域で大