流体力学理论天文学宇宙物理学恳谈会.ppt

流体力学（格子法）宇宙気体力学への応用 松田卓也 神戸大学?理学部?地球惑星科学科 目次 宇宙気体力学 数値計算法の歴史と現状 近接連星系内のガスの流れ 降着円盤 恒星気象学 宇宙気体力学 宇宙はほとんどガス(気体) 中性ガス プラズマ???磁場 時間依存性 静的 恒星の内部構造論 定常 降着円盤、星風 動的 爆発 近接連星系内のガス流 降着円盤 解析的モデル アルファ円盤モデル ADAF 数値シミュレーション 星風降着 解析的モデル Bondi-Hoyle-Lyttletonの理論 数値シミュレーション 対象天体 降着円盤 激変星 Cataclysmic variables 新星 Novae 小質量X線星 Low mass X-ray binaries 星風降着 共生星 Symbiotic stars 大質量X線星 High mass X-ray binaries 激変星の想像図 基礎方程式系 Boltzmann方程式 気体運動論 Chapman-Enskog展開 Navier-Stokes方程式 粘性あり Euler方程式 粘性なし 衝撃波、接触不連続面あり Euler方程式 Euler方程式の性質 非粘性の圧縮性流体 双曲型方程式 楕円型方程式 Poisson方程式 放物型方程式 熱伝導方程式 時間依存 非線形 衝撃波、接触不連続面などの不連続 保存系か非保存系か 保存系 Rankine-Hugoniot条件をみたす 衝撃波の位置が正確 非保存系 保存量の保存性を調べる 陽的解法と陰的解法 陽的解法 右辺は既知量のみ アルゴリズムが簡単 Courant条件で時間刻みが限定される 陰的解法 右辺にも未知量 行列の反転 Courant条件なし しかし時間刻みの制限はある 数値計算法の歴史1 Von Neumann Richtmayer(1950) 人工粘性 Lax-Friedrichs法(1954) 非常になめらか Lax-Wendroff法(1964) 2次精度 Taylor展開に基づく方法 MacCormack法(1969) Beam-Warming法(80年代) 数値計算法の歴史2 風上差分法 Courant-Isaacson-Rees法(1952) 1次精度 Godunovの定理(1959) 次の2つしか満足しない 線形性 単調性 1次精度以上 Euler方程式は非線形 数値計算法の歴史3 天文業界はガラパゴス島であった 進化に取り残されて、独自に発展 今はグローバルスタンダードの波が 天文業界3大手法 Total Variation Diminishing(TVD) PPM SPH 私のアルゴリズム遍歴 1次元ラグランジュ法+von Neuman人工粘性 Lax-Wendroff法 オイラー+ラグランジュ混合法 Fluid in Cell法(1次精度) MacCormack法 Beam-Warming法(隠的法) Osher法(風上差分、1次、2次?隠的、陽的) Van Leer Flux-Splitting(MHD) SFS法(嶋、城之内) 風上差分法 セル境界でRiemann問題を解く Riemann問題の解法 逐次近似で厳密に Godunov法,PLM(2次近似), PPM(3次近似) 近似解法 Osher Solomon(1982) Harten, Lax, van Leer(1983) Roe(1986) 線形Riemann問題 風上差分法 セル境界でRiemann問題を解く Riemann問題の解法 逐次近似で厳密に Godunov法,PLM(2次近似), PPM(3次近似) 近似解法 Osher Solomon(1982) Harten, Lax, van Leer(1983) Roe(1986) Shock tube 問題 Shock tube　問題の解 1次精度と2次精度 Van Leer と　Roe 近接連星系内のガスの流れ 数値実験の歴史 Prendergast (1960) Prendergast Taam (1974) ビーム法 Sorensen, Matsuda Sakurai (1974, 1975) Fluid in Cell Method (FLIC法) Lin Pringle (1976) 粘着粒子法 降着円盤上の渦状衝撃波 数値実験で発見さる(Sorensen et al.から11年後) Sawada, Matsuda Hachisu, M.N. 219 (1986a) Sawada, Matsuda Hachisu, M.N. 221 (1986b) Sawada,