-Discontinuous Galerkin methods for the Helmholtz equation with large wave number

Authors:
Xiaobing Feng and Haijun Wu

Journal:
Math. Comp. **80** (2011), 1997-2024

MSC (2010):
Primary 65N12, 65N15, 65N30, 78A40

Published electronically:
February 25, 2011

MathSciNet review:
2813347

Full-text PDF

Abstract | References | Similar Articles | Additional Information

Abstract: In this paper we develop and analyze some interior penalty -discontinuous Galerkin (-DG) methods for the Helmholtz equation with first order absorbing boundary condition in two and three dimensions. The proposed -DG methods are defined using a sesquilinear form which is not only mesh-dependent (or -dependent) but also degree-dependent (or -dependent). In addition, the sesquilinear form contains penalty terms which not only penalize the jumps of the function values across the element edges but also the jumps of the first order tangential derivatives as well as jumps of all normal derivatives up to order . Furthermore, to ensure the stability, the penalty parameters are taken as complex numbers with positive imaginary parts, so essentially and practically no constraint is imposed on the penalty parameters. It is proved that the proposed -discontinuous Galerkin methods are stable (hence, well-posed) without any mesh constraint. For each fixed wave number , sub-optimal order (with respect to and ) error estimates in the broken -norm and the -norm are derived without any mesh constraint. The error estimates as well as the stability estimates are improved to optimal order under the mesh condition by utilizing these stability and error estimates and using a stability-error iterative procedure, where is some constant independent of , , , and the penalty parameters. To overcome the difficulty caused by strong indefiniteness (and non-Hermitian nature) of the Helmholtz problems in the stability analysis for numerical solutions, our main ideas for stability analysis are to make use of a local version of the Rellich identity (for the Laplacian) and to mimic the stability analysis for the PDE solutions given in [19, 20, 33],

which enable us to derive stability estimates and error bounds with explicit dependence on the mesh size , the polynomial degree , the wave number , as well as all the penalty parameters for the numerical solutions.

