渦巻銀河シミュレーションにおけるPython 2D Jean不安定性

渦巻銀河シミュレーションの約1000のスナップショットがあります（分子雲の摂動t = 0からt> 1Gyrまで。スナップショット時間のいずれかの位置がジーンズの不安定性の基準を満たすかどうかを判断することです。

シミュレーションは2Dであり、スナップショットスクリプトは位置、速度、表面密度readを読み取ることができますが、温度は読み取ることができません。シミュレーションにはガスも含まれていません。代わりに、シミュレーションは、個々の星にガスに対応する追加の質量を与えることによってモデル化されました。私はオンラインで見たところ、いくつかのサイトではジーンズの質量が2dでc / G * ∑に等しいと述べています。しかし、音の速さcを決定する方法がわかりません。また、私はジーンズの質量にはあまり興味がありません。そのため、星を作るために任意の位置が崩壊する可能性があるかどうか、つまり、ジーンズの不安定性を満たす必要があります。$_{s}$$^4$$^2$$_{s}$

最初に警告、次にハッキング：

1）このシミュレーションを実行しなかった場合、学部レベル以上の訓練された理論的な天文学者ではない場合、およびこのシミュレーションに必要な種類のデータがない場合は、信頼してはいけないという質問に答える必要があります。まったく結果。

奇妙な数値効果（設定された方法が原因でシミュレーションで発生するものであり、実際には発生しない）、背景の物理学、およびこの種のものに当てはまる仮定と近似の負荷を理解するには、何年もの研究が必要です。

私たちはしばしば「ゴミを入れて、ゴミを出して」と言います。これから得られる答えは、実際の空間よりもシミュレーションについて多くを語っているので、おそらくあまり目立たないでしょう。

2）とはいえ、これを単に楽しみのために、またはスキルを習得するために実行している場合は、すばらしいです！再生し続けます！また、いくつかの数値を取得するために行うことができる非常に基本的な仮定（かなり間違っています）は、空間内のすべてのガスが約10Kであるということです。