Numerical Solution of the Two-Phase Stefan Problem in the Enthalpy Formulation with Smoothing the Coefficients

To simulate heat transfer processes with phase transitions, the classical enthalpy model of Stefan is used, accompanied by phase transformations of the medium with absorption and release of latent heat of a change in the state of aggregation. The paper introduces a solution to the two-phase Stefan problem for a one-dimensional quasilinear second-order parabolic equation with discontinuous coefficients. A method for smearing the Dirac delta function using the smoothing of discontinuous coefficients by smooth functions is proposed. The method is based on the use of the integral of errors and the Gaussian normal distribution with an automated selection of the value of the interval of their smoothing by the desired function (temperature). The discontinuous coefficients are replaced by bounded smooth temperature functions. For the numerical solution, the finite difference method and the finite element method with an automated selection of the smearing and smoothing parameters for the coefficients at each time layer are used. The results of numerical calculations are compared with the solution of Stefan’s two-phase self-similar problem --- with a mathematical model of the formation of the ice cover of the reservoir. Numerical simulation of the thawing effect of installing additional piles on the existing pile field is carried out. The temperature on the day surface of the base of the structure is set with account for the amplitude of air temperature fluctuations, taken from the data of the Yakutsk meteorological station. The study presents the results of numerical calculations for concrete piles installed in the summer in large-diameter drilled wells using cement-sand mortars with a temperature of 25 °С. The distributions of soil temperature are obtained for different points in time

This work was supported by the Russian Science Foundation (project RSF no. 19-11-00230)

References

[1] Samarskii A.A., Moiseyenko B.D. An economic continuous calculation scheme for the Stefan multidimensional problem. USSR Comput. Math. Math. Phys., 1965, vol. 5, iss. 5, pp. 43--58. DOI: https://doi.org/10.1016/0041-5553(65)90004-2

[2] Budak B.M., Sobol’eva E.N., Uspenskii A.B. A difference method with coefficient smoothing for the solution of Stefan problems. USSR Comput. Math. Math. Phys., 1965, vol. 5, iss. 5, pp. 59--76. DOI: https://doi.org/10.1016/0041-5553(65)90005-4

[3] Danilyuk I.I. On the Stefan problem. Russian Math. Surveys, 1985, vol. 40, no. 5, pp. 157--223. DOI: https://doi.org/10.1070/RM1985v040n05ABEH003684

[4] Vasilyev V.I. Chislennoe integrirovanie differentsial’nykh uravneniy s nelokal’nymi granichnymi usloviyami. [Numerical integration of differential equations with non-local boundary conditions] Yakutsk, YaF SO AN SSSR Publ., 1985.

[5] Vasilyev V.I., Maksimov A.M., Petrov E.E., et al. Teplomassoperenos v promerzayushchikh i protaivayushchikh gruntakh [Heat and mass transfer in frozen and melted soils]. Moscow, Nauka Publ., 1996.

[6] Samarskiy A.A., Vabishchevich P.N. Vychislitel’naya teploperedacha [Computational heat and mass transfer]. Moscow, Librokom Publ., 2009.

[7] Krylov D.A., Sidnyaev N.I., Fedotov A.A. Mathematical modelling of temperature distribution. Matematicheskoe modelirovanie, 2013, vol. 25, no. 7, pp. 3--27 (in Russ.).

[8] Vasilyeva M.V., Stepanov S.P., Sirditov I.K. Reduced dimension model for heat transfer of ground heat exchange in permafrost. J. Phys.: Conf. Ser., 2017, vol. 937, art. 012056. DOI: https://doi.org/10.1088/1742-6596/937/1/012056

[9] Stepanov S., Vasilyeva M., Vasil’ev V. Generalized multiscale discontinuous Galerkin method for solving the heat problem with phase change. J. Comput. Appl. Math., 2018, vol. 340, pp. 645--652. DOI: https://doi.org/10.1016/j.cam.2017.12.004

[10] Sivtsev P.V., Kolesov A.E., Sirditov I.K., et al. The numerical solution of thermo-poroelastoplasticity problems. AIP Conf. Proc., 2016, vol. 1773, art. 110010. DOI: https://doi.org/10.1063/1.4965014

[11] Stepanov S.P., Vasilyeva M.V., Vasilyev V.I. Numerical simulation of the convective heat transfer on high-performance computing systems. AIP Conf. Proc., 2016, vol. 1773, art. 110011. DOI: https://doi.org/10.1063/1.4965015

[12] Vasilyev V.I., Vasilyeva M.V., Sirditov I.K., et al. Mathematical modeling of temperature regime of soils of foundation on permafrost. Herald of the Bauman Moscow State Technical University, Series Natural Sciences, 2017, no. 1 (70), pp. 142--159 (in Russ.). DOI: https://doi.org/10.18698/1812-3368-2017-1-142-159

[13] Vasilyeva M., Stepanov S., Spiridonov D., et al. Multiscale Finite Element Method for heat transfer problem during artificial ground freezing. J. Comput. Appl. Math., 2020, vol. 371, art. 112605. DOI: https://doi.org/10.1016/j.cam.2019.112605

[14] Vasil’ev V., Vasilyeva M. An accurate approximation of the two-phase Stefan problem with coefficient smoothing. Mathematics, 2020, vol. 8, iss. 11, art. 1924. DOI: https://doi.org/10.3390/math8111924

[15] Samarskiy A.A. Vvedenie v teoriyu raznostnykh skhem [Introduction into theory of difference schemes]. Moscow, Nauka Publ., 1971.