焓-多孔介质法在固液相变界面研究中的应用

doi:10.3969/j.issn.1007-5461.2026.03.001

量子电子学报 ›› 2026, Vol. 43 ›› Issue (3): 327-336.doi: 10.3969/j.issn.1007-5461.2026.03.001

• 综述 • 下一篇

焓-多孔介质法在固液相变界面研究中的应用

王瑞仙 1,2,3, 杨明亮 1,2,3, 何異 1,2,3, 孙贵华 2,3, 张庆礼 2,3*

1 中国科学技术大学研究生院科学岛分院, 安徽合肥 230031; 2 中国科学院合肥物质科学研究院安徽光学精密机械研究所, 光子器件与材料安徽省重点实验室, 安徽合肥 230031; 3 先进激光技术安徽省实验室, 安徽合肥 230037

收稿日期:2023-08-31 修回日期:2023-10-12 出版日期:2026-05-28 发布日期:2026-05-28
通讯作者: E-mail: zql@aiofm.ac.cn E-mail:zql@aiofm.ac.cn
作者简介:王瑞仙 ( 1997 - ), 女, 云南曲靖人, 博士生, 主要从事晶体生长数值模拟方面的研究。E-mail: rxwang2021@163.com
基金资助:
国家重点研发计划 (2022YFB3605700), 国家自然科学基金 (52272011, 11875248), 安徽省实验室重点基金 (AHL20220ZR04)

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

WANG Ruixian 1,2,3 , YANG Mingliang 1,2,3 , HE Yi 1,2,3 , SUN Guihua 2,3 , ZHANG Qingli 2,3*

1 Science Island Branch of Graduate School, University of Science and Technology of China, Hefei 230031, China; 2 Anhui Provincial Key Laboratory of Photonic Devices and Materials, Anhui Institute of Optics and Fine Mechanics, HFIPS, Chinese Academy of Sciences, Hefei 230031, China; 3 Advanced Laser Technology Laboratory of Anhui Province, Hefei 230037, China

Received:2023-08-31 Revised:2023-10-12 Published:2026-05-28 Online:2026-05-28

摘要/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

中图分类号:

O414.13

王瑞仙, 杨明亮, 何異, 孙贵华, 张庆礼, . 焓-多孔介质法在固液相变界面研究中的应用[J]. 量子电子学报, 2026, 43(3): 327-336.

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.

参考文献

[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

焓-多孔介质法在固液相变界面研究中的应用

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

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 9

编辑推荐

Metrics

本文评价

[1]	张瑞, 梅大江, ∗, 石小兔, 马荣国, 张庆礼, ∗, 窦仁勤, 刘文鹏, . YAG 晶体的位错研究进展[J]. 量子电子学报, 2022, 39(5): 687-706.
[2]	江俊勤∗, 涂菁. 多势阱中电子的束缚态与能带形成研究[J]. 量子电子学报, 2021, 38(3): 316-326.
[3]	徐民, 龚巧瑞, 李善明, 杭寅∗. 钛宝石激光晶体研究进展[J]. 量子电子学报, 2021, 38(2): 148-159.
[4]	韩卫民, 倪友保∗, 吴海信, 王振友, 黄昌保. 新型长波红外材料PbGa6Te10 的生长[J]. 量子电子学报, 2021, 38(2): 172-179.
[5]	李壮, 李春霄, 姚吉勇, ∗, 吴以成, . BaGa4Se7 和BaGa2GeSe6 晶体研究进展[J]. 量子电子学报, 2021, 38(2): 185-191.
[6]	张中晗, 戴云, 王阳啸, 张振, 武安华, 苏良碧, ∗. 单晶光纤的生长技术与应用研究[J]. 量子电子学报, 2021, 38(2): 192-213.
[7]	文林芳方建兴. 塞曼能级电磁诱导透明暗态的演化[J]. J4, 2015, 32(6): 706-712.
[8]	曹文华徐平刘颂豪. 基于分布增益非线性光纤环镜的弱脉冲放大与压缩[J]. J4, 2012, 29(1): 80-88.
[9]	王迪万松明张庆礼孙敦陆殷绍唐尤静林王媛媛. LiB₃O₅晶体及其自助溶生长溶液的高温拉曼光谱研究[J]. J4, 2011, 28(2): 241-246.