Search

A mathematical analysis is established for the weak Galerkin finite element methods for the Poisson equation with Dirichlet boundary value when the curved elements are involved on the interior edges of the finite element partition or/and on the boundary of the whole domain in two dimensions. The optimal orders of error estimates for the weak Galerkin approximations in both the $H^1$-norm and the $L^2$-norm are established. Numerical results are reported to demonstrate the performance of the weak Galerkin methods on general curved polygonal partitions., Comment: 25 pages, 7 figures, 3 tables

This article proposes and analyzes the generalized weak Galerkin ({\rm g}WG) finite element method for the second order elliptic problem. A generalized discrete weak gradient operator is introduced in the weak Galerkin framework so that the {\rm g}WG methods would not only allow arbitrary combinations of piecewise polynomials defined in the interior and on the boundary of each local finite element, but also work on general polytopal partitions. Error estimates are established for the corresponding numerical functions in the energy norm and the usual $L^2$ norm. A series of numerical experiments are presented to demonstrate the performance of the newly proposed {\rm g}WG method., Comment: 30 pages, 19 tables

This paper introduces an extension of the Morley element for approximating solutions to biharmonic equations. Traditionally limited to piecewise quadratic polynomials on triangular elements, the extension leverages weak Galerkin finite element methods to accommodate higher degrees of polynomials and the flexibility of general polytopal elements. By utilizing the Schur complement of the weak Galerkin method, the extension allows for fewest local degrees of freedom while maintaining sufficient accuracy and stability for the numerical solutions. The numerical scheme incorporates locally constructed weak tangential derivatives and weak second order partial derivatives, resulting in an accurate approximation of the biharmonic equation. Optimal order error estimates in both a discrete $H^2$ norm and the usual $L^2$ norm are established to assess the accuracy of the numerical approximation. Additionally, numerical results are presented to validate the developed theory and demonstrate the effectiveness of the proposed extension., Comment: 22 pages, 4 tables, 4 figures

A new $L^p$-primal-dual weak Galerkin method ($L^p$-PDWG) with $p>1$ is proposed for the first-order transport problems. The existence and uniqueness of the $L^p$-PDWG numerical solutions is established. In addition, the $L^p$-PDWG method offers a numerical solution which retains mass conservation locally on each element. An optimal order error estimate is established for the primal variable. A series of numerical results are presented to verify the efficiency and accuracy of the proposed $L^p$-PDWG scheme., 26 pages, 7 figures, 10 tables. arXiv admin note: text overlap with arXiv:1906.07336

This paper introduces an extension of the well-known Morley element for the biharmonic equation, extending its application from triangular elements to general polytopal elements using the weak Galerkin finite element methods. By leveraging the Schur complement of the weak Galerkin method, this extension not only preserves the same degrees of freedom as the Morley element on triangular elements but also expands its applicability to general polytopal elements. The numerical scheme is devised by locally constructing weak tangential derivatives and weak second-order partial derivatives. Error estimates for the numerical approximation are established in both the energy norm and the $L^2$ norm. A series of numerical experiments are conducted to validate the theoretical developments., 21 pages, 4 tables, 1 figure

This article introduces a weak Galerkin (WG) finite element method for quad-curl problems in three dimensions. It is proved that the proposed WG method is stable and accurate in an optimal order of error estimates for the exact solution in discrete norms. In addition, an $L^2$ error estimate in an optimal order except the lowest orders $k=1, 2$ is derived for the WG solution. Some numerical experiments are conducted to verify the efficiency and accuracy of our WG method and furthermore a superconvergence has been observed from the numerical results., 16 pages, 2 tables

This paper presents a new $L^p$-primal-dual weak Galerkin (PDWG) finite element method for the div-curl system with the normal boundary condition for $p>1$. Two crucial features for the proposed $L^p$-PDWG finite element scheme are as follows: (1) it offers an accurate and reliable numerical solution to the div-curl system under the low $W^{\alpha, p}$-regularity ($\alpha>0$) assumption for the exact solution; (2) it offers an effective approximation of the normal harmonic vector fields on domains with complex topology. An optimal order error estimate is established in the $L^q$-norm for the primal variable where $\frac{1}{p}+\frac{1}{q}=1$. A series of numerical experiments are presented to demonstrate the performance of the proposed $L^p$-PDWG algorithm., Comment: 22 pages, 2 figures, 8 tables. arXiv admin note: text overlap with arXiv:2101.03466

