MathLabs
定理証明済み

ラプラス・ルンゲ・レンツベクトルの保存

内容

逆二乗力r¨=−GMr3r\ddot{\mathbf r} = -\dfrac{GM}{r^{3}}\mathbf rのもとでの運動において、ベクトルA=r˙×h−GMr^\mathbf A = \dot{\mathbf r}\times\mathbf h - GM\hat{\mathbf r}は時間的に一定であり、大きさは∣A∣=GMe|\mathbf A| = GMe、向きは焦点から近点へ向かう。

なぜ正しいのか?

エネルギーと角運動量だけではケプラー軌道の大きさと形は決まるが、平面内での向きは決まらない。ラプラス・ルンゲ・レンツベクトルは近点の位置を固定する追加の保存量であり、その存在自体が1/r21/r^2力に特有である――大半の中心力にはこのようなベクトルは存在せず、それがまさに(下記の制限三体問題のように)大半の中心力軌道が固定された閉曲線上に留まらない理由である。

証明の概略

微分するとdAdt=r¨×h+r˙×h˙−GMdr^dt\dfrac{d\mathbf A}{dt} = \ddot{\mathbf r}\times\mathbf h + \dot{\mathbf r}\times\dot{\mathbf h} - GM\dfrac{d\hat{\mathbf r}}{dt}。h\mathbf hは一定(定理1)なので中間項は消える。r¨=−GMr2r^\ddot{\mathbf r}=-\dfrac{GM}{r^{2}}\hat{\mathbf r}とh=r×r˙\mathbf h=\mathbf r\times\dot{\mathbf r}を用い、ベクトル三重積の恒等式r^×(r×r˙)=r(r^⋅r˙)−r˙(r^⋅r)\hat{\mathbf r}\times(\mathbf r\times\dot{\mathbf r}) = \mathbf r(\hat{\mathbf r}\cdot\dot{\mathbf r}) - \dot{\mathbf r}(\hat{\mathbf r}\cdot\mathbf r)よりr¨×h=−GMr2[rr˙ r^−rr˙]=GM(r˙r−r˙rr^)\ddot{\mathbf r}\times\mathbf h = -\dfrac{GM}{r^{2}}\big[r\dot r\,\hat{\mathbf r} - r\dot{\mathbf r}\big] = GM\left(\dfrac{\dot{\mathbf r}}{r} - \dfrac{\dot r}{r}\hat{\mathbf r}\right)を得る。ここでr˙=r^⋅r˙\dot r=\hat{\mathbf r}\cdot\dot{\mathbf r}は距離rrのスカラー変化率(r^⋅r=r\hat{\mathbf r}\cdot\mathbf r=rを使用)。

一方、r^=r/r\hat{\mathbf r}=\mathbf r/rを直接微分するとdr^dt=r˙r−r r˙r2=r˙r−r˙rr^\dfrac{d\hat{\mathbf r}}{dt} = \dfrac{\dot{\mathbf r}}{r} - \dfrac{\mathbf r\,\dot r}{r^{2}} = \dfrac{\dot{\mathbf r}}{r} - \dfrac{\dot r}{r}\hat{\mathbf r}となり、これは上で求めたr¨×h/GM\ddot{\mathbf r}\times\mathbf h/GMの括弧内の式とまったく同じである。よってr¨×h=GMdr^dt\ddot{\mathbf r}\times\mathbf h = GM\dfrac{d\hat{\mathbf r}}{dt}が恒等的に成り立つ。

したがってdAdt=GMdr^dt−GMdr^dt=0\dfrac{d\mathbf A}{dt} = GM\dfrac{d\hat{\mathbf r}}{dt} - GM\dfrac{d\hat{\mathbf r}}{dt} = \mathbf 0となり、A\mathbf Aが一定であることが証明された。近点(そこではr˙\dot{\mathbf r}はr^\hat{\mathbf r}に垂直な純接線方向で速さvp=h/rpv_p=h/r_p)でA\mathbf Aを評価すると、r˙×h\dot{\mathbf r}\times\mathbf hは大きさh2/rph^{2}/r_pで−r^-\hat{\mathbf r}方向を向くためA=(h2rp−GM)r^\mathbf A = \left(\dfrac{h^{2}}{r_p}-GM\right)\hat{\mathbf r}となり、軌道方程式からrp=p/(1+e)=h2GM(1+e)r_p=p/(1+e)=\dfrac{h^{2}}{GM(1+e)}を用いるとA=GMe r^p\mathbf A = GMe\,\hat{\mathbf r}_pとなって、∣A∣=GMe|\mathbf A|=GMeが近点方向を正しく向くことが確認できる。■\blacksquare

この定理を使うトピック

ステップごとの証明

この定理のステップごとの証明はまだありません。

参考文献

  1. Carl D. Murray, Stanley F. Dermott (1999). Solar System Dynamics
  2. Alain Chenciner, Richard Montgomery (2000). A remarkable periodic solution of the three-body problem in the case of equal masses
  3. NASA Science (2024). Webb's Orbit at Sun-Earth Lagrange Point 2 (L2)
  4. NASA Science (2024). Basics of Spaceflight: A Gravity Assist Primer