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

Author:

Feng Xiaobing,Wu Haijun

Abstract

In this paper we develop and analyze some interior penalty h p hp -discontinuous Galerkin ( h p hp -DG) methods for the Helmholtz equation with first order absorbing boundary condition in two and three dimensions. The proposed h p hp -DG methods are defined using a sesquilinear form which is not only mesh-dependent (or h h -dependent) but also degree-dependent (or p p -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 p p . 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 h p hp -discontinuous Galerkin methods are stable (hence, well-posed) without any mesh constraint. For each fixed wave number k k , sub-optimal order (with respect to h h and p p ) error estimates in the broken H 1 H^1 -norm and the L 2 L^2 -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 k 3 h 2 p 2 C 0 k^3h^2p^{-2}\le C_0 by utilizing these stability and error estimates and using a stability-error iterative procedure, where C 0 C_0 is some constant independent of k k , h h , p p , 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 h h , the polynomial degree p p , the wave number k k , as well as all the penalty parameters for the numerical solutions.

Publisher

American Mathematical Society (AMS)

Subject

Applied Mathematics,Computational Mathematics,Algebra and Number Theory

Reference48 articles.

1. Dispersive and dissipative properties of discontinuous Galerkin finite element methods for the second-order wave equation;Ainsworth, M.;J. Sci. Comput.,2006

2. A discontinuous finite element formulation for Helmholtz equation;Alvarez, Gustavo Benitez;Comput. Methods Appl. Mech. Engrg.,2006

3. An interior penalty finite element method with discontinuous elements;Arnold, Douglas N.;SIAM J. Numer. Anal.,1982

4. Unified analysis of discontinuous Galerkin methods for elliptic problems;Arnold, Douglas N.;SIAM J. Numer. Anal.,2001

5. A scattering problem for the Helmholtz equation;Aziz, A. K.,1979

同舟云学术

1.学者识别学者识别

2.学术分析学术分析

3.人才评估人才评估

"同舟云学术"是以全球学者为主线,采集、加工和组织学术论文而形成的新型学术文献查询和分析系统,可以对全球学者进行文献检索和人才价值评估。用户可以通过关注某些学科领域的顶尖人物而持续追踪该领域的学科进展和研究前沿。经过近期的数据扩容,当前同舟云学术共收录了国内外主流学术期刊6万余种,收集的期刊论文及会议论文总量共计约1.5亿篇,并以每天添加12000余篇中外论文的速度递增。我们也可以为用户提供个性化、定制化的学者数据。欢迎来电咨询!咨询电话:010-8811{复制后删除}0370

www.globalauthorid.com

TOP

Copyright © 2019-2024 北京同舟云网络信息技术有限公司
京公网安备11010802033243号  京ICP备18003416号-3