Simulation of heat transfer with considering permafrost thawing in 3D media

Abstract

An approach to mathematical modeling of heat transfer with a permafrost algorithm in 3D media based on the idea of localizing the phase transition area is considered. The paper presents a problem statement for a non-stationary heat transfer and a description of a numerical method based on a predictor-corrector scheme. For a better understanding of the proposed splitting method, the accuracy order of approximation considering inhomogeneous right-hand side was studied. The phase changes in the numerical implementation of permafrost thawing is considered in the temperature range and requires recalculation of coefficients values of the heat equation at each iteration step with respect to time. A brief description of the parallel algorithm based on a 3D decomposition method and the parallel sweep method is presented. A study of the parallel algorithm implementations using a high-performance computing system of the Siberian Supercomputer Center of the SB RAS was performed. The results of the permafrost algorithm on models with wellbores are also presented.