A Weak Galerkin Mixed Finite Element Method for Second Order Elliptic Equations on 2D Curved Domains

A Weak Galerkin Mixed Finite Element Method for Second Order Elliptic Equations on 2D Curved Domains

Year:    2022

Author:    Yi Liu, Wenbin Chen, Yanqiu Wang

Communications in Computational Physics, Vol. 32 (2022), Iss. 4 : pp. 1094–1128

Abstract

This article concerns the weak Galerkin mixed finite element method (WG-MFEM) for second order elliptic equations on 2D domains with curved boundary. The Neumann boundary condition is considered since it becomes the essential boundary condition in this case. It is well-known that the discrepancy between the curved physical domain and the polygonal approximation domain leads to a loss of accuracy for discretization with polynomial order $α>1.$ The purpose of this paper is two-fold. First, we present a detailed error analysis of the original WG-MFEM for solving problems on curved domains, which exhibits an $O(h^{1/2})$ convergence for all $α ≥ 1.$ It is a little surprising to see that even the lowest-order WG-MFEM $(α = 1)$ experiences a loss of accuracy. This is different from known results for the finite element method (FEM) or the mixed FEM, and appears to be a combined effect of the WG-MFEM design and the fact that the outward normal vector on the polygonal approximation domain is different from the one on the curved domain. Second, we propose a remedy to bring the approximation rate back to optimal by employing two techniques. One is a specially designed boundary correction technique. The other is to take full advantage of the nice feature that weak Galerkin discretization can be defined on polygonal meshes, which allows the curved boundary to be better approximated by multiple short edges without increasing the total number of mesh elements. Rigorous analysis shows that a combination of the above two techniques renders optimal convergence for all $α.$ Numerical results further confirm this conclusion.

You do not have full access to this article.

Already a Subscriber? Sign in as an individual or via your institution

Journal Article Details

Publisher Name:    Global Science Press

Language:    English

DOI:    https://doi.org/10.4208/cicp.OA-2022-0106

Communications in Computational Physics, Vol. 32 (2022), Iss. 4 : pp. 1094–1128

Published online:    2022-01

AMS Subject Headings:    Global Science Press

Copyright:    COPYRIGHT: © Global Science Press

Pages:    35

Keywords:    Weak Galerkin method polygonal mesh curved domain mixed formulation.

Author Details

Yi Liu

Wenbin Chen

Yanqiu Wang