skip to main content

Title: A discrete maximum principle for the weak Galerkin finite element method on nonuniform rectangular partitions

This article establishes a discrete maximum principle (DMP) for the approximate solution of convection–diffusion–reaction problems obtained from the weak Galerkin (WG) finite element method on nonuniform rectangular partitions. The DMP analysis is based on a simplified formulation of the WG involving only the approximating functions defined on the boundary of each element. The simplified weak Galerkin (SWG) method has a reduced computational complexity over the usual WG, and indeed provides a discretization scheme different from the WG when the reaction terms are present. An application of the SWG on uniform rectangular partitions yields some 5‐ and 7‐point finite difference schemes for the second order elliptic equation. Numerical experiments are presented to verify the DMP and the accuracy of the scheme, particularly the finite difference scheme.

more » « less
Author(s) / Creator(s):
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Numerical Methods for Partial Differential Equations
Page Range / eLocation ID:
p. 552-578
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    This article introduces a particular weak Galerkin (WG) element on rectangular/cuboid partitions that uses th order polynomial for weak finite element functions and th order polynomials for weak derivatives. This WG element is highly accurate with convergence two orders higher than the optimal order in an energy norm and the norm. The superconvergence is verified analytically and numerically. Furthermore, the usual stabilizer in the standard weak Galerkin formulation is no longer needed for this element.

    more » « less
  2. Abstract

    The weak Galerkin (WG) finite element method is an effective and flexible general numerical technique for solving partial differential equations. A simple WG finite element method is introduced for second‐order elliptic problems. First we have proved that stabilizers are no longer needed for this WG element. Then we have proved the supercloseness of order two for the WG finite element solution. The numerical results confirm the theory.

    more » « less
  3. It is well known that convergence rate of finite element approximation is suboptimal in theL2norm for solving biharmonic equations whenP2orQ2element is used. The goal of this paper is to derive a weak Galerkin (WG)P2element with theL2optimal convergence rate by assuming the exact solution sufficiently smooth. In addition, our new WG finite element method can be applied to general mesh such as hybrid mesh, polygonal mesh or mesh with hanging node. The numerical experiments have been conducted on different meshes including hybrid meshes with mixed of pentagon and rectangle and mixed of hexagon and triangle.

    more » « less
  4. Abstract

    An adaptive modified weak Galerkin method (AmWG) for an elliptic problem is studied in this article, in addition to its convergence and optimality. The modified weak Galerkin bilinear form is simplified without the need of the skeletal variable, and the approximation space is chosen as the discontinuous polynomial space as in the discontinuous Galerkin method. Upon a reliable residual‐baseda posteriorierror estimator, an adaptive algorithm is proposed together with its convergence and quasi‐optimality proved for the lowest order case. The primary tool is to bridge the connection between the modified weak Galerkin method and the Crouzeix–Raviart nonconforming finite element. Unlike the traditional convergence analysis for methods with a discontinuous polynomial approximation space, the convergence of AmWG is penalty parameter free. Numerical results are presented to support the theoretical results.

    more » « less
  5. In this paper we propose and analyze finite element discontinuous Galerkin methods for the one- and two-dimensional stochastic Maxwell equations with multiplicative noise. The discrete energy law of the semi-discrete DG methods were studied. Optimal error estimate of the semi-discrete method is obtained for the one-dimensional case, and the two-dimensional case on both rectangular meshes and triangular meshes under certain mesh assumptions. Strong Taylor 2.0 scheme is used as the temporal discretization. Both one- and two-dimensional numerical results are presented to validate the theoretical analysis results. 
    more » « less