現在私はN体問題を解いて銀河のシミュレーションをしようとしているのですが、肝心の銀河の初期値を決め方で悩んでいます。
色々検索したところこちらのページに様々なモデルのポテンシャルと密度分布と重力と釣り合う速度が書かれているのを見つけたのでこれを参考にしようと思いました。
>http://www.astro.utu.fi/~cflynn/galdyn/lecture4. …
しかし粒子がちょうど系を飛び出さない速度は1/2mv^2+Φ(r)=0になるはずなので計算できるのですが、肝心の粒子の位置を決める方法が全くわかりません。
一様球の場合直感的に想像できるので簡単に分かったのですが、それ以外の密度に偏りのあるモデルのになると全く分からなくなります。
Plummer modelの計算をしていると思われるコードを見つけたのですがどうしてこうなるのか分かりません。正しいのかどうかも確かめようもありません。
>http://www.artcompsci.org/kali/vol/plummer/ch02. …
(二つめの緑枠内の”radius = 1.0 / sqrt( rand ** (-2.0/3.0) - 1.0)”の行です)
前置きが長くなりましたが、上記のページに載っているような銀河モデルの初期位置を決める式の導出の方法を教えて頂けないでしょうか。
一様球以外の銀河も作ってみたいのでよろしくお願いします。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・14歳の自分に衝撃の事実を告げてください
- ・架空の映画のネタバレレビュー
- ・「お昼の放送」の思い出
- ・昨日見た夢を教えて下さい
- ・ちょっと先の未来クイズ第4問
- ・【大喜利】【投稿~10/21(月)】買ったばかりの自転車を分解してひと言
- ・メモのコツを教えてください!
- ・CDの保有枚数を教えてください
- ・ホテルを選ぶとき、これだけは譲れない条件TOP3は?
- ・家・車以外で、人生で一番奮発した買い物
- ・人生最悪の忘れ物
- ・【コナン30周年】嘘でしょ!?と思った○○周年を教えて【ハルヒ20周年】
- ・ハマっている「お菓子」を教えて!
- ・最近、いつ泣きましたか?
- ・夏が終わったと感じる瞬間って、どんな時?
- ・10秒目をつむったら…
- ・人生のプチ美学を教えてください!!
- ・あなたの習慣について教えてください!!
- ・牛、豚、鶏、どれか一つ食べられなくなるとしたら?
- ・都道府県穴埋めゲーム
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
フェルミエネルギー
-
量子力学的に光子のスピンは何...
-
ラザフォードの実験について
-
置換演算子
-
シュレーディンガーの人間
-
酢酸ビニルの乳化重合について。
-
この世に、完全な球体は存在し...
-
量子力学における同種粒子の個別性
-
スピンとパリティについて・・・
-
一番浴びている素粒子は?
-
ブラッグの式とデビソン-ガーマ...
-
ダイラタンシー現象と液状化現...
-
量子の観測問題と放射線被爆に...
-
普段の生活の中での位置エネルギー
-
位置エネルギー U
-
「U = mgh」の「U」は何の略な...
-
なぜ、イオンのM殻は18個まで入...
-
高校物理の力学の質問
-
XPS(ESCA):Ga2p1/2,3/2、この1...
-
電子の移動度と正孔の移動度に...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
おすすめ情報
内容が伝わりにくいようなので整理します。
要するに上記ページのρ(r)の密度分布に従う乱数の作り方を教えてください。
aは球の半径、Mは質量の総和です。