**1.**M. Ainsworth, P. Monk, and W. Muniz,*Dispersive and dissipative properties of discontinuous Galerkin finite element methods for the second-order wave equation*, J. Sci. Comput.**27**(2006), no. 1-3, 5–40. MR**2285764**, 10.1007/s10915-005-9044-x**2.**Gustavo Benitez Alvarez, Abimael Fernando Dourado Loula, Eduardo Gomes Dutra do Carmo, and Fernando Alves Rochinha,*A discontinuous finite element formulation for Helmholtz equation*, Comput. Methods Appl. Mech. Engrg.**195**(2006), no. 33-36, 4018–4035. MR**2229833**, 10.1016/j.cma.2005.07.013**3.**Douglas N. Arnold,*An interior penalty finite element method with discontinuous elements*, SIAM J. Numer. Anal.**19**(1982), no. 4, 742–760. MR**664882**, 10.1137/0719052**4.**Douglas N. Arnold, Franco Brezzi, Bernardo Cockburn, and L. Donatella Marini,*Unified analysis of discontinuous Galerkin methods for elliptic problems*, SIAM J. Numer. Anal.**39**(2001/02), no. 5, 1749–1779. MR**1885715**, 10.1137/S0036142901384162**5.**A. K. Aziz and R. B. Kellogg,*A scattering problem for the Helmholtz equation*, Advances in computer methods for partial differential equations, III (Proc. Third IMACS Internat. Sympos., Lehigh Univ., Bethlehem, Pa., 1979), IMACS, New Brunswick, N.J., 1979, pp. 93–95. MR**603460****6.**A. K. Aziz and A. Werschulz,*On the numerical solutions of Helmholtz’s equation by the finite element method*, SIAM J. Numer. Anal.**17**(1980), no. 5, 681–686. MR**588754**, 10.1137/0717058**7.**I. Babuška and B. Q. Guo,*Approximation properties of the ℎ-𝑝 version of the finite element method*, Comput. Methods Appl. Mech. Engrg.**133**(1996), no. 3-4, 319–346. MR**1399640**, 10.1016/0045-7825(95)00946-9**8.**I. Babuška and Manil Suri,*The ℎ-𝑝 version of the finite element method with quasi-uniform meshes*, RAIRO Modél. Math. Anal. Numér.**21**(1987), no. 2, 199–238 (English, with French summary). MR**896241****9.**Garth A. Baker,*Finite element methods for elliptic equations using nonconforming elements*, Math. Comp.**31**(1977), no. 137, 45–59. MR**0431742**, 10.1090/S0025-5718-1977-0431742-5**10.**Gang Bao,*Finite element approximation of time harmonic waves in periodic structures*, SIAM J. Numer. Anal.**32**(1995), no. 4, 1155–1169. MR**1342287**, 10.1137/0732053**11.**Susanne C. Brenner and L. Ridgway Scott,*The mathematical theory of finite element methods*, Texts in Applied Mathematics, vol. 15, Springer-Verlag, New York, 1994. MR**1278258****12.**Erik Burman and Alexandre Ern,*Continuous interior penalty ℎ𝑝-finite element methods for advection and advection-diffusion equations*, Math. Comp.**76**(2007), no. 259, 1119–1140 (electronic). MR**2299768**, 10.1090/S0025-5718-07-01951-5**13.**C. L. Chang,*A least-squares finite element method for the Helmholtz equation*, Comput. Methods Appl. Mech. Engrg.**83**(1990), no. 1, 1–7. MR**1078694**, 10.1016/0045-7825(90)90121-2**14.**Eric T. Chung and Björn Engquist,*Optimal discontinuous Galerkin methods for wave propagation*, SIAM J. Numer. Anal.**44**(2006), no. 5, 2131–2158 (electronic). MR**2263043**, 10.1137/050641193**15.**Philippe G. Ciarlet,*The finite element method for elliptic problems*, North-Holland Publishing Co., Amsterdam-New York-Oxford, 1978. Studies in Mathematics and its Applications, Vol. 4. MR**0520174****16.**Bernardo Cockburn, George E. Karniadakis, and Chi-Wang Shu (eds.),*Discontinuous Galerkin methods*, Lecture Notes in Computational Science and Engineering, vol. 11, Springer-Verlag, Berlin, 2000. Theory, computation and applications; Papers from the 1st International Symposium held in Newport, RI, May 24–26, 1999. MR**1842160****17.**Bernardo Cockburn and Chi-Wang Shu,*The local discontinuous Galerkin method for time-dependent convection-diffusion systems*, SIAM J. Numer. Anal.**35**(1998), no. 6, 2440–2463 (electronic). MR**1655854**, 10.1137/S0036142997316712**18.**David Colton and Peter Monk,*The numerical solution of the three-dimensional inverse scattering problem for time harmonic acoustic waves*, SIAM J. Sci. Statist. Comput.**8**(1987), no. 3, 278–291. MR**883771**, 10.1137/0908035**19.**P. Cummings.*Analysis of Finite Element Based Numerical Methods for Acoustic Waves, Elastic Waves and Fluid-Solid Interactions in the Frequency Domain*.

PhD thesis, The University Tennessee, 2001.**20.**Peter Cummings and Xiaobing Feng,*Sharp regularity coefficient estimates for complex-valued acoustic and elastic Helmholtz equations*, Math. Models Methods Appl. Sci.**16**(2006), no. 1, 139–160. MR**2194984**, 10.1142/S021820250600108X**21.**Jim Douglas Jr. and Todd Dupont,*Interior penalty procedures for elliptic and parabolic Galerkin methods*, Computing methods in applied sciences (Second Internat. Sympos., Versailles, 1975) Springer, Berlin, 1976, pp. 207–216. Lecture Notes in Phys., Vol. 58. MR**0440955****22.**Jim Douglas Jr., Dongwoo Sheen, and Juan E. Santos,*Approximation of scalar waves in the space-frequency domain*, Math. Models Methods Appl. Sci.**4**(1994), no. 4, 509–531. MR**1291136**, 10.1142/S0218202594000297**23.**Björn Engquist and Andrew Majda,*Radiation boundary conditions for acoustic and elastic wave calculations*, Comm. Pure Appl. Math.**32**(1979), no. 3, 314–358. MR**517938**, 10.1002/cpa.3160320303**24.**Björn Engquist and Olof Runborg,*Computational high frequency wave propagation*, Acta Numer.**12**(2003), 181–266. MR**2249156**, 10.1017/S0962492902000119**25.**E. J. Kubatko, J. J. Westerink, and C. Dawson.

