The interior penalty virtual element method for fourth-order singular perturbation problems

Fang Feng,Yue Yu

2023-12-18

Abstract:This paper is dedicated to the numerical solution of a fourth-order singular perturbation problem using the interior penalty virtual element method (IPVEM) proposed in [42]. The study introduces modifications to the jumps and averages in the penalty term, as well as presents an automated mesh-dependent selection of the penalty parameter. Drawing inspiration from the modified Morley finite element methods, we leverage the conforming interpolation technique to handle the lower part of the bilinear form. Through our analysis, we establish optimal convergence in the energy norm and provide a rigorous proof of uniform convergence concerning the perturbation parameter in the lowest-order case.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is about the numerical solution of fourth - order singular perturbation problems. Specifically, the author focuses on how to effectively solve the fourth - order partial differential equation in the following form when the perturbation parameter \(\varepsilon\) approaches 0:

\[
\begin{cases}
\varepsilon^2 \Delta^2 u - \Delta u = f & \text{in } \Omega, \\
u = \frac{\partial u}{\partial n} = 0 & \text{on } \partial \Omega,
\end{cases}
\]

where \(\Omega\) is a bounded polygonal region, \(f\in L^2(\Omega)\), \(\Delta\) is the standard Laplace operator, \(\varepsilon\) is a small positive real number (\(0 < \varepsilon\leq1\)), and \(n\) is the unit outer normal vector of the boundary \(\partial\Omega\).

### Main Challenges
1. **When \(\varepsilon\) is relatively large**: At this time, the problem can be approximated as a bi - harmonic equation, and the standard fourth - order finite element method can be used to solve it.
2. **When \(\varepsilon\rightarrow0\)**: The equation degenerates into a Poisson equation, and the traditional \(H^2\)-conforming finite element methods (such as Morley elements) perform poorly or even diverge when dealing with this limit situation. Therefore, a new method needs to be found to ensure the stability and convergence of the numerical solution.

### Solutions
To meet the above challenges, the author introduced the Interior Penalty Virtual Element Method (IPVEM). The main features of this method include:
- **Modify the jump terms and average terms**: Adjust the jump and average in the penalty terms to better adapt to the characteristics of singular perturbation problems.
- **Adaptive mesh - dependent selection**: A method for automatically selecting the penalty parameter is proposed, so that the method can maintain stability under different mesh scales.
- **Elliptic projection**: Use the conforming interpolation technique to deal with the low - order part of the bilinear form, thereby simplifying the implementation and improving the computational efficiency.

Through these improvements, IPVEM shows optimal convergence in the energy norm and has uniform convergence with respect to the perturbation parameter in the lowest - order case.

### Summary
The core objective of this paper is to develop a robust numerical method to solve fourth - order singular perturbation problems, especially when the perturbation parameter approaches zero, to ensure the stability and high - precision of the numerical solution. By introducing the Interior Penalty Virtual Element Method and its related improvements, the author has successfully solved this problem and provided strict theoretical analysis and numerical verification.