In the present paper we initiate the challenging task of building a mathematically sound theory for adaptive virtual element methods (AVEMs). Among the realm of polygonal meshes, we restrict our analysis to triangular meshes with hanging nodes in two dimensions-the simplest meshes with a systematic refinement procedure that preserves shape regularity and optimal complexity. A major challenge in the a posteriori error analysis of AVEMs is the presence of the stabilization term, which is of the same order as the residual-type error estimator but prevents the equivalence of the latter with the energy error. Under the assumption that any chain of recursively created hanging nodes has uniformly bounded length, we show that the stabilization term can be made arbitrarily small relative to the error estimator provided the stabilization parameter of the scheme is sufficiently large. This quantitative estimate leads to stabilization-free upper and lower a posteriori bounds for the energy error. This novel and crucial property of VEMs hinges on the largest subspace of continuous piecewise linear functions and the delicate interplay between its coarser scales and the finer ones of the VEM space. An important consequence for piecewise constant data is a contraction property between consecutive loops of AVEMs, which we also prove. Our results apply to H1-conforming (lowest order) VEMs of any kind, including the classical and enhanced VEMs.
Beirao Da Veiga, L., Canuto, C., Nochetto, R., Vacca, G., Verani, M. (2023). Adaptive VEM: Stabilization-Free A Posteriori Error Analysis and Contraction Property. SIAM JOURNAL ON NUMERICAL ANALYSIS, 61(2), 457-494 [10.1137/21M1458740].
Adaptive VEM: Stabilization-Free A Posteriori Error Analysis and Contraction Property
Beirao Da Veiga L.;
2023
Abstract
In the present paper we initiate the challenging task of building a mathematically sound theory for adaptive virtual element methods (AVEMs). Among the realm of polygonal meshes, we restrict our analysis to triangular meshes with hanging nodes in two dimensions-the simplest meshes with a systematic refinement procedure that preserves shape regularity and optimal complexity. A major challenge in the a posteriori error analysis of AVEMs is the presence of the stabilization term, which is of the same order as the residual-type error estimator but prevents the equivalence of the latter with the energy error. Under the assumption that any chain of recursively created hanging nodes has uniformly bounded length, we show that the stabilization term can be made arbitrarily small relative to the error estimator provided the stabilization parameter of the scheme is sufficiently large. This quantitative estimate leads to stabilization-free upper and lower a posteriori bounds for the energy error. This novel and crucial property of VEMs hinges on the largest subspace of continuous piecewise linear functions and the delicate interplay between its coarser scales and the finer ones of the VEM space. An important consequence for piecewise constant data is a contraction property between consecutive loops of AVEMs, which we also prove. Our results apply to H1-conforming (lowest order) VEMs of any kind, including the classical and enhanced VEMs.File | Dimensione | Formato | |
---|---|---|---|
BeiraodaVeiga-2022-Arxiv-PrePrint.pdf
accesso aperto
Descrizione: Preprint Submitted Version
Tipologia di allegato:
Submitted Version (Pre-print)
Licenza:
Creative Commons
Dimensione
3.3 MB
Formato
Adobe PDF
|
3.3 MB | Adobe PDF | Visualizza/Apri |
BeiraodaVeiga-2023-SINUM-VoR.pdf
Solo gestori archivio
Descrizione: Published Version
Tipologia di allegato:
Publisher’s Version (Version of Record, VoR)
Licenza:
Tutti i diritti riservati
Dimensione
6.16 MB
Formato
Adobe PDF
|
6.16 MB | Adobe PDF | Visualizza/Apri Richiedi una copia |
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.