Abstract:We propose a new stable variational formulation for the quad-div problem in three dimensions and prove its well-posedness. Using this weak form, we develop and analyze the $\boldsymbol{H}(\operatorname{grad-div})$-conforming virtual element method of arbitrary approximation orders on polyhedral meshes. Three families of $\boldsymbol{H}(\operatorname{grad-div})$-conforming virtual elements are constructed based on the structure of a de Rham sub-complex with enhanced smoothness, resulting in an exact discrete virtual element complex. In the lowest-order case, the simplest element has only one degree of freedom at each vertex and face, respectively. We rigorously prove the interpolation error estimates, the stability of discrete bilinear forms, the well-posedness of discrete formulation and the optimal error estimates. Some numerical examples are shown to verify the theoretical results.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to propose a new stable variational formulation for the quad - div problem in three - dimensional space, and develop and analyze the H(grad - div)-conforming virtual element method (VEM) of arbitrary approximation order based on this. Specifically, the authors aim to solve the following quad - div problem:
\[
(\nabla \text{div})^2 u = f \quad \text{in } \Omega,
\]
\[
\text{curl} \, u = 0 \quad \text{in } \Omega,
\]
\[
u \cdot n = 0 \quad \text{on } \partial \Omega,
\]
\[
\text{div} \, u = 0 \quad \text{on } \partial \Omega,
\]
where \( u \) and \( f \) are vector fields defined on the bounded domain \( \Omega\subset \mathbb{R}^3 \).
### Main research objectives
1. **Propose a new variational formulation**: Propose a new stable variational formulation for the above quad - div problem and prove its well - posedness.
2. **Develop the virtual element method**: Based on the proposed variational formulation, develop the H(grad - div)-conforming virtual element method applicable to polyhedral meshes, including constructing three different types of H(grad - div)-conforming virtual elements.
3. **Theoretical analysis and numerical verification**: Rigorously prove the interpolation error estimate, the stability of the discrete bilinear form, the well - posedness of the discrete scheme and the optimal error estimate, and verify the theoretical results through numerical examples.
### Background and significance
- **Application background**: The quad - div operator originates from linear elasticity mechanics and represents shear strain energy. It is closely related to other fourth - order elliptic operators such as the biharmonic operator \( \Delta^2 \) and the quad - curl operator, and has important application value in the fields of materials science and engineering.
- **Research status**: Although the research on the biharmonic operator and the quad - curl operator is relatively extensive, the research on the quad - div problem is relatively scarce. Therefore, the work in this paper fills this gap in this field and provides new numerical methods and theoretical support.
### Method innovation points
- **Utilize the discrete de Rham complex**: By constructing a de Rham sub - complex with enhanced smoothness, ensure that the designed virtual element space maintains the compatibility of differential operators.
- **Introduce auxiliary variables**: In order to solve the ill - posed problem of the original variational problem, introduce auxiliary variables \( \phi \) and \( p \), and propose a new equivalent variational formulation, thereby ensuring the well - posedness of the problem.
- **Efficient computation**: Through reasonable selection of degrees of freedom and definition of projection operators, the computation in the case of high - order approximation is still efficient and accurate.
In summary, this paper aims to provide a new numerical solution method for the three - dimensional quad - div problem and verify its effectiveness and superiority through strict theoretical analysis and numerical experiments.