Discretization of div-curl systems by weak Galerkin finite element methods on polyhedral partitions (Q333219)

The authors present a new numerical scheme for div-curl systems defined in connected domains with heterogeneous media by using the weak Galerkin finite element method. Two types of boundary value problems are considered. The first type is the normal boundary condition whereas the second one is the tangential boundary condition. A new variational formulation is developed for the normal boundary value problem by using the Helmholtz decomposition which avoids the computation of functions in the harmonic fields. Both boundary value problems are reduced to a general saddle-point problem involving the curl and divergence operators, for which the weak Galerkin finite element method is employed and analyzed. One of the main results is that the numerical technique lies in the discretization of the divergence operator applied to vector fields in heterogeneous media. Some error estimates with optimal order are proved for the corresponding finite element approximations in several discrete Sobolev norms.