双曲型方程式-千叶大学.ppt

* * * * * * * * * * * * * * * * * * * * * * * * * . * . * * * * 近似リーマン解法の多次元化 近似リーマン解法の多次元化 多次元の特性の理論に基づく多次元解法 Euler方程式でも容易ではない MHD方程式では想像を絶する 磁場による波動の指向性 磁場のソレノイダル性 近似リーマン解法の多次元化 近似リーマン解法の多次元化 多次元の特性の理論に基づく多次元解法 Euler方程式でも容易ではない MHD方程式では想像を絶する 磁場による波動の指向性 磁場のソレノイダル性 1次元数値解法の利用 Split法 Unsplit法 数値的な磁場発散の生成 近似リーマン解法の多次元化 数値的な磁場発散の影響 非物理的な磁気力が解全体に影響 数値的な磁場発散の処理は必須！ （補正なし） （補正あり） 磁場発散の数値処理 プロジェクション法 ソレノイダルベクトル場への射影 　　課題：連立一次方程式の計算コスト 移流拡散法 数値的な磁場発散の移流、拡散 　　課題：磁場発散の停留、蓄積 Constrained-Transport（CT）法 ソレノイダル条件を維持する離散化 　　課題：高安定化、高次精度化 磁場発散の数値処理 プロジェクション法 [Brackbill, Barnes, 1980] ソレノイダル条件を満足する最小補正ベクトル場 各ステップの計算後に連立一次方程式の計算 チェッカーボード現象 磁場発散の数値処理 境界プロジェクション法 [Miyoshi, Kusano, 2011] 各ステップの計算前に連立一次方程式の計算 数値流束の段階で非物理的磁気力を排除 磁場発散の数値処理 移流拡散法（8-wave法） [Powell, 1994] 磁場発散はエントロピー波で移流 非保存型解法 流れのよどみ点での磁場発散の蓄積 拡散項は実効的ではない 磁場発散の数値処理 移流拡散法（9-wave法） [Dedner+, 2002] 磁場発散は追加された固有値で等方的に移流 固有値は流れと直接的には無関係 保存型解法 非保存（Powell型のソース項）への拡張も可能 ：波動方程式 ：拡散方程式 磁場発散の数値処理 CT法 [Evans, Hawley, 1988] Field-CT法 Flux-CT法 [Balsara, Spricer, 1999] 磁場発散の数値処理 CT法 HLL-Flux-CT法 [Miyoshi, Kusano, 2011] 1次元近似リーマン解法とコンシステント 電場の微分をHLL数値流束で評価 磁場発散の数値処理 CT法 HLL-Flux-CT法 [Miyoshi, Kusano, 2011] 1次元近似リーマン解法とコンシステント 　[Gardiner, Stone, 2005] 電場の微分をHLL数値流束で評価 磁場発散の数値処理 数値実験 HLLD近似リーマン解法 2次MUSCL+minmod制限関数 2次Runge-Kutta-TVD法 ヤコビ法（連立一次方程式） Orszag-Tang渦問題 Field loop移流問題 爆発風問題 Field loop |B|2 Orszag-Tang渦問題 磁場発散の数値処理 Projection Face-projection 8-wave 9-wave flux-CT HLL-flux-CT Field loop問題 Field loop移流問題 Projection 8-wave 9-wave flux-CT v0=0 v0=0 v0=0 v0=0 v0=0 v0=0 |B|2 |B|2 |B|2 |B|2 |B|2 |B|2 |B|2 |B|2 |B|2 |B|2 |B|2 |B|2 Face-projection HLL-flux-CT 爆発風問題 磁場発散の数値処理 Projection 8-wave 9-wave flux-CT P P P P P P Face-projection HLL-flux-CT y=0.3 y=0.3 y=0.3 y=0.3 y=0.3 y=0.3 磁場発散の数値処理 連立一次方程式の処理の手抜き エネルギー補正: literation= 10 literation = 1000 flux-CT + E-fix Face-projection HLL-flux-CT + E-fix literation = 1000 literation= 10 横方向速度（接線速度）を一定と仮定 HLLD近似リーマン解法の展開 横方向速度（接線速度）を一定と仮定