非凸多边形障碍域上Helmholtz方程线性CIP-FEM的波数显式拟最优误差分析

《Mathematics》：Wavenumber-Explicit Quasi-Optimal Error Analysis of Linear CIP-FEM for the Helmholtz Equation on Nonconvex Polygonal Obstacle Domains

　　

研究人员研究了二维非凸多边形障碍域上带混合Dirichlet-阻抗边界条件的Helmholtz方程的线性连续内罚有限元方法(CIP-FEM)。凹角将全局正则性降至H²(Ω)以下，要求对CIP稳定项进行对角敏感的处理。利用波数显式正则-奇异分解，研究人员在降正则性下建立了CIP格式的精确相容性。还证明了截断角奇异函数的Scott–Zhang拟插值的关键阶惩罚半范估计|I_hS_j|_J ? h^α_j。将该估计与端点Scott–Zhang逼近以及正则和奇异分量的波数显式界相结合，得到inf_vh∈Vh‖u?v_h‖_h,k ? (kh + k^α?1/2h^α)‖f‖_0,Ω。随后，在充分分辨条件k²h≤δ下，Schatz型对偶论证得到相应的拟最优误差估计‖u?u_h‖_h,k ? (kh + k^α?1/2h^α)‖f‖_0,Ω，其中δ与k和h无关。离散唯一性与该条件无关地被证明。数值实验支持所预测的角奇异性行为；对于所测问题，在固定惩罚值γ=0.1下，CIP-FEM在高波数预渐近区的若干相对较粗网格上比标准FEM误差更小。本分析未单独识别加性污染误差项。

研究背景与问题。高频Helmholtz方程的有限元逼近中，波长随波数k增大而变短；标准有限元法(FEM)的误差可能包含不受kh≤C控制的污染分量，因此需要波数显式误差界。经典波数显式框架通常基于H²或更高正则性；然而当计算域是带凹角的非凸多边形障碍域时，凹角使全局正则性低于H²(Ω)，标准Galerkin正交性和插值论证不再直接适用，CIP惩罚项也需要在分数正则性下重新验证相容性。已有CIP-FEM分析多针对凸域或正则性足够好的情形，未系统处理物理凹角带来的低正则性。为此，研究人员在《Mathematics》发表的论文中，针对二维混合Dirichlet-阻抗边界条件的非凸多边形障碍域，建立了线性CIP-FEM的波数显式拟最优误差分析。该研究利用波数显式正则-奇异分解，在低正则性下证明了CIP格式的精确相容性，并给出了角敏感的关键阶惩罚估计；结合Scott–Zhang拟插值和Schatz型对偶论证，在充分分辨条件k²h≤δ下得到拟最优误差估计，且离散唯一性不需要该分辨条件。数值实验支持理论预测的角奇异性行为，并显示固定正惩罚γ=0.1在部分较粗网格上优于标准FEM。

关键技术方法（不超过250字）。本研究不涉及生物样本队列；数值实验采用截断Bessel基准、正方形障碍(α=2/3)和正六边形障碍(α=3/4)。主要技术方法包括：①波数显式正则-奇异分解，分离正则分量与截断角奇异函数并获得波数显式系数界；②保持Dirichlet边界条件的Scott–Zhang拟插值；③角奇异函数CIP跳跃临界阶估计，结合角点附近局部H¹稳定性和二进环带求和；④Schatz型对偶论证与网格依赖范数控制；⑤离散唯一性由同质CIP形式虚部独立证明。

研究结果。

6.1 基准验证。研究人员采用Chaumont-Frelet-Nicaise的截断Bessel基准，构造满足齐次Dirichlet与阻抗边界条件的精确解，比较拟最优插值、标准FEM和CIP-FEM的误差。数值结果显示收敛速率在加密过程中逐渐下降，最细三层网格上的最小二乘拟合斜率高于理论角指数2/3，说明尚未完全进入渐近区，因此该实验仅作为实现验证基准。

6.2 多边形障碍域上的角效应。对正方形障碍，凹角开度对应α=2/3；对正六边形障碍，对应α=3/4。研究人员在选定凹角处构造带径向截断的精确解，统计全局、角区域和过渡区域的H¹半范误差。正方形障碍角区域最细三层的拟合斜率接近2/3，过渡区域仍接近一阶；正六边形障碍角区域呈现向3/4过渡的非单调行为，过渡区域误差仍保持近一阶。结果支持预测的角奇异性控制收敛率。

6.3 高波数预渐近行为。回到正方形障碍，固定γ=0.1的CIP-FEM与γ=0的标准FEM比较相对H¹半范误差。在较粗网格上CIP-FEM误差小于标准FEM，最细网格上优势消失。波数从50增至100时，CIP-FEM误差增幅小于标准FEM。惩罚参数敏感性测试显示，在固定预渐近网格上多个正惩罚值均优于标准FEM，且γ从0.01增至0.1后结果几乎不变。该比较仅为描述性观察，不推断独立污染误差项。

讨论与结论。作者指出，本文分析未单独识别加性污染误差项；k²h≤δ只是充分条件而非必要条件。数值上观察到的优势是经验性的，FEM与CIP-FEM曲线交叉点位置可能依赖惩罚参数、几何、网格族和解，不应解释为最优惩罚选择。分析限于线性元、常系数、拟均匀贴体剖分、正一致有界惩罚参数和角指数α_j<1；拟均匀网格不能在凹角集中加密，因此未反映分级或