[论文解读] Convergence of the Point Integral method for the Poisson equation with Dirichlet boundary on point cloud

Zuoqiang Shi, Jian Sun|arXiv (Cornell University)|Dec 16, 2013

Numerical methods in engineering参考文献 13被引用 4

一句话总结

本文证明了在点云上求解带狄利克雷边界条件的泊松方程时，点积分法（PIM）的收敛性。通过采用罗宾型积分公式近似问题，并利用基于核的积分算子，该方法在适当的网格和核参数条件下实现了最优收敛速率。

ABSTRACT

The Poisson equation on manifolds plays an fundamental role in many applications. Recently, we proposed a novel numerical method called the Point Integral method (PIM) to solve the Poisson equations on manifolds from point clouds. In this paper, we prove the convergence of the point integral method for solving the Poisson equation with the Dirichlet boundary condition.

研究动机与目标

• 建立点积分法（PIM）在点云上求解带狄利克雷边界条件的泊松方程的收敛性。
• 解决将积分近似应用于狄利克雷问题的挑战，因为边界条件不提供法向导数。
• 使用罗宾问题作为狄利克雷问题的正则化近似，从而能够使用积分公式。
• 推导离散PIM解在网格尺寸 $ h $、核参数 $ t $ 和边界正则性方面的误差界。
• 证明在适当的参数缩放下，解的离散 $ L^2 $ 和 $ H^1 $ 范数收敛于其连续对应量。

提出的方法

• 通过基于核的积分算子近似流形上的拉普拉斯-贝尔特拉米算子，其中涉及 $ R_t(\bx,\by) $ 和 $ \bar{R}_t(\bx,\by) $，后者通过径向基函数 $ R $ 定义。
• 利用恒等式 $ -\int_{\mathcal{M}} \Delta_{\mathcal{M}} u(\by) \bar{R}_t(\bx,\by) d\mu_\by \approx \frac{1}{t} \int_{\mathcal{M}} R_t(\bx,\by)(u(\bx)-u(\by)) d\mu_\by - 2 \int_{\partial\mathcal{M}} \bar{R}_t(\bx,\by) \frac{\partial u}{\partial \bn}(\by) d\tau_\by $，将泊松方程转化为积分方程。
• 使用罗宾型边界条件 $ u + \beta \frac{\partial u}{\partial \bn} = b $ 作为狄利克雷问题的正则化近似，其中 $ \beta $ 较小。
• 使用点云 $ \{\bp_i\} $ 及其权重 $ V_i $ 离散化积分方程，以近似流形测度 $ \mu $ 和边界测度 $ \tau $。
• 定义离散解向量 $ \bu $，并基于积分近似建立线性系统，其中核权重依赖于 $ |\bp_i - \bp_j|^2 $。
• 通过核衰减和网格正则性，建立连续与离散 $ L^2 $ 和 $ H^1 $ 范数之间差值的误差估计。

实验结果

研究问题

• RQ1 能否严格证明点积分法（PIM）在点云上求解带狄利克雷边界条件的泊松方程时具有收敛性？
• RQ2 PIM 的近似误差如何依赖于网格尺寸 $ h $、核参数 $ t $ 和边界正则性？
• RQ3 在什么 $ h $ 与 $ t $ 的最优比例下，可确保离散解收敛于真实解？
• RQ4 在基于核的积分方法背景下，狄利克雷问题能否通过罗宾型公式有效近似？
• RQ5 PIM 解的离散 $ L^2 $ 和 $ H^1 $ 范数与精确解的连续范数之间有何关系？

主要发现

• PIM 解的离散 $ L^2 $ 范数满足 $ \left| \|u_{t,h}\|_{L^2(\mathcal{M})}^2 - \sum_{i=1}^n u_i^2 V_i \right| \leq \frac{Ch}{t^{1/2}} \left( \|u_{t,h}\|_{H^1(\mathcal{M})}^2 + t^{3/2} \|f\|_{\infty}^2 \right) $，其中 $ \frac{h}{t^{1/2}} \leq \frac{1}{2} $。
• 边界上的离散 $ L^2 $ 范数满足 $ \left| \|u_{t,h}\|_{L^2(\partial\mathcal{M})}^2 - \sum_{l \in I_S} u_l^2 A_l \right| \leq \frac{Ch}{t} \left( \|u_{t,h}\|_{H^1(\mathcal{M})}^2 + t^{3/2} \|f\|_{\infty}^2 \right) $。
• 在条件 $ \frac{h}{t^{1/2}} \leq \frac{1}{2} $ 下，离散的类似 $ H^1 $ 范数被限制为 $ \left( \sum u_i^2 V_i \right)^{1/2} + t^{1/4} \left( \sum u_l^2 A_l \right)^{1/2} \leq C \|u_{t,h}\|_{H^1(\mathcal{M})} + C \sqrt{h} \, t^{3/4} \|f\|_{\infty} $。
• 当 $ h \to 0 $ 且 $ t \to 0 $ 时，若 $ \frac{h}{t^{1/2}} \to 0 $，则方法在 $ H^1(\mathcal{M}) $-范数下收敛，确保最优误差衰减。
• 误差估计在核参数 $ t $ 上是一致的，且收敛速率由 $ h $、$ t $ 和解的正则性之间的相互作用所控制。
• 分析证实，PIM 是一种在点云上求解狄利克雷问题的稳定且收敛的方法，即使在高维环境空间中亦成立。