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

現在、太陽黒点数の推移と、太陽磁極の逆転メカニズム、を追っています

5月度その24:太陽磁極の逆転を追う ➡ MHDシミュレータ ➡ リーマン問題・Bx初期設定値を振ってみる!

太陽磁極の逆転を追う ➡ MHDシミュレータ ➡ リーマン問題・Bx初期設定値を振ってみる!

 

 という訳で、Bx初期値はどのような働きをしているのか?を知る為に、これを振ってみる。

 

 まず、最初はオリジナル ! Brio & Wu (1988) モデルで:

if ( x(i) < 0.d0 ) then
! Brio & Wu (1988)
ro0 = 1.d0
pr0 = 1.d0
vx0 = 0.d0
vy0 = 0.d0
vz0 = 0.d0
bx0 = 0.75d0 ! 正の固定値を設定
by0 = 1.d0 ! プラス1
bz0 = 0.d0

else ! ゼロ未満で切り替え
! Brio & Wu (1988)
ro0 = 0.125d0
pr0 = 0.1d0
vx0 = 0.d0
vy0 = 0.d0
vz0 = 0.d0
bx0 = 0.75d0 ! 正の固定値を設定
by0 =-1.d0 ! マイナス1で、ステップ入力
bz0 = 0.d0

endif

である。一次元Xの長さは300で中央がXゼロで左右対象。

Byの走行結果は:

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

である。 尚、Bxは0.75固定で、Bzはゼロ。

 

ここでBx = 0.0 としてみると:

if ( x(i) < 0.d0 ) then
! Brio & Wu (1988)
ro0 = 1.d0
pr0 = 1.d0
vx0 = 0.d0
vy0 = 0.d0
vz0 = 0.d0
bx0 = 0.d0 ! 0.75d0
by0 = 1.d0
bz0 = 0.d0

else

! Brio & Wu (1988)
ro0 = 0.125d0
pr0 = 0.1d0
vx0 = 0.d0
vy0 = 0.d0
vz0 = 0.d0
bx0 = 0.d0 ! 0.75d0
by0 =-1.d0
bz0 = 0.d0

endif

得られたByは:

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

 であり、シャープ性に欠けてくる。 尚、Bx、Bzは共にゼロ。

 

 Bxの初期値依存性については、もう少し調べてみないと分からない。

かなかなリーマンから抜け出せないけれど、基本中の基本との事なので、ジックリ理解しないとイケナイと考えております。

 

 

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