[论文解读] A Brief Introduction to the Adomian Decomposition Method, with Applications in Astronomy and Astrophysics

M. K. Mak, Chun Sing Leung+1 more|arXiv (Cornell University)|Feb 21, 2021

Fractional Differential Equations Solutions参考文献 12被引用 6

一句话总结

本文对Adomian分解法（ADM）进行了教学性介绍，这是一种强大的解析技术，可用于求解非线性常微分方程和偏微分方程，无需离散化或线性化。该方法在求解一阶方程（线性、Bernoulli型、Riccati型、Abel型）、二阶与五阶常微分方程，以及Fisher-Kolmogorov型偏微分方程时，与解析解完全一致。该方法成功应用于天体物理学中的关键问题，包括开普勒方程、Lane-Emden方程和Schwarzschild方程。

ABSTRACT

The Adomian Decomposition Method (ADM) is a very effective approach for solving broad classes of nonlinear partial and ordinary differential equations, with important applications in different fields of applied mathematics, engineering, physics and biology. It is the goal of the present paper to provide a clear and pedagogical introduction to the Adomian Decomposition Method and to some of its applications. In particular, we focus our attention to a number of standard first-order ordinary differential equations (the linear, Bernoulli, Riccati, and Abel) with arbitrary coefficients, and present in detail the Adomian method for obtaining their solutions. In each case we compare the Adomian solution with the exact solution of some particular differential equations, and we show their complete equivalence. The second order and the fifth order ordinary differential equations are also considered. An important extension of the standard ADM, the Laplace-Adomian Decomposition Method is also introduced through the investigation of the solutions of a specific second order nonlinear differential equation. We also present the applications of the method to the Fisher-Kolmogorov second order partial nonlinear differential equation, which plays an important role in the description of many physical processes, as well as three important applications in astronomy and astrophysics, related to the determination of the solutions of the Kepler equation, of the Lane-Emden equation, and of the general relativistic equation describing the motion of massive particles in the spherically symmetric and static Schwarzschild geometry.

研究动机与目标

• 为不熟悉该技术的研究人员提供Adomian分解法（ADM）的清晰、分步式教学介绍。
• 展示该方法在求解广泛非线性常微分方程与偏微分方程方面的有效性。
• 通过展示ADM解与标准方程已知解析解之间的完全一致，验证ADM的有效性。
• 将该方法扩展至高阶常微分方程，并引入拉普拉斯-Adomian分解法以提升收敛性。
• 展示ADM在求解天文学与天体物理学中关键方程（包括开普勒方程、Lane-Emden方程与广义相对论运动方程）中的实用性。

提出的方法

• 将非线性微分方程的解分解为无穷多项函数的级数形式。
• 使用Adomian多项式表示方程中的非线性项，通过幂级数展开递归计算。
• 应用递归算法，利用Adomian多项式与线性算子逐项计算级数解的每一项。
• 对于拉普拉斯-Adomian变体，先对原方程应用拉普拉斯变换，再应用ADM以增强收敛性。
• 通过迭代计算级数中的各项直至收敛，无需线性化或摄动法。
• 通过将ADM级数解与基准方程的已知精确解进行比较，验证结果的正确性。

实验结果

研究问题

• RQ1 Adomian分解法能否为标准非线性一阶常微分方程产生精确或高度精确的解析解？
• RQ2 ADM在求解高阶非线性常微分方程（如二阶与五阶方程）时表现如何？
• RQ3 拉普拉斯-Adomian分解法能否提升非线性二阶常微分方程的收敛性与精度？
• RQ4 ADM在求解具有物理意义的偏微分方程（如Fisher-Kolmogorov方程）方面可应用到何种程度？
• RQ5 ADM能否有效求解天体物理学中的基本方程，包括开普勒方程、Lane-Emden方程与Schwarzschild几何中的相对论运动方程？

主要发现

• 对于线性常微分方程 $ \frac{dy}{dx} + 2xy = 4x^3 $，ADM解与解析解完全一致，证实了该方法的准确性。
• 对于Bernoulli方程 $ \frac{dy}{dx} - 2xy = -4x^3 y^2 $，ADM得到的级数解收敛于精确解。
• Riccati方程 $ \frac{dy}{dx} = 2e^x - e^{-x} y^2 $ 的ADM解与已知解析解完全相同。
• Abel方程 $ \frac{dy}{dx} = x + 3xy + 3xy^2 + xy^3 $ 的ADM解与精确解一致，证明其适用于更高阶非线性项。
• ADM成功求解了二阶常微分方程 $ \frac{d^2y}{dx^2} + 4\frac{dy}{dx} + 3y = 3 $，其级数解收敛于精确解。
• 该方法被扩展至五阶常微分方程 $ \frac{d^5y}{dx^5} - 3\frac{d^4y}{dx^4} - 5\frac{d^3y}{dx^3} + 15\frac{d^2y}{dx^2} + 4\frac{dy}{dx} - 12y = 0 $，并得到收敛的级数解。
• 拉普拉斯-Adomian分解法有效求解了非线性二阶常微分方程，显著提升了收敛速度。
• ADM为初始条件为 $ c(0,x,y,z) = x^2 + y^2 + z^2 $ 的Fisher-Kolmogorov偏微分方程提供了解析解，展示了其在非线性PDE中的实用性。
• ADM对开普勒方程、Lane-Emden方程以及Schwarzschild几何中的相对论运动方程均给出了高精度解，证实了其在天体物理建模中的价值。