Special Topic: Anti-icing and De-icing Technology for Aeroengines

Numerical simulation of melting and shape evolution characteristics of non-spherical ice crystals

  • Fuhao ZHONG ,
  • Xiufang LIU ,
  • Jiajun CHEN ,
  • Bo HAN ,
  • Zhou FANG ,
  • Yu HOU
Expand
  • 1.School of Energy and Power Engineering,Xi’an Jiaotong University,Xi’an 710049,China
    2.Key Laboratory of Cryogenic Technology and Equipment of Ministry of Education,Xi’an Jiaotong University,Xi’an 710049,China

Received date: 2025-10-17

  Revised date: 2025-12-02

  Accepted date: 2025-12-25

  Online published: 2026-01-09

Supported by

National Natural Science Foundation of China(52476018);National Science and Technology Major Project of China (J2019-Ⅲ-0010-0054)

Abstract

To investigate the melting characteristics and shape evolution of non-spherical ice crystals in hot airflow environment, a melting model for non-spherical ice crystals was developed based on a decoupling strategy for flow and phase-change heat transfer, enabling a solution to the complex heat transfer problem involving gas-liquid-solid coupling. The results demonstrate that the proposed model accurately predicts the shape evolution process during ice crystal melting, and the predicted melting time agrees well with experimental results, with a deviation of ±15%. The shape evolution of non-spherical ice crystals undergoes three distinct stages: an initial warming stage where the shape remains stable, a partial water coverage stage where liquid water partially covers the ice core surface, and a complete water coverage stage where the ice core is fully enveloped. The initial aspect ratio is identified as the key factor governing the shape evolution during melting. A larger initial aspect ratio results in a longer time required for the ice crystal to evolve into a spherical shape. Furthermore, an empirical correlation between dimensionless sphericity and the melting ratio was established, overcoming the limitations of the traditional linear approximation model for high-aspect-ratio (initial aspect ratio λ > 2) ice crystals and significantly improving computational accuracy.

Cite this article

Fuhao ZHONG , Xiufang LIU , Jiajun CHEN , Bo HAN , Zhou FANG , Yu HOU . Numerical simulation of melting and shape evolution characteristics of non-spherical ice crystals[J]. ACTA AERONAUTICAET ASTRONAUTICA SINICA, 2026 , 47(11) : 632933 -632933 . DOI: 10.7527/S1000-6893.2025.32933

References

[1] TETTEH E, LOTH E, NEUTEBOOM M O, et al. In-flight gas turbine engine icing: Review[J]. AIAA Journal202260(10): 5610-5632.
[2] 黄平, 卜雪琴, 刘一鸣, 等. 混合相/冰晶条件下的结冰研究综述[J]. 航空学报202243(5): 025178.
  HUANG P, BU X Q, LIU Y M, et al. Mixed phase/glaciated ice accretion: Review[J]. Acta Aeronautica et Astronautica Sinica202243(5): 025178 (in Chinese).
[3] CAO Y H, TAN W Y, WU Z L. Aircraft icing: An ongoing threat to aviation safety[J]. Aerospace Science and Technology201875: 353-385.
[4] 苏龙伟, 申世才, 田晓平, 等. 航空发动机冰晶结冰研究进展[J]. 航空发动机202349(5): 89-99.
  SU L W, SHEN S C, TIAN X P, et al. Research progress on ice crystal icing in aeroengine[J]. Aeroengine202349(5): 89-99 (in Chinese).
[5] 卜雪琴, 黄平, 林贵平, 等. 多步法冰晶积冰数值模拟[J]. 航空学报202445(14): 129308.
  BU X Q, HUANG P, LIN G P, et al. Multi-step method for numerical simulation of ice crystal icing[J]. Acta Aeronautica et Astronautica Sinica202445(14): 129308 (in Chinese).
[6] 袁庆浩, 樊江, 白广忱. 航空发动机内部冰晶结冰研究综述[J]. 推进技术201839(12): 2641-2650.
  YUAN Q H, FAN J, BAI G C. Review of ice crystal icing in aero-engines[J]. Journal of Propulsion Technology201839(12): 2641-2650 (in Chinese).
[7] 刘宗辉, 卜雪琴, 包佳仪, 等. 部分融化冰晶撞击表面临界破碎速度实验[J]. 航空学报202546(14): 131558.
  LIU Z H, BU X Q, BAO J Y, et al. Experiment on critical fragmentation velocity of partially melted ice particles impacting a rigid surface[J]. Acta Aeronautica et Astronautica Sinica202546(14): 131558 (in Chinese).
