Stefano Bonetti,Michele Botti,Ilario Mazzieri,Paola F. Antonietti
Abstract:We present and analyze a high-order discontinuous Galerkin method for the space discretization of the wave propagation model in thermo-poroelastic media. The proposed scheme supports general polytopal grids. Stability analysis and $hp$-version error estimates in suitable energy norms are derived for the semi-discrete problem. The fully-discrete scheme is then obtained based on employing an implicit Newmark-$\beta$ time integration scheme. A wide set of numerical simulations is reported, both for the verification of the theoretical estimates and for examples of physical interest. A comparison with the results of the poroelastic model is provided too, highlighting the differences between the predictive capabilities of the two models.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the numerical modeling and analysis of wave propagation phenomena in thermo - poroelastic media. Specifically, the authors propose and analyze a high - order discontinuous Galerkin method for spatial discretization to simulate wave propagation phenomena in thermo - poroelastic media. This method supports general polyhedral meshes, and a stability analysis and hp - version error estimate for the semi - discrete problem are carried out. In addition, a fully - discrete scheme is obtained by using the implicit Newmark - β time - integration scheme. To verify the validity of the theoretical estimates and demonstrate the physical application value of this model, a wide range of numerical simulation experiments are reported in the paper. At the same time, a comparison with the results of the poroelastic model is also provided, highlighting the differences between the prediction capabilities of the two models.
### Main contributions
1. **Accurate model derivation and mathematical formulation**: The mathematical model of the fully - dynamic thermo - poroelastic problem is derived in detail.
2. **PolyDG discretization analysis**: A comprehensive analysis of PolyDG discretization is carried out, and stability and error estimates are established.
3. **In - depth numerical study**: The convergence of the method is verified through numerical experiments, and the applicability of the model in complex media is demonstrated, especially for wave propagation phenomena considering the temperature effect.
### Mathematical model
The model consists of three equations: mass conservation, momentum conservation, and energy conservation. It is assumed that the medium is a linear thermo - poroelastic medium, and inertial terms are included in all three equations. The poroelastic part is expressed in the two - displacement form, that is, the unknowns are solid displacement and filtration displacement, and the pressure field is obtained by post - processing.
### Discretization method
- **Spatial discretization**: The PolyDG finite - element method is used for spatial discretization. This method has arbitrary - order accuracy and high geometric flexibility on polyhedral meshes, and is suitable for dealing with complex geometries, local low - regularity, and layered heterogeneous materials.
- **Time discretization**: The implicit Newmark - β time - integration scheme is used.
### Numerical experiments
- **Convergence test**: The convergence of the method is verified, including errors in the L2 norm and dG norm.
- **Physical test**: The numerical performance is evaluated through a series of physically meaningful test cases, and a comparison with the results of the poroelastic model is made.
### Conclusion
By proposing a high - order discontinuous Galerkin method, this paper successfully solves the numerical modeling problem of wave propagation phenomena in thermo - poroelastic media, and verifies the effectiveness and accuracy of the method through extensive numerical experiments.