A Numerical Study on the Weak Galerkin Method for the Helmholtz Equation

被引:57
作者
Mu, Lin [1 ]
Wang, Junping [2 ]
Ye, Xiu [3 ]
Zhao, Shan [4 ]
机构
[1] Michigan State Univ, Dept Math, E Lansing, MI 48824 USA
[2] Natl Sci Fdn, Div Math Sci, Arlington, VA 22230 USA
[3] Univ Arkansas, Dept Math & Stat, Little Rock, AR 72204 USA
[4] Univ Alabama, Dept Math, Tuscaloosa, AL 35487 USA
基金
美国国家科学基金会;
关键词
Galerkin finite element methods; discrete gradient; Helmholtz equation; large wave numbers; weak Galerkin; FINITE-ELEMENT-METHOD; DISCONTINUOUS GALERKIN; WAVE-NUMBER; APPROXIMATION; SCATTERING; VERSION;
D O I
10.4208/cicp.251112.211013a
中图分类号
O4 [物理学];
学科分类号
0702 ;
摘要
A weak Galerkin (WG) method is introduced and numerically tested for the Helmholtz equation. This method is flexible by using discontinuous piecewise polynomials and retains the mass conservation property. At the same time, the WG finite element formulation is symmetric and parameter free. 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. Challenging problems with high wave numbers are also examined. Our numerical experiments indicate that the 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.
引用
收藏
页码:1461 / 1479
页数:19
相关论文
共 40 条
[1]   Discrete dispersion relation for hp-version finite element approximation at high wave number [J].
Ainsworth, M .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 2004, 42 (02) :553-575
[2]   DISPERSIVE AND DISSIPATIVE BEHAVIOR OF THE SPECTRAL ELEMENT METHOD [J].
Ainsworth, Mark ;
Wajid, Hafiz Abdul .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 2009, 47 (05) :3910-3937
[3]   Unified analysis of discontinuous Galerkin methods for elliptic problems [J].
Arnold, DN ;
Brezzi, F ;
Cockburn, B ;
Marini, LD .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 2002, 39 (05) :1749-1779
[4]   A Generalized Finite Element Method for solving the Helmholtz equation in two dimensions with minimal pollution [J].
Babuska, I ;
Ihlenburg, F ;
Paik, ET ;
Sauter, SA .
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 1995, 128 (3-4) :325-359
[5]  
Babuska IM, 2000, SIAM REV, V42, P451
[6]   Numerical solution of the Helmholtz equation with high wavenumbers [J].
Bao, G ;
Wei, GW ;
Zhao, S .
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, 2004, 59 (03) :389-408
[7]  
Brezzi F., 1991, Mixed and Hybrid Finite Element Methods, V15
[8]   Using plane waves as base functions for solving time harmonic equations with the ultra weak variational formulation [J].
Cessenat, O ;
Després, B .
JOURNAL OF COMPUTATIONAL ACOUSTICS, 2003, 11 (02) :227-238
[9]   Application of an ultra weak variational formulation of elliptic PDEs to the two-dimensional Helmholtz problem [J].
Cessenat, O ;
Despres, B .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 1998, 35 (01) :255-299
[10]   A Galerkin boundary element method for high frequency scattering by convex polygons [J].
Chandler-Wilde, S. N. ;
Langdon, S. .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 2007, 45 (02) :610-640