Search
Close
Search
Advanced Search
Search Menu
AI Discovery Assistant

Abstract

Nondivergence-free discretizations for the incompressible Stokes problem may suffer from a lack of pressure-robustness characterized by large discretizations errors due to irrotational forces in the momentum balance. This paper argues that also divergence-free virtual element methods on polygonal meshes are not really pressure-robust as long as the right-hand side is not discretized in a careful manner. To be able to evaluate the right-hand side for the test functions, some explicit interpolation of the virtual test functions is needed that can be evaluated pointwise everywhere. The standard discretization via an L2-best approximation does not preserve the divergence, and so destroys the orthogonality between divergence-free test functions and possibly eminent gradient forces in the right-hand side. To repair this orthogonality and restore pressure-robustness, another divergence-preserving reconstruction is suggested based on Raviart–Thomas approximations on local subtriangulations of the polygons. All findings are proven theoretically and are demonstrated numerically in two dimensions. The construction is also interesting for hybrid high-order methods on polygonal or polyhedral meshes.

pressure-robustness, incompressible Stokes equations, virtual element methods, polyhedral meshes, divergence-preserving reconstruction
Issue Section:
Article
© The Author(s) 2020. Published by Oxford University Press on behalf of the Institute of Mathematics and its Applications. All rights reserved.
This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://dbpia.nl.go.kr/journals/pages/open_access/funder_policies/chorus/standard_publication_model)
You do not currently have access to this article.
Download all slides