Convergence and stability of a high-order leap-frog based discontinuous Galerkin method for the Maxwell equations on non-conforming meshes

Hassan Fahs (1) , Stephane Lanteri (1)

1 NACHOS - Numerical modeling and high performance computing for evolution problems in complex domains and heterogeneous media

In this work, we discuss the formulation, stability, convergence and numerical validation of a high-order leap-frog based non-dissipative discontinuous Galerkin time-domain (DGTD) method for solving Maxwell's equations on non-conforming simplicial meshes. This DGTD method makes use of a nodal polynomial interpolation method for the approximation of the electromagnetic field within a simplex, and of a centered scheme for the calculation of the numerical flux at an interface between neighboring elements. Moreover, the interpolation degree is defined at the element level and the mesh is refined locally in a non-conforming way resulting in arbitrary level hanging nodes. The method is proved to be stable and conserves a discrete analog of the electromagnetic energy for metallic cavities. The convergence of the semi-discrete approximation to Maxwell's equations is established rigorously and bounds on the global divergence error are provided. Numerical experiments with high-order elements show the potential of the method.

Mots clés: time-domain Maxwell's equations; discontinuous Galerkin method; leap-frog time scheme; non-conforming meshes

Domaines: Analyse numérique [math.NA]; Equations aux dérivées partielles [math.AP]; Modélisation et simulation

Hassan Fahs, Stephane Lanteri. Convergence and stability of a high-order leap-frog based discontinuous Galerkin method for the Maxwell equations on non-conforming meshes. [Research Report] RR-6699, INRIA. 2008, pp.27. ⟨inria-00332277⟩