In this paper, we introduce a simple local flux recovery for $ \mathcal{Q}_k $ finite element of a scalar coefficient diffusion equation on quadtree meshes, with no restriction on the irregularities of hanging nodes. The construction requires no specific ad hoc tweaking for hanging nodes on $ l $-irregular ($ l\geq 2 $) meshes thanks to the adoption of virtual element families. The rectangular elements with hanging nodes are treated as polygons as in the flux recovery context. An efficient a posteriori error estimator is then constructed based on the recovered flux, and its reliability is proved under common assumptions, both of which are further verified in numerics.
Citation: Shuhao Cao. A simple virtual element-based flux recovery on quadtree[J]. Electronic Research Archive, 2021, 29(6): 3629-3647. doi: 10.3934/era.2021054
In this paper, we introduce a simple local flux recovery for $ \mathcal{Q}_k $ finite element of a scalar coefficient diffusion equation on quadtree meshes, with no restriction on the irregularities of hanging nodes. The construction requires no specific ad hoc tweaking for hanging nodes on $ l $-irregular ($ l\geq 2 $) meshes thanks to the adoption of virtual element families. The rectangular elements with hanging nodes are treated as polygons as in the flux recovery context. An efficient a posteriori error estimator is then constructed based on the recovered flux, and its reliability is proved under common assumptions, both of which are further verified in numerics.
| [1] | MFEM: A modular finite element library. Computers & Mathematics with Applications (2021) 81: 42-74. |
| [2] |
W. Bangerth, R. Hartmann and G. Kanschat, deal.II - a general purpose object oriented finite element library, ACM Trans. Math. Software, 33 (2007), Art. 24, 27 pp. doi: 10.1145/1268776.1268779
|
| [3] |
Asymptotically exact a posteriori error estimators, part ii: General unstructured grids. SIAM J. Numer. Anal. (2003) 41: 2313-2332.
|
| [4] |
Basic principles of virtual element methods. Mathematical Models and Methods in Applied Sciences (2013) 23: 199-214.
|
| [5] |
Serendipity face and edge VEM spaces. Atti Accad. Naz. Lincei Rend. Lincei Mat. Appl. (2017) 28: 143-180.
|
| [6] |
Adaptive finite element methods for elliptic equations with non-smooth coefficient. Numer. Math. (2000) 85: 579-608.
|
| [7] |
Orthogonal polynomials in badly shaped polygonal elements for the virtual element method. Finite Elem. Anal. Des. (2017) 129: 14-31.
|
| [8] |
Basic principles of mixed virtual element methods. ESAIM Math. Model. Numer. Anal. (2014) 48: 1227-1240.
|
| [9] |
On the Poisson equation with intersecting interfaces. Applicable Anal. (1974) 4: 101-129.
|
| [10] |
A recovery-based a posteriori error estimator for H(curl) interface problems. Comput. Methods in Appl. Mech. Eng. (2015) 296: 169-195.
|
| [11] |
Recovery-based error estimators for interface problems: conforming linear elements. SIAM J. Numer. Anal. (2009) 47: 2132-2156.
|
| [12] |
A posteriori error estimates for the virtual element method. Numer. Math. (2017) 137: 857-893.
|
| [13] |
Anisotropic error estimates of the linear nonconforming virtual element methods. SIAM J. Numer. Anal. (2019) 57: 1058-1081.
|
| [14] | Hanging nodes in the unifying theory of a posteriori finite element error control. J. Comput. Math. (2009) 27: 215-236. |
| [15] |
J. Červený, V. Dobrev and T. Kolev, Nonconforming mesh refinement for high-order finite elements, SIAM J. Sci. Comput., 41 (2019), C367-C392. doi: 10.1137/18M1193992
|
| [16] | L. Chen, iFEM: An Innovative Finite Element Methods Package in MATLAB, Technical report, 2008, URLhttps://github.com/lyc102/ifem. |
| [17] |
On the efficiency of adaptive finite element methods for elliptic problems with discontinuous coefficients. SIAM J. Sci. Comput. (2002) 24: 443-462.
|
| [18] |
A simple and effective gradient recovery scheme and a posteriori error estimator for the virtual element method (VEM). Comput. Methods Appl. Mech. Engrg. (2019) 347: 21-58.
|
| [19] | F. Dassi, J. Gedicke and L. Mascotto, Adaptive virtual element methods with equilibrated fluxes, arXiv preprint, arXiv: 2004.11220. |
| [20] |
Toward a universal hp adaptive finite element strategy, part 1. constrained approximation and data structure. Comput. Methods Appl. Mech. Engrg. (1989) 77: 79-112.
|
| [21] |
An easy treatment of hanging nodes in $hp$-finite elements. Finite Elem. Anal. Des. (2016) 121: 101-117.
|
| [22] |
Flux reconstruction and a posteriori error estimation for discontinuous Galerkin methods on general nonmatching grids. C. R. Math. Acad. Sci. Paris (2009) 347: 441-444.
|
| [23] |
V. Girault and P.-A. Raviart, Finite Element Methods for Navier-Stokes Equations: Theory and Algorithms, Springer-Verlag, 1986. doi: 10.1007/978-3-642-61623-5
|
| [24] |
Superconvergent gradient recovery for virtual element methods. Math. Models Methods Appl. Sci. (2019) 29: 2007-2031.
|
| [25] |
Ill-conditioning in the virtual element method: Stabilizations and bases. Numer. Methods Partial Differential Equations (2018) 34: 1258-1281.
|
| [26] |
Arbitrary-level hanging nodes and automatic adaptivity in the $hp$-FEM. Math. Comput. Simulation (2008) 77: 117-132.
|
| [27] |
Error estimates for some quasi-interpolation operators. M2AN Math. Model. Numer. Anal. (1999) 33: 695-713.
|
| [28] |
The superconvergent patch recovery and a posteriori error estimates. part 1: The recovery technique. Internat. J. Numer. Methods Engrg. (1992) 33: 1331-1364.
|
Figures(6)
Shuhao Cao. A simple virtual element-based flux recovery on quadtree[J]. Electronic Research Archive, 2021, 29(6): 3629-3647. doi: 10.3934/era.2021054
DownLoad: