A Weak Galerkin Finite Element Method for the Maxwell Equations

This paper introduces a numerical scheme for the time-harmonic Maxwell equations by using weak Galerkin (WG) finite element methods. The WG finite element method is based on two operators: discrete weak curl and discrete weak gradient, with appropriately defined stabilizations that enforce a weak continuity of the approximating functions. This WG method is highly flexible by allowing the use of discontinuous approximating functions on arbitrary shape of polyhedra and, at the same time, is parameter free. Optimal-order of convergence is established for the WG approximations in various discrete norms which are either $$H^1$$H1-like or $$L^2$$L2 and $$L^2$$L2-like. An effective implementation of the WG method is developed through variable reduction by following a Schur-complement approach, yielding a system of linear equations involving unknowns associated with element boundaries only. Numerical results are presented to confirm the theory of convergence.

[1]  Junping Wang,et al.  An efficient numerical scheme for the biharmonic equation by weak Galerkin finite element methods on polygonal or polyhedral meshes , 2013, Comput. Math. Appl..

[2]  Susanne C. Brenner,et al.  Hodge decomposition for divergence-free vector fields and two-dimensional Maxwell's equations , 2012, Math. Comput..

[3]  Ronald H. W. Hoppe,et al.  Finite element methods for Maxwell's equations , 2005, Math. Comput..

[4]  Junping Wang,et al.  Weak Galerkin finite element methods for the biharmonic equation on polytopal meshes , 2013, 1303.0927.

[5]  Thomas Rylander,et al.  Computational Electromagnetics , 2005, Electronics, Power Electronics, Optoelectronics, Microwaves, Electromagnetics, and Radar.

[6]  Ilaria Perugia,et al.  Mixed discontinuous Galerkin approximation of the Maxwell operator: The indefinite case , 2005 .

[7]  Junping Wang,et al.  A weak Galerkin finite element method for second-order elliptic problems , 2011, J. Comput. Appl. Math..

[8]  Carretera de Valencia,et al.  The finite element method in electromagnetics , 2000 .

[9]  Susanne C. Brenner,et al.  A Locally Divergence-Free Interior Penalty Method for Two-Dimensional Curl-Curl Problems , 2008, SIAM J. Numer. Anal..

[10]  Junping Wang,et al.  A weak Galerkin finite element method for the stokes equations , 2013, Adv. Comput. Math..

[11]  Ilaria Perugia,et al.  The hp-local discontinuous Galerkin method for low-frequency time-harmonic Maxwell equations , 2003, Math. Comput..

[12]  Wei Cai,et al.  Computational Methods for Electromagnetic Phenomena: Electrostatics in solvation , 2013 .

[13]  Junping Wang,et al.  A weak Galerkin mixed finite element method for second order elliptic problems , 2012, Math. Comput..

[14]  D. Schötzau,et al.  Stabilized interior penalty methods for the time-harmonic Maxwell equations , 2002 .

[15]  L. Demkowicz,et al.  hp-adaptive finite elements in electromagnetics , 1999 .

[16]  J. Nédélec.  Mixed finite elements in ℝ3 , 1980 .

[17]  D. Schötzau,et al.  hp -DGFEM for Maxwell’s equations , 2003 .