[论文解读] Lane-Emden Equation: perturbation method

Z. F. Seidov|arXiv (Cornell University)|Feb 5, 2004

Mathematics and Applications参考文献 3被引用 4

一句话总结

本文提出一种数值摄动方法，用于求解任意多级指数 $ n $ 的Lane-Emden方程，将解析解的适用范围扩展至已知情况 $ n = 0, 1, 5 $ 之外。通过将解按 $ \delta = n - n_0 $ 的幂级数展开，高精度计算了无量纲半径 $ X $、质量 $ \mu $ 以及中心密度与平均密度之比 $ \rho_c/\rho_m $，并与已知数值数据对比验证，结果表明该方法在 $ 0 < n < 5 $ 范围内具有适用性。

ABSTRACT

The perturbation method is applied to numerical solution of the Lane-Emden Equation of arbitrary index n, and the global parameters of polytropes are found as function of polytropic index n.

研究动机与目标

• 开发一种用于求解任意多级指数 $ n $ 的Lane-Emden方程（LEE）的数值摄动方法，突破传统可解析求解的三种情况（$ n = 0, 1, 5 $）的限制。
• 以高精度计算全局结构参数——无量纲半径 $ X $、质量 $ \mu $，以及中心密度与平均密度之比 $ \rho_c/\rho_m $，作为 $ n $ 的函数。
• 通过与 Jabbar（1984）、Horedt（1987, 1990）及其他研究者的文献数据对比，对摄动方法进行数值验证。
• 评估该方法在摄动项与未扰动解量级相近时的有效性与局限性，特别是在 $ n \approx 5 $ 附近区域。

提出的方法

• 将Lane-Emden解按 $ \delta = n - n_0 $ 的幂级数展开，至二阶项：$ y = y_0 + \delta y_1 + \delta^2 y_2 $，其中 $ n_0 \in \{0, 1, 5\} $。
• 通过非线性项 $ y^n $ 的摄动展开，推导出 $ y_0 $、$ y_1 $ 和 $ y_2 $ 的耦合二阶常微分方程，得到方程（6）至（8）。
• 引入辅助变量 $ z_i = -x^2 y_i' $，将系统转化为六阶一阶常微分方程组，利用 Mathematica 中的 NDSolve 实现稳定数值积分。
• 在 $ x = 0 $ 处采用高阶幂级数展开（最高至 $ x^{24} $）作为初值，其来源于Lane-Emden解的已知泰勒级数。
• 在中间点 $ x = X_0 $ 处应用匹配技术，传播解并利用 $ \mu = -X^2 y'(X) $ 确定 $ X $、$ \mu $ 和 $ \rho_c/\rho_m $。
• 通过将计算得到的 $ X(n) $、$ \mu(n) $ 和 $ \rho_c/\rho_m $ 与 Jabbar（1984）和 Horedt（1987, 1990）的已发表值对比，验证结果，显示四位至六位有效数字的一致性。

实验结果

研究问题

• RQ1 该数值摄动方法能否在 $ n \in (0,5) $ 范围内，特别是在 $ n \approx 5 $ 附近，高精度计算多级指数的全局参数？
• RQ2 当 $ n $ 接近 0、1 和 5 时，二阶摄动修正项 $ y_1 $ 和 $ y_2 $ 在预测 $ X(n) $ 和 $ \mu(n) $ 时的精度如何？
• RQ3 当摄动项 $ y_1 $、$ y_2 $ 的大小与未扰动解 $ y_0 $ 相当时，该摄动方法的有效性在多大程度上仍成立？
• RQ4 计算得到的 $ X(n) $、$ \mu(n) $ 和 $ \rho_c/\rho_m $ 与 Jabbar（1984）和 Horedt（1987, 1990）的既定数值结果相比如何？
• RQ5 在 $ x $ 区间较大时，特别是 $ n \approx 5 $ 时，该摄动方案的数值稳定性和收敛性表现如何？

主要发现

• 当 $ n = 3.1 $ 时，计算得到的半径 $ X(n=3.1) = 7.3073764235 $ 与 Jabbar 的结果一致至六位有效数字，与 Horedt 的结果一致至五位有效数字。
• 当 $ n = 2.9 $ 时，$ X(n=2.9) = 6.5273727715 $ 与 Jabbar 的结果一致至六位有效数字，与 Horedt 的结果一致至五位有效数字。
• 质量参数 $ \mu(n=3.1) = 1.9897152342 $ 与 Horedt 的估计值 $ 1.9897028553 $ 一致，差异小于 $ 10^{-5} $。
• 当 $ n = 3.6 $ 时，$ X(n=3.6) = 10.2979713253 $ 与 Jabbar 的值一致至六位有效数字，与 Horedt 的值一致至五位有效数字。
• 当 $ n = 4.6 $ 时，$ X(n=4.6) = 40.0619871651 $ 与 Jabbar 的值 $ 40.4132296343 $ 相差约 $ 0.86\% $，且 $ \mu(n=4.6) = 1.7306833002 $ 与 Horedt 的结果 $ 1.7307484954 $ 一致至 $ 3.8 \times 10^{-5} $。
• 该方法成功计算了 $ n \in (0,5) $ 范围内的 $ \mu(n) $ 和 $ X(n) $，对大多数测试值，$ \mu $ 的相对误差低于 $ 10^{-4} $，证实其在天体物理建模中的实际应用价值。