Abstract:We consider a mixed finite element method for a biharmonic equation with clamped boundary conditions based on biorthogonal systems with weakly imposed Dirichlet boundary condition. We show that the weak imposition of the boundary condition arising from a natural minimisation formulation allows to get an optimal a priori error estimate for the finite element scheme improving the existing error estimate for such a formulation without weakly imposed Dirichlet boundary condition. We also briefly outline the algebraic formulation arising from the finite element method.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to improve the a priori error estimate of the finite element method by weakly imposing Dirichlet boundary conditions in the bi - harmonic equation with clamped boundary conditions. Specifically, the author proposes a mixed finite element method, which uses the Nitsche technique to weakly impose Dirichlet boundary conditions, thereby obtaining the optimal convergence rate and improving the sub - optimal error estimate of existing methods in the energy norm.
### Background and Motivation
The bi - harmonic equation often appears in the modeling of thin plates and beams, strain - gradient elasticity, phase separation of binary mixtures, and fluid - flow problems. When dealing with such fourth - order elliptic problems, the traditional finite element method leads to sub - optimal error estimates because it needs to impose Dirichlet boundary conditions. Specifically, for a finite element method of order \(k\), its convergence rate in the energy norm is \(h^{k - 1/2}\), rather than the expected \(h^k\).
### Main Contributions
1. **Weak Imposition of Dirichlet Boundary Conditions**:
- The author proposes a method of weakly imposing Dirichlet boundary conditions using the Nitsche technique, thereby avoiding the loss of accuracy caused by imposing boundary conditions.
- This method does not need to introduce a penalty parameter, which simplifies the implementation process.
2. **Optimal Error Estimate**:
- By weakly imposing boundary conditions, the author proves that the new finite element method can achieve the optimal a priori error estimate, that is, the convergence rate is \(h^k\) in the energy norm.
3. **Efficient Algorithm**:
- The author constructs an efficient finite element method using a biorthogonal system, so that all auxiliary variables (such as vorticity \(\phi_h\) and Lagrange multiplier \(p_h\)) can be statically condensed out of the system, resulting in a sparse linear system.
### Mathematical Model and Method
1. **Mixed Finite Element Method**:
- The bi - harmonic equation \(\Delta^2 u = f\) on the region \(\Omega\), with clamped boundary conditions \(u=\frac{\partial u}{\partial n} = 0\).
- By introducing a new unknown \(\phi=\Delta u\), the bi - harmonic equation is transformed into a three - field problem and solved using the mixed finite element method.
2. **Weak Imposition of Dirichlet Boundary Conditions**:
- Using the Nitsche technique, the Dirichlet boundary conditions are weakly imposed in the form of a minimization problem.
- The specific form is:
\[
\int_\Omega \phi \psi \, dx+\langle u, v\rangle_{1/2, \Gamma}=\int_\Omega f v \, dx+\langle g_D, v\rangle_{1/2, \Gamma}
\]
- Here, \(\langle u, v\rangle_{1/2, \Gamma}\) represents the inner product on the boundary.
3. **Numerical Analysis**:
- By introducing an appropriate mesh - dependent norm, the continuity and coercivity of the bilinear form are proved.
- Using the inf - sup condition, the well - posedness of the discrete saddle - point problem is proved.
### Conclusion
The mixed finite element method proposed by the author successfully improves the a priori error estimate of the bi - harmonic equation under clamped boundary conditions by weakly imposing Dirichlet boundary conditions and achieves the optimal convergence rate. In addition, this method is also efficient and easy to implement, providing a new and effective way to solve fourth - order elliptic problems.