期刊
JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS
卷 386, 期 -, 页码 -出版社
ELSEVIER
DOI: 10.1016/j.cam.2020.113228
关键词
Weak Galerkin spectral element; Penalty; Weak derivatives; Piola transform; Error estimates
资金
- National Natural Science Foundation of China [11871455, 11971016]
- 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.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据