[8] 齐海峰, 常士楠, 任战鹏, 等. 基于沿程冷却法的冰晶沿程参数变化模拟[J]. 航空学报202647(4): 132227.
  QI H F, CHANG S N, REN Z P, et al. Simulation of ice crystal parameters along path change based on path cooling method [J]. Acta Aeronautica et Astronautica Sinica202647(4): 132227 (in Chinese).
[9] 黄平, 卜雪琴, 林贵平, 等. 冰晶粒子运动过程中的相变特性[J]. 航空动力学报202237(7): 1379-1391.
  HUANG P, BU X Q, LIN G P, et al. Phase transition characteristics of ice crystal particles in motion[J]. Journal of Aerospace Power202237(7): 1379-1391 (in Chinese).
[10] 马乙楗, 柴得林, 王强, 等. 翼面结冰过程中的冰晶运动相变与黏附特性[J]. 航空学报202344(1): 627817.
  MA Y J, CHAI D L, WANG Q, et al. Phase change and adhesion characteristics of ice crystal movements in wing icing[J]. Acta Aeronautica et Astronautica Sinica202344(1): 627817 (in Chinese).
[11] XIAO H X, ZHANG F, HE Q S, et al. Classification of ice crystal habits observed from airborne cloud particle imager by deep transfer learning[J]. Earth and Space Science20196(10): 1877-1886.
[12] LAWSON R P, ANGUS L J, HEYMSFIELD A J. Cloud particle measurements in thunderstorm anvils and possible weather threat to aviation[J]. Journal of Aircraft199835(1): 113-121.
[13] BAILEY I H, MACKLIN W C. The surface configuration and internal structure of artificial hailstones[J]. Quarterly Journal of the Royal Meteorological Society196894(399): 1-11.
[14] YAN S H, PALACIOS J. A physics-based sphericity, drag coefficient and Nusselt number model of a partially-melted ice crystal[J]. International Journal of Thermal Sciences2024197: 108798.
[15] ZHONG F H, LIU X F, CHEN J J, et al. Experimental study on the melting characteristics of a suspended ice crystal based on new identification methods of ice-water interface and melting termination[J]. International Journal of Heat and Fluid Flow2024110: 109629.
[16] KINTEA D M, HAUK T, ROISMAN I V, et al. Shape evolution of a melting nonspherical particle[J]. Physical Review E, Statistical, Nonlinear, and Soft Matter Physics201592(3): 033012.
[17] 钟富豪, 刘秀芳, 刘森云, 等. 微小冰晶融化特性的实验与数值模拟研究[J]. 中南大学学报(自然科学版)202556(6): 2477-2486.
  ZHONG F H, LIU X F, LIU S Y, et al. Experimental and numerical investigation on melting characteristics of small ice crystals[J]. Journal of Central South University (Science and Technology)202556(6): 2477-2486 (in Chinese).
[18] ZHONG F H, WEI Z, CHEN J J, et al. Dynamic behavior and heat transfer characteristics of non-spherical ice crystals in high-temperature air flow[J]. Transactions of Nanjing University of Aeronautics and Astronautics202340(2): 205-215.
[19] YANG R, VAN DEN HAM T, VERZICCO R, et al. Circular objects do not melt the slowest in water[J]. Physical Review Fluids20249(8): 083501.
[20] YANG R, HOWLAND C J, LIU H R, et al. Shape effect on solid melting in flowing liquid[J]. Journal of Fluid Mechanics2024980: R1.
[21] YANG R, CHONG K L, LIU H R, et al. Abrupt transition from slow to fast melting of ice[J]. Physical Review Fluids20227(8): 083503.
[22] ZHU P F, ZHANG J C, HAN B B, et al. Three-dimensional numerical simulation of ice crystal melting in jet engine[J]. Journal of Thermal Science201928(5): 984-992.
[23] HAUK T, BONACCURSO E, VILLEDIEU P, et al. Theoretical and experimental investigation of the melting process of ice particles[J]. Journal of Thermophysics and Heat Transfer201630(4): 946-954.
[24] HAUK T, ROISMAN I V, TROPEA C D. Investigation of the melting behaviour of ice particles in an acoustic levitator[C]∥ 11th AIAA/ASME Joint Thermophysics and Heat Transfer Conference. Reston: AIAA, 2014.
[25] BAUMERT A, BANSMER S, TRONTIN P, et al. Experimental and numerical investigations on aircraft icing at mixed phase conditions[J]. International Journal of Heat and Mass Transfer2018123: 957-978.
[26] WRIGHT W. Further refinement of the LEWICE SLD model[C]∥ 44th AIAA Aerospace Sciences Meeting and Exhibit. Reston: AIAA, 2006.
[27] NILAMDEEN S, HABASHI W, AUBé M, et al. FENSAP-ICE: Modeling of water droplets and ice crystals[C]∥ 1st AIAA Atmospheric and Space Environments Conference. Reston: AIAA, 2009.
[28] WRIGHT W, JORGENSON P, VERES J. Mixed phase modeling in GlennICE with application to engine icing[C]∥ AIAA Atmospheric and Space Environments Conference. Reston: AIAA, 2010.
[29] QI H F, CHANG S N, YANG Y L. Numerical study of mixed phase ice accumulation in aero-engine inlet system[J]. Applied Thermal Engineering2023231: 120909.
[30] 张丽芬, 赵建辉, 余邦拓, 等. 旋转叶片通道内冰晶运动和融化数值研究[J]. 西北工业大学学报202442(6): 987-995.
  ZHANG L F, ZHAO J H, YU B T, et al. Numerical study of ice crystal movement and melting in rotating blade channels[J]. Journal of Northwestern Polytechnical University202442(6): 987-995 (in Chinese).
