Abstract:This paper presents a numerical method for the simulation of elastic solid materials coupled to fluid inclusions. The application is motivated by the modeling of vascularized tissues and by problems in medical imaging which target the estimation of effective (i.e., macroscale) material properties, taking into account the influence of microscale dynamics, such as fluid flow in the microvasculature. The method is based on the recently proposed Reduced Lagrange Multipliers framework. In particular, the interface between solid and fluid domains is not resolved within the computational mesh for the elastic material but discretized independently, imposing the coupling condition via non-matching Lagrange multipliers. Exploiting the multiscale properties of the problem, the resulting Lagrange multipliers space is reduced to a lower-dimensional characteristic set. We present the details of the stability analysis of the resulting method considering a non-standard boundary condition that enforces a local deformation on the solid-fluid boundary. The method is validated with several numerical examples.
What problem does this paper attempt to address?
### Problems the paper attempts to solve
The paper aims to solve the non - matching coupling problem in multi - scale elasticity problems, especially the interaction between solid materials and fluid inclusions. Specifically, the paper focuses on how to consider the influence of fluid flow in the microvascular network when simulating the mechanical properties of vascularized tissues (such as the liver). Due to geometric complexity and difficulties in handling the fluid - solid interface, traditional fluid - solid coupling models are difficult to achieve efficient and accurate simulations. Therefore, this paper proposes a method based on the simplified Lagrange multiplier framework to handle the mismatched interface conditions between the solid and fluid domains, and verifies the effectiveness and stability of the method through numerical experiments.
### Key point summary
1. **Problem background**:
- **Application fields**: Bio - tissue imaging (such as multi - parameter MRI, magnetic resonance elastography), which requires estimating effective (i.e., macroscopic) material properties.
- **Challenges**: Fully resolved fluid - solid coupling models are computationally expensive and difficult to implement due to geometric complexity and difficulties in handling multi - physics fields.
2. **Method overview**:
- **Simplified Lagrange multiplier framework**: By using Lagrange multipliers on the mismatched interface to impose coupling conditions, the need for precise meshing of the fluid and solid domains is avoided.
- **Multi - scale characteristics**: Utilize the multi - scale characteristics of the problem to reduce the Lagrange multiplier space to a low - dimensional feature set.
- **Local deformation boundary conditions**: A new local deformation Dirichlet boundary condition is introduced, making the coupling conditions independent of the solid deformation at the macroscopic scale.
3. **Theoretical analysis**:
- **Stability analysis**: Analyze the stability of the method in detail, considering a non - standard boundary condition that imposes local deformation on the solid - fluid interface.
- **Numerical verification**: Verify the effectiveness and stability of the method through multiple numerical experiments.
4. **Contributions**:
- **Expand existing methods**: Extend the simplified Lagrange multiplier framework to multi - scale elasticity problems, providing a natural way to handle local deformation boundary conditions.
- **Theoretical and practical applications**: Discuss the reduced - order formula from both theoretical and practical application aspects, and conduct detailed numerical verification in different examples.
### Formula summary
- **Cauchy stress tensor**:
\[
\sigma(u):=\mu\nabla u+\lambda\text{div}(u)I
\]
- **Average trace operator**:
\[
(A_0f)(s): = \frac{1}{|D(s)|}\int_{D(s)}f\,dD(s)
\]
- **Extension operator**:
\[
(E_0w)(x):=(w\circ\Pi)(x)
\]
- **Weak form**:
\[
(\sigma(u),\nabla(v))_\Omega+\langle[\![\sigma(u)]\!]\cdot n,v\rangle_\Gamma=(f,v)_\Omega,\quad\forall v\in H^1_0(\Omega)^d
\]
- **Lagrange multiplier weak form**:
\[
\langle Au,v\rangle+\langle B^T\lambda,v\rangle=\langle f,v\rangle,\quad\forall v\in V_\Omega
\]
\[
\langle Bu,q\rangle=\langle g,q\rangle,\quad\forall q\in Q_\Gamma\setminus Q_0^\Gamma
\]
- **Reduced - dimension Lagrange multiplier space**:
\[
\text{Span}\{\phi_1,\ldots,\phi_N\}
\]
### Conclusion