[论文解读] 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

一句话总结

本文为由两个均匀、不可压缩组分组成的嵌套、刚性旋转椭球平衡构型，其界面由共焦扁椭球面分隔，提出了近似解析解。通过在共焦参数上的微扰展开，推导出在小椭圆率和小共焦参数条件下，旋转速率与压强分布的闭式表达式，采用两种界面压强模型——恒定或与柱坐标半径呈二次关系——结果与数值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.

研究动机与目标

• 推导由共焦扁椭球面分隔的两个刚性旋转、均匀、不可压缩椭球组分的平衡构型的近似解析解。
• 研究在共焦性与共椭圆性被放松时，此类嵌套系统在刚性旋转下实现力学平衡的条件。
• 探讨两种不同的界面压强模型——恒定压强与与柱坐标半径呈二次关系——对旋转速率与结构参数的影响。
• 在小椭圆率极限下，提供实用的低阶解析公式，适用于缓慢旋转的恒星与行星内部。
• 将解析方法与数值自洽场（SCF）解进行对比，验证其在小共焦参数条件下的准确性。

提出的方法

• 在共焦参数 $ c $ 上使用微扰展开，将其视为小量并保留至一阶，而非在椭圆率上展开。
• 在旋转参考系中应用伯努利方程，结合离心势与引力势，以保证流体静力平衡。
• 利用依赖于椭圆率 $ \epsilon $ 的椭球谐系数 $ A_0, A_1, A_3 $ 推导内部引力势 $ \Psi_{\text{int.}} $。
• 在两组分界面处施加压强平衡条件，考虑两种情形：恒定界面压强与与柱坐标半径 $ R $ 呈二次依赖的压强。
• 通过伯努利方程与边界条件求解旋转速率 $ \Omega_1, \Omega_2 $ 与中心压强 $ p_c $，在小-$ c $ 极限下得到闭式表达式。
• 将结果与数值SCF模拟进行对比，发现在 $ |c| \ll 1 $ 条件下具有良好一致性，证实了微扰方法的有效性。

实验结果

研究问题

• RQ1 在何种条件下，两个刚性旋转、均匀、不可压缩的椭球组分在共焦扁椭球几何构型下可实现力学平衡？
• RQ2 对界面压强的不同假设（恒定 vs. 与 $ R $ 呈二次关系）如何影响最终的旋转速率与结构参数？
• 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 $ 的一阶下有效。
• 对于相对旋转（V型解），推导出嵌入组分 $ \Omega_1 $ 与主组分 $ \Omega_2 $ 的旋转速率，其值依赖于 $ \epsilon_1, \epsilon_2, q, \alpha $，且 $ \Omega_1 \neq \Omega_2 $，系统在指定压强模型下仍保持平衡。
• 主组分的旋转速度分布 $ \Omega_i(R) $ 对 $ R $ 呈二次依赖，伴有微小的四次校正项，表明即使在解析模型中也存在对完美刚性旋转的偏离。
• 该模型提供了一个简单的 F90 程序（附录B），可用于基于推导公式计算旋转速率与压强，特别适用于小椭圆率的缓慢旋转恒星与行星。