[论文解读] Estimate of the Convergence Rate of Finite Element Solutions to Elliptic Equations of Second Order with Discontinuous Coefficients

Jinchao Xu|arXiv (Cornell University)|Nov 17, 2013

Advanced Mathematical Modeling in Engineering参考文献 1被引用 29

一句话总结

该论文为二维二阶椭圆问题在系数不连续情况下应用线性有限元方法建立了最优误差估计。通过基于变点泰勒展开和加权积分估计的创新分析方法，证明了收敛速率 ‖u−uₕ‖₁,Ω ≤ Ch|ln h|¹ᐟ²‖u‖₂,Ω₁₊Ω₂，表明不连续系数仅轻微降低收敛速度，与光滑情况相比影响甚微。

ABSTRACT

In this paper, we consider elliptic boundary value problems with discontinuous coefficients and obtain the asymptotic optimal error estimate $\|u-u_k\|_{1,Ω}\leqslant Ch|\ln h|^{1/2}\|u\|_{2,Ω_1+Ω_2}$ for triangle linear elements.

研究动机与目标

• 分析具有不连续系数的二阶椭圆PDE的线性有限元方法的收敛行为。
• 量化系数不连续性对二维有限元误差估计的影响。
• 推导一个渐近最优的误差界，该界考虑了分隔子区域 Ω₁ 和 Ω₂ 的分段光滑界面 S 的存在。
• 建立一个在系数跳跃导致全局 H² 正则性缺失的情况下仍保持近乎最优的收敛速率。

提出的方法

• 采用带不连续系数 B 和低阶项 σ 的模型问题的变分形式，通过双线性形式 a(u,v) 定义。
• 在尊重界面 S 的拟均匀三角剖分上使用连续分段线性函数的有限元空间 Sₕ。
• 采用基于变基点泰勒展开的新颖误差分析技术，以处理与界面 S 相交的非规则单元。
• 引入变量变换和加权积分估计，以控制梯度在界面附近的性质。
• 利用两个关键引理：一个用于奇异积分的 L² 估计（各向异性核），另一个用于通过测度依赖的索博列夫型嵌入在子集上估计 L² 范数。
• 通过选择 ε = 1/(2|ln h|) 优化误差界，以最小化对 h 和 |ln h| 的依赖，从而得出最终的收敛速率。

实验结果

研究问题

• RQ1 不连续系数的存在如何影响二维线性有限元方法的收敛速率？
• RQ2 当精确解因系数跳跃而缺乏全局 H² 正则性时，是否可以为 H¹-合条件有限元推导出最优误差估计？
• RQ3 当界面 S 为分段光滑且系数 B 在 S 上不连续时，误差对网格尺寸 h 的精确依赖关系是什么？
• RQ4 当系数不连续时，标准有限元误差估计在多大程度上失效或退化？

主要发现

• 有限元解 uₕ 在 H¹-范数下以速率 ‖u−uₕ‖₁,Ω ≤ Ch|ln h|¹ᐟ²‖u‖₂,Ω₁₊Ω₂ 收敛于精确解 u。
• 对数因子 |ln h|¹ᐟ² 源于对界面 S 附近非规则单元的分析，且在给定假设下被证明为最优。
• 尽管存在系数 B 的不连续性，误差估计仍保持渐近最优，表明与光滑系数情况相比，收敛速度仅轻微退化。
• 该分析证明，即使 u 不属于 H²(Ω)，H¹-误差仍被有界于常数乘以 h|ln h|¹ᐟ² 再乘以 u 在 Ω₁ 和 Ω₂ 上的 H² 半范数。
• 该方法通过将误差分解为规则和非规则单元，并对每一类应用定制化估计，成功处理了全局 H² 正则性的缺失。