Application of the enthalpy‑porosity method in studies on solid‑liquid phase change interfaces

doi:10.3969/j.issn.1007-5461.2026.03.001

Abstract

Abstract: Solid-liquid phase transition is a fundamental problem in material science research, and numerical simulation of solid-liquid interfaces is of great significance for the growth of large-size and high-quality crystals. The enthalpy-porosity method uses a unified control equation to solve problems in the liquid and solid computational domains, without the need to track the liquid-solid interface separately, which can solve phase transition problems with complex interfaces. This paper comprehensively introduces the basic principle and numerical methods of the enthalpy-porosity method to deal with solidliquid interfaces. By giving several examples of experimental and calculated results, the effectiveness of the method is evaluated, and its application in the field of crystal growth is introduced. It is shown that the enthalpy-porosity method has the characteristics of simplicity and reliability, providing an effective numerical calculation method for the studies of solid-liquid phase transitions in many fields.

Key words: solid-liquid interface, enthalpy-porosity method, crystal growth, numerical simulation

CLC Number:

O414.13

WANG Ruixian , , YANG Mingliang , , HE Yi , , SUN Guihua , , ZHANG Qingli , . Application of the enthalpy‑porosity method in studies on solid‑liquid phase change interfaces[J]. Chinese Journal of Quantum Electronics, 2026, 43(3): 327-336.

References

