[论文解读] The DJ method for exact solutions of Laplace equation

Muhammad Yaseen, Muhammad Samraiz+1 more|arXiv (Cornell University)|Aug 16, 2012

Fractional Differential Equations Solutions参考文献 9被引用 5

一句话总结

本文将Daftardar-Gejji与Jafari（DJ）迭代法应用于拉普拉斯方程在狄利克雷与诺伊曼边界条件下的精确解析解求解。该方法通过连续逼近构造解，无需使用阿多米安多项式或拉格朗日乘子，所有测试的物理模型均得到高精度的闭式精确解，且计算量极小。

ABSTRACT

In this paper, the iterative method developed by Daftardar-Gejji and Jafari (DJ method) is employed for analytic treatment of Laplace equation with Dirichlet and Neumann boundary conditions. The method is demonstrated by several physical models of Laplace equation. The obtained results show that the present approach is highly accurate and requires reduced amount of calculations compared with the existing iterative methods.

研究动机与目标

• 展示DJ方法在求解拉普拉斯方程时获得精确解析解的有效性。
• 将DJ方法的精度与计算效率与现有迭代技术（如阿多米安分解法与变分迭代法）进行比较。
• 确立该方法在无需限制性假设或复杂多项式计算的情况下处理线性与非线性项的能力。
• 求解拉普拉斯方程的四个物理模型，其中两个为狄利克雷边界条件，两个为诺伊曼边界条件，并验证其精确解。

提出的方法

• DJ方法将拉普拉斯方程重新表述为包含非线性算子与已知函数的泛函方程。
• 通过递归逼近构造级数解：$ u_0 = f $，$ u_1 = N(u_0) $，且 $ u_{m+1} = N(\sum_{j=0}^{m} u_j) - N(\sum_{j=0}^{m-1} u_j) $。
• 该方法避免了在非线性项中使用阿多米安多项式，且无需像变分迭代法那样确定拉格朗日乘子。
• 解以无限级数形式获得：$ u = \sum_{i=0}^{\infty} u_i $，当级数终止或与已知函数匹配时，该级数收敛至精确解。
• 通过将边界条件应用于级数解，可确定未知函数（如 $ g(x) $），并通过泰勒展开中的系数匹配，得到闭式解。
• 该方法应用于四个模型：两个狄利克雷边界条件与两个诺伊曼边界条件，每个模型均得到闭式精确解。

实验结果

研究问题

• RQ1 DJ方法能否为带有狄利克雷边界条件的拉普拉斯方程产生精确解析解？
• RQ2 DJ方法能否在无需复杂多项式分解的情况下高效求解带有诺伊曼边界条件的拉普拉斯方程？
• RQ3 与阿多米安分解法和变分迭代法相比，DJ方法在精度与计算量方面表现如何？
• RQ4 DJ方法是否能在不施加非线性项限制性假设的情况下，通过连续逼近收敛至精确解？

主要发现

• DJ方法在所有四个测试的拉普拉斯方程模型中均成功获得高精度的闭式精确解，且计算量极小。
• 在例1中，通过DJ方法精确导出解 $ u(x,y) = \sinh(x)\cos(y) $，其中 $ g(x) = 0 $。
• 在例2中，通过系数匹配确定 $ g(x) = \cosh(x) $，获得解 $ u(x,y) = \cosh(x)\sin(y) $。
• 例3得到精确解 $ u(x,y) = \cos(2x)\cosh(2y) + C_0 $，其中 $ g(x) = \cos(2x) + C_0 $，由 $ g''(x) = -4\cos(2x) $ 推导得出。
• 例4得到精确解 $ u(x,y) = \cos(x)\sinh(y) + C_0 $，其中 $ g(x) = C_0 $，由 $ g''(x) = 0 $ 验证。
• 该方法仅通过少数几项即实现对精确解的收敛，证实其在具有混合边界条件的线性偏微分方程中具有高效性与可靠性。