In this article, the authors present a new $L^p$- primal-dual weak Galerkin method ($L^p$-PDWG) for convection-diffusion equations with $p>1$. The existence and uniqueness of the numerical solution is discussed, and an optimal-order error estimate is derived in the $L^q$-norm for the primal variable, where $\frac 1p+\frac 1q=1$. Furthermore, error estimates are established for the numerical approximation of the dual variable in the standard $W^{m,p}$ norm, $0\le m\le 2$. Numerical results are presented to demonstrate the efficiency and accuracy of the proposed $L^p$-PDWG method., 22 pages, 10 tables, original research article

This article presents a new primal-dual weak Galerkin method for second order elliptic equations in non-divergence form. The new method is devised as a constrained $L^p$-optimization problem with constraints that mimic the second order elliptic equation by using the discrete weak Hessian locally on each element. An equivalent min-max characterization is derived to show the existence and uniqueness of the numerical solution. Optimal order error estimates are established for the numerical solution under the discrete $W^{2,p}$ norm, as well as the standard $W^{1,p}$ and $L^p$ norms. An equivalent characterization of the optimization problem in term of a system of fixed-point equations via the proximity operator is presented. An iterative algorithm is designed based on the fixed-point equations to solve the optimization problems. Implementation of the iterative algorithm is studied and convergence of the iterative algorithm is established. Numerical experiments for both smooth and non-smooth coefficients problems are presented to verify the theoretical findings., 36 pages

In this work, semi-discrete and fully-discrete error estimates are derived for the Biot's consolidation model described using a three-field finite element formulation. The fields include displacements, total stress and pressure. The model is implemented using a backward Euler discretization in time for the fully-discrete scheme and validated for benchmark examples. Computational experiments presented verifies the convergence orders for the lowest order finite elements with discontinuous and continuous finite element appropriation for the total stress., 18 pages, 2 figures

The goal of this article is to clarify some misunderstandings and inappropriate claims made in [6] regarding the relation between the weak Galerkin (WG) finite element method and the hybridizable discontinuous Galerkin (HDG). In this paper, the authors offered their understandings and interpretations on the weak Galerkin finite element method by describing the basics of the WG method and how WG can be applied to a model PDE problem in various variational forms. In the authors' view, WG-FEM and HDG methods are based on different philosophies and therefore represent different methodologies in numerical PDEs, though they share something in common in their roots. A theory and an example are given to show that the primal WG-FEM is not equivalent to the existing HDG [9].

A new primal-dual weak Galerkin (PD-WG) finite element method was developed and analyzed in this article for first-order linear convection equations in non-divergence form. The PD-WG method results in a symmetric discrete system involving not only the original equation for the primal variable, but also the dual/adjoint equation for the dual variable (also known as Lagrangian multiplier). Optimal order of error estimates in various discrete Sobolev norms are derived for the numerical solutions arising from the PD-WG scheme. Numerical results are produced and reported to illustrate the accuracy and effectiveness of the new PD-WG method., Comment: 35 pages, 10 figures, 38 tables. arXiv admin note: text overlap with arXiv:1906.07336

The authors propose and analyze a well-posed numerical scheme for a type of ill-posed elliptic Cauchy problem by using a constrained minimization approach combined with the weak Galerkin finite element method. The resulting Euler-Lagrange formulation yields a system of equations involving the original equation for the primal variable and its adjoint for the dual variable, and is thus an example of the primal-dual weak Galerkin finite element method. This new primal-dual weak Galerkin algorithm is consistent in the sense that the system is symmetric, well-posed, and is satisfied by the exact solution. A certain stability and error estimates were derived in discrete Sobolev norms, including one in a weak $L^2$ topology. Some numerical results are reported to illustrate and validate the theory developed in the paper., 27 pages, 4 figures

In this article a simplified weak Galerkin finite element method is developed for the Dirichlet boundary value problem of convection-diffusion-reaction equations. The simplified weak Galerkin method utilizes only the degrees of freedom on the boundary of each element and, hence, has significantly reduced computational complexity over the regular weak Galerkin finite element method. A stability and some optimal order error estimates in the $H^1$ and $L^2$ norms are established for the corresponding numerical solutions. Numerical results are presented to verify the theory error estimates and a superconvergence phenomena on rectangular partitions.

