Search

In this paper, the authors constructed an auxiliary space multigrid preconditioner for the weak Galerkin finite element method for second-order diffusion equations, discretized on simplicial 2D or 3D meshes. The idea of the auxiliary space multigrid preconditioner is to use an auxiliary space as a "coarse" space in the multigrid algorithm, where the discrete problem in the auxiliary space can be easily solved by an existing solver. In this construction, the authors conveniently use the $H^1$ conforming piecewise linear finite element space as an auxiliary space. The main technical difficulty is to build the connection between the weak Galerkin discrete space and the $H^1$ conforming piecewise linear finite element space. The authors successfully constructed such an auxiliary space multigrid preconditioner for the weak Galerkin method, as well as a reduced system of the weak Galerkin method involving only the degrees of freedom on edges/faces. The preconditioned systems are proved to have condition numbers independent of the mesh size. Numerical experiments are conducted to support the theoretical results., Comment: 19 pages, 2 figures, 5 tables

A family of weak Galerkin finite element discretization is developed for solving the coupled Darcy-Stokes equation. The equation in consideration admits the Beaver-Joseph-Saffman condition on the interface. By using the weak Galerkin approach, in the discrete space we are able to impose the normal continuity of velocity explicitly. Or in other words, strong coupling is achieved in the discrete space. Different choices of weak Galerkin finite element spaces are discussed, and error estimates are given.

A weak Galerkin discretization of the boundary value problem of a general anisotropic diffusion problem is studied for preservation of the maximum principle. It is shown that the direct application of the $M$-matrix theory to the stiffness matrix of the weak Galerkin discretization leads to a strong mesh condition requiring all of the mesh dihedral angles to be strictly acute (a constant-order away from 90 degrees). To avoid this difficulty, a reduced system is considered and shown to satisfy the discrete maximum principle under weaker mesh conditions. The discrete maximum principle is then established for the full weak Galerkin approximation using the relations between the degrees of freedom located on elements and edges. Sufficient mesh conditions for both piecewise constant and general anisotropic diffusion matrices are obtained. These conditions provide a guideline for practical mesh generation for preservation of the maximum principle. Numerical examples are presented.

This paper focuses on interior penalty discontinuous Galerkin methods for second order elliptic equations on very general polygonal or polyhedral meshes. The mesh can be composed of any polygons or polyhedra which satisfies certain shape regularity conditions characterized in a recent paper by two of the authors in [17]. Such general meshes have important application in computational sciences. The usual $H^1$ conforming finite element methods on such meshes are either very complicated or impossible to implement in practical computation. However, the interior penalty discontinuous Galerkin method provides a simple and effective alternative approach which is efficient and robust. This article provides a mathematical foundation for the use of interior penalty discontinuous Galerkin methods in general meshes., Comment: 12 pages, research results

This article introduces and analyzes a weak Galerkin mixed finite element method for solving the biharmonic equation. The weak Galerkin method, first introduced by two of the authors (J. Wang and X. Ye) in an earlier publication for second order elliptic problems, is based on the concept of discrete weak gradients. The method allows the use of completely discrete finite element functions on partitions of arbitrary polygon or polyhedron. In this article, the idea of weak Galerkin method is applied to discretize the Ciarlet-Raviart mixed formulation for the biharmonic equation. In particular, an a priori error estimation is given for the corresponding finite element approximations. The error analysis essentially follows the framework of Babuska, Osborn, and Pitkaranta and uses specially designed mesh-dependent norms. The proof is technically tedious due to the discontinuous nature of the weak Galerkin finite element functions. Some computational results are presented to demonstrate the efficiency of the method., Comment: 23 pages, 4 tables

The weak Galerkin finite element method is a novel numerical method that was first proposed and analyzed by Wang and Ye for general second order elliptic problems on triangular meshes. The goal of this paper is to conduct a computational investigation for the weak Galerkin method for various model problems with more general finite element partitions. The numerical results confirm the theory established by Wang and Ye. The results also indicate that the weak Galerkin method is efficient, robust, and reliable in scientific computing., Comment: 19 pages