重力井戸シミュレーター 戻る
対話型シミュレーター

重力井戸シミュレーター

クリックで重力源を配置し、時空のゴム膜アナロジーでグリッドが歪む様子を体験。テスト粒子が軌道を描き、一般相対論の直感的イメージを掴もう。

クリックで重力源配置 右クリックで削除 ニュートン重力 + Verlet積分
重力ポテンシャル井戸
パラメータ
新規質量
重力定数 G
グリッド密度
アクション
プリセット
表示設定
グリッド表示
パーティクルトレイル
力ベクトル
統計
計算結果
0
重力源
0
粒子数
0
全運動E
60
FPS
50
新規質量
1.0
G
理論・主要公式

$$F = -\frac{G m_1 m_2}{r^2}$$

万有引力の法則(N):$G = 6.674 \times 10^{-11}$ N·m²/kg²、$r$ は質点間距離(m)。

$$U = -\frac{G m_1 m_2}{r}$$

重力ポテンシャルエネルギー(J):$r \to \infty$ でゼロ、近づくほど負に深くなる。

$$v_{esc} = \sqrt{\frac{2GM}{r}}$$

脱出速度(m/s):半径 $r$ の位置から引力圏を脱出するのに必要な最低速度。

重力井戸シミュレーターとは

🙋
このシミュレーターでグリッドが歪んで見えるのは、何を表しているんですか?
🎓
これは、アインシュタインの「ゴム膜アナロジー」を可視化したものだよ。大まかに言うと、質量が置かれると、時空という舞台そのものが歪む、という一般相対論の考え方を2次元のシートで表現しているんだ。上の「重力源の質量」スライダーを動かしてみて。質量が大きいほど、グリッドが深く沈み込むのがわかるよね。これが「重力井戸」だ。
🙋
え、そうなんですか?でも、画面をクリックして動かしている粒子は、なぜあんなふうにクルクル回ったり、吸い込まれたりするんですか?
🎓
粒子は、歪んだグリッド(時空)の上を、できるだけまっすぐに進もうとするんだ。でも、舞台自体が歪んでいるから、その軌道は曲がって見える。これが重力による「落下」や「公転」の正体だ。例えば、太陽の周りを地球が回るのも同じ理屈だね。「粒子の初速度」スライダーを変えると、軌道が楕円になったり、放物線になったりするのが観察できるよ。
🙋
「ソフトニング」ってパラメータもありますね。これは何のためにあるんですか?
🎓
いいところに気が付いたね!これは実務(シミュレーション)で非常に重要なパラメータなんだ。重力は距離がゼロに近づくと、理論上は無限大に発散して計算が破綻してしまう。それを防ぐために、分母に小さな値を足して計算を安定させているんだ。右端の「ソフトニング」スライダーを小さくしすぎると、粒子が重力源に激突して計算が吹っ飛ぶ様子を体験できるよ。現場では、星の衝突シミュレーションなどでこうした工夫が必須なんだ。

よくある質問

複数の重力源がある場合、粒子は各重力源からの引力のベクトル和を受けて運動します。軌道は単純な楕円ではなく、複雑なカオス的な動きになることがあります。例えば、2つの重力源の間を8の字に回る軌道や、一方に捕獲されてから放出されるスイングバイのような現象も観察できます。
グリッドの各頂点の高さ(変形の深さ)は、その位置における重力ポテンシャルの合計値を視覚化したものです。重力源の質量が大きいほど、また重力源に近いほど深く沈みます。これは一般相対論で言う時空の歪みを、ゴム膜アナロジーで直感的に理解するための表現です。
現バージョンでは、粒子をクリック&ドラッグして放すことで、放した方向と速度で初速度を与えられます。ドラッグ距離が長いほど初速が大きくなります。また、設定パネルから数値で直接初速度ベクトル(x, y成分)を入力することも可能です。
重力源や粒子の数が多すぎると計算負荷が高まります。まずは不要な重力源を右クリックで削除してください。また、設定メニューで「グリッド解像度」を下げる(例:高→中)と描画負荷が軽減されます。それでも重い場合は、ブラウザのタブをリロードして初期状態に戻してください。

実世界での応用

