2005世界物理年-quasar.cc.osaka.ppt

輻射輸送問題三題噺～最近の進展～１　輻射場の非等方性について（復習）２　速度依存変動エディントン因子を用いた球対称一般相対論的輻射流３　有限の厚みをもつ平行平板における非相対論的輻射輸送の厳密解 １　輻射場の非等方性について（復習） ０　準備 　　輻射輸送方程式 １　静的平行平板大気 　　周縁減光効果とピーキング効果 ２　静的球対称大気 　　膨張大気で先鋭化するピーキング ３　球対称風 　　高速流における希薄化　 ０．準備　輻射の流れ 平均自由行程進むと、物質粒子によって吸収や散乱を受ける。 ０．準備　輻射輸送方程式 原理的には、輻射輸送方程式を解けば、輻射輸送の問題はまぎれなく解ける。 ただし、7つの独立変数（r、l、t、ν）をもった微分積分方程式である。こんなの解きたくない！ 相対論：（座標）静止系／実験室系と（流体）静止系／共動系を区別しなければならない。　 ０．準備　モーメント定式化 輻射の“非等方性”はあまり強くないと仮定し、光線の角度依存性は弱いとして、輻射輸送方程式を角度方向に展開し、角度について積分して、0次のモーメント、1次のモーメント、などと呼ばれる一群の方程式セットを得ることができる。 方程式系を閉じるために別の関係式が必要。 ０．準備　エディントン近似 輻射場が等方的な場合に成り立つ関係： 　　　 P＝E/3　（一般にはPij=δijE/3） 　　このエディントン近似でモーメント式を閉じる。 この関係は常に成り立つとは限らない。 天体の表面近傍など輻射場が光学的に薄くなる領域では、輻射場の非等方性が無視できなくなる。（ピーキング効果） 亜光速に加速される流れで速度勾配が非常に大きい領域でも、共動系でさえ輻射場が非等方的になる。 １　静的平行平板大気状況：無限半平面 １　静的平行平板大気基礎方程式（ｚ） 輻射輸送 0次モーメント 1次モーメント Eddington近似 輻射平衡（R.E.）：j – cκabsE=0 １　静的平行平板大気基礎方程式（τ） 光学的深さτ導入 輻射輸送 0次モーメント 1次モーメント Eddington近似 １　静的平行平板大気Milne-Eddington解 境界条件 モーメント 輻射強度 表面での強度： 　周縁減光効果 １　静的平行平板大気Milne-Eddington解 表面での強度： 　周縁減光効果 １　静的平行平板大気Milne-Eddington解 表面での強度： 　周縁減光効果 １　静的平行平板大気Eddington近似の精度 得られた輻射強度を、表面で再度積分 ２　静的球対称大気状況：希薄大気 ２　静的球対称大気基礎方程式（ｒ） 輻射輸送 0次モーメント 1次モーメント Eddington近似 輻射平衡：J=S ２　静的球対称大気Milne-Eddington解 境界条件 モーメント 源泉関数＝一定 ２　静的球対称大気Milne-Eddington解 ２　静的球対称大気Milne-Eddington解 輻射輸送方程式 　視線s方向 輻射強度 表面での強度 ２　静的球対称大気状況：膨張希薄大気 表面での強度： 　周縁減光効果 τ0＝（κ+σ）ρR ２　静的球対称大気ピーキング効果はより強い ２　静的球対称大気過去の研究 Transfer (Analytical) Kosirev 1934 Chandrasekhar 1934 Chandrasekhar 1945 Chapman 1966 Unno Kondo 1976, 1977 Masaki Unno 1978 Unno 1989 Transfer (Numerical) Hummer Rybicki 1971 Schmid-Burgk 1975 Hundt et al. 1975 ３　球対称風状況：膨張希薄大気 ３　球対称風基礎方程式 球対称流 速度一定 　τはだらだらと減少 ３　球対称風過去の研究 (Relativistic) Radiation Hydrodynamics Lindquist 1966 Castor 1972 Ruggles Bath 1979 Mihalas 1980 Quinn Paczynski 1985 Paczynski Proszynski 1986 Turolla et al. 1986 Paczynski 1990 Nobili et al. 1993 Nobili et al. 1994 King Pound 2003 低速近似（v/cの1次まで） 速度一定を仮定していて、ダイナミックスを解いていない 共動系で拡散近似を使