Abstract
This paper aims to explore a robust and efficient discrete scheme for poroelastic-elastic coupled problems. Based on staggered grids, a discrete scheme using finite difference methods is constructed in a form similar to domain decomposition. On this basis, another discrete scheme with a uniform form over the entire domain is derived. These two discrete schemes serve different purposes and are shown to be equivalent. The stability is easily derived by establishing a discrete variational formulation. To prove the convergence, appropriate discrete interpolations are introduced. It is proved that the discrete scheme has second-order superconvergence. Then, combined with the inf-sup condition, the first-order uniform convergence is obtained. This means the discrete scheme has great potential to overcome Poisson locking and pressure oscillations. Some numerical experiments are also carried out, and the results show the robustness and efficiency of the discrete scheme.
Similar content being viewed by others
Data Availability
Enquiries about data availability should be directed to the authors.
References
Anaya, V., De, W., Zoa, G.-V., Bryan, M., David, R.-B.: Rotation-based mixed formulations for an elasticity-poroelasticity interface problem. SIAM J. Scientif. Comput. 42(1), B225–B249 (2020)
Babuška, I., Suri, M.: Locking effects in the finite element approximation of elasticity problems. Numerische Mathematik 62(1), 439–463 (1992)
Babuška, I., Suri, M.: On locking and robustness in the finite element method. SIAM J. Numer. Anal. 29(5), 1261–1293 (1992)
Biot, M.-A.: Consolidation settlement under a rectangular load distribution. J. Appl. Phys. 12(5), 426–430 (1941)
Biot, M.-A.: General theory of three-dimensional consolidation. J. Appl. Phys. 12(2), 155–164 (1941)
Ewing, R.E., Iliev, O.P., Lazarov, R.D., Naumovich, A.: On convergence of certain finite volume difference discretizations for 1D poroelasticity interface problems. Numer. Methods Partial. Diff. Eq.: Int. J. 23(3), 652–671 (2007)
Girault, V., Pencheva, G.V., Wheeler, M.F.: Domain decomposition for linear elasticity with DG jumps and mortars. Computer Methods Appl. Mech. Eng. 198(21–26), 1751–1765 (2009)
Girault, V, Wheeler, M F., Almani, T, Dana, S: A priori error estimates for a discretized poro-elastic–elastic system solved by a fixed-stress algorithm. Oil Gas Sci. Technol.–Revue d’IFP Energies nouvelles 74, 24 (2019)
Haga, J.B., Osnes, H., Langtangen, H.P.: On the causes of pressure oscillations in low-permeable and low-compressible porous media. Int. J. Numer. Anal. Methods Geomech. 36(12), 1507–1522 (2012)
Lebedev, V.I.: Difference analogues of orthogonal decompositions, basic differential operators and some boundary problems of mathematical physics. I. USSR Comput. Math. Math. Phys. 4(3), 69–92 (1964)
Mikelić, Andro W., Mary, F.: On the interface law between a deformable porous medium containing a viscous fluid and an elastic body. Math. Models Methods Appl. Sci. 22, 1250031 (2012)
Monk, P., Süli, E.: A convergence analysis of Yee’s scheme on nonuniform grids. SIAM J. Numer. Anal. 31(2), 393–412 (1994)
Oyarzúa, R., Ruiz-Baier, R.: Locking-free finite element methods for poroelasticity. SIAM J. Numer. Anal. 54(5), 2951–2973 (2016)
Phillips, P.J., Wheeler, M.F.: A coupling of mixed and discontinuous Galerkin finite-element methods for poroelasticity. Comput. Geosci. 12(4), 417–435 (2008)
Rui, H., Li, X.: Stability and superconvergence of MAC scheme for Stokes equations on nonuniform grids. SIAM J. Numer. Anal. 55(3), 1135–1158 (2017)
Rui, H., Liu, W.: A two-grid block-centered finite difference method for Darcy-Forchheimer flow in porous media. SIAM J. Numer. Anal. 53(4), 1941–1962 (2015)
Rui, H., Pan, H.: A block-centered finite difference method for the Darcy-Forchheimer model. SIAM J. Numer. Anal. 50(5), 2612–2631 (2012)
Rui, H., Sun, M.: A locking-free finite difference method on staggered grids for linear elasticity problems. Computers Math. Appl. 76(6), 1301–1320 (2018)
Rui, H., Sun, Y.: A MAC scheme for coupled stokes-darcy equations on non-uniform grids. J. Scientif. Comput. 82(3), 1–29 (2020)
Serepilli, M.: Classical and higher order interface conditions in poroelasticity. Annals Solid Struct. Mech. 11(1), 1–10 (2019)
Welch, J Eddie, Harlow, Francis Harvey, Shannon, John P, Daly, Bart J: The MAC method-a computing technique for solving viscous, incompressible, transient fluid-flow problems involving free surfaces. Technical report, Los Alamos National Lab.(LANL), Los Alamos, NM (United States), (1965)
Wheeler, M.-F.: A multipoint flux mixed finite element method. SIAM J. Numer. Anal. 44(5), 2082–2106 (2006)
Yee, K.: Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media. IEEE Trans. Antennas Propagation 14(3), 302–307 (1966)
Yi, S.-Y.: A study of two modes of locking in poroelasticity. SIAM J. Numer. Anal. 55(4), 1915–1936 (2017)
Funding
The work is supported by the National Natural Science Foundation of China Grant No.12131014.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors have not disclosed any competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Yang, B., Rui, H. A Robust Discrete Scheme based on Staggered Grids for Poroelastic-Elastic Coupled Problems. J Sci Comput 95, 25 (2023). https://doi.org/10.1007/s10915-023-02149-x
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1007/s10915-023-02149-x
Keywords
- Poroelasticity-elasticity coupling
- Interface problem
- Staggered grids
- Finite difference
- Superconvergence
- Uniform convergence