Iteratively coupled mixed finite element solver for thermo-hydro-mechanical modeling of permafrost thaw

• Published in: Results in Applied Mathematics
• Main Authors: Naren Vohra, Malgorzata Peszynska
• Format: Article in Journal/Newspaper
• Language: English
• Published: Elsevier 2024
• Subjects: Thermo-hydro-mechanical models; Heterogeneous permafrost; Heat conduction with phase change; Poroelasticity; Mixed finite elements; Nonlinear constitutive relationships; Mathematics; QA1-939; permafrost
• Online Access: https://doi.org/10.1016/j.rinam.2024.100439; https://doaj.org/article/d09af7ac45ed4890ac0897c9a8df3d86

In this paper we consider computational challenges associated with thermo-hydro-mechanical models for simulation of subsidence due to permafrost thaw. The model we outline couples heat conduction with phase change and thermal advection to Biot’s poroelasticity equations with attention paid to the dependence of the constitutive parameters on temperature. Our numerical scheme uses the lowest order mixed finite elements for discretization of thermal and hydrological flow, and Galerkin finite elements for mechanics, and uses an implicit–explicit time stepping. We set up an iterative solver that solves the thermal subproblem followed by the hydro-mechanical subproblem, and demonstrate its robustness in practical heterogeneous permafrost scenarios. We also identify the challenges associated with the roughness of the dependence of mechanical parameters on the temperature.