[論文レビュー] Nested spheroidal figures of equilibrium I. Approximate solutions for rigid rotations

J-M Huré|arXiv (Cornell University)|Feb 22, 2022

Astrophysics and Star Formation Studies参考文献 18被引用数 4

ひとこと要約

本稿では、共焦点の扁平な回転楕円体面で分離された2つの均一で非圧縮性の成分から成る、ネストされた剛体回転楕円体的平衡状態について、近似的な解析解を提示する。共焦点パラメータにおける摂動展開を用い、小楕円率および小共焦点パラメータの範囲で、定常圧または円柱座標 $ R $ の2次関数としての界面圧力モデルの下で、回転速度および圧力分布の閉形式表現を導出する。数値的自己無撞着場（SCF）解との強い一致が得られている。

ABSTRACT

We discuss the equilibrium conditions for a body made of two homogeneous components separated by oblate spheroidal surfaces and in relative motion. While exact solutions are not permitted for rigid rotation (unless a specific ambient pressure), approximations can be obtained for configurations involving a small confocal parameter. The problem then admits two families of solutions, depending on the pressure along the common interface (constant or quadratic with the cylindrical radius). We give in both cases the pressure and the rotation rates as a function of the fractional radius, ellipticities and mass-density jump. Various degrees of flattening are allowed but there are severe limitations for global rotation, as already known from classical theory (e.g. impossibility of confocal and coelliptical solutions, gradient of ellipticity outward). States of relative rotation are much less constrained, but these require a mass-density jump. This analytical approach compares successfully with the numerical solutions obtained from the Self-Consistent-Field method. Practical formula are derived in the limit of small ellipticities appropriate for slowly-rotating star/planet interiors.

研究の動機と目的

• 共焦点の扁平な回転楕円体面で分離された2つの剛体回転する均一で非圧縮性の楕円体的成分の平衡配置について、近似的な解析的解を導出すること。
• 共焦点性および共楕円率が緩和された場合に、このようなネスト構造が剛体回転下で機械的平衡に達する条件を調査すること。
• 定常界面圧力と円柱座標 $ R $ の2次関数としての界面圧力という2つの異なる界面圧力モデルを検討し、それらが回転速度および構造的パラメータに与える影響を調査すること。
• 小楕円率の極限において、実用的で低次の解析的公式を提供し、ゆっくり回転する星間物質および惑星内部に適用可能であること。
• 自己無撞着場（SCF）数値解と照合し、小共焦点パラメータ条件下で解析的手法の正確性を検証すること。

提案手法

• 共焦点パラメータ $ c $ を小量として扱い、楕円率ではなく $ c $ の線形の摂動展開を用いる。
• 回転系におけるベルヌーイ方程式を適用し、遠心力と重力ポテンシャルを組み合わせて静水的平衡を満たす。
• 楕円率 $ \epsilon $ に依存する楕円体調和係数 $ A_0, A_1, A_3 $ を用いて、内部重力ポテンシャル $ \Psi_{\text{int.}} $ を導出する。
• 2つの成分間の界面で圧力バランスを課し、2つのケースを検討する：定常界面圧力、および円柱座標 $ R $ の2次関数としての圧力依存性。
• ベルヌーイ方程式および境界条件を用いて、回転速度 $ \Omega_1, \Omega_2 $ と中心圧力 $ p_c $ を解き、小 $ c $ における閉形式表現を導出する。
• 数値的SCFシミュレーションとの照合により、$ |c| \ll 1 $ 時に良好な一致が得られ、摂動的手法の有効性が確認された。

実験結果

リサーチクエスチョン

• RQ1 共焦点の扁平な回転楕円体面を有する2つの剛体回転する均一で非圧縮性の楕円体的成分が、どのような条件下で機械的平衡に達できるか？
• RQ2 界面圧力の仮定（定常 vs. $ R $ の2次関数）の違いが、得られる回転速度および構造的パラメータにどのように影響を与えるか？
• RQ3 小楕円率および小共焦点パラメータの極限において、回転速度および圧力の解析的近似式をどのように導出できるか？
• RQ4 同様の構造に対して、解析的解と数値的SCF解はどの程度一致するか？
• RQ5 安定で剛体回転するネストされた楕円体的系において、質量密度のジャンプおよび部分半径の制限は何か？

主な発見

• 解析モデルは $ |c| \ll 1 $ の配置に対して、数値的SCF解を良好に再現し、摂動的手法の有効性が検証された。
• 小楕円率の下で、主成分の回転速度は $ \Omega_2 \approx \sqrt{ \frac{2\pi G \rho_2}{3} } \left( 1 + \text{補正項} \right) $ で近似可能であり、補正項は $ \epsilon_2^2 $, $ \epsilon_1^2 $, および質量密度ジャンプ $ \alpha $ に依存する。
• グローバル回転（タイプC解）を実現するための質量密度ジャンプ $ \alpha $ は $ \alpha_C \approx 1 + \frac{3}{2} \left( \frac{\epsilon_2^2 - \epsilon_1^2}{1 - \epsilon_2^2} \right) q^2 $ で与えられ、$ \epsilon^2 $ および $ c $ の1次にまで有効である。
• 相対回転（タイプV解）では、埋め込まれた楕円体の回転速度 $ \Omega_1 $ と主成分の $ \Omega_2 $ が $ \epsilon_1, \epsilon_2, q, \alpha $ の関数として導出され、$ \Omega_1 \neq \Omega_2 $ であり、指定された圧力モデル下でも系は平衡を保つ。
• 主成分の回転速度プロファイル $ \Omega_i(R) $ は $ R $ に対して2次関数的であり、わずかな4次補正項を含むため、完全な剛体回転からわずかに逸脱していることが示された。
• 導出された公式を用いた簡単なF90コード（付録B）が提供されており、特に小楕円率のゆっくり回転する星や惑星の計算に有用である。