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

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.