Optimal error estimates for a discontinuous Galerkin method on curved boundaries with polygonal meshes

技术摘要：基于多边形网格的曲线边界间断伽辽金方法的优化误差估计

问题陈述
本文研究了在具有任意光滑曲线边界（$\Omega$）的定义域上，二位线性平流-扩散-反应方程的数值解问题。一个主要的挑战在于，当使用多边形网格（$\Omega_h$）近似物理定义域时，会发生收敛阶数的损失。当标准的间断伽辽金（DG）方法应用于此类多边形近似，并将边界条件直接施加在计算边界（$\partial\Omega_h$）上时，通常会出现次优收敛的问题。由于曲线物理边界与分片线性计算边界之间的不匹配所引入的几何误差，会污染解，使得精度无论使用何种多项式次数 $N$，最高只能达到二阶。

方法论
作者提出了并分析了 DG–ROD（间断伽辽金–离岸数据重构） 方法。该方法扩展了先前在有限元（FE）和有限体积（FV）方法中开发的策略。

• 网格与空间： 物理定义域使用标准的多边形（三角形）网格进行离散化。该方法利用破碎索伯列夫空间（broken Sobolev spaces）和间断分片多项式空间（$S_{h,N}$）。
• 边界条件处理： DG–ROD 方法并非将预设的边界条件直接施加在多边形边界 $\partial\Omega_h$ 上，而是采用多项式重构。对于每个边界单元，通过求解一个约束最小化问题来确定一个局部多项式，该多项式在满足特定物理边界 $\partial\Omega$ 上点处的精确边界条件的同时，使到数值解的距离最小化。
• 变分形式： 该方法结合了用于扩散项的对称内罚（SIP）方案和用于平流-反应项的上风方案。重构后的边界数据（$g_{ROD}$）被直接纳入原变分形式中，从而在弱形式内有效地修正了几何不匹配。
• 实现： 作者采用了全局形式，即重构过程与 DG 求解器在单个步骤中同时完成，而非采用迭代交替的方法。

主要贡献与理论分析
本文对 DG–ROD 方法进行了严谨的理论分析，建立了针对凸定义域和非凸定义域的存在性、唯一性以及最优误差估计。

1. 适定性（Well-Posedness）： 作者证明了双线性型的有界性，并通过 Banach–Nečas–Babuška 定理建立了 inf-sup 条件，确保了在网格尺寸 $h$ 足够小且惩罚参数 $\eta$ 足够大时，离散解的存在性与唯一性。
2. 凸定义域： 对于凸定义域，作者证明了伽辽金正交性成立。在解具有适当正则性假设的前提下，他们推导出：
   • 在 DG 范数（能量范数）下具有 $O(h^N)$ 的最优收敛阶。
   • 利用对偶理论，在 $L^2$ 范数下具有 $O(h^{N+1})$ 的最优收敛阶。
3. 非凸定义域： 在非凸设定下，几何近似破坏了标准的伽辽金正交性。作者引入了一个由将解扩展到更大光滑定义域 $\tilde{\Omega}$ 而产生的残差项。他们推导出的误差估计包含一个取决于几何不匹配的残差项。
   • 理论分析表明，对于 $N > 2$ 的非凸定义域，除非计算定义域以更高的精度（阶数 $q > 2$）近似物理定义域，或者解具有满足特定条件的正则扩展，否则可能无法实现 $L^2$ 范数下的最优收敛。
   • 分析强调，正交性的丧失是非凸情形下的主要技术挑战。

数值结果
通过三个基准几何形状（圆盘——凸、环形——非凸、玫瑰形定义域——非凸）的数值实验验证了理论发现。

• 对比： 结果将经典的 DG 方法（在 $\partial\Omega_h$ 上施加条件）与 DG–ROD 方法进行了比较。
• 性能： 经典 DG 方法始终表现出精度退化，其收敛速率受限于约二阶，且与多项式次数 $N$ 无关。
• DG–ROD 有效性： 对于所有测试的几何形状及各种平流-反应系数，DG–ROD 方法成功恢复了最优收敛速率（在 $L^2$ 范数下为 $O(h^{N+1})$），特别是在 $N=2, 3, 4$ 时。值得注意的是，即使在非凸的玫瑰形定义域中，该方法也实现了最优速率，作者将其归功于使用了允许正则扩展的光滑制造解，从而有效地绕过了针对一般非凸情形所推导出的理论最坏情况限制。

意义与主张
本文声称提供了对 DG–ROD 方法的首次严谨理论分析，将理论框架从有限元扩展到了间断伽基本方法。

• 恢复精度： 其主要意义在于证明了，只要通过 ROD 技术重构边界数据，就可以在简单的多边形网格上恢复最优收敛速率，而无需使用复杂的曲线单元或等参数映射。
• 鲁棒性： 数值结果证实了该方法在处理具有曲线边界的凸及非凸定义域时的鲁棒性。
• 对局限性的审慎说明： 作者承认，虽然针对非凸定义域的理论误差估计表明，在没有特定正则性条件下 $N > 2$ 时可能存在次优性，但通过使用光滑解的数值基准测试显示，该方法在实践中表现最优。他们指出，向非线性方程、时间相关问题以及其他边界条件（Neumann, Robin）的扩展仍是未来研究中开放的挑战。这项工作也暗示该方法可推广至三维问题，并引用了经典 Lagrange 有限元的相关工作。