Abstract:In the present work we introduce a novel refinement algorithm for two-dimensional elliptic partial differential equations discretized with Virtual Element Method (VEM). The algorithm improves the numerical solution accuracy and the mesh quality through a controlled refinement strategy applied to the generic polygonal elements of the domain tessellation. The numerical results show that the outlined strategy proves to be versatile and applicable to any two-dimensional problem where polygonal meshes offer advantages. In particular, we focus on the simulation of flow in fractured media, specifically using the Discrete Fracture Network (DFN) model. A residual a-posteriori error estimator tailored for the DFN case is employed. We chose this particular application to emphasize the effectiveness of the algorithm in handling complex geometries. All the numerical tests demonstrate optimal convergence rates for all the tested VEM orders.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is that when traditional triangular mesh generation methods are computationally expensive and impractical in dealing with complex geometric structures, how to improve the accuracy of numerical solutions and mesh quality through improved polygon mesh generation and refinement algorithms. Specifically, for two - dimensional elliptic partial differential equations (PDEs), the article proposes a new refinement algorithm based on the virtual element method (VEM). This algorithm aims to optimize arbitrary polygon elements through control strategies to improve the accuracy of numerical solutions and mesh quality. ### Main problems and goals 1. **Improve the accuracy of numerical solutions**: By introducing a new refinement algorithm, ensure that more accurate numerical solutions can be obtained in complex geometric structures. 2. **Improve mesh quality**: By optimizing the quality of polygon meshes, make the simulation results more reliable. 3. **Apply to complex geometric structures**: Especially apply to fluid flow simulation in discrete fracture network (DFN) models, emphasizing the effectiveness of the algorithm in dealing with complex geometric structures. 4. **Ensure convergence and optimality**: Verify the convergence and optimality of the algorithm through a series of numerical tests. ### Research background In the fields of science and engineering, many applications require solving partial differential equations (PDEs) on highly complex geometric structures. For these complex geometric structures, generating high - quality triangular meshes is usually computationally expensive or impractical. Therefore, in recent years, the method of general polygon meshes has been adopted. In this context, there is an urgent need for a numerical method that can handle general polygon tessellations, and the virtual element method (VEM) is one of them. In addition, a posteriori error estimation is crucial for ensuring the reliability of simulation tools. By developing adaptive algorithms, the ability to control approximation errors can be used to reduce computational costs. The refinement algorithm proposed in this paper is especially suitable for two - dimensional second - order elliptic PDEs and extends previous research results. ### Key contributions 1. **Introduce a new refinement strategy**: By combining the maximum momentum direction (MAX - MOMENTUM) and smooth direction (SMOOTH - DIRECTION) algorithms, achieve the optimal segmentation of polygon elements. 2. **Improve mesh quality**: By checking conditions (CHECK - QUALITY), ensure that the newly generated sub - polygons have higher quality and fewer aligned edges. 3. **Extension operation**: Through the extension operation, ensure that all polygon units containing marked edges can meet the quality check conditions, thereby further improving the convergence speed and the accuracy of numerical solutions. ### Application examples The article pays special attention to fluid flow simulation in discrete fracture network (DFN) models. The DFN model naturally generates an initial discretization mesh, which contains a large number of aligned edges and poorly - shaped polygon units. By applying the new refinement algorithm, the mesh quality and the accuracy of numerical solutions are significantly improved. In conclusion, this research proposes a new polygon mesh generation and refinement algorithm, which can effectively improve the accuracy of numerical solutions and mesh quality in complex geometric structures and has broad application prospects.

Effective polygonal mesh generation and refinement for VEM

A new quality preserving polygonal mesh refinement algorithm for Virtual Element Methods

Novel Methodology for Viscous-Layer Meshing by the Boundary Element Method

Virtual finite element and hyperbolic problems: the PAMPA algorithm

The Mixed Virtual Element Method on curved edges in two dimensions

Mesh Optimization for the Virtual Element Method: How Small Can an Agglomerated Mesh Become?

Vem++, a c++ library to handle and play with the Virtual Element Method

Machine learning based refinement strategies for polyhedral grids with applications to virtual element and polyhedral discontinuous Galerkin methods

A Priori and a Posteriori Error Estimates for a Virtual Element Method for the Non-Self-adjoint Steklov Eigenvalue Problem

A Virtual Element Method for 2D linear elastic fracture analysis

A posteriori error estimates for the virtual element method

Virtual element method for nonlinear Sobolev equation on polygonal meshes

The role of mesh quality and mesh quality indicators in the Virtual Element Method

Mesh Quality Agglomeration algorithm for the Virtual Element Method applied to Discrete Fracture Networks

3D Adaptive VEM with stabilization-free a posteriori error bounds

Virtual element method for modeling the deformation of multiphase composites

Nonconforming Virtual Element basis functions for space-time Discontinuous Galerkin schemes on unstructured Voronoi meshes

Virtual elements on agglomerated finite elements to increase the critical time step in elastodynamic simulations

General Order Virtual Element Methods for Neumann Boundary Optimal Control Problems in Saddle Point Formulation

Extended virtual element method for two-dimensional linear elastic fracture

Unconditionally optimal error estimates of linearized Crank-Nicolson Virtual element methods for quasilinear parabolic problems on general polygonal meshes