A non-uniform basis order for the discontinuous Galerkin method of the 3D dissipative wave equation with perfectly matched layer

  • Authors:

  • T. Lähivaara;T. Huttunen

  • Affiliations:

  • Department of Physics and Mathematics, University of Eastern Finland, P.O. Box 1627, FI-70211 Kuopio, Finland;Department of Physics and Mathematics, University of Eastern Finland, P.O. Box 1627, FI-70211 Kuopio, Finland and Kuava Ltd., P.O. Box 1188, FI-70211 Kuopio, Finland

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2010

Quantified Score

Hi-index 31.45

Visualization

Abstract

In this study a discontinuous Galerkin method (DG) for solving the three-dimensional time-dependent dissipative wave equation is investigated. In the case of unbounded problems, the perfectly matching layer (PML) is used to truncate the computational domain. The aim of this work is to investigate a simple selection method for choosing the basis order for elements in the computational mesh in order to obtain a predetermined error level. The selection method studied here relies on the error estimates provided by Ainsworth [M. Ainsworth, Dispersive and dissipative behaviour of high order discontinuous Galerkin finite element methods, Journal of Computational Physics 198(1) (2004) 106-130]. The performance of the non-uniform basis is examined using numerical experiments. In the simulated model problems, a feasible method choosing the basis order for arbitrary sized elements is achieved. In simulations, the effect of dissipation and the choices of the PML parameters on the performance of the DG method are also investigated.