Simulation of the process of infiltration into fractured porous soil in permafrost

The article provides mathematical modeling of the complex multiphysical problem relevant for the territories of the Far North and the Arctic. The relevance of this task is characterized by importance of the seepage process in the formation and thawing of the permafrost layer. Modern applications for the most part require consideration of complex geometries, as well as a large number of different processes and their mutual relationship. The multiphysical model consists of the Richards equation to describe the seepage process, the double porosity model to describe natural soil fracturing, the Stefan task to describe the temperature regime of the soil in permafrost zone conditions. The computational algorithm is based on finite-element space approximation on triangulated Delone meshes and using of a time splitting scheme using linearization from a previous time layer.