降着円盤磁気流体シミュレータの開発(2) 輻射磁気流体ソルバー

Slides:



Advertisements
Similar presentations
エリスワームホール時空における ダスト流解とそのシャドウ Yamaguchi University Takayuki Ohgami, Nobuyuki Sakai ブラックホール地平面勉強会 10 月 4,5 日 湯田温泉.
Advertisements

相対論的場の理論における 散逸モードの微視的同定 斎藤陽平( KEK ) 共同研究者:藤井宏次、板倉数記、森松治.
宇宙ジェット形成シミュレー ションの 可視化 宇宙物理学研究室 木村佳史 03S2015Z. 発表の流れ 1. 本研究の概要・目的・動機 2. モデルの仮定・設定と基礎方程式 3. シンクロトロン放射 1. 放射係数 2. 吸収係数 4. 輻射輸送方程式 5. 結果 6. まとめと今後の発展.
磁気流体降着円盤から噴出する磁気タワージェットの数値実験
Computational Fluid Dynamics(CFD) 岡永 博夫
佐野孝好 (大阪大学レーザーエネルギー学研究センター)
自己重力多体系の 1次元シミュレーション 物理学科4年 宇宙物理学研究室  丸山典宏.
宇宙年齢10億年以前におけるSMBHの存在 遠方宇宙の観測で宇宙10億歳(z~6)未満で10億M⦿程度以上の活動銀河核中のSMBHの存在を確認 赤方偏移 z SMBH質量 [M⦿] URAS J ~2×109 M⦿ 宇宙7.5億歳(z~7)
(Fri) Astrophysics Laboratory MATSUO Kei
第5回 黒体放射とその応用 東京大学教養学部前期課程 2013年冬学期 宇宙科学II 松原英雄(JAXA宇宙研)
観測的宇宙物理学を講ずる。特に,基礎方程式の導出とその解法,観測との比較について詳細な講義を行う。 ***前半の講義内容 ***
第5回 黒体放射とその応用 東京大学教養学部前期課程 2012年冬学期 宇宙科学II 松原英雄(JAXA宇宙研)
ブラックホール磁気圏研究会 ( 大阪市立大学)
加藤真理子1、藤本正樹2、井田茂1 1) 東京工業大学 2) JAXA/ISAS
渦位(Potential Vorticity)と角運動量の保存
スパッタ製膜における 膜厚分布の圧力依存性
原始惑星系円盤の形成と進化の理論 1. 導入:円盤の形成と進化とは? 2. 自己重力円盤の進化 3. 円盤内での固体物質の輸送
観測的宇宙物理学を講ずる。特に,基礎方程式の導出とその解法,観測との比較について詳細な講義を行う。 ***前半の講義内容 ***
相対論的輻射流体力学における 速度依存変動エディントン因子 Velocity-Dependent Eddington Factor in Relativistic Photohydrodynamics 福江 純@大阪教育大学.
SHIGENOBU HIROSE AND JULIAN H. KROLIK
山崎祐司(神戸大) 粒子の物質中でのふるまい.
輻射優勢円盤のMHD数値実験 千葉大学宇宙物理学研究室 M2 松尾 圭 Thu.
ブラックホール降着・噴出流の 輻射磁気流体シミュレーション 大須賀 健 (国立天文台) -新型ジェットと母銀河への影響- 輻射加速 磁気収束
HLLD法に基づく 磁気流体方程式の差分解法
反応性流体力学特論  -燃焼流れの力学- 燃焼の流体力学 4/22,13 燃焼の熱力学 5/13.
日本物理学会年次大会・総合パネル討論「現代プラズマ科学の 最前線:学際連携によるプラズマ理工学のさらなる展開」
降着円盤からの 相対論的輻射流体風 Relativistic Radiation Hydrodynamical Winds from Accretion Disks with Velocity-Dependent Eddington Factor Plane-Parallel Case 福江 純、秋月千鶴@大阪教育大学.
ブラックホール降着流の大局的 3次元磁気流体数値実験
輻射圧駆動風の臨界点ついて 非相対論的領域 Radiatively Driven Spherical Wind Critical Points and Curves Nonrelativistic Regime 福江 純@大阪教育大学.
数値相対論の展望        柴田 大 (東大総合文化:1月から京大基研).
原子核物理学 第4講 原子核の液滴模型.
天文・宇宙分野2 松元亮治 「次世代スーパーコンピュータでせまる物質と宇宙の起源と構造」
速度依存変動エディントン因子を用いた 相対論的鉛直輻射流 Velocity-Dependent Eddington Factor in Relativistic Photohydrodynamics Plane-Parallel Case 福江 純、秋月千鶴@大阪教育大学.
速度勾配依存 変動エディントン因子 Velocity-Gradient-Dependent Relativistic Variable Eddington Factor Plane-Parallel Case 福江 純@大阪教育大学.
1次元電子系の有効フェルミオン模型と GWG法の発展
黒体輻射とプランクの輻射式 1. プランクの輻射式  2. エネルギー量子 プランクの定数(作用量子)h 3. 光量子 4. 固体の比熱.
ブラックホール降着流・噴出流による フィードバック過程の 磁気流体シミュレーション
明るいハードステートに対応する 光学的に薄い降着円盤モデル
ブラックホール周辺の 磁場構造について 大阪市立大学 孝森 洋介 共同研究者 石原秀樹,木村匡志,中尾憲一(阪市大),柳哲文(京大基研)
磁気リコネクション (Craig-Henton解の安定性) ~シミュレーションサマースクール@千葉大より~
YT2003 論文紹介 荻原弘尭.
流体の粘性項を 気体分子運動論の助けを借りて、 直感的に理解する方法
基礎宇宙物理学 II 電磁流体力学入門 第1回 天体活動現象入門 2011年4月8日.
燃焼の流体力学 4/22 燃焼の熱力学 5/13 燃焼流れの数値解析 5/22
磁気回転不安定性によるブラックホール降着流の角運動量輸送機構の 解明
Taniguchi Lab. meeting 2004/10/15 Shigefumi TOKUDA
ON THE VERTICAL STRUCTURE OF RADIATION-DOMINATED ACCRETION DISKS N. J
量子力学の復習(水素原子の波動関数) 光の吸収と放出(ラビ振動)
LHC計画が目指す物理とは × 1:ヒッグス粒子の発見 2:標準理論を越える新しい物理の発見 未発見!
渦位(Potential Vorticity)と角運動量の保存
応用課題 8.太陽風磁気圏相互作用 ベクトル化とベクトル並列化(MPI)の3次元グローバルMHDコードを用いて、SUNワークステーションとベクトル並列型のスーパーコンピュータ Fujitsu VPP5000で太陽風と地球磁気圏相互作用のシミュレーションを行い、惑星間磁場(IMF)が北向きと南向きの場合の磁気圏構造を調べる。図形処理として、PostScript言語を用いた断面図や3次元磁力線の描画、VRMLを用いた3次元可視化を実行する。
連続体とは 連続体(continuum) 密度*が連続関数として定義できる場合
新潟大学集中講義 ープラズマ物理学特論ー (天体電磁流体力学入門) 2004年1月19日ー1月21日
Iida, O., N. Kasagi and Y. Nagano
原子核の殻構造の相対論的記述 n n n σ ω n σ ω n 柴田研究室 石倉 徹也 1.Introduction n n
速度依存変動エディントン因子を用いた 相対論的鉛直&球対称輻射流 Velocity-Dependent Eddington Factor in Relativistic Photohydrodynamics 福江 純、秋月千鶴@大阪教育大学.
シミュレーションサマースクール課題 降着円盤とジェット
ブラックホール近傍からの高エネルギー輻射について
九州大学 猿渡元彬 共同研究者 橋本正章 (九州大学)、江里口良治(東京大学)、固武慶 (国立 天文台)、山田章一(早稲田理工)
大阪市立大学 宇宙物理(重力)研究室 D2 孝森 洋介
インフレーション宇宙における 大域的磁場の生成
河川工学 -洪水流(洪水波の伝播)- 昼間コース 選択一群 2単位 朝位
Massive Gravityの基礎と宇宙論
惑星と太陽風 の相互作用 惑星物理学研究室 4年 深田 佳成 The Interaction of The Solar
熱量 Q:熱量 [ cal ] or [J] m:質量 [g] or [kg] c:比熱 [cal/(g・K)] or [J/(kg・K)]
川島 朋尚 (国立天文台)、朝比奈 雄太 (国立天文台)、工藤祐己 (千葉大) supervised by 松本 洋介 (千葉大)
大型ヘリカル装置における実座標を用いた 粒子軌道追跡モンテカルロコードの開発
How shall we do “Numerical Simulation”?
Massive Gravityの基礎と宇宙論
(昨年度のオープンコースウェア) 10/17 組み合わせと確率 10/24 確率変数と確率分布 10/31 代表的な確率分布
Presentation transcript:

降着円盤磁気流体シミュレータの開発(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 初期条件: 物質は一様 輻射エネルギーのみ原点付近で大 初期条件: 物質は一様 輻射エネルギーのみ左側の領域で大