4.7 Article

A penalized weak Galerkin spectral element method for second order elliptic equations

出版社

ELSEVIER
DOI: 10.1016/j.cam.2020.113228

关键词

Weak Galerkin spectral element; Penalty; Weak derivatives; Piola transform; Error estimates

资金

  1. National Natural Science Foundation of China [11871455, 11971016]
  2. National Key R&D Program of China [2018YFB0204404]

向作者/读者索取更多资源

The proposed weak Galerkin spectral element method introduces weak functions and weak derivatives to solve second order partial differential equations, which is validated through numerical experiments for both typical square and L-shaped domains.
A weak Galerkin spectral element method is proposed to solve second order partial differential equations. Following the idea of the weak Galerkin finite element method, this method introduces for the unknown solution a weak function consisting of one component on the elements together with one component on the element interfaces, and replaces derivatives in the standard variational form with weak derivatives defined on the space of weak functions on each element. As in a classic spectral element method, approximation spaces for weak functions on each triangular or parallelogram element are defined from orthogonal polynomials on the reference domain through a one-to-one mapping, and approximation spaces for weak derivatives are then established via the Piola transform after an insightful investigation of the weak gradient and the discrete weak gradient. To eliminate the effect of the possible nullity of the discrete weak gradient and guarantee the wellposedness of the resulting algebraic system, a penalty term defined on the edges is supplemented into the Galerkin approximation scheme. Error estimates for both the source problem and the eigen-value problem on meshes consisting of affine families of triangles and quadrilaterals are obtained in the sequel, which are optimal in the mesh size and suboptimal by one-half order with respect to the polynomial degree. Numerical experiments for the eigen-value problems are performed on both the typical square domain and L-shaped domain with triangular meshes and quadrilateral meshes, which illustrate the effectiveness and high accuracy of our penalized weak Galerkin spectral element method. (C) 2020 Elsevier B.V. All rights reserved.

作者

我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。

评论

主要评分

4.7
评分不足

次要评分

新颖性
-
重要性
-
科学严谨性
-
评价这篇论文

推荐

暂无数据
暂无数据