Abstract:We study a fictitious domain approach with Lagrange multipliers to discretize Stokes equations on a mesh that does not fit the boundaries. A mixed finite element method is used for fluid flow. Several stabilization terms are added to improve the approximation of the normal trace of the stress tensor and to avoid the inf-sup conditions between the spaces of the velocity and the Lagrange multipliers. We generalize first an approach based on eXtended Finite Element Method due to Haslinger-Renard involving a Barbosa-Hughes stabilization and a robust reconstruction on the badly cut elements. Secondly, we adapt the approach due to Burman-Hansbo involving a stabilization only on the Lagrange multiplier. Multiple choices for the finite elements for velocity, pressure and multiplier are considered. Additional stabilization on pressure (Brezzi-Pitkäranta, Interior Penalty) is added, if needed. We prove the stability and the optimal convergence of several variants of these methods under appropriate assumptions. Finally, we perform numerical tests to illustrate the capabilities of the methods.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: in the context of fluid - structure interaction, especially when simulating granular flow, how to stably and efficiently discretize the Stokes equations on a grid that does not fit the boundary. Specifically, the author has studied a fictitious domain approach based on Lagrange multipliers for discretizing the Stokes equations on a grid that is not aligned with the boundary. To improve the approximation of the normal trace of the stress tensor and avoid the inf - sup condition problem between the velocity and Lagrange multiplier spaces, several stabilization terms are introduced in the paper.
### Main problem description
1. **Discretization of the Stokes equations**:
- Researchers are concerned with how to discretize the Stokes equations on a grid that does not fit the boundary in the context of fluid - structure interaction, especially when simulating the movement of rigid or elastic particles in a fluid.
- The Stokes equations are as follows:
\[
- 2\text{div}D(u)+\nabla p = f\quad\text{in }F,\tag{1}
\]
\[
\text{div}u = 0\quad\text{in }F,\tag{2}
\]
\[
u = g\quad\text{on }\Gamma,\tag{3}
\]
\[
u = g_{\text{wall}}\quad\text{on }\Gamma_{\text{wall}},\tag{4}
\]
where \(D(u)=\frac{1}{2}(\nabla u+(\nabla u)^T)\), assuming the viscosity is 1 for simplicity.
2. **Introduction of Lagrange multipliers**:
- Introduce the force \(\lambda\) exerted by the fluid on the solid on the interface \(\Gamma\):
\[
\lambda=- 2D(u)n+pn\quad\text{on }\Gamma,\tag{5}
\]
where \(n\) is the unit normal vector pointing outward from the fluid domain \(F\), and \(\lambda\) is regarded as the Lagrange multiplier associated with the Dirichlet condition (3).
3. **Stabilization techniques**:
- The paper proposes several stabilization techniques to ensure the stability of discretization, including Barbosa - Hughes stabilization, local treatment in the Haslinger - Renard method, Lagrange multiplier stabilization in the Burman - Hansbo method, etc.
- For different types of finite - element spaces (such as P1 - P1, P1 - P0, Taylor - Hood, etc.), the corresponding stabilization methods are studied, and the stability and optimal convergence of these methods are proved.
### Solution
By introducing appropriate stabilization terms and using the extended finite - element method (XFEM), the author proposes a method that can effectively solve the Stokes equations on a grid that does not fit the boundary. This method not only avoids the need for remeshing but also can handle complex fluid - structure interaction problems, especially when the interface moves over time.
In summary, the main objective of this paper is to develop a robust and efficient numerical method that can stably solve the Stokes equations on a grid that does not fit the boundary and apply it to fluid - structure interaction problems.