An over-penalized weak Galerkin method for parabolic interface problems with time-dependent coefficients

Abstract

Based on the idea of classical discontinuous Galerkin and weak Galerkin finite element methods, we introduce an over-penalized term in the new scheme as a part of stabilization for solving parabolic interface problems. From the double-valued functions defined on interior edges of elements, it is natural to generate jumps of the over-penalized term. An over-penalized weak Galerkin (OPWG) finite element method can be applied very well to interface problems with general imperfect interface. Importantly, the diffusion coefficients of the interface problems depend on both temporal and spacial variables, not only limited in space as usual. With the use of (Pk(K),Pk−1(e),[Pk−1(K)]d) elements, semi-discrete and fully discrete schemes with backward Euler approximation in time are presented, and then the semi-discrete one is analyzed to be unconditionally stable. To analyze error estimates of semi-discrete and fully discrete schemes directly by error equation, we can just have optimal convergence order in energy norm. By virtue of the introduction of an elliptic projection operator, optimal error estimates of those schemes in L2 norm can be proved. Numerical examples are given to validate the efficiency and optimal convergence orders of the new schemes.