1.2NAJan 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.
3.3NAJun 13
A Staggered Discontinuous Galerkin Method for linear elasticity problem on Polytopal MeshesLong Chen, Xuehai Huang, Ruishu Wang et al.
This paper develops a novel staggered discontinuous Galerkin (SDG) method for linear elasticity based on the Hellinger-Reissner variational principle. We construct symmetric stress spaces with normal continuity across element boundaries on arbitrary polytopal meshes, while approximating the displacement field using piecewise polynomial functions defined on the same meshes. The method is locking-free and satisfies a local balance of linear momentum and angular momentum. We present a comprehensive theoretical analysis, including proofs of stability and error estimates. The formulation admits a hybridizable structure, which significantly simplifies the numerical implementation. Numerical experiments validate the theoretical results and demonstrate the effectiveness of the proposed approach.
1.2NASep 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.
1.2NAAug 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.