Exponential convergence of the hp Virtual Element Method with corner singularities

Lorenzo Mascotto,Lourenço Beirão da Veiga,Alexey Chernov,Alessandro Russo

DOI: https://doi.org/10.48550/arXiv.1611.10165

2016-11-30

Abstract:In the present work, we analyze the $hp$ version of Virtual Element methods for the 2D Poisson problem. We prove exponential convergence of the energy error employing sequences of polygonal meshes geometrically refined, thus extending the classical choices for the decomposition in the $hp$ Finite Element framework to very general decomposition of the domain. A new stabilization for the discrete bilinear form with explicit bounds in $h$ and $p$ is introduced. Numerical experiments validate the theoretical results. We also exhibit a numerical comparison between $hp$ Virtual Elements and $hp$ Finite Elements.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: when dealing with two - dimensional Poisson problems with corner singularities, by using geometrically refined polygonal meshes and appropriate local accuracy distributions, it is proved that the hp - version virtual element method (VEM) can achieve exponential convergence of the energy error. Specifically, the authors hope to extend the decomposition choices in the classical hp finite element framework to very general domain decompositions, and for this purpose, they introduce a new stabilization method for the discrete bilinear form, which has an explicit p - dependent bound. In addition, they also conduct numerical experiments to verify the theoretical results and compare the performance between hp - virtual elements and hp - finite elements.

### Main research contents
1. **Generalize hp - VEM to corner singularity problems**:
- The authors apply hp - VEM to two - dimensional Poisson problems with corner singularities.
- Use geometrically refined polygonal meshes to deal with non - smooth solutions.

2. **Introduce a new stabilization method**:
- Propose a new stabilization method applicable to the hp - version of VEM.
- Explicitly give the local accuracy bound of the new stabilization term with respect to the H1 semi - norm.

3. **Theoretical analysis and numerical verification**:
- Prove the exponential convergence of the energy error on geometrically refined polygonal meshes.
- Conduct numerical experiments to verify the validity of the theoretical results.

4. **Comparison with hp - finite elements**:
- Show the comparison of the numerical performance between hp - virtual elements and hp - finite elements.

### Key points for solving the problem
- **Geometrically refined polygonal meshes**: In order to deal with corner singularity problems, geometrically refined polygonal meshes are adopted, which enables the method to capture the singular behavior of the solution more accurately.
- **New stabilization method**: A new stabilization method is introduced to ensure the stability and convergence of the method.
- **Exponential convergence**: Through theoretical analysis and numerical experiments, the exponential convergence of the energy error is proved, which is one of the main contributions of this research.

In conclusion, this paper successfully solves the two - dimensional Poisson problem with corner singularities by improving and extending the virtual element method, and shows its superior performance on geometrically refined polygonal meshes.