Charles M. Elliott,Achilleas Mavrakis
Abstract:We consider the surface Stokes equation with Lagrange multiplier and approach it numerically. Using a Taylor-Hood surface finite element method, along with an appropriate estimation for the additional Lagrange multiplier, we establish optimal convergence results for the velocity in $H^1$ and $L^2$, and in $L^2$ for the two pressures. Furthermore, we present a new inf-sup condition to help with the stability and convergence results. This formulation offers the advantage of simplified implementation without approximating geometric quantities, although it requires a higher approximation of the parameterized surface. In addition, we provide numerical simulations that confirm the established error bounds and perform a comparative analysis against the penalty approach.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is about the numerical solution of the surface Stokes problem, especially the discretization and numerical simulation of the surface Stokes equations with Lagrange multipliers. Specifically, the author focuses on how to effectively numerically solve the Stokes equations defined on the surface while maintaining numerical stability and convergence.
### Specific description of the problem
1. **Mathematical model**:
- The object of study is a steady - state generalized surface Stokes problem defined on a two - dimensional smooth closed surface \(\Gamma\). The given system of equations is:
\[
\begin{cases}
- 2\mu \operatorname{div}_{\Gamma}(E(u))+u + \lambda n_{\Gamma}+\nabla_{\Gamma}p = f, \\
\operatorname{div}_{\Gamma}u = g, \\
u\cdot n_{\Gamma}= 0,
\end{cases}
\]
where \(u\) is the velocity field, \(p\) is the surface pressure, \(\lambda\) is the Lagrange multiplier, \(\mu\) is the viscosity coefficient, \(f\) is the given force field, \(g\) is the given source term, \(n_{\Gamma}\) is the unit normal vector of the surface, and \(E(u)\) is the strain rate tensor (deformation tensor).
2. **Numerical method**:
- The author adopts the Taylor - Hood surface finite element method (SFEM) to discretize the above problem. To ensure the stability and convergence of the numerical method, they introduce an additional Lagrange multiplier \(\lambda_h\) to weakly impose the tangential constraint conditions.
- Through this method, the author is able to establish optimal \(H^1\) and \(L^2\) convergence results, and no additional approximation of geometric quantities (such as mean curvature or Gaussian curvature) is required.
3. **Innovation points**:
- A new inf - sup condition is proposed to assist in the stability and convergence analysis.
- Numerical simulation results are provided to verify the proposed error bounds and a comparative analysis with the penalty function method is carried out.
### Summary
The core problem of this paper is to develop an efficient and stable numerical method to solve the surface Stokes equations with Lagrange multipliers. By introducing an additional Lagrange multiplier and using the Taylor - Hood surface finite element method, the author has successfully achieved this goal and provided theoretical analysis and numerical verification.