The Gautschi time stepping scheme for edge finite element discretizations of the Maxwell equations

  • Authors:
  • M. A. Botchev;D. Harutyunyan;J. J. W. van der Vegt

  • Affiliations:
  • Department of Applied Mathematics, University of Twente, Enschede, The Netherlands;Department of Applied Mathematics, University of Twente, Enschede, The Netherlands;Department of Applied Mathematics, University of Twente, Enschede, The Netherlands

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

Quantified Score

Hi-index 31.45

Visualization

Abstract

For the time integration of edge finite element discretizations of the three-dimensional Maxwell equations, we consider the Gautschi cosine scheme where the action of the matrix function is approximated by a Krylov subspace method. First, for the space-discretized edge finite element Maxwell equations, the dispersion error of this scheme is analyzed in detail and compared to that of two conventional schemes. Second, we show that the scheme can be implemented in such a way that a higher accuracy can be achieved within less computational time (as compared to other implicit schemes). We also analyzed the error made in the Krylov subspace matrix function evaluations. Although the new scheme is unconditionally stable, it is explicit in structure: as an explicit scheme, it requires only the solution of linear systems with the mass matrix.