Dispersion analysis of plane wave discontinuous Galerkin methods

SUMMARY The plane wave DG (PWDG) method for the Helmholtz equation was introduced and analyzed in [GITTELSON, C., HIPTMAIR, R., AND PERUGIA, I. Plane wave discontinuous Galerkin methods: analysis of the h-version. Math. Model. Numer. Anal. 43 (2009), 297–331] as a generalization of the so-called ultra-weak variational formulation, see [O. CESSENAT AND B. DESPReS, Application of an ultra weak variational formulation of elliptic PDEs to the two-dimensional Helmholtz equation, SIAM J. Numer. Anal., 35 (1998), pp. 255–299]. The method relies on Trefftz-type local trial spaces spanned by plane waves of different directions, and links cells of the mesh through numerical fluxes in the spirit of DG methods. We conduct a partly empirical dispersion analysis of the method in a discrete translation-invariant setting by studying the mismatch of wave numbers of discrete and continuous plane waves traveling in the same direction. We find agreement of the wave numbers for directions represented in the local trial spaces. For other directions, the PWDG methods turn out to incur both phase and amplitude errors. This manifests itself as a pollution effect haunting the h-version of the method. Our dispersion analysis allows a quantitative prediction of the strength of this effect and its dependence on the wave number and number of plane waves. Copyright © 2014 John Wiley & Sons, Ltd.

[1]  Mark Ainsworth,et al.  Dispersive and dissipative behaviour of high order discontinuous Galerkin finite element methods , 2004 .

[2]  I. Babuska,et al.  The Partition of Unity Method , 1997 .

[3]  Peter Monk,et al.  Improvements for the ultra weak variational formulation , 2013 .

[4]  J. Kaipio,et al.  Computational aspects of the ultra-weak variational formulation , 2002 .

[5]  Charbel Farhat,et al.  Convergence Analysis of a Discontinuous Galerkin Method with Plane Waves and Lagrange Multipliers for the Solution of Helmholtz Problems , 2009, SIAM J. Numer. Anal..

[6]  U. Hetmaniuk,et al.  Multiple-stencil dispersion analysis of the Lagrange multipliers in a discontinuous Galerkin method for the Helmholtz equation , 2003 .

[7]  I. Babuska,et al.  Dispersion and pollution of the FEM solution for the Helmholtz equation in one, two and three dimensions , 1999 .

[8]  Stefan A. Sauter,et al.  Is the Pollution Effect of the FEM Avoidable for the Helmholtz Equation Considering High Wave Numbers? , 1997, SIAM Rev..

[9]  R. J. Astley,et al.  Modelling of short wave diffraction problems using approximating systems of plane waves , 2002 .

[10]  R. J. Astley,et al.  The partition of unity finite element method for short wave acoustic propagation on non‐uniform potential flows , 2006 .

[11]  Charbel Farhat,et al.  A discontinuous Galerkin method with Lagrange multipliers for the solution of Helmholtz problems in the mid-frequency regime , 2003 .

[12]  Ralf Hiptmair,et al.  PLANE WAVE DISCONTINUOUS GALERKIN METHODS: ANALYSIS OF THE h-VERSION ∗, ∗∗ , 2009 .

[13]  R. Hiptmair,et al.  Trefftz discontinuous Galerkin methods for acoustic scattering on locally refined meshes , 2014 .

[14]  I. Babuska,et al.  Finite element solution of the Helmholtz equation with high wave number Part I: The h-version of the FEM☆ , 1995 .

[15]  Dispersion analysis in scalar generalized finite element method , 2008, 2008 IEEE Antennas and Propagation Society International Symposium.

[16]  R. J. Astley,et al.  Special short wave elements for flow acoustics , 2005 .

[17]  Ralf Hiptmair,et al.  STABILITY RESULTS FOR THE TIME-HARMONIC MAXWELL EQUATIONS WITH IMPEDANCE BOUNDARY CONDITIONS , 2011 .

[18]  Ralf Hiptmair,et al.  Plane Wave Discontinuous Galerkin Methods for the 2D Helmholtz Equation: Analysis of the p-Version , 2011, SIAM J. Numer. Anal..

[19]  Bruno Després,et al.  Using Plane Waves as Base Functions for Solving Time Harmonic Equations with the Ultra Weak Variational Formulation , 2003 .

[20]  Peter Monk,et al.  Solving Maxwell's equations using the ultra weak variational formulation , 2007, J. Comput. Phys..

[21]  Ivo Babuška,et al.  The generalized finite element method for Helmholtz equation: Theory, computation, and open problems , 2006 .

[22]  Mark Ainsworth,et al.  Discrete Dispersion Relation for hp-Version Finite Element Approximation at High Wave Number , 2004, SIAM J. Numer. Anal..

[23]  P Bettess,et al.  Short–wave scattering: problems and techniques , 2004, Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[24]  Mark Ainsworth,et al.  Dispersive and Dissipative Properties of Discontinuous Galerkin Finite Element Methods for the Second-Order Wave Equation , 2006, J. Sci. Comput..

[25]  T. Strouboulis,et al.  Partition of unity method for Helmholtz equation: q-convergence for plane-wave and wave-band local bases , 2006 .

[26]  Pierre Ladevèze,et al.  THE MULTISCALE VTCR APPROACH APPLIED TO ACOUSTICS PROBLEMS , 2008 .

[27]  M. Ainsworth.  Dispersive properties of high–order Nédélec/edge element approximation of the time–harmonic Maxwell equations , 2004, Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[28]  Charbel Farhat,et al.  Three‐dimensional discontinuous Galerkin elements with plane waves and Lagrange multipliers for the solution of mid‐frequency Helmholtz problems , 2006 .

[29]  O. Cessenat,et al.  Application of an Ultra Weak Variational Formulation of Elliptic PDEs to the Two-Dimensional Helmholtz Problem , 1998 .

[30]  Charbel Farhat,et al.  Higher‐order extensions of a discontinuous Galerkin method for mid‐frequency Helmholtz problems , 2004 .