Abstract:This manuscript develops edge-averaged virtual element (EAVE) methodologies
to address convection-diffusion problems effectively in the
convection-dominated regime. It introduces a variant of EAVE that ensures
monotonicity (producing an $M$-matrix) on Voronoi polygonal meshes, provided
their duals are Delaunay triangulations with acute angles. Furthermore, the
study outlines a comprehensive framework for EAVE methodologies, introducing
another variant that integrates with the stiffness matrix derived from the
lowest-order virtual element method for the Poisson equation. Numerical
experiments confirm the theoretical advantages of the monotonicity property and
demonstrate an optimal convergence rate across various mesh configurations.
What problem does this paper attempt to address?
This paper attempts to solve the problem that the standard finite element method (FEM) cannot provide accurate numerical solutions in convection - dominated convection - diffusion problems. Specifically, when the convection term is much larger than the diffusion term (i.e., \(\alpha(x) \ll |\beta(x)|\)), the traditional FEM will produce spurious oscillations, thereby reducing the quality of the solution.
To solve this problem, the paper introduces the edge - averaged virtual element method (EAVE) and specifically develops an EAVE variant that ensures monotonicity and can generate M - matrices on Voronoi polygonal meshes, provided that its dual mesh is a Delaunay triangulation with acute angles. In addition, the paper also proposes an EAVE framework combined with the lowest - order virtual element method (the stiffness matrix for the Poisson equation). The theoretical advantages of these methods, including the optimal convergence rate and stability under different mesh configurations, are verified through numerical experiments.
In short, the main objectives of this paper are:
1. **Develop effective numerical methods**: For convection - dominated convection - diffusion problems, develop numerical methods that can avoid spurious oscillations and maintain the stability and accuracy of the solution.
2. **Ensure monotonicity**: By designing specific discrete schemes, ensure that the numerical solution satisfies the discrete maximum principle (DMP), thereby eliminating spurious oscillations.
3. **Applicable to polygonal meshes**: Extend traditional methods to make them applicable to a wider range of polygonal meshes, not just simple triangular or quadrilateral meshes.
In terms of formulas, the convection - diffusion equation has the following form:
\[
-\nabla \cdot (\alpha(x) \nabla u+\beta(x) u)=f \quad \text{in } \Omega,
\]
and the boundary condition is:
\[
u = 0 \quad \text{on } \partial \Omega.
\]
Here, \(\alpha(x)\) is the diffusion coefficient, \(\beta(x)\) is the convection velocity, and \(f\) is the source term. In the convection - dominated case, \(\alpha(x) \ll |\beta(x)|\), which makes the problem more complex and requires special numerical processing methods to ensure the stability and accuracy of the solution.