L. Beirão da Veiga,D. Mora,A. Silgado
Abstract:In this work we design a novel $C^1$-conforming virtual element method of arbitrary order $k \geq 2$, to solve the biharmonic problem on a domain with curved boundary and internal curved interfaces in two dimensions. By introducing a suitable stabilizing form, we develop a rigorous interpolation, stability, and convergence analysis obtaining optimal error estimates in the energy norm. Finally, we validate the theoretical findings through numerical experiments.
What problem does this paper attempt to address?
### What problem does this paper attempt to solve?
This paper aims to design a new C1 - conforming Virtual Element Method (VEM) for solving two - dimensional biharmonic problems with curved boundaries. Specifically, this method can handle fourth - order partial differential equations (PDEs) on regions with curved boundaries and internal curved interfaces. By introducing an appropriate stabilization form, the authors develop strict interpolation, stability and convergence analyses, and obtain optimal error estimates in the energy norm. Finally, the theoretical results are verified by numerical experiments.
#### Key Point Summary:
1. **Research Background**:
- Fourth - order partial differential equations (such as the biharmonic equation) are widely present in physical phenomena, for example, plate - bending problems, fluid - flow problems, Cahn - Hilliard phase - field models, etc.
- Standard numerical methods face the problem of geometric approximation error when dealing with regions with curved boundaries, which will affect the analysis and optimal convergence rate, especially in high - order schemes.
2. **Proposed Method**:
- A new type of C1 - conforming Virtual Element Method (VEM) is designed, which is suitable for biharmonic problems of any order \( k \geq 2 \).
- This method can be used on polygonal meshes with curved boundaries, avoiding the influence of geometric errors.
- A new stabilization form is introduced to ensure the H2 - consistency of the discrete space and handle high - order terms related to the curved boundary.
3. **Technical Contributions**:
- Strict interpolation, stability and convergence analyses are provided, and the optimal error estimates are proved.
- Methods for calculating discrete bilinear forms and load terms are developed.
- The correctness and effectiveness of the theoretical results are verified by numerical experiments.
4. **Application Scenarios**:
- This method is especially suitable for fields requiring high precision and complex geometric shapes, such as aerodynamics, biomedicine, geophysics, etc.
### Formula Presentation
- **Variational Form of the Biharmonic Problem**:
Given \( f \in L^2(\Omega) \), find \( u: \Omega \to \mathbb{R} \) such that:
\[
\begin{cases}
\Delta^2 u = f & \text{in } \Omega, \\
u = \partial_n u = 0 & \text{in } \Gamma.
\end{cases}
\]
Its variational form is: find \( u \in V := H^2_0(\Omega) \) such that
\[
a(u, v) = (f, v)_{0, \Omega} \quad \forall v \in V,
\]
where the bilinear form \( a(u, v) \) is defined as:
\[
a(u, v) := (D^2 u, D^2 v)_{0, \Omega}.
\]
- **Discrete Bilinear Form**:
\[
a_h^E(v_h, w_h) := a_E(\Pi_{D,k}^E v_h, \Pi_{D,k}^E w_h) + S_E((I - \Pi_{D,k}^E)v_h, (I - \Pi_{D,k}^E)w_h),
\]
where \( S_E \) is a symmetric stabilization bilinear form, defined as:
\[
S_E(\eta_h, \sigma_h) := h_E^{-2} \sum_{j = 1}^{(k - 3)(k - 2)/2} D_j^\circ(\eta_h) D_j^\circ(\sigma_h) + 3N_E \sum_{i = 1}^{N_E} h_{v_i}^{-2} D_{\partial, v_i}(\eta_h) D_{\partial, v_i}(\sigma_h) + h_E N_e \sum_{i = 1}^{N_e} \left( (\partial_{t_E}(\partial_{n_E} \right)