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の走行結果は:
である。 尚、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は:
であり、シャープ性に欠けてくる。 尚、Bx、Bzは共にゼロ。
* Bxの初期値依存性については、もう少し調べてみないと分からない。
かなかなリーマンから抜け出せないけれど、基本中の基本との事なので、ジックリ理解しないとイケナイと考えております。