Weak Galerkin finite element methods for elasticity and coupled flow problems

Includes bibliographical references.2020 Summer.We present novel stabilizer-free weak Galerkin finite element methods for linear elasticity and coupled Stokes-Darcy flow with a comprehensive treatment of theoretical results and the numerical methods for each. Weak Galerkin finite element methods take a discontinuous approximation space and bind degrees of freedom together through the discrete weak gradient, which involves solving a small symmetric positive-definite linear system on every element of the mesh. We introduce notation and analysis using a general framework that highlights properties that unify many existing weak Galerkin methods. This framework makes analysis for the methods much more straightforward. The method for linear elasticity on quadrilateral and hexahedral meshes uses piecewise constant vectors to approximate the displacement on each cell, and it uses the Raviart-Thomas space for the discrete weak gradient. We use the Schur complement to simplify the solution of the global linear system and increase computational efficiency further. We prove first-order convergence in the L2 norm, verify our analysis with numerical experiments, and compare to another weak Galerkin approach for this problem. The method for coupled Stokes-Darcy flow uses an extensible multinumerics approach on quadrilateral meshes. The Darcy flow discretization uses a weak Galerkin finite element method with piecewise constants approximating pressure and the Arbogast-Correa space for the weak gradient. The Stokes domain discretization uses the classical Bernardi-Raugel pair. We prove first-order convergence in the energy norm and verify our analysis with numerical experiments. All algorithms implemented in this dissertation are publicly available as part of James Liu's DarcyLite and Darcy+ packages and as part of the deal.II library

Addressing Integration Error for Polygonal Finite Elements Through Polynomial Projections: A Patch Test Connection

Polygonal finite elements generally do not pass the patch test as a result of quadrature error in the evaluation of weak form integrals. In this work, we examine the consequences of lack of polynomial consistency and show that it can lead to a deterioration of convergence of the finite element solutions. We propose a general remedy, inspired by techniques in the recent literature of mimetic finite differences, for restoring consistency and thereby ensuring the satisfaction of the patch test and recovering optimal rates of convergence. The proposed approach, based on polynomial projections of the basis functions, allows for the use of moderate number of integration points and brings the computational cost of polygonal finite elements closer to that of the commonly used linear triangles and bilinear quadrilaterals. Numerical studies of a two-dimensional scalar diffusion problem accompany the theoretical considerations