Álvaro Pé de la Riva,Francisco J. Gaspar,Xiaozhe Hu,James Adler,Carmen Rodrigo,Ludmil Zikatanov
Abstract:In this work, we present a new stabilization method aimed at removing spurious oscillations in the pressure approximation of Biot's model for poroelasticity with low permeabilities and/or small time steps. We consider different finite-element discretizations and illustrate how not only does such a stabilized scheme provide numerical solutions that are free of non-physical oscillations, but it also allows one to iterate the fluid and mechanics problems in a fashion similar to the well-known fixed-stress split method. The resulting solution method is convergent without the necessity for additional terms to stabilize the iteration. Finally, we present numerical results illustrating the robust behavior of both the stabilization and iterative solver with respect to the physical and discretization parameters of the model.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is: in the case of low permeability and/or small time steps, the non - physical oscillation problem in the pressure approximation in the Biot model. Specifically, the author proposes a new stabilization method, aiming to eliminate these non - physical oscillations and ensure the stability of the numerical solution.
### Problem Background
The Biot model describes the coupling relationship between fluid flow and mechanical deformation in porous media and is widely used in multiple fields, such as geothermal energy extraction, CO2 storage, hydraulic fracturing, and cancer research. However, in some cases (for example, low permeability or small time steps), the standard finite - element discretization method will lead to the appearance of non - physical oscillation phenomena in the pressure field. These oscillations not only affect the accuracy of the numerical solution but may also lead to unreliable calculation results.
### Paper Objectives
To solve the above - mentioned problems, the paper proposes a new stabilization method. This method eliminates non - physical oscillations by adding specific stabilization terms. In addition, this method also allows the fluid and mechanical problems to be solved iteratively in a manner similar to the fixed - stress splitting method without the need for additional stabilization terms to ensure iterative convergence.
### Main Contributions
1. **Elimination of Non - physical Oscillations**: The new method can provide oscillation - free pressure approximation in the case of low permeability and/or small time steps.
2. **Iterative Coupled Solution**: This method allows the fluid and mechanical problems to be solved iteratively in an alternating manner, similar to the fixed - stress splitting method, but without the need for additional stabilization terms.
3. **Parameter Robustness**: The proposed method has good robustness to changes in physical and discretization parameters and is suitable for a wider range of practical applications.
### Mathematical Formulas
To better understand this method, the following are the mathematical formulas related to stabilization:
- Time - derivative term of fluid pressure:
\[
\bar{\partial}_t p^n_h=\frac{p^n_h - p^{n - 1}_h}{\tau}
\]
- Stabilization term:
\[
L(\bar{\partial}_t p^n_h, q_h)_0 - L(\bar{\partial}_t p^n_h, q_h)
\]
where \((\cdot,\cdot)_0\) is the mass - lumped approximation of the \(L^2(\Omega)\) inner product.
Through these improvements, the paper successfully solves the common non - physical oscillation problem in the Biot model and provides an efficient and stable numerical solution method.