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

• Stepanov Sergei P., cepe2a@inbox.ru Ammosov North-Eastern Federal University, Institute of Mathematics and Informatics, 42 Kulakovsky Street, Yakutsk 677000, Russia
• Grigoriev Aleksandr V., re5itsme@gmail.com Ammosov North-Eastern Federal University, Institute of Mathematics and Informatics, 42 Kulakovsky Street, Yakutsk 677000, Russia
• Afanasyeva Nadezhda M., afanasieva.nm@gmail.com Ammosov North-Eastern Federal University, Institute of Mathematics and Informatics, 42 Kulakovsky Street, Yakutsk 677000, Russia
Keywords: the Richards equation, the Stefan problem, double porosity, fractured porous media

### Abstract

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.

### References

[1]
Samarskii A. A. and Moiseenko B. D., “Through calculation scheme for the multidimensional Stefan problem,” Comput. Math. Math. Phys., 5, No. 5, 816–827 (1965).

[2]
Samarskii A. A., Vabishchevich P. N., Iliev O. P., and Churbanov A. G., “Numerical simulation of convection/diffusion phase change problems: a review,” Int. J. Heat Mass Transfer, 36, No. 17, 4095–4106 (1993).

[3]
Samarskii A. A. and Vabishchevich P. N., Mathematical Modelling, vol. 1, Computational Heat Transfer, Wiley (1995).

[4]
Cheng Q., Sun Y., Jones S. B. et al., “In situ measured and simulated seasonal freeze-thaw cycle: A 2-year comparative study between layered and homogeneous field soil profiles,” J. Hydrology, 519, 1466–1473 (2014).

[5]
Vasil’ev V. I., Vasilyeva M. V., Stepanov S. P. et al., “Mathematical modeling of the temperature regime of soils of basement foundation on permafrost [in Russian],” Vestn. Mosk. Gos. Tekhn. Univ. im. N. E. Baumana, Ser. Estestv. Nauki, No. 1, 142–159 (2017).

[6]
Vabishchevich P. N., Varlamov S. P., Vasiliev V. I. et al., “Numerical simulation of the temperature dynamics of railway foundation material in permafrost,” Math. Models Comput. Simul., 9, No. 3, 292–304 (2017).

[7]
Ivanov V. A., “Numerical research on effect of thermal insulation on a gas pipeline’s performance in Far North environment,” Mat. Zamet. SVFU, 24, No. 4, 96–108 (2017). doi: https://doi.org/10.25587/SVFU.2018.4.11320.

[8]
Masters I., Pao W. K. S., and Lewis R. W., “Coupling temperature to a double–porosity model of deformable porous media,” Int. J. Numer. Methods Engin., 49, No. 3, 421–438 (2000).

[9]
Helmig R., Multiphase Flow and Transport Processes in the Subsurface, A Contribution to the Modeling of Hydrosystems, Springer-Verl., Berlin (1997).

[10]
Bai M. and Roegiers J.-C., “Fluid flow and heat flow in deformable fractured porous media,” Int. J. Eng. Sci., 32, No. 10, 1615–1633 (1994).

[11]
Vasilyeva M. V. and Prokopiev G. A., “Numerical solution to the problem of two-phase filtration with heterogeneous coefficients by the finite element method,” Mat. Zamet. SVFU, 24, No. 2, 46–62 (2017). doi: https://doi.org/10.25587/SVFU.2017.2.9245.

[12]
Langtangen H. P., “A FEniCS tutorial,” in: Automated Solution of Differential Equations by the Finite Element Method, pp. 1–73, Springer-Verl., Berlin; Heidelberg (2012).

[13]
Samarskii A. A., The Theory of Difference Schemes, Marcel Dekker, New York; Basel (2001).

[14]
Barenblatt G. I., Zheltov I. P., and Kochina I. N., “Basic concepts in the theory of seepage of homogeneous liquids in fissured rocks,” J. Appl. Math. Mech., 24, No. 5, 1286–1303 (1960).

[15]
Rathfelder K. and Abriola L. M., “Mass conservative numerical solutions of the head-based Richards equation,” Water Resources Res., 30, No. 9, 2579–2586 (1994).

[16]
Ross P. J., “Efficient numerical methods for infiltration using Richards equation,” Water Resources Res., 26, No. 2, 279–290 (1990).

[17]
Celia M. A., Bouloutas E. T., and Zarba R. L., “A general mass-conservative numerical solution for the unsaturated flow equation,” Water Resources Res., 26, No. 7, 1483–1496 (1990).

[18]
Van Genuchten M. T., “A closed-form equation for predicting the hydraulic conductivity of unsaturated soils 1,” Soil Sci. Soc. America J., 44, No. 5, 892–898 (1980).

[19]
Vasiliev V. I., Maksimov A. M., Petrov E. E., and Tsypkin G. G., Heat and Mass Transfer in Freezing and Thawing Soils, Nauka, Moscow (1996).

[20]
Droniou J., Eymard R., and Guichard C., “Uniform-in-time convergence of numerical schemes for Richards’ and Stefan’s models,” in: Finite Volumes for Complex Applications VII: Methods and Theoretical Aspects, pp. 247–254, Springer, Cham (2014).
How to Cite
Stepanov, S., Grigoriev, A. and Afanasyeva, N. (2020) “Simulation of the process of infiltration into fractured porous soil in permafrost”, Mathematical notes of NEFU, 27(2), pp. 105-117. doi: https://doi.org/10.25587/SVFU.2020.15.67.007.
Issue
Section
Mathematical Modeling