CutVEM: Conforming virtual element method on embedded domains with shape-agnostic element agglomeration

TL;DR

CutVEM introduces a novel element agglomeration algorithm for the virtual element method, improving robustness on embedded domains with minimal impact on accuracy and maintaining optimal convergence rates.

Contribution

The paper presents a shape-agnostic, eigenvalue-based element agglomeration algorithm that enhances VEM stability on cut cell meshes for embedded domain simulations.

Findings

Significantly improved condition numbers of stiffness matrices.

Retains optimal convergence rates in heat conduction problems.

Demonstrates robustness across various embedded interface configurations.

Abstract

The virtual element method (VEM) is a stabilized Galerkin method that is robust and accurate on general polygonal meshes. This feature makes it an appealing candidate for simulations involving meshes with embedded interfaces and evolving geometries. However, similar to the finite element method, in such scenarios the VEM can also yield poorly conditioned stiffness matrices due to meshes having cut cells. With the objective of developing an embedded domain method, we propose a novel element agglomeration algorithm for the VEM to address this issue. The agglomeration algorithm renders the VEM robust over planar polygonal meshes, particularly on finite element meshes cut by immersed geometries. The algorithm relies on the element stability ratio, which we define using the extreme eigenvalues of the element stiffness matrix. The resulting element agglomeration criterion is free from nebulous polygon quality metrics and is defined independently of polygon shapes. The algorithm proceeds iteratively and element-wise to maximize the minimum element stability ratio, even at the expense of degrading elements with better ratios. The resulting method, which we label as CutVEM, retains node locations of cut elements unchanged, and yields discretizations that conform to embedded interfaces. This, in turn, facilitates straightforward imposition of boundary conditions and interfacial constraints. Through detailed numerical experiments that sample varied element-interface intersections, we demonstrate that CutVEM enjoys dramatically improved condition numbers of global stiffness matrices over the VEM. Furthermore, simulations of prototypical heat conduction problems with Dirichlet and Neumann boundary conditions on domains with immersed geometries show that element agglomeration does not noticeably degrade solution accuracy and that CutVEM retains the VEM's optimal convergence rate.