宇宙ジェット形成シミュレーションの可視化宇宙物理学研究室木村佳史 03S2015Z 発表の流れ本研究の概要・目的・動機モデルの仮定・設定と基礎方程式シンクロトロン放射 1. 2. 3. 1. 2. 4. 5. 6. 放射係数吸収係数輻射輸送方程式結果まとめと今後の発展本研究の概要・目的  磁気流体シミュレーションの結果をもとにして、宇宙ジェットの可視化をする。  放射強度の計算による実際のジェットの写真のようなイメージを作る。電波銀河（電波と可視光）銀河核からのジェット本研究の動機等密度面のアニメーションシンクロトロン放射で可視化する理由。これまでは密度の図などで可視化されてきたが、実際にそれが見える訳ではない。それを実際に見えるようなイメージにするため、放射光を考える。ジェットの密度分布と磁力線シミュレーションモデルの仮定・設定①    鉛直磁場に貫かれた降着円盤を考える。理想磁気流体（非粘性・磁気拡散なし）を仮定。中心重力源による外場ｇが存在。磁場B 初期状態では円盤の外に高温ガスコロナ（温度一様・回転無し）鉛直方向のみ降着円盤（鉛直方向を軸として回転）中心重力源（例えばブラックホール）シミュレーションモデルの仮定・設定②  円盤密度は重力と遠心力と圧力勾配力のつり合いより決める。 n3 11 n としてシミュレー P  圧力はポリトロープを仮定。ションをした。  高温ガスコロナの密度は静水圧平衡（圧力勾配力と重力のつり合い）より決める。 k B  cTc と仮定。  圧力は pc  m 基礎方程式     質量保存の式（連続の式）運動方程式（外力 g ）     ( v )  0 t  B   B   v      v   v   p   g  4  t  誘導方程式 B     v  B t エネルギー式       v  Pv  0 t ε：内部エネルギー γ：比熱比 5 3   とした。 e 1 2 v  2 p    1  シンクロトロン放射①    光速近くまで加速された電子が磁場中で曲げられるときに電磁波を発生する。 B e 運動方程式 d q   mv   v × B dt c プラズマ（多数の陽子電子）を考えたときの、１つの電子が放射する放射エネルギー（全振動数について積分したもの） 2q 4 B2 2  2 2 P sin  2 3 α：速度と磁場のなす角 3me c 1 v  たくさんの電子があるとしてい  2 1   るので平均化して１／２とした。 c シンクロトロン放射②  電子は熱的な分布をしていると仮定する。  エネルギーＥを持つ電子の個数 N E   KE exp E kBT  2 マクスウェルの速度分布より N ( p)  Kp2 exp   E kBT  相対論的な電子を考えているので   また m  0 N  E  dE  2K  kBT  e 3 K  E  pc  2me  kBT  3 シンクロトロン放射③   放射係数 1 j  4   mc （シンクロトロン自己吸収）キルヒホッフの法則 B T ：黒体輻射強度レイリージーンズ極限 N  E  P  , E  dE 2 吸収係数  2 j c   2 2 k BT  B T   jv h kT 輻射輸送方程式 dI   I  j ds  ｓは光線の通る道筋光線放射のみの場合   0  放射と吸収の両方がある場合の2つの場合で、シンクロトロン放射から得られた放射・吸収係数を用いて、放射強度の図を描く。振動数ν＝１０１０Hz （電波）とした。  可視化  3次元ビジュアリゼーションソフトＡＶＳ（Application Visualization System）を使用。真横から見た時放射のみの場合のジェット放射と吸収を含めたジェット放射のみの場合スクリーンを傾けて見たジェット真横から見たジェット放射と吸収の場合スクリーンを傾けて見たジェット真横から見たジェット温度との比較  温度が高いところは放射が強い。同じ時刻の放射強度と温度の図まとめと今後の発展   今回は放射強度の図を描くことにより、宇宙ジェットの可視化をした。ドップラー効果や屈折や他の吸収要因を加えてみることもできそう。終参考：スクリーンと積分方法スクリーンの原点の座標 r0  r0 , 0 , z0  x  xex  yey  スクリーン上の座標  光線が通る点の座標 r  r,  , z   スクリーンからの距離ｓスクリーンの法線ベクトルｎと平行に積分していく。  r  r, , z   r0  x  sn スクリーンの直交基底はｘ方向の基底のｚ座標が０になるようにとった。