[论文解读] Solution of boundary value and eigenvalue problems for second order elliptic operators in the plane using pseudoanalytic formal powers

Raúl Castillo-Pérez, Vladislav V. Kravchenko+1 more|arXiv (Cornell University)|Feb 5, 2010

Algebraic and Geometric Analysis参考文献 13被引用 15

一句话总结

本文提出伪解析形式幂法（MPFP）以求解平面上二阶椭圆型微分算子的边值问题与特征值问题。通过利用伪解析函数理论，并基于生成序列的递归积分构造完备的零空间解系，该方法实现了高精度、无网格的配点法求解，表现出快速收敛性，并在性能上与有限元法具有竞争力。

ABSTRACT

We propose a method for solving boundary value and eigenvalue problems for the elliptic operator D=divpgrad+q in the plane using pseudoanalytic function theory and in particular pseudoanalytic formal powers. Under certain conditions on the coefficients p and q with the aid of pseudoanalytic function theory a complete system of null solutions of the operator can be constructed following a simple algorithm consisting in recursive integration. This system of solutions is used for solving boundary value and spectral problems for the operator D in bounded simply connected domains. We study theoretical and numerical aspects of the method.

研究动机与目标

• 开发一种用于求解平面上二阶椭圆型微分算子边值与特征值问题的无网格数值方法。
• 利用伪解析函数理论，为算子 $ D = \operatorname{div} p \operatorname{grad} + q $ 构造完备的零空间解系。
• 提供对传统基于网格划分方法（如有限元法）的替代方案，尤其适用于分离变量法失效的区域。
• 通过数值实验展示该方法的效率、收敛速率与精度。
• 通过完备解系近似特征值，采用行列式求根法，将方法推广至特征值问题。

提出的方法

• 该方法通过伪解析形式幂构造椭圆型算子 $ D = \operatorname{div} p \operatorname{grad} + q $ 的完备解系，其源自Vekua方程。
• 需要一个特定解 $ u_0 $ 满足 $ D u = 0 $，使得 $ f = p^{1/2} u_0 $ 光滑、非零且在正交坐标系下可分离为 $ f = S(s)T(t) $。
• 利用近期文献中的算法，递归构造伪解析形式幂的生成序列，从而实现解系的显式构造。
• 边值问题的解通过这些形式幂的线性组合近似，系数由边界上的配点法确定。
• 对于特征值问题，方法寻找使矩阵 $ U = (u_k(z_j)) $ 的行列式为零的 $ \lambda $，其中 $ u_k $ 为形式幂，$ z_j \in \partial\Omega $。
• 该方法推广了拉普拉斯方程的调和多项式方法，当 $ p \equiv 1, q \equiv 0, u_0 \equiv 1 $ 时，退化为 $ \operatorname{Re}(z-z_0)^n, \operatorname{Im}(z-z_0)^n $。

实验结果

研究问题

• RQ1 在有界平面区域中，当系数满足可分离条件时，能否显式构造算子 $ D = \operatorname{div} p \operatorname{grad} + q $ 的完备零空间解系？
• RQ2 与传统基于网格划分的方法（如有限元法）相比，MPFP 的收敛速率与精度如何？
• RQ3 MPFP 能否有效适应求解椭圆型算子（尤其是亥姆霍兹方程）的特征值问题？
• RQ4 函数 $ f = p^{1/2} u_0 $ 的可分离性在形式幂构造与解系完备性中起何作用？
• RQ5 在实际应用中，边界配点矩阵 $ U $ 的行列式在多大程度上可可靠用于计算特征值？

主要发现

• MPFP 在求解边值问题时表现出高精度，最大绝对误差从 $ N = 4 $ 时的 0.079 降至 $ N = 32 $ 时的 $ 2.8 \times 10^{-10} $，表明收敛速度极快。
• 在单位圆盘中求解亥姆霍兹方程时，仅使用 $ N = 21 $ 即可将前五个特征值计算至四位小数精度，使用 $ N = 23 $ 可得六个特征值的相同精度。
• 与有限元法相比，该方法性能具有竞争力，尤其在可显式构造形式幂的问题中表现更优。
• 在给定 $ p $ 与 $ q $ 的条件下，伪解析形式幂系的完备性在形式上等价于解析函数理论中 $ (z - z_0)^n $ 的完备性。
• 该方法推广了经典调和多项式法，当 $ p \equiv 1 $, $ q \equiv 0 $, 且 $ u_0 \equiv 1 $ 时，可恢复标准调和多项式。
• 尽管在复值 $ u_0 $ 情况下存在理论挑战，数值实验表明该方法在特征值问题中依然有效，$ U $ 的行列式虽快速衰减，但仍能提供高精度的特征值近似。