-discontinuous Galerkin methods for advection dominated problems in shallow water flow.*Comput. Methods Appl. Mech. Engrg.*, 196:437-451, 2006.**26.**Xiaobing Feng and Ohannes A. Karakashian,*Fully discrete dynamic mesh discontinuous Galerkin methods for the Cahn-Hilliard equation of phase transition*, Math. Comp.**76**(2007), no. 259, 1093–1117 (electronic). MR**2299767**, 10.1090/S0025-5718-07-01985-0**27.**Xiaobing Feng and Haijun Wu,*Discontinuous Galerkin methods for the Helmholtz equation with large wave number*, SIAM J. Numer. Anal.**47**(2009), no. 4, 2872–2896. MR**2551150**, 10.1137/080737538**28.**Emmanuil H. Georgoulis and Endre Süli,*Optimal error estimates for the ℎ𝑝-version interior penalty discontinuous Galerkin finite element method*, IMA J. Numer. Anal.**25**(2005), no. 1, 205–220. MR**2110241**, 10.1093/imanum/drh014**29.**C. I. Goldstein,*The finite element method with nonuniform mesh sizes applied to the exterior Helmholtz problem*, Numer. Math.**38**(1981/82), no. 1, 61–82. MR**634753**, 10.1007/BF01395809**30.**Benqi Guo,*Approximation theory for the 𝑝-version of the finite element method in three dimensions. I. Approximabilities of singular functions in the framework of the Jacobi-weighted Besov and Sobolev spaces*, SIAM J. Numer. Anal.**44**(2006), no. 1, 246–269 (electronic). MR**2217381**, 10.1137/040614803**31.**Benqi Guo and Weiwei Sun,*The optimal convergence of the ℎ-𝑝 version of the finite element method with quasi-uniform meshes*, SIAM J. Numer. Anal.**45**(2007), no. 2, 698–730. MR**2300293**, 10.1137/05063756X**32.**Isaac Harari and Thomas J. R. Hughes,*Analysis of continuous formulations underlying the computation of time-harmonic acoustics in exterior domains*, Comput. Methods Appl. Mech. Engrg.**97**(1992), no. 1, 103–124. MR**1182436**, 10.1016/0045-7825(92)90109-W**33.**U. Hetmaniuk,*Stability estimates for a class of Helmholtz problems*, Commun. Math. Sci.**5**(2007), no. 3, 665–678. MR**2352336****34.**Paul Houston, Christoph Schwab, and Endre Süli,*Discontinuous ℎ𝑝-finite element methods for advection-diffusion-reaction problems*, SIAM J. Numer. Anal.**39**(2002), no. 6, 2133–2163. MR**1897953**, 10.1137/S0036142900374111**35.**Paul Houston, Max Jensen, and Endre Süli,*ℎ𝑝-discontinuous Galerkin finite element methods with least-squares stabilization*, Proceedings of the Fifth International Conference on Spectral and High Order Methods (ICOSAHOM-01) (Uppsala), 2002, pp. 3–25. MR**1910549**, 10.1023/A:1015180009979**36.**Frank Ihlenburg and Ivo Babuška,*Finite element solution of the Helmholtz equation with high wave number. II. The ℎ-𝑝 version of the FEM*, SIAM J. Numer. Anal.**34**(1997), no. 1, 315–358. MR**1445739**, 10.1137/S0036142994272337**37.**J. M. Melenk and S. Sauter,*Convergence analysis for finite element discretizations of the Helmholtz equation with Dirichlet-to-Neumann boundary conditions*, Math. Comp.**79**(2010), no. 272, 1871–1914. MR**2684350**, 10.1090/S0025-5718-10-02362-8**38.**J.M. Melenk and S. Sauter.

Wave-Number Explicit Convergence Analysis for Galerkin Discretizations of the Helmholtz Equation (extended version).*Preprint 09-2009, University of Zurich*, (2009).**39.**I. Perugia.

