Abstract:The numerical approximation of dynamic poroelasticity, modeling flow in deformable porous media, by a family of continuous space-time finite element methods is investigated. Equal order approximation in space without any further stabilization is used for the displacement and pore pressure variable. Optimal order $L^\infty(L^2)$ error estimates are proved and numerically confirmed.
What problem does this paper attempt to address?
This paper attempts to address some key challenges in the numerical approximation of dynamic poroelasticity problems, especially the error estimation and stability issues encountered when using equal - order finite element methods (FEM). Specifically, the goals of the paper are:
1. **Prove optimal - order error estimates**: The author aims to prove that in the case of using equal - order approximation in space (i.e., approximating the displacement \(u\) and pore pressure \(p\) with polynomials of the same order), the continuous space - time finite element method can achieve optimal - order \(L^\infty(L^2)\) error estimates.
2. **Improve computational efficiency**: By using equal - order techniques, the numerical implementation can be simplified and computational efficiency can be improved without the need for additional stabilization terms.
3. **Verify theoretical results**: Verify the proposed theoretical results through numerical experiments to ensure their effectiveness and reliability in practical applications.
### Main contributions of the paper
- **Optimal - order error estimates**: The paper proves that in the dynamic poroelastic model, when using the equal - order finite element method for spatial discretization, optimal - order error estimates can be obtained. This solves the sub - optimal error estimation problem that occurred in previous studies.
- **No need for additional stabilization**: Through an improved energy technique, the paper shows that without introducing additional stabilization terms, the equal - order finite element method can still maintain good stability and convergence.
- **Numerical verification**: Through a series of numerical experiments, the paper verifies the correctness of the theoretical analysis and shows the effectiveness of the equal - order finite element method in practical applications.
### Mathematical model and equations
The dynamic poroelastic model studied in the paper can be represented by the following system of equations:
\[
\begin{cases}
\partial_t u - v = 0, \\
\rho \partial_t v - \nabla \cdot (C \varepsilon(u)) + \alpha \nabla p = \rho f, \\
c_0 \partial_t p + \alpha \nabla \cdot \partial_t u - \nabla \cdot (K \nabla p) = g,
\end{cases}
\]
where:
- \(u\) is the displacement field,
- \(v\) is the velocity field,
- \(p\) is the pore pressure,
- \(\varepsilon(u)=\frac{\nabla u + (\nabla u)^T}{2}\) is the symmetric gradient,
- \(\rho\), \(\alpha\), \(c_0\) are positive constant coefficients,
- \(C\) and \(K\) are symmetric positive definite tensors.
### Conclusion
Through rigorous mathematical derivations and numerical experiments, the paper proves that when using the equal - order finite element method for spatial discretization in dynamic poroelastic problems, optimal - order error estimates can be achieved and no additional stabilization terms are required. This result provides theoretical basis and technical support for the efficient numerical solution of dynamic poroelastic problems.