NAJan 21, 2016
A Weak Galerkin Finite Element Scheme for solving the stationary Stokes EquationsRuishu Wang, Xiaoshen Wang, Qilong Zhai et al.
A weak Galerkin (WG) finite element method for solving the stationary Stokes equations in two- or three- dimensional spaces by using discontinuous piecewise polynomials is developed and analyzed. The variational form we considered is based on two gradient operators which is different from the usual gradient-divergence operators. The WG method is highly flexible by allowing the use of discontinuous functions on arbitrary polygons or polyhedra with certain shape regularity. Optimal-order error estimates are established for the corresponding WG finite element solutions in various norms. Numerical results are presented to illustrate the theoretical analysis of the new WG finite element scheme for Stokes problems.
NADec 31, 2016
Supercloseness analysis and polynomial preserving recovery for a class of weak Galerkin methodRuishu Wang, Ran Zhang, Xu Zhang et al.
In this paper, we analyze convergence and supercloseness properties of a class of weak Galerkin (WG) finite element methods for solving second-order elliptic problems. It is shown that the WG solution is superclose to the Lagrange type interpolation using Lobatto points. This supercloseness behavior is obtained through some newly designed stabilization terms. A post-processing technique using the polynomial preserving recovery (PPR) is introduced for WG approximation. Superconvergence analysis is carried out for the PPR approximation. Numerical examples are provided to verify our theoretical results.
NASep 29, 2015
A weak Galerkin finite element scheme with boundary continuity for second-order elliptic problemsQilong Zhai, Xiu Ye, Ruishu Wang et al.
A new weak Galerkin (WG) finite element method for solving the second-order elliptic problems on polygonal meshes by using polynomials of boundary continuity is introduced and analyzed. The WG method is utilizing weak functions and their weak derivatives which can be approximated by polynomials in different combination of polynomial spaces. Different combination gives rise to different weak Galerkin finite element methods, which makes WG methods highly flexible and efficient in practical computation. This paper explores the possibility of certain combination of polynomial spaces that minimize the degree of freedom in the numerical scheme, yet without losing the accuracy of the numerical approximation. Error estimates of optimal order are established for the corresponding WG approximations in both a discrete $H^1$ norm and the standard $L^2$ norm. In addition, the paper also presents some numerical experiments to demonstrate the power of the WG method. The numerical results show a great promise of the robustness, reliability, flexibility and accuracy of the WG method.
NAAug 16, 2015
A Locking-Free Weak Galerkin Finite Element Method for Elasticity Problems in the Primal FormulationChunmei Wang, Junping Wang, Ruishu Wang et al.
This paper presents an arbitrary order locking-free numerical scheme for linear elasticity on general polygonal/polyhedral partitions by using weak Galerkin (WG) finite element methods. Like other WG methods, the key idea for the linear elasticity is to introduce discrete weak strain and stress tensors which are defined and computed by solving inexpensive local problems on each element. Such local problems are derived from weak formulations of the corresponding differential operators through integration by parts. Locking-free error estimates of optimal order are derived in a discrete $H^1$-norm and the usual $L^2$-norm for the approximate displacement when the exact solution is smooth. Numerical results are presented to demonstrate the efficiency, accuracy, and the locking-free property of the weak Galerkin finite element method.