[1] Stefan J.Uber einige probleme der theorie der warmeleitung[J].Sber Akad Wiss Wien, 1889, 98(0):473-484
[2]Sharma A, Tyagi V V, Chen C R, et al.Review on thermal energy storage with phase change materials and applications[J].Renewable and Sustainable Energy Reviews, 2009, 13(2):318-345
[3]Martínez-Tomás C, Mu?oz V.CdTe crystal growth process by the Bridgman method: numerical simulation[J].Journal of Crystal Growth, 2001, 222(3):435-451
[4] Sun X Q, Zhang Q, Liu Y J, et al.Numerical calculation and simulation on melting and solidification time periods of a phase change material[J]. Advanced Materials Research, 2012, 455-456(0):374-381
[5] Lu T, Jiang P X.Numerical simulation of natural convection for melting phase change in porous media[J].Journal of Engineering Thermophysics, 2005, 26(0):167-170
[6] Zhang B B, Xing Y M, Sheng Q.Numerical simulation of two phase change materials melting process[J].Advanced Materials Research, 2014, 1049-1050(0):94-100
[7] Mao Q J, Zhang Y M.Thermal energy storage performance of a three-PCM cascade tank in a high-temperature packed bed system[J].Renewable Energy, 2020, 152(0):110-119
[8]Jaluria Y.Fluid flow phenomena in materials processing-The 2000 freeman scholar lecture[J].Journal of Fluids Engineering, 2001, 123(2):173-210
[9]Dutil Y, Rousse D R, Salah N B, et al.A review on phase-change materials: Mathematical modeling and simulations[J].Renewable and Sustainable Energy Reviews, 2011, 15(1):112-130
[10]Gupta S C.A moving grid numerical scheme for multi-dimensional solidification with transition temperature range[J].Computer Methods in Applied Mechanics and Engineering, 2000, 189(2):525-544
[11]Tenchev R T, Mackenzie J A, Scanlon T J, et al.Finite element moving mesh analysis of phase change problems with natural convection[J].International Journal of Heat and Fluid Flow, 2005, 26(4):597-612
[12] Voller V R, Cross M, Markatos N C.An enthalpy method for convection/diffusion phase change[J].International Journal for Numerical Methods in Engineering, 197, 24(1):271-284
[13]Voller V R, Prakash C.A fixed grid numerical modeling methodology for convection-diffusion mushy region phase-change problem[J].ASME International Journal of Heat and Mass Transfer, 1987, 30(8):1709-1719
[14] Khattari Y, Rhafiki T E, Choab N, et al.Apparent heat capacity method to investigate heat transfer in a composite phase change material[J].Journal of Energy Storage, 2020, 28(0):101239-101239
[15] Yao W A, Yao H X, Zuo C.Precise integration boundary element method for solving dual phase change problems based on the effective heat capacity model[J].Engineering Analysis with Boundary Elements, 2019, 108(0):411-421
[16] Ma Z W, Zhang P.Modeling the heat transfer characteristics of flow melting of phase change material slurries in the circular tubes[J].International Journal of Heat and Mass Transfer, 2013, 64(0):874-881
[17]Pham Q T.Modelling heat and mass transfer in frozen foods: a review[J].International Journal of Refrigeration, 2006, 29(6):876-888
[18]Yang B, Fu X R, Yang W, et al.Effect of melt and mold temperatures on the solidification behavior of HDPE during gas-assisted injection molding: an enthalpy transformation approach[J].Macromolecular Materials and Engineering, 2009, 294(5):336-344
[19]Liang S P, Yang B, Fu X R, et al.A simple method for forecast of cooling time of high-density polyethylene during gas-assisted injection molding[J].Journal of Applied Polymer Science, 2010, 117(2):729-735
[20] Danaila I, Moglan R, Hecht F, et al.A Newton method with adaptive finite elements for solving phase-change problems with natural convection[J].Journal of Computational Physics, 2014, 274(0):826-840
[21]Salzmann C, Radaelli P G, Hallbrucker A, et al.The preparation and structures of hydrogen ordered phases of ice[J].Science, 2006, 311(5768):1758-1761
[22]Jana S, Ray S, Durst F.A numerical method to compute solidification and melting processes[J].Applied Mathematical Modelling, 2007, 31(1):93-119
[23]Udaykumar H S, Mittal R, Shyy W.Computation of solid–liquid phase fronts in the sharp interface limit on fixed grids[J].Journal of Computational Physics, 1999, 153(2):535-574
[24] Furmański P, Banaszek J.Modelling of the mushy zone permeability for solidifi-cation of binary alloys[J].Materials Science Forum, 2006, 508(0):411-418
[25] ANSYS Inc.ANSYS Fluent 14.0 User' s Guide[J].ANSYS Fluent 14.0 User' s Guide, 2009, 0(0):0-0
[26] Brent A D, Voller V R, Reid K J.Enthalpy-porosity technique for modeling convection-diffusion phase change: Application to the melting of a pure metal[J].Numerical Heat Transfer, 1988, 13(0):297-318
[27]刘夷平, 黄为民, 王经.利用焓-多孔介质法对垂直Bridgman生长CdTe的数值模拟[J].材料研究学报, 2006, 20(3):225-230
[28]何宜谦, 杨海天.基于Sigmoid函数光滑化的等效热容和有限元法求解相变传热问题[J].应用基础与工程科学学报, 2011, 19(5):817-829
[29]许波, 王恪典, 李雪芝, 等.焊接数值计算中几种相变潜热处理方法[J].热加工工艺, 2020, 49(3):7-10+6
[30]潘艾刚, 王俊彪, 张贤杰.基于等效热容法和焓法的相变传热数值分析[J].计算机仿真, 2014, 31(2):315-319
[31]张尚中, 曾忠, 张永祥, 等.Czochralski法晶体生长全局数值模拟[J].重庆交通大学学报自然科学版, 2009, 28(S1):355-357+365
[32] Younsi Z, Naji H.A numerical investigation of melting phase change process via the enthalpyporosity approach: Application to hydrated salts[J].International Communications in Heat and Mass Transfer, 2017, 86(0):12-24
[33] Iten M, Liu S, Shukla A.Experimental Validation of an Air-PCM Storage Unit Comparing the Effective Heat Capacity and Enthalpy Methods through CFD Simulations[J].Energy, 2018, 155(0):495-503
[34] Patankar S V.Numerical Heat Transfer and Fluid Flow[M]. New York: McGraw-Hill, 1980, p.285.
[35] Kulkarni S, Lin C X, Ghenai C.Numerical modeling of solidification process using the enthalpy-porosity technique[C]. ASME Summer Heat Transfer Conference, San Francisco, California, USA, 2005: 401-407.
[36] Duggirala R K, Lin C X, Ghenai C.Experimental Investigation of Fluid Flow and Phase moving interface during the solidification of binary mixture in trapezoidal cavity[C]. ASME Heat Transfer/Fluid Engineering conference, Charlotte, North Carolina, 2004: 11-15.
[37]Saleel C A.Numerical study on melting and heat transfer characteristics of paraffin wax SiC paraffin using enthalpy-porosity model[J].Journal of Thermal Analysis and Calorimetry, 2022, 147(19):10497-10508
[38] Shokouhmand H, Kamkari B.Experimental investigation on melting heat transfer characteristics of lauric acid in a rectangular thermal storage unit[J].Experimental Thermal and Fluid Science, 2013, 50(0):201-212
[39]刘永才, Simon B.石榴石晶体生长过程的生长动力学模型[J].人工晶体学报, 1999, 28(1):8-16
[40] Ding J L, Liu L J.Real-time prediction of crystal/melt interface shape during Czochralski crystal growth[J].CrystEngComm, 2018, 20(0):6925-6931
[41] Carrhà S, Fogliani S, Masi M, et al.Melt-solid interface shape in LEC GaAs crystals: comparison between calculated and experimentally observed shapes[J].Journal of Crystal Growth, 199, 166(1-4):641-645
[42] Wu X B, Geng X, Guo Z Y.Fundamental study of crystal/melt interface shape change in Czochralski crystal growth[J].Journal of Crystal Growth, 1996, 169(4):786-794
[43] Shiraishi Y, Maeda S, Nakamura K.Prediction of solid-liquid interface shape during CZ Si crystal growth using experimental and global simulation[J].Journal of Crystal Growth, 2004, 266(1-3):28-33
[44] Zhu Y Z, Ma D C, Long S W, et al.In-situ detection of growth striations by crystallization electromotive force measurement during Czochralski crystal growth[J].Journal of Crystal Growth, 2017, 475(0):70-76
[45]李有荣, 阮登芳, 彭岚, 等.硅单晶CzochralskiI法生长全局数值模拟传热与流动特性[J].材料研究学报, 2004, 18(2):102-108