中文

求解高波数亥姆霍兹方程的可杂交不连续伽辽金方法(第一部分:线性情形)

数值分析 2020-05-01 v1 数值分析

摘要

本文讨论了带阻抗边界条件的高频亥姆霍兹方程线性可杂交不连续伽辽金方法(HDG)的若干方面。首先,推导了HDG对精确解uu及其负梯度q=u\mathbf{q}=-\nabla u的逼近误差估计,显式依赖于波数kk。在k3h2k^3h^2足够小且惩罚参数τk\tau\eqsim k(其中hh为网格尺寸)的条件下,证明了kuuhL2(Ω)+qqhL2(Ω)=O(k2h2+k4h3)k\Vert u - u_h \Vert_{L^2(\Omega)} + \Vert \mathbf{q} -\mathbf{q}_h \Vert_{L^2(\Omega)} = O(k^2h^2+k^4h^3)。注意qh\mathbf{q}_h的收敛阶是满阶且污染误差为O(k4h3)O(k^4h^3),改进了已有结果。其次,利用椭圆问题HDG方法的标准后处理过程,得到分片二次函数uhu_h^*使得kuuhL2(Ω)=O(k3h3+k4h3)k\Vert u-u_h^*\Vert_{L^2(\Omega)}=O(k^3h^3+k^4h^3)。注意该后处理仅改进插值误差(从O(k2h2)O(k^2h^2)O(k3h3)O(k^3h^3))而污染误差O(k4h3)O(k^4h^3)不变。第三,频散分析和大量数值试验表明,通过选取适当惩罚参数,在一维情形可完全消除污染效应,在二维情形可大幅减弱。高波数亥姆霍兹方程高阶HDG方法的预渐近误差分析在第二部分研究。

关键词

引用

@article{arxiv.2004.14553,
  title  = {Hybridizable Discontinuous Galerkin Methods for Helmholtz Equation with High Wave Number. Part I: Linear case},
  author = {Bingxin Zhu and Haijun Wu},
  journal= {arXiv preprint arXiv:2004.14553},
  year   = {2020}
}

备注

19 pages, 16 figures