Stephan B. Lunowa,Barbara Wohlmuth
Abstract:Biot's consolidation model is a classical model for the evolution of deformable porous media saturated by a fluid and has various interdisciplinary applications. While numerical solution methods to solve poroelasticity by typical schemes such as finite differences, finite volumes or finite elements have been intensely studied, lattice Boltzmann methods for poroelasticity have not been developed yet. In this work, we propose a novel semi-implicit coupling of lattice Boltzmann methods to solve Biot's consolidation model in two dimensions. To this end, we use a single-relaxation-time lattice Boltzmann method for reaction-diffusion equations to solve the Darcy flow and combine it with a recent pseudo-time multi-relaxation-time lattice Boltzmann scheme for quasi-static linear elasticity by Boolakee, Geier and De Lorenzis (2023, DOI: <a class="link-https link-external" data-doi="10.1016/j.cma.2022.115756" href="https://doi.org/10.1016/j.cma.2022.115756" rel="external noopener nofollow">https://doi.org/10.1016/j.cma.2022.115756</a>). The numerical results demonstrate that naive coupling schemes lead to instabilities when the poroelastic system is strongly coupled. However, the newly developed centered coupling scheme using fully explicit and semi-implicit contributions is stable and accurate in all considered cases, even for the Biot--Willis coefficient being one. Furthermore, the numerical results for Terzaghi's consolidation problem and a two-dimensional extension thereof highlight that the scheme is even able to capture discontinuous solutions arising from instantaneous loading.
What problem does this paper attempt to address?
This paper aims to solve the numerical solution problem of the Biot consolidation model in porous elastic media. Specifically, the paper proposes a novel lattice Boltzmann method (LBM) to solve Biot's consolidation model, especially in the two - dimensional case. The following are the main problems and solutions in the paper:
### Problems the paper attempts to solve
1. **Lack of lattice Boltzmann method applicable to porous elasticity**
- Although traditional numerical methods such as the finite difference method, the finite volume method, and the finite element method have been widely used in porous elasticity problems, no lattice Boltzmann method specifically for porous elasticity has been developed yet.
2. **Stability problems of strongly - coupled systems**
- When the Biot - Willis coefficient is close to 1, the traditional simple coupling scheme will lead to numerical instability. The paper proposes a new central coupling scheme, which ensures the stability and accuracy of the numerical method through the combination of fully explicit and semi - implicit contributions.
3. **Handling discontinuous solutions caused by instantaneous loading**
- The paper shows that the proposed scheme can capture the discontinuous solutions caused by instantaneous loading, which is difficult to achieve with traditional methods.
### Main contributions
1. **Proposing a new lattice Boltzmann method**
- Use the single - relaxation - time lattice Boltzmann method to solve Darcy flow, and combine the recent pseudo - time multi - relaxation - time lattice Boltzmann scheme to solve the quasi - static linear elasticity problem.
2. **Solving the stability problem of strongly - coupled systems**
- By introducing a central coupling scheme with fully explicit and semi - implicit contributions, the stability and accuracy of the numerical method under all Biot - Willis coefficients are ensured.
3. **Verifying the effectiveness of the method**
- Through numerical results, especially the Terzaghi consolidation problem and its two - dimensional extension, it is proved that this method can capture the discontinuous solutions caused by instantaneous loading.
### Mathematical models and methods
1. **Basic equations of the Biot consolidation model**
- Force balance equation:
\[
-\nabla \cdot \sigma+\alpha \nabla p = f\quad \text{in }Q
\]
- Stress - strain relationship:
\[
\sigma=\lambda(\nabla \cdot \eta)I+\mu(\nabla \eta+\nabla \eta^{\top})\quad \text{in }Q
\]
- Fluid conservation equation:
\[
\partial_{t}(c_{0}p+\alpha \nabla \cdot \eta)+\nabla \cdot u = s\quad \text{in }Q
\]
- Darcy flow equation:
\[
u = -\kappa(\nabla p - g)\quad \text{in }Q
\]
2. **Dimensionless system**
- Introduce characteristic length \(L\), time \(T\), displacement \(H\) and pressure \(P = c_{0}^{-1}\), and define dimensionless variables:
\[
\hat{x}=\frac{x}{L},\quad \hat{t}=\frac{t}{T},\quad \hat{\eta}=\frac{\eta}{H},\quad \hat{\sigma}=\frac{c_{0}L\sigma}{H},\quad \hat{p}=c_{0}p,\quad \hat{u}=\frac{T u}{L}
\]
- Dimensionless equations:
\[
-\hat{\nabla}\cdot \hat{\sigma}+\varepsilon \alpha \h