Time-integration methods for finite element discretisations of the second-order Maxwell equation

  • Authors:

  • D. SáRmáNy;M. A. Botchev;J. J. W. Van Der Vegt

  • Affiliations:

  • School of Computing, University of Leeds, LS2 9JT, Leeds, United Kingdom;Department of Applied Mathematics, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands;Department of Applied Mathematics, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands

  • Venue:
  • Computers & Mathematics with Applications
  • Year:
  • 2013

Quantified Score

Hi-index 0.09

Visualization

Abstract

This article deals with time integration for the second-order Maxwell equations with possibly non-zero conductivity in the context of the discontinuous Galerkin finite element method (DG-FEM) and the H(curl)-conforming FEM. For the spatial discretisation, hierarchic H(curl)-conforming basis functions are used up to polynomial order p=3 over tetrahedral meshes, meaning fourth-order convergence rate. A high-order polynomial basis often warrants the use of high-order time-integration schemes, but many well-known high-order schemes may suffer from a severe time-step stability restriction owing to the conductivity term. We investigate several possible time-integration methods from the point of view of accuracy, stability and computational work. We also carry out a numerical Fourier analysis to study the dispersion and dissipation properties of the semi-discrete DG-FEM scheme as well as the fully-discrete schemes with several of the time-integration methods. The dispersion and dissipation properties of the spatial discretisation and those of the time-integration methods are investigated separately, providing additional insight into the two discretisation steps.