日本地球惑星科学連合2022年大会

講演情報

[J] 口頭発表

セッション記号 P (宇宙惑星科学) » P-EM 太陽地球系科学・宇宙電磁気学・宇宙環境

[P-EM16] 宇宙プラズマ理論・シミュレーション

2022年5月22日(日) 10:45 〜 12:15 105 (幕張メッセ国際会議場)

コンビーナ:天野 孝伸(東京大学 地球惑星科学専攻)、コンビーナ:三宅 洋平(神戸大学計算科学教育センター)、梅田 隆行(名古屋大学 宇宙地球環境研究所)、コンビーナ:中村 匡(福井県立大学)、座長:梅田 隆行(名古屋大学 宇宙地球環境研究所)、三好 隆博(広島大学大学院先進理工系科学研究科)

11:30 〜 11:45

[PEM16-10] 高階微分項を用いた陽的FDTD法の数値分散の低減

*関戸 晴宇1,2梅田 隆行2三好 由純2 (1.名古屋大学、2.名大ISEE)


FDTD(Finite-Difference Time-Domain)法は、Yee(1966)によって開発された、電磁界の時間発展を解く数値計算手法であり、空間と時間ともに2次精度の差分でMaxwell方程式を近似することで求められる。不連続波形で数値振動が発生するほか、連続波形でも傾きが大きいと数値振動や振幅の減衰が起こるといった欠点が存在する。これらの欠点を改善するため、Petropoulos(1994)は空間微分項の差分精度を4次精度としたFDTD(2,4)法を提案したが、数値不安定が生じないクーラン数Cの上限値が小さくなるという問題を抱えている。本研究では電界および磁界のそれぞれの時間発展方程式に高階空間微分項を付加することにより、新たな陽的時間発展式を導出した。その結果、従来の手法では数値不安定が生じ、計算することができなかった範囲においても計算が可能となることを確認した。また、従来の手法で計算が可能であった範囲においては、位相速度の誤差を抑えることができた。