[论文解读] Exact solutions of Laplace equation by differential transform method

M. Jamil Amir, Muhammad Yaseen+1 more|arXiv (Cornell University)|Dec 22, 2013

Differential Equations and Numerical Methods被引用 4

一句话总结

本文将二维微分变换法（DTM）应用于求解具有狄利克雷和诺伊曼边界条件的拉普拉斯方程的四个模型，以解析方式获得精确的闭式解。该方法通过极少的迭代次数，高效地生成已知函数的级数展开形式的精确解，相较于其他迭代方法显著减少了计算量，在所有测试案例中均表现出高度的准确性与快速收敛性。

ABSTRACT

In this paper, we solve Laplace equation analytically by using differential transform method. For this purpose, we consider four models with two Dirichlet and two Neumann boundary conditions and obtain the corresponding exact solutions. The obtained results show the simplicity of the method and massive reduction in calculations when one compares it with other iterative methods, available in literature. It is worth mentioning that here only a few number of iterations are required to reach the closed form solutions as series expansions of some known functions.

研究动机与目标

• 开发并应用二维微分变换法（DTM）求解具有精确解析解的拉普拉斯方程模型。
• 展示该方法在计算复杂度方面相较于现有迭代技术（如ADM、HAM和VIM）的高效性。
• 求解四个物理模型——其中两个具有狄利克雷边界条件，两个具有诺伊曼边界条件——且无需线性化或离散化处理。
• 表明仅需少数几次迭代即可获得与已知特殊函数匹配的闭式级数解。
• 通过在所有测试案例中精确恢复已知解，验证该方法的准确性和收敛速度。

提出的方法

• 通过使用二阶偏导数的微分变换，将拉普拉斯方程转化为递推关系，应用二维微分变换法。
• 该方法基于如下定义：$ U(m,n) = \frac{1}{m!n!} \left[ \frac{\partial^{m+n} u(x,y)}{\partial x^m \partial y^n} \right]_{(x_0,y_0)} $，其中 $ u(x,y) = \sum_{m=0}^\infty \sum_{n=0}^\infty U(m,n) (x-x_0)^m (y-y_0)^n $。
• 边界条件被转换为谱域中的初始条件，从而实现从已知函数展开式递归计算 $ U(m,n) $。
• 从拉普拉斯方程 $ u_{xx} + u_{yy} = 0 $ 推导出递推关系 $ (m+1)(m+2)U(m+2,n) + (n+1)(n+2)U(m,n+2) = 0 $。
• 通过逆变换重构解，并将级数模式与已知函数（如 $ \sinh x $、$ \cos y $、$ \cosh x \sin y $）进行匹配，以识别精确的闭式表达。
• 该方法避免了数值离散化，通过将边界数据转换为级数域中的初始系数条件，同时处理狄利克雷和诺伊曼条件。

实验结果

研究问题

• RQ1 二维微分变换法（DTM）能否为具有狄利克雷边界条件的拉普拉斯方程模型产生精确的解析解？
• RQ2 DTM能否在无需线性化或离散化的情况下，高效求解具有诺伊曼边界条件的拉普拉斯方程模型？
• RQ3 与ADM、HAM和VIM等其他迭代方法相比，DTM在求解拉普拉斯方程时的计算复杂度如何？
• RQ4 DTM在恢复已知精确解时，其收敛速度有多快，所需迭代次数是否极少？
• RQ5 DTM能否将精确解重构为 $ \sinh x $、$ \cos y $、$ \cosh x \sin y $ 和 $ \cos 2x \cosh 2y $ 等初等函数的级数展开形式？

主要发现

• DTM成功恢复了第一个狄利克雷问题的精确解 $ u(x,y) = \sinh x \sin y $，其级数展开与已知函数形式完全一致。
• 对于第二个狄利克雷问题，该方法得到了精确解 $ u(x,y) = \cosh x \sin y $，并通过级数重构得到验证。
• 在第一个诺伊曼问题中，DTM生成了精确解 $ u(x,y) = \cos 2x \cosh 2y $，其系数与推导出的递推关系及已知函数级数完全匹配。
• 对于第二个诺伊曼问题，该方法成功恢复了精确解 $ u(x,y) = \cos x \sinh y $，表明其在不同边界类型下均具有一致性。
• 该方法仅需少数几次迭代即可实现快速收敛，且所需项数极少，即可与已知特殊函数匹配。
• 与其它迭代方法相比，DTM显著降低了计算工作量，这体现在可直接推导出闭式解，而无需迭代逼近。