天体物理学・宇宙探査:探査機を他の惑星に送る際の軌道設計(スイングバイ)に、このような重力場中の運動計算が不可欠です。シミュレーターで初速度を変えて放物線軌道を作る操作は、その基礎にあたります。

宇宙望遠鏡の位置決め:ハッブル宇宙望遠鏡やジェイムズ・ウェッブ宇宙望遠鏡は、地球の重力井戸の中の特定の点(ラグランジュ点)に配置され、安定して観測を行います。重力場の理解が軌道決定に役立ちます。

高精度測位システム(GPS):GPS衛星は地球の重力井戸の外側を高速で回っています。一般相対論効果(時空の歪みによる時間の進み方の違い)を補正しなければ、位置情報は数分で数kmもずれてしまいます。

数値相対論シミュレーション:ブラックホール連星の合体など、強い重力場中の極限現象を計算する「数値相対論」の分野では、時空の歪みを格子で表現し、その上で物理量を計算する手法が発展しています。このシミュレーターはそのごく初歩的な概念体験と言えます。

よくある誤解と注意点

まず、このシミュレーターは2次元のアナロジーだという点を押さえよう。実際の時空の歪みは3次元空間+時間の4次元で起きており、この「ゴム膜」のイメージはあくまで理解の助け。特に、粒子が「落ちる」のは膜の「下」方向ではなく、歪んだ幾何学自体の中を進む結果だというのが本質だよ。

次に、「ソフトニング」パラメータの設定。これを0.1や0.01といった小さな値にすると、粒子が重力源に激突して計算が発散しやすい。実務では、シミュレーションの対象スケールに合わせて調整する。例えば、銀河衝突シミュレーションでは星同士の距離が非常に大きいので、ソフトニングは比較的大きく取って計算を安定させる。逆に、小さい系を精密に計算する時は小さくするけど、その分タイムステップ(時間の刻み幅)も小さくする必要がある。このシミュレーターで「粒子の初速度」を大きくした時に軌道が暴れるのは、タイムステップが固定だから起こる現象で、これも実務では注意すべき点だ。

最後に、「Verlet積分」の特性。この手法はエネルギー保存性が良い反面、速度は直接計算しない(位置の差分から求める)ので、速度に依存する力(空気抵抗など)を後で追加する時は一工夫必要になる。このシミュレーターは純粋な重力だけだからシンプルに使えているけど、応用する時は頭の片隅に入れておこう。

使い方ガイド

  1. スライダー「質量」(sl-mass)で重力源の質量を設定します。太陽質量相当の10⁸kg~10³⁰kgの範囲で調整可能です。
  2. 「重力定数」(sl-g)でニュートン定数G=6.67×10⁻¹¹を基準値とし、効果を強調するため10~1000倍にスケール調整できます。
  3. シミュレータ上の任意位置をクリックして重力源を配置すると、周囲のグリッド(tog-grid)が即座に歪み、テスト粒子(tog-trail)が曲線軌道を描きます。
  4. 「グリッド解像度」(sl-grid)でメッシュ分割数を変更し、時空歪みの詳細度を調整します。

具体的な計算例

質量M=1.989×10³⁰kg(太陽質量)、ソフトニングパラメータh=0.15au、初期速度v=42.1km/s(地球の公転速度)でシミュレーション実行時、Verlet積分により計算タイムステップΔt=0.001年で約365ステップ進めると、テスト粒子は�楕円軌道(軌道長半径a≈1au)を描きます。重力加速度a=GM/r²から、1au地点ではa≈5.9mm/s²となり、グリッド変形は半径r²に反比例する強度で視覚化されます。

実務での注意点

  1. Verlet積分はエネルギー保存性が優れているため、数万ステップの長期軌道計算でも軌道減衰が最小限に抑えられます。ただしソフトニング係数h(内部パラメータ)を過度に小さくするとスケール不安定性が発生するため、グリッド間隔の5~10%程度に設定してください。
  2. 複数の重力源を配置する場合、各源間の距離がソフトニング長より小さいと特異性が抑制され、物理的な接近では発散が防止されます。
  3. 時空グリッド変形の曲率は局所的な質量密度に依存するため、高密度領域では黒体放射計算のシュバルツシルト半径rs=2GM/c²の概念的理解に役立ちます。標準ソルバー(odeint等)との結果検証時には相対誤差1%以内を目安にしてください。