Papers › An efficient preconditioner for mixed-dimensional contact poromechanics based on the...
An efficient preconditioner for mixed-dimensional contact poromechanics based on the fixed stress splitting scheme
Yury Zabegaev, Inga Berre, Eirik Keilegavlen, Kundan Kumar
The archive published only this paper's code-link row. Authors, date and abstract are from arXiv's metadata (CC0), read from the Kaggle arXiv metadata snapshot of 2026-09-12 where its title matched the archive's; the title is the archive's.
Numerical simulation of fracture contact poromechanics is essential for various applications, including CO2 sequestration, geothermal energy production and underground gas storage. Modeling this problem accurately presents significant challenges due to the complex physics involved in strongly coupled poromechanics and frictional contact mechanics of fractures. The robustness and efficiency of the simulation heavily depends on a preconditioner for the linear solver, which addresses the Jacobian matrices arising from Newton's method in fully implicit time-stepping schemes. Developing an effective preconditioner is difficult because it must decouple three interdependent subproblems: momentum balance, fluid mass balance, and contact mechanics. The challenge is further compounded by the saddle-point structure of the contact mechanics problem, a result of the Augmented Lagrange formulation, which hinders the direct application of the well-established fixed stress approximation to decouple the poromechanics subproblem. In this work, we propose a preconditioner hat combines nested Schur complement approximations with a linear transformation, which addresses the singular nature of the contact mechanics subproblem. This approach extends the fixed stress scheme to both the matrix and fracture subdomains. We investigate analytically how the contact mechanics subproblem affects the convergence of the proposed fixed stress-based iterative scheme and demonstrate how it can be translated into a practical preconditioner. The scalability and robustness of the method are validated through a series of numerical experiments.
Code
Repository list and official/mentioned flags are the archive's, frozen 2025-07-28. Reachability, where shown, is from one Syntology probe window (2026-09-16 to 2026-09-18); repositories not probed show nothing. GitHub stars are not tracked.
Code Syntology ran Syntology
Not run by Syntology. Nothing on this page verifies that the listed code works.
Results from the paper archive 2025-07-28
No leaderboard rows for this paper in the archive.
Report a problem or propose a change · a person checks every report against the paper or source before anything changes; decisions are listed on /corrections