A note on the discontinuous Galerkin approximation of the Helmholtz equation.*Lecture Notes*, ETH Zürich, 2006.**40.**Jie Shen and Li-Lian Wang,*Analysis of a spectral-Galerkin approximation to the Helmholtz equation in exterior domains*, SIAM J. Numer. Anal.**45**(2007), no. 5, 1954–1978. MR**2346366**, 10.1137/060665737**41.**Nikolaos A. Kampanis, Vassilios A. Dougalis, and John A. Ekaterinaris (eds.),*Effective computational methods for wave propagation*, Numerical Insights, vol. 5, Chapman & Hall/CRC, Boca Raton, FL, 2008. MR**2406004****42.**Béatrice Rivière, Mary F. Wheeler, and Vivette Girault,*Improved energy estimates for interior penalty, constrained and discontinuous Galerkin methods for elliptic problems. I*, Comput. Geosci.**3**(1999), no. 3-4, 337–360 (2000). MR**1750076**, 10.1023/A:1011591328604**43.**Ch. Schwab,*𝑝- and ℎ𝑝-finite element methods*, Numerical Mathematics and Scientific Computation, The Clarendon Press, Oxford University Press, New York, 1998. Theory and applications in solid and fluid mechanics. MR**1695813****44.**Pavel Šolín, Karel Segeth, and Ivo Doležel,*Higher-order finite element methods*, Studies in Advanced Mathematics, Chapman & Hall/CRC, Boca Raton, FL, 2004. With 1 CD-ROM (Windows, Macintosh, UNIX and LINUX). MR**2000261****45.**B. Stamm and T. P. Wihler.*-optimal discontinuous Galerkin methods for linear elliptic problems*. CMCS-REPORT-2007-006.**46.**T. Warburton and J. S. Hesthaven,*On the constants in ℎ𝑝-finite element trace inverse inequalities*, Comput. Methods Appl. Mech. Engrg.**192**(2003), no. 25, 2765–2773. MR**1986022**, 10.1016/S0045-7825(03)00294-9**47.**Mary Fanett Wheeler,*An elliptic collocation-finite element method with interior penalties*, SIAM J. Numer. Anal.**15**(1978), no. 1, 152–161. MR**0471383****48.**O. C. Zienkiewicz,*Achievements and some unsolved problems of the finite element method*, Internat. J. Numer. Methods Engrg.**47**(2000), no. 1-3, 9–28. Richard H. Gallagher Memorial Issue. MR**1744287**, 10.1002/(SICI)1097-0207(20000110/30)47:1/3<9::AID-NME793>3.0.CO;2-P

Retrieve articles in *Mathematics of Computation*
with MSC (2010):
65N12,
65N15,
65N30,
78A40

Retrieve articles in all journals with MSC (2010): 65N12, 65N15, 65N30, 78A40

Additional Information

**Xiaobing Feng**

Affiliation:
Department of Mathematics, The University of Tennessee, Knoxville, Tennessee 37996

Email:
xfeng@math.utk.edu

**Haijun Wu**

Affiliation:
Department of Mathematics, Nanjing University, Nanjing, Jiangsu 210093, People’s Republic of China

Email:
hjw@nju.edu.cn

DOI:
https://doi.org/10.1090/S0025-5718-2011-02475-0

Keywords:
Helmholtz equation,
time harmonic waves,
absorbing boundary conditions,
$hp$-discontinuous Galerkin methods,
error estimates

Received by editor(s):
July 20, 2009

Received by editor(s) in revised form:
August 22, 2010

Published electronically:
February 25, 2011

Additional Notes:
The work of the first author was partially supported by the NSF grants DMS-0410266 and DMS-0710831.

The work of the second author was partially supported by the national basic research program of China under grant 2005CB321701 and by the program for the new century outstanding talents in universities of China. Part of this joint work was carried out during the second author’s recent visit to the University of Tennessee, the author would like to thank the host institution for its hospitality and financial support of the visit.

Article copyright:
© Copyright 2011
American Mathematical Society

The copyright for this article reverts to public domain 28 years after publication.