[31] 钟富豪, 刘森云, 刘秀芳, 等. 机翼混合相冰晶结冰现象的数值研究[J]. 西安交通大学学报202458(10): 168-177.
  ZHONG F H, LIU S Y, LIU X F, et al. Numerical study on mixed-phase ice crystal icing on airfoils[J]. Journal of Xi’an Jiaotong University202458(10): 168-177 (in Chinese).
[32] AGUILAR B, TRONTIN P, REITTER L, et al. Ice crystal drag model extension to snowflakes: Experimental and numerical investigations[J]. AIAA Journal202260(12): 6633-6646.
[33] K?BSCHALL K, TRAUT B, ROISMAN I V, et al. Melting of fractal snowflakes: Experiments and modeling[J]. International Journal of Heat and Mass Transfer2023212: 124254.
[34] TRONTIN P, BLANCHARD G, VILLEDIEU P. A comprehensive numerical model for mixed-phase and glaciated icing conditions[C]∥ 8th AIAA Atmospheric and Space Environments Conference. Reston: AIAA, 2016.
[35] 陶文铨. 传热学[M]. 5版. 北京: 高等教育出版社, 2019: 514-516.
  TAO W Q. Heat transfer[M]. 5th ed. Beijing: Higher Education Press, 2019: 514-516 (in Chinese).
[36] 魏震, 刘秀芳, 钟富豪, 等. 微小冰晶粒子融化特性可视化实验[J]. 航空学报202344(): 729301.
  WEI Z, LIU X F, ZHONG F H, et al. Visual experimental investigation on melting characteristics of minuscule ice crystal particles[J]. Acta Aeronautica et Astronautica Sinica202344(Sup2): 729301 (in Chinese).
[37] YANG X, MCGILVRAY M, GILLESPIE D R H. Modelling the particle trajectory and melting behaviour of non-spherical ice crystal particles[J]. International Journal of Multiphase Flow2022148: 103949.
[38] WHITHAM G B. Linear and nonlinear waves[M]. Hoboken: John Wiley & Sons, Inc., 1999: 19-42.
[39] OSHER S, FEDKIW R P. Level set methods: An overview and some recent results[J]. Journal of Computational Physics2001169(2): 463-502.
[40] DEL R??O O I, NEUMANN A W. Axisymmetric drop shape analysis: Computational methods for the measurement of interfacial properties from the shape and dimensions of pendant and sessile drops[J]. Journal of Colloid and Interface Science1997196(2): 136-147.
[41] BUTCHER J C. A history of Runge-Kutta methods[J]. Applied Numerical Mathematics199620(3): 247-260.
[42] HUANG L Y, LIU Z L, LIU Y M, et al. Effect of contact angle on water droplet freezing process on a cold flat surface[J]. Experimental Thermal and Fluid Science201240: 74-80.
[43] THIéVENAZ V, JOSSERAND C, SéON T. Retraction and freezing of a water film on ice[J]. Physical Review Fluids20205(4): 041601.
[44] YANG X, MCGILVRAY M, GILLESPIE D. Numerical modeling and parametric study of the melting behavior of ice crystal particles[J]. AIAA Journal202159(11): 4660-4668.
Outlines

/