なぜ地球磁極は逆転するのか?

太陽黒点数/オゾン全数/エルニーニョ/太陽活動と米国日本の地磁気変動を追います!

5月度その20:太陽磁極の逆転を追う ➡ MHDシミュレータ ➡ 一次元リーマン問題を変形してみる!

太陽磁極の逆転を追う ➡ MHDシミュレータ ➡ 一次元リーマン問題を変形してみる!

 

 という訳で、そろそろ応用段階に入ろうか?と。 そこでまずビギナーとしては、一次元リーマン問題を変形してみよう。

即ち、オリジナルはステップ応答であるが、これをパルス応答にしてみる。 私としては、まずはこの辺りから踏み出さないと、fortranを理解するという点も含めて。

以下、パルス応答させたリーマン解の gif で、Pressure(圧力)、Density(密度)、By(磁束のY成分)です。

 

 以下、圧力 Pressure:

f:id:yoshihide-sugiura:20190527044024g:plain

 以下、密度 Density:

f:id:yoshihide-sugiura:20190527044049g:plain

 以下、磁束Y成分 By:

f:id:yoshihide-sugiura:20190527044108g:plain 

 これで正しいのかどうか?が、スグには分からない。 ステップ応答の結果と見比べてみよう。 

MHDでは、自分の結果が正しいのかどうかの確認検証をどうするのか?が常につきまとう、と思われます。

 

 

最後までお読み頂きまして、誠にありがとう御座いました。