• 研究課題をさがす
  • 研究者をさがす
  • KAKENの使い方
  1. 課題ページに戻る

2001 年度 実績報告書

モンテカルロ直接数値シミュレーションによる降着円盤の研究

研究課題

研究課題/領域番号 13640241
研究機関神戸大学

研究代表者

松田 卓也  神戸大学, 理学部, 教授 (20026206)

研究分担者 中川 義次  神戸大学, 理学部, 教授 (30172282)
キーワード降着円盤 / モンテカルロ直接法 / 数値シュミレーション
研究概要

モンテカルロ直接法(Direct Simulation Monte Carlo method=DSMC)を用いて、近接連星系の降着円盤の数値シミュレーションを行った。
降着円盤においては、気体からの輻射熱冷却が重要な役割を果たす。本研究では、それを近似的に扱うために、二種類の冷却法を試みた。1)粒子間の非弾性衝突、2)粒子の熱速度の強制的減少。いずれの方法でも冷却を模擬できることが分かったが、1の方式ではガスがグラニュラーになるので、2の方が適当であることが分かった。
質量比が1:1の近接連星系の場合を計算し、計算領域は質量降着星まわりの空間でL1点を含む3次元領域とした。L1点から気体粒子をタイムステップごとに一定数、温度ゼロで放出して、進化を追跡した。その結果、数公転周期の後に薄い降着円盤が得られた。放出粒子数の数により、降着円盤の形が変わることが分かった。放出粒子数が少ない場合は、降着円盤はほぼ円形をなしている。しかし放出粒子数が多いと、L1流は降着円盤に接するような形になる。つまり円盤の下限がL1流で制限された形になる。この結果は実は藤原たちによる有限体積法の結果と異なっている。その場合は降着円盤の縁が幾何学的に厚く、L1流は降着円盤に貫入するかたちになっている。DSMC計算では、冷却効果が十分に効いているようで、円盤が薄く、L1流は貫入しない。どちらが正しいかは、輻射冷却をきちんと考慮した計算をするまでは分からない。
もうひとつの相違点は、DSMCでは有限体積法で見られたような顕著な渦状衝撃波が見られなかった点である。その原因は、DSMCでは降着円盤の周辺部で粒子数が少なくなり、流れの連続性を保証できないからであろう。この点は改良を必要とする。

  • 研究成果

    (4件)

すべて その他

すべて 文献書誌 (4件)

  • [文献書誌] 林 英二: "Is angular momentum in an accretion disk transported inwards?"Progress of Theoretical Physics. 105. 531-536 (2001)

  • [文献書誌] 藤原秀和: "Numerical Simulation of interaction between an L1 stream and an accretion disk in a close binary system"Progress of Theoretical Physics. 106. 729-749 (2001)

  • [文献書誌] 松田卓也: "Direct Monte Carlo simulation of accretion disc"Proceedings of "The Physics of Cataclismic Variables and Related Objects 2001". (2002)

  • [文献書誌] 松田卓也: "Application of DSMC method to astrophysical flows"Proceedings of IAU Symposium 208 "Astrophysical Supercomputing using Particle Simulations". (2002)

URL: 

公開日: 2003-04-03   更新日: 2016-04-21  

サービス概要 検索マニュアル よくある質問 お知らせ 利用規程 科研費による研究の帰属

Powered by NII kakenhi