|
われわれはここ数年間、星風降着流の数値シミュレーションをさまざまな側面から研究してきた。星風降着流は天体に一様な風が吹き付けるとして近似する。まず流れが軸対称の場合、つまり2次元の計算を行った。つぎに平面対称、ないしは薄板近似の2次元計算を行い、フリップ・フロップ不安定を見つけた。本研究の目的のひとつは、その不安定が3次元の場合にも存在するかどうかを確認することである。ガスが等温的と仮定した場合の3次元計算を行い、ながれは薄板近似の場合ほどは不安定でないことを見つけた。ということは、フリップ・フロップ不安定は計算が薄板2次元近似の場合のみ、起こるらしいことを示唆している。
これらの計算はすべて、空間に格子を配置して離散化する有限体積法を用いてきた。天体物理学で多用されるSPH法とよばれる粒子法がある。われわれはSPH法の研究を行い、それが衝撃波を捉えるには不十分な方法であることを確認した。そこでその欠点を回避するために、SPH法にリーマン問題の手法を併用する犬塚の手法を研究した。これはまだテスト段階で、それが有用かどうかの結論を得るにいたっていない。
つぎにわれわれは、軸対称の星風降着流に磁場の効果をいれる研究を行っている。計算手法はMHDのRoe法である。この研究もまだとりかかったばかりで、完全な成果は得ていない。
|
