A new numerical method for div-curl systems with low regularity assumptions

This paper presents a new numerical method for div-curl systems with the normal boundary condition by using a finite element technique known as primal-dual weak Galerkin (PDWG). The PDWG finite element scheme for the div-curl system has two prominent features in that it offers not only an accurate and reliable numerical solution to the div-curl system under the low Hα-regularity (α>0) assumption for the true solution, but also an effective approximation of the normal harmonic vector fields on domains with complex topology. Seven numerical experiments are conducted and the results are presented to demonstrate the performance of the PDWG algorithm, including one example on the computation of discrete normal harmonic vector fields.