Comparing DG and Nedelec finite element discretisations of the second-order time-domain Maxwell equation

D. Sarmany, Mikhail A. Bochev, Jacobus J.W. van der Vegt, J.G. Verwer

    Research output: Book/ReportReportProfessional

    47 Downloads (Pure)

    Abstract

    This article compares the discontinuous Galerkin finite element method (DG-FEM) with the $H(\mathrm{curl})$-conforming FEM in the discretisation of the second-order time-domain Maxwell equations with possibly nonzero conductivity term. While DG-FEM suffers from an increased number of degrees of freedom compared with $H(\mathrm{curl})$-conforming FEM, it has the advantage of a purely block-diagonal mass matrix. This means that, as long as an explicit time-integration scheme is used, it is no longer necessary to solve a linear system at each time step -- a clear advantage over $H(\mathrm{curl})$-conforming FEM. It is known that DG-FEM generally favours high-order methods whereas $H(\mathrm{curl})$-conforming FEM is more suitable for low-order ones. The novelty we provide in this work is a direct comparison of the performance of the two methods when hierarchic $H(\mathrm{curl})$-conforming basis functions are used up to polynomial order $p=3$. The motivation behind this choice of basis functions is its growing importance in the development of $p$- and $hp$-adaptive FEMs. The fact that we allow for nonzero conductivity requires special attention with regards to the time-integration methods applied to the semi-discrete systems. High-order polynomial basis warrants the use of high-order time-integration schemes, but existing high-order schemes may suffer from a too severe time-step stability restriction as result of the conductivity term. We investigate several alternatives from the point of view of accuracy, stability and computational work. Finally, we carry out a numerical Fourier analysis to study the dispersion and issipation properties of the semi-discrete DG-FEM scheme and several of the time-integration methods. It is instructive in our approach that 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.
    Original languageUndefined
    Place of PublicationEnschede
    PublisherUniversity of Twente, Faculty of Mathematical Sciences
    Number of pages23
    Publication statusPublished - 1 Dec 2009

    Publication series

    NameMemorandum / Department of Applied Mathematics
    PublisherDepartment of Applied Mathematics, University of Twente
    No.1912
    ISSN (Print)1874-4850
    ISSN (Electronic)1874-4850

    Keywords

    • second-order Maxwell wave equation
    • numerical time integration
    • Discontinuous Galerkin finite element method
    • MSC-65L06
    • METIS-264211
    • MSC-65M60
    • IR-68865
    • H(curl) conforming finite element methods
    • MSC-65M20
    • EWI-16951

    Cite this