[論文レビュー] Gravitational potential and energy of homogeneous rectangular parallelepiped

Z. F. Seidov, P. I. Skvirsky|arXiv (Cornell University)|Feb 28, 2000

Dynamics and Control of Mechanical Systems参考文献 2被引用数 7

ひとこと要約

本稿では、記号的積分を用いて、均一な長方形平行六面体（RP）の重力ポテンシャルおよびエネルギーの正確な解析的式を提示する。Mathematicaを用いた記号的計算により、初等関数を用いた閉形式の解が得られ、任意の点（内部・外部を問わず）での正確な計算が可能となる。主な結果として、立方体の中心および頂点におけるポテンシャル、および立方体、細長い棒、平板に対する正確なポテンシャルエネルギーの公式が得られる。

ABSTRACT

Gravitational potential and gravitational energy are presented in analytical form for homogeneous rectangular parallelepiped.

研究の動機と目的

• 均一な長方形平行六面体（RP）の重力ポテンシャルについて、初等関数を用いた閉形式の正確な解析的式を導出すること。
• RPの体積上でのポテンシャルの積分を用いて重力ポテンシャルエネルギーを計算し、正確な公式を導出すること。
• 立方体、細長い棒、薄板などの極限状態を検討することで、一般解の妥当性を検証し、スケーリングの性質を調査すること。
• 楕円体やRPを含むさまざまな形状におけるWUM比（ポテンシャルエネルギー / （中心でのポテンシャル × 質量））を比較し、普遍的な傾向を明らかにすること。

提案手法

• 固定された積分限界を用いて、RPの体積上でのニュートン的重力ポテンシャルを三重積分として定式化する。
• 記号的計算（Mathematica）を用い、積分を順にz、次にy、最後にxの順に逐次評価し、対称的な最終式を得る。
• ポテンシャルの核心的要素として、関数 $ v(x',y',z') = y'z'\ln[x' + \sqrt{x'^2 + y'^2 + z'^2}] - \frac{x'^2}{2}\arctan\left(\frac{y'z'}{x'\sqrt{x'^2 + y'^2 + z'^2}}\right) $ を用いる。
• ポテンシャルエネルギーを、体積上でのポテンシャルの三重積分により導出し、対称性を活用して簡略化し、極限状態による検証も行う。
• スケーリング不変性を適用して、形状間での普遍的比較が可能な次元なしの形で結果を表現する。
• 既知の極限状態（線分、長方形、楕円体）を用いて結果を検証し、WUM比を楕円体の既知の値 $ 2/5 $ と比較する。

実験結果

リサーチクエスチョン

• RQ1 空間内の任意の点における、均一な長方形平行六面体の重力ポテンシャルの正確な解析的式は何か？
• RQ2 RPのポテンシャルエネルギーはその寸法および密度にどのように依存するか。また、閉形式で表現できるか？
• RQ3 均一な立方体の対称的中心点（中心、頂点、面の中心）における重力ポテンシャルの値は何か？
• RQ4 ポテンシャルエネルギーと中心ポテンシャルおよび質量の積の比（WUM比）は、さまざまなRPの幾何形状においてどのように変化するか。また、楕円体の場合と比較するとどうなるか？
• RQ5 得られたRPの公式は、細長い棒や薄板といった極限状態において、既知の結果に簡略化できるか？

主な発見

• 中心における均一なRPの重力ポテンシャルは、$ U(0,0,0) = 24\ln\left(\frac{1+\sqrt{3}}{\sqrt{2}}\right) - 2\pi \approx 9.52017 $ であり、単位は $ a^2 G \rho $ である。
• 立方体の頂点におけるポテンシャルは、中心ポテンシャルのちょうど半分に等しい：$ U(1,1,1) = \frac{1}{2} U(0,0,0) $。
• 辺の長さが $ 2a $ の均一な立方体のポテンシャルエネルギーは、一般式（31）から導かれて $ W_C = 30.117 G \rho^2 a^5 $ である。
• 細長い棒（$ a \gg b = c $）の場合、ポテンシャルエネルギーは $ W_{\text{stick}} = \frac{32}{3} G \rho^2 a b^4 \ln\left(\frac{a}{b}\right) $ に比例する。
• 薄い正方形の平板（$ a \to 0 $、$ b = c $）の場合、ポテンシャルエネルギーは $ W_S = 47.5714 G \rho^2 b^3 a^2 $ であり、WUM比は約 0.421673 である。
• WUM比（ポテンシャルエネルギー / （中心ポテンシャル × 質量））は、立方体では約 0.39544 に最小値を示し、細長い棒では 0.5 に近づき、薄板では約 0.421673 に近づく。これは、楕円体の $ WUM = 2/5 = 0.4 $ とは異なり、非単調な挙動を示しており、特徴的である。