TY - JOUR

T1 - Numerical modeling of non-stationary heat problems in a two-phase medium

AU - Gladkikh, V. S.

AU - Ilin, V. P.

AU - Petukhov, A. V.

AU - Krylov, A. M.

N1 - Funding Information: The work was carried out using the resources of the Siberian Supercomputer Center (CKP SSCC ICM&MG SB RAS), supported by RFBR grant 18-01-00295 in theoretical part and RSF grant 19-11-00048 in computational experiments part. Publisher Copyright: © 2021 Institute of Physics Publishing. All rights reserved. Copyright: Copyright 2021 Elsevier B.V., All rights reserved.

PY - 2021/1/4

Y1 - 2021/1/4

N2 - The simulation of thermal fields and melting permafrost is an important step for engineering construction in many regions of Russia. Water-ice phase transitions in a 3-D complicated computation domain must be taken into account. The enthalpy statement with an implicit finite volume method was used to create a good discretization and a local mesh refinement to focus an area near the oil well. System of linear algebraic equations is solved using the iterative conjugate residual or conjugate gradient method in Krylov subspaces with the incomplete factorization algorithm in the Eisenstat modification as a preconditioner. A special high-parallel version of matrix generation and solving a system of the linear algebraic equations code have been developed and then efficiency has been estimated. The performance results of the developed code and temperature fields for different wells configuration (1 well and 4 wells placed as squared) during 5 years have been presented. The results validation was based on the comparison with previous papers and other algorithms. The results of simulations are close those to presented by other authors.

AB - The simulation of thermal fields and melting permafrost is an important step for engineering construction in many regions of Russia. Water-ice phase transitions in a 3-D complicated computation domain must be taken into account. The enthalpy statement with an implicit finite volume method was used to create a good discretization and a local mesh refinement to focus an area near the oil well. System of linear algebraic equations is solved using the iterative conjugate residual or conjugate gradient method in Krylov subspaces with the incomplete factorization algorithm in the Eisenstat modification as a preconditioner. A special high-parallel version of matrix generation and solving a system of the linear algebraic equations code have been developed and then efficiency has been estimated. The performance results of the developed code and temperature fields for different wells configuration (1 well and 4 wells placed as squared) during 5 years have been presented. The results validation was based on the comparison with previous papers and other algorithms. The results of simulations are close those to presented by other authors.

UR - http://www.scopus.com/inward/record.url?scp=85100783076&partnerID=8YFLogxK

U2 - 10.1088/1742-6596/1715/1/012002

DO - 10.1088/1742-6596/1715/1/012002

M3 - Conference article

AN - SCOPUS:85100783076

VL - 1715

JO - Journal of Physics: Conference Series

JF - Journal of Physics: Conference Series

SN - 1742-6588

IS - 1

M1 - 012002

T2 - International Conference on Marchuk Scientific Readings 2020, MSR 2020

Y2 - 19 October 2020 through 23 October 2020

ER -