[论文解读] IP-DGFEM method for the $p(x)$- Laplacian

Leandro M. Del Pezzo, Ariel L. Lombardi+1 more|arXiv (Cornell University)|Sep 10, 2010

Advanced Numerical Methods in Computational Mathematics参考文献 13被引用 1

一句话总结

本文提出了一种针对变指数 p(x)-Laplacian 的内部惩罚不连续伽辽金有限元方法（IP-DGFEM），其中 p(x) 为对数 Hölder 连续函数，且满足 1 < p₁ ≤ p₂ < ∞。该方法证明了离散极小化器收敛于真实解，并在 p₁ 接近 1 时的 1D 测试中表现出优于协调伽辽金方法的性能，尤其适用于图像处理应用。

ABSTRACT

In this paper we construct an Interior Penalty Discontinuous Galerkin method to approximate the minimizer of a variational problem related to the $p(x)-$Laplacian. The function $p:\Omega o [p_1,p_2]$ is log Holder continuous and $1<p_1\leq p_2<\infty$. We prove that the minimizers of the discrete functional converge to the solution. We also make some numerical experiments in dimension one to compare this method with the Conforming Galerkin Method, in the case where $p_1$ is close to one. This example is motivated by its applications to image processing.

研究动机与目标

• 开发一种针对变指数 p(x) 的稳健数值方法，其中 p(x) 为对数 Hölder 连续函数。
• 建立离散极小化器收敛于变分问题真实解的理论结果。
• 在一维情形下将 IP-DGFEM 与协调伽辽金方法进行比较，尤其关注当 p₁ 接近 1 时的表现。
• 评估该方法在图像处理相关应用中的性能，其中 p(x) 接近 1 时可建模为保边正则化。

提出的方法

• 构建一种内部惩罚不连续伽辽金方法，用于离散与 p(x)-Laplacian 相关的变分泛函。
• 通过在单元界面施加罚项，以弱形式强制连续性，从而允许使用不连续的有限元空间。
• 在有限维不连续伽辽金空间上最小化离散泛函，罚参数的选择确保了稳定性与收敛性。
• 收敛性分析依赖于 p(x) 的对数 Hölder 连续性以及离散泛函的强制性。
• 在一维情形下采用分段多项式逼近进行数值实验，比较 IP-DGFEM 与协调伽辽金方法的性能。
• 测试案例聚焦于 p₁ 接近 1 的参数区间，模拟保边图像复原问题。

实验结果

研究问题

• RQ1 在 p(x) 满足对数 Hölder 连续性的假设下，IP-DGFEM 是否收敛于 p(x)-Laplacian 变分问题的真实极小化器？
• RQ2 当 p₁ 接近 1 时，IP-DGFEM 在数值表现上是否优于协调伽辽金方法？
• RQ3 IP-DGFEM 是否能有效处理 p(x)-Laplacian 的非线性与变指数结构？
• RQ4 采用内部惩罚的不连续有限元空间对解的精度与稳定性有何影响？

主要发现

• 在假设 p(x) 为对数 Hölder 连续的条件下，IP-DGFEM 的离散极小化器收敛于 p(x)-Laplacian 问题的真实解。
• 在一维测试中，当 p₁ 接近 1 时，IP-DGFEM 表现出优于协调伽辽金方法的稳定性和精度。
• 该方法能有效捕捉 p(x)-Laplacian 在图像处理相关参数区间的特性，其中 p(x) 接近 1 时可实现边缘保持。
• 数值实验验证了理论收敛性，并凸显了 IP-DGFEM 在变指数问题中的鲁棒性。