[论文解读] A Numerical Study on the Weak Galerkin Method for the Helmholtz Equation with Large Wave Numbers

Lin Mu, Junping Wang+2 more|arXiv (Cornell University)|Nov 2, 2011

Advanced Numerical Methods in Computational Mathematics参考文献 40被引用 15

一句话总结

本文研究了用于求解大波数 Helmholtz 方程的弱伽辽金（WG）有限元方法，展示了其鲁棒性、高精度和易于实现的特性。该方法采用不连续多项式逼近并以弱形式定义梯度，在高波数下仍能实现稳定收敛，当 $ kh = 0.25 $ 时，相对 $ H^1 $ 误差几乎不增加，表明有效控制了污染误差。

ABSTRACT

Weak Galerkin (WG) refers to general finite element methods for partial differential equations in which differential operators are approximated by weak forms through the usual integration by parts. In particular, WG methods allow the use of discontinuous finite element functions in the algorithm design. One of such examples was recently introduced by Wang and Ye for solving second order elliptic problems. The goal of this paper is to apply the WG method of Wang and Ye to the Helmholtz equation with high wave numbers. Several test scenarios are designed for a numerical investigation on the accuracy, convergence, and robustness of the WG method in both inhomogeneous and homogeneous media over convex and non-convex domains. Our numerical experiments indicate that weak Galerkin is a finite element technique that is easy to implement, and provides very accurate and robust numerical solutions for the Helmholtz problem with high wave numbers.

研究动机与目标

• 评估弱伽辽金（WG）有限元方法在求解高波数 Helmholtz 方程时的性能。
• 评估 WG 方法在凸域与非凸域中，于均匀与非均匀介质下的精度、收敛性与鲁棒性。
• 将 WG 方法与现有不连续伽辽金方法进行比较，尤其关注参数敏感性与污染误差方面的表现。
• 研究网格加密与多项式阶次对高波数问题求解精度的影响。

提出的方法

• WG 方法对解使用不连续分段多项式函数，并通过分部积分构造梯度算子的弱形式。
• 为每个单元构造离散弱梯度，实现在不需跨单元边界连续性的前提下局部逼近微分算子。
• 变分形式在每个单元上以残差意义强制满足 Helmholtz 方程，同时采用罗宾型吸收边界条件以近似 Sommerfeld 辐射条件。
• 与标准不连续伽辽金格式不同，该方法无需罚参数，从而简化实现并减少参数调优。
• 数值实验采用分片常数和分片线性单元，测试在不同波数与网格尺寸下的收敛性。
• 通过计算相对 $ H^1 $ 误差来评估收敛速率与污染效应，保持 $ kh $ 恒定以检验波数鲁棒性。

实验结果

研究问题

• RQ1 弱伽辽金方法在大波数下是否能保持 Helmholtz 方程的精度与收敛性？
• RQ2 与罚函数型不连续伽辽金方法相比，WG 方法在误差与参数依赖性方面表现如何？
• RQ3 WG 方法在高波数问题中在多大程度上缓解了污染效应？
• RQ4 网格加密与多项式阶次如何影响 WG 方法对振荡型 Helmholtz 解的求解精度？

主要发现

• 当 $ k = 100 $ 时，WG 方法在 $ kh = 0.25 $ 条件下实现了稳定的收敛，相对 $ H^1 $ 误差保持较低水平，表明污染效应极小。
• 在 $ kh = 0.25 $ 条件下，收敛曲线几乎平坦且斜率很小，表明该方法对波数增加具有鲁棒性。
• 在 $ k = 100 $ 且 $ h = 1/200 $ 时，WG 方法与精确解高度一致，随着网格细化，可视化收敛效果持续改善。
• 在 $ k = 100 $ 且 $ h = 1/60 $ 条件下，使用分片线性单元时，WG 方法能准确捕捉解的振幅与高频振荡特性。
• 该方法无需罚参数，相较于标准 DG 方法实现更简单，同时达到相当的精度水平。
• 即使在 $ k = 240 $ 时，当 $ kh $ 保持恒定，误差随 $ k $ 单调增加，但其增长斜率仍很小，表明污染误差得到良好控制。