渦状銀河の数値モデルに現れるWiggle instability の原因を探るため、強い衝撃波を伴う流れの数値解について詳細に調べた。中でも体積粘性の役割について検討した。薄い衝撃波遷移領域では体積粘性が圧力テンソルの主要項を生み出す。体積粘性は極めて小さく、遷移領域の外では無視できるが、解像度が有限であるために数値解では遷移領域が人為的に広がる。体積粘性による圧力を適正に捉えるためには、粘性の量を人為的に増さなればならない。人為的に増やすと、解像度を高くした極限で古典的な跳躍をふくむ解に収束する。また数値格子の回転やガリレイ変換に対して不変である望ましい解が得られる
|