A discrete divergence-free weak Galerkin finite element method is developed for the Stokes equations based on a weak Galerkin (WG) method introduced in the reference [15]. Discrete divergence-free bases are constructed explicitly for the lowest order weak Galerkin elements in two and three dimensional spaces. These basis functions can be derived on general meshes of arbitrary shape of polygons and polyhedrons. With the divergence-free basis derived, the discrete divergence-free WG scheme can eliminate the pressure variable from the system and reduces a saddle point problem to a symmetric and positive definite system with many fewer unknowns. Numerical results are presented to demonstrate the robustness and accuracy of this discrete divergence-free WG method., 12 pages

In this paper, the authors devise a new discretization scheme for div-curl systems defined in connected domains with heterogeneous media by using the weak Galerkin finite element method. Two types of boundary value problems are considered in the algorithm development: (1) normal boundary condition, and (2) tangential boundary condition. A new variational formulation is developed for the normal boundary value problem by using the Helmholtz decomposition which avoids the computation of functions in the harmonic fields. Both boundary value problems are reduced to a general saddle-point problem involving the curl and divergence operators, for which the weak Galerkin finite element method is devised and analyzed. The novelty of the technique lies in the discretization of the divergence operator applied to vector fields with heterogeneous media. Error estimates of optimal order are established for the corresponding finite element approximations in various discrete Sobolev norms., Comment: 27 pages

This paper presents a hybridized formulation for the weak Galerkin finite element method for the biharmonic equation. The hybridized weak Galerkin scheme is based on the use of a Lagrange multiplier defined on the element boundaries. The Lagrange multiplier is verified to provide a numerical approximation for certain derivatives of the exact solution. An optimal order error estimate is established for the numerical approximations arising from the hybridized weak Galerkin finite element method. The paper also derives a computational algorithm (Schur complement) by eliminating all the unknown variables on each element, yielding a significantly reduced system of linear equations for unknowns on the boundary of each element., Comment: 17 pages

A new weak Galerkin (WG) method is introduced and analyzed for the second order elliptic equation formulated as a system of two first order linear equations. This method, called WG-MFEM, is designed by using discontinuous piecewise polynomials on finite element partitions with arbitrary shape of polygons/polyhedra. The WG-MFEM is capable of providing very accurate numerical approximations for both the primary and flux variables. Allowing the use of discontinuous approximating functions on arbitrary shape of polygons/polyhedra makes the method highly flexible in practical computation. Optimal order error estimates in both discrete $H^1$ and $L^2$ norms are established for the corresponding weak Galerkin mixed finite element solutions., 26 pages

A newly developed weak Galerkin method is proposed to solve parabolic equations. This method allows the usage of totally discontinuous functions in approximation space and preserves the energy conservation law. Both continuous and discontinuous time weak Galerkin finite element schemes are developed and analyzed. Optimal order error estimates in both H^1 and L^2 norms are established. Numerical tests are performed and reported., Comment: 18 pages, 6 tables with computational results

A C^0-weak Galerkin (WG) method is introduced and analyzed for solving the biharmonic equation in 2D and 3D. A weak Laplacian is defined for C^0 functions in the new weak formulation. This WG finite element formulation is symmetric, positive definite and parameter free. Optimal order error estimates are established in both a discrete H^2 norm and the L^2 norm, for the weak Galerkin finite element solution. Numerical results are presented to confirm the theory. As a technical tool, a refined Scott-Zhang interpolation operator is constructed to assist the corresponding error estimate. This refined interpolation preserves the volume mass of order (k+1-d) and the surface mass of order (k+2-d) for the P_{k+2} finite element functions in d-dimensional space., Comment: 21 pages

Weak Galerkin (WG) refers to general finite element methods for partial differential equations in which differential operators are approximated by weak forms through the usual integration by parts. In particular, WG methods allow the use of discontinuous finite element functions in the algorithm design. One of such examples was recently introduced by Wang and Ye for solving second order elliptic problems. The goal of this paper is to apply the WG method of Wang and Ye to the Helmholtz equation with high wave numbers. Several test scenarios are designed for a numerical investigation on the accuracy, convergence, and robustness of the WG method in both inhomogeneous and homogeneous media over convex and non-convex domains. Our numerical experiments indicate that weak Galerkin is a finite element technique that is easy to implement, and provides very accurate and robust numerical solutions for the Helmholtz problem with high wave numbers., 20 pages, many figures from computation