基于沿程冷却法的冰晶沿程参数变化模拟

doi:10.7527/S1000-6893.2025.32227

流体力学与飞行力学

本期目录 | 过刊浏览 | 高级检索

基于沿程冷却法的冰晶沿程参数变化模拟

齐海峰¹, 常士楠¹(), 任战鹏², 杨映麟¹

^1.北京航空航天大学航空科学与工程学院，北京 100191
^2.中国飞机强度研究所强度与结构完整性全国重点实验室，西安 710065

收稿日期:2025-05-13 修回日期:2025-06-08 接受日期:2025-07-11 出版日期:2025-07-16 发布日期:2025-07-15
通讯作者: 常士楠 E-mail:sn_chang@buaa.edu.cn
基金资助:
国家自然科学基金(12172031);强度与结构完整性国家重点实验室开放基金(LSSIKFJJ202405002)

Simulation of ice crystal parameters along path change based on path cooling method

Haifeng QI¹, Shinan CHANG¹(), Zhanpeng REN², Yinglin YANG¹

^1.School of Aeronautic Science and Engineering，Beijing University，Beijing 100191，China
^2.National Key Laboratory of Strength and Structural Integrity，Aircraft Strength Research Institute of China，Xi’an 710065，China

Received:2025-05-13 Revised:2025-06-08 Accepted:2025-07-11 Online:2025-07-16 Published:2025-07-15
Contact: Shinan CHANG E-mail:sn_chang@buaa.edu.cn
Supported by:
National Natural Science Foundation of China(12172031);National Key Laboratory of Strength and Structural Integrity Open Fund(LSSIKFJJ202405002)

摘要/Abstract

摘要：

为了解决冰晶生成的技术难题，采用沿程冷却法对水滴相变成冰晶的过程进行了数值模拟，以确定冰晶生成的具体位置、范围，即进行冰晶结冰模拟试验的试验段区间。冰晶在随气流运动过程中会与周围空气发生热、质传递，其直径、温度和速度等参数变化显著。基于沿程冷却的方法，加入相变判断准则，研究了不同条件下冰晶沿程参数变化的规律。研究结果表明，冰风洞稳定段进口气流速度增加时，试验段内部的风速大幅增加。进口气流速度为10、15、20 m/s的条件下，试验段内部的结冰比例相对较高，均达到了90%以上。在进口气流温度为233.15 K和223.15 K时，水滴在x=2.3 m附近已完全冻结成冰晶。随着水滴直径的增加，试验段内部沿程不同截面上的总水含量逐渐降低，水滴开始结冰的位置沿运动方向向后推移。

关键词: 水滴冻结, 沿程冷却法, 冰风洞, 传热传质, 冰晶生成

Abstract:

In order to solve the technical problems of ice crystal formation， the process of water droplet phase turning into ice crystals was numerically simulated by the path cooling method to determine the specific location and range of ice crystal formation， that is， the test section range for conducting the ice crystal icing simulation test. In the process of moving with the air flow， the ice crystal will have heat/mass transfer with the surrounding air， and its parameters such as diameter， temperature and velocity change significantly. Based on the path cooling method and phase transition criteria， the variation of ice crystal parameters under different conditions is studied. The results show that when the inlet wind speed of the stable section increases， the wind speed inside the test section increases significantly. Under the conditions of inlet wind speed of 10， 15， 20 m/s， the icing proportion inside the test section is relatively high， reaching more than 90%. At ambient temperatures of 233.15 K and 223.15 K， the water droplets have completely frozen into ice crystals near x = 2.3 m. Furthermore， with the increase of the diameter of the water drop， the total water content in different sections along the test section gradually decreases， and the position where the water drops begins to freeze increases along the direction of movement.

Key words: water droplet icing, path cooling method, ice wind tunnel, heat and mass transfer, ice crystal formation

中图分类号:

V211.3

齐海峰, 常士楠, 任战鹏, 杨映麟. 基于沿程冷却法的冰晶沿程参数变化模拟[J]. 航空学报, 2026, 47(4): 132227.

Haifeng QI, Shinan CHANG, Zhanpeng REN, Yinglin YANG. Simulation of ice crystal parameters along path change based on path cooling method[J]. Acta Aeronautica et Astronautica Sinica, 2026, 47(4): 132227.

图/表 16

图 1

图 2

图 3

图 4

表 1

冰风洞出口参数的结果对比

参数	类别	数值
$u o u t$ /（m·s^-1）	试验值	35.8
	模拟值	35.2
	相对误差	-1.676%
RH $o u t$ /%	试验值	57.00
	模拟值	55.91
	相对误差	-1.912%
LWC/（g·m^-3）	试验值	0.130
	模拟值	0.128
	相对误差	-1.539%
MVD/μm	试验值	22.80
	模拟值	21.95
	相对误差	-3.728%
T_out/K	试验值	271.55
	模拟值	272.23
	相对误差	0.250%

表 1

图 5

图 6

图 7

图 8

图 9

图 10

图 11

图 12

图 13

图 14

图 15

参考文献 29

[1]	钟富豪，刘森云，刘秀芳，等. 机翼混合相冰晶结冰现象的数值研究［J］. 西安交通大学学报， 2024， 58（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 University， 2024， 58（10）： 168-177 （in Chinese）.
[2]	黄平，卜雪琴，刘一鸣，等. 混合相/冰晶条件下的结冰研究综述［J］. 航空学报， 2022， 43（5）： 025178.
	HUANG P， BU X Q， LIU Y M， et al. Mixed phase/glaciated ice accretion： Review［J］. Acta Aeronautica et Astronautica Sinica， 2022， 43（5）： 025178 （in Chinese）.
[3]	MASON J， STRAPP W， CHOW P. The ice particle threat to engines in flight［C］∥ 44th AIAA Aerospace Sciences Meeting and Exhibit. Reston： AIAA， 2006.
[4]	苏龙伟，申世才，田晓平，等. 航空发动机冰晶结冰研究进展［J］. 航空发动机， 2023， 49（5）： 89-99.
	SU L W， SHEN S C， TIAN X P， et al. Research progress on ice crystal icing in aeroengine［J］. Aeroengine， 2023， 49（5）： 89-99 （in Chinese）.
[5]	QI H F， CHANG S N， YANG Y L， et al. Numerical investigations on aero-engine icing characteristics at mixed phase conditions［J］. Applied Thermal Engineering， 2024， 238： 122044.
[6]	宋建宇，张森，涂杰，等. 航空发动机冰晶结冰适航审定要点及符合性方法研究［J］. 燃气涡轮试验与研究， 2022， 35（3）： 8-15.
	SONG J Y， ZHANG S， TU J， et al. Study on airworthiness certification key points and compliance method of aero-engine ice crystal icing［J］. Gas Turbine Experiment and Research， 2022， 35（3）： 8-15 （in Chinese）.
[7]	马乙楗，柴得林，易贤，等. 考虑侵蚀效应的冰晶/混合相结冰计算方法［J］. 航空学报， 2023， 44（15）： 528609.
	MA Y J， CHAI D L， YI X， et al. Calculation method for ice crystal/mixed phase icing considering ice crystal erosion［J］. Acta Aeronautica et Astronautica Sinica， 2023， 44（15）： 528609 （in Chinese）.
[8]	魏震，刘秀芳，钟富豪，等. 微小冰晶粒子融化特性可视化实验［J］. 航空学报， 2023， 44（）： 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 Sinica， 2023， 44（S2）： 729301 （in Chinese）.
[9]	丁军亮，赵利利，杨涛，等. 自然结冰飞行试验技术综述［J］. 航空学报， 2023， 44（17）： 028270.
	DING J L， ZHAO L L， YANG T， et al. Flight test technology of natural icing［J］. Acta Aeronautica et Astronautica Sinica， 2023， 44（17）： 028270 （in Chinese）.
[10]	胡应交，徐峰，杨志军. 航空发动机整机防结冰试验能力综述［J］. 航空学报， 2023， 44（）： 729449.
	HU Y J， XU F， YANG Z J. Overview of aero-engine ice testing capability［J］. Acta Aeronautica et Astronautica Sinica， 2023， 44（S2）： 729449 （in Chinese）.
[11]	李斯，束珺，张志强，等. 冰风洞过冷大水滴云雾水滴质量分布模拟［J］. 南京航空航天大学学报， 2023， 55（1）： 146-153.
	LI S， SHU J， ZHANG Z Q， et al. Study on simulation for droplet mass distribution of supercooled large droplet cloud in icing wind tunnel［J］. Journal of Nanjing University of Aeronautics & Astronautics， 2023， 55（1）： 146-153 （in Chinese）.
[12]	WANG B Y， ZHANG L T， HUANG X Y， et al. Simultaneous measurements of morphology， size distribution and concentration of ice crystals in icing wind tunnel with HACPI［J］. Experimental Thermal and Fluid Science， 2024， 159： 111266.
[13]	WANG Q， CHEN N L， WANG Y B， et al. Scallop ice shape characteristics of swept wing based on large-scale icing wind tunnel experiment［J］. Chinese Journal of Aeronautics， 2023， 36（12）： 214-230.
[14]	MA L S， SONG H， CHEN Y， et al. A numerical simulation of the distribution and the variation law of the liquid water content in icing wind tunnel［J］. Applied Thermal Engineering， 2024， 236： 121539.
[15]	蔡英磊，董威. 冰风洞试验中水滴的传热传质计算研究［J］. 航空发动机， 2013， 39（3）： 36-40.
	CAI Y L， DONG W. Calculation and analysis of heat and mass transfer for water droplet in icing tunnel test［J］. Aeroengine， 2013， 39（3）： 36-40 （in Chinese）.
[16]	柯鹏，杨春信，吴江浩，等. 结冰实验段空气-水混合流动过程的数值模拟［J］. 航空动力学报， 2009， 24（7）： 1449-1456.
	KE P， YANG C X， WU J H， et al. Numerical simulation of mixing flow of air and water droplet in the test duct of icing research tunnel［J］. Journal of Aerospace Power， 2009， 24（7）： 1449-1456 （in Chinese）.
[17]	易贤，马洪林，王开春，等. 结冰风洞液滴运动及传质传热特性分析［J］. 四川大学学报（工程科学版）， 2012， 44（）： 132-135.
	YI X， MA H L， WANG K C， et al. Analysis of water droplet movement and heat/mass transfer in an icing wind tunnel［J］. Advanced Engineering Sciences， 2012， 44（S2）： 132-135 （in Chinese）.
[18]	WILLBANKS C E， SCHULZ R J. Analytical study of icing simulation for turbine engines in altitude test cells［J］. Journal of Aircraft， 1975， 12（12）： 960-967.
[19]	SCHULZ R. A model for predicting mixed-phase flow in ground test facilities［C］∥37th Aerospace Sciences Meeting and Exhibit. Reston： AIAA， 1999.
[20]	李晓峰，常士楠，冷梦尧. 冰风洞内水滴参数沿程变化的数值模拟［J］. 计算机仿真， 2018， 35（10）： 66-71.
	LI X F， CHANG S N， LENG M Y. Numerical simulation of water droplet parameters variation along the icing wind tunnel［J］. Computer Simulation， 2018， 35（10）： 66-71 （in Chinese）.
[21]	何鹏，田骏丹，田震，等. 航空发动机试车台空气雾化加湿系统设计与研究［J］. 科学技术与工程， 2024， 24（11）： 4495-4500.
	HE P， TIAN J D， TIAN Z， et al. Design and research of air atomization humidification system for an aero-engine test bench［J］. Science Technology and Engineering， 2024， 24（11）： 4495-4500 （in Chinese）.
[22]	田小江，孙玥，吴锋，等. 航空发动机结冰试验中水滴特性参数分析［J］. 燃气涡轮试验与研究， 2022， 35（1）： 41-44， 57.
	TIAN X J， SUN Y， WU F， et al. Characteristics of droplet in the icing test of aero-engine［J］. Gas Turbine Experiment and Research， 2022， 35（1）： 41-44， 57 （in Chinese）.
[23]	孙澄川，卢静，陈东，等. 悬浮液冷喷涂的液滴蒸发及加速过程研究［J］. 材料研究与应用， 2023， 17（1）： 87-101.
	SUN C C， LU J， CHEN D， et al. Investigation on evaporation and acceleration process of droplets during suspension cold spraying［J］. Materials Research and Application， 2023， 17（1）： 87-101 （in Chinese）.
[24]	LENG M Y， CHANG S N， WU M L， et al. Simulation of air-droplet mixed phase flow in icing wind-tunnel ［C］∥AIP Conference Proceedings， 2013， 1547（1）： 659-670.
[25]	WILLIS K D， ORME M E. Experiments on the dynamics of droplet collisions in a vacuum［J］. Experiments in Fluids， 2000， 29（4）： 347-358.
[26]	DENG H Y， CHANG S N， SONG M J. The optimization of simulated icing environment by adjusting the arrangement of nozzles in an atomization equipment for the anti-icing and deicing of aircrafts［J］. International Journal of Heat and Mass Transfer， 2020， 155： 119720.
[27]	SAZHIN S S. Advanced models of fuel droplet heating and evaporation［J］. Progress in Energy and Combustion Science， 2006， 32（2）： 162-214.
[28]	潘艳秋，杜宇杰，俞路. 小型喷雾冷却装置两相区传热特性模拟及优化［J］. 天津大学学报（自然科学与工程技术版）， 2023， 56（6）： 564-572.
	PAN Y Q， DU Y J， YU L. Numerical simulation and optimization on two-phase heat transfer characteristics of small spray cooling device［J］. Journal of Tianjin University （Science and Technology）， 2023， 56（6）： 564-572 （in Chinese）.
[29]	DAVISON C R， MACLEOD J D， CHALMERS J L. Droplet evaporation model for determining liquid water content in engine icing tunnels and examination of the factors affecting liquid water content［C］∥9th AIAA Atmospheric and Space Environments Conference. Reston： AIAA， 2017.

E-mail：hkxb@buaa.edu.cn

关于我们

期刊社服务

专业学科

封面文章

友情链接

主管单位：中国科学技术协会主办单位：中国航空学会北京航空航天大学

基于沿程冷却法的冰晶沿程参数变化模拟

Simulation of ice crystal parameters along path change based on path cooling method

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

图/表 16

参考文献 29

相关文章 15

编辑推荐

Metrics

本文评价

[1]	刘欣乐, 姜亚楠, 辛荣提, 李庆辉, 蔡晋生. 超疏水电热复合分区防冰策略[J]. 航空学报, 2025, 46(9): 130784-130784.
[2]	易贤, 任靖豪, 赖庆仁, 刘宇, 王强. 大型飞机全尺寸多段翼结冰特性计算和试验[J]. 航空学报, 2025, 46(5): 531575-531575.
[3]	杨秋明, 朱永峰, 陈华伟, 刘晓林. 镂空图案化电加热组件防冰性能试验研究[J]. 航空学报, 2023, 44(S2): 729334-729334.
[4]	许岭松, 吴渊, 朱东宇, 张付昆, 刘昱. FL⁃61结冰风洞翼型俯仰振荡机构研制[J]. 航空学报, 2023, 44(S2): 729296-729296.
[5]	容易, 刘辉, 于子文, 朱平平, 彭越, 于佳晖. 可重复使用运载火箭返回段低温流体行为特性[J]. 航空学报, 2023, 44(23): 629563-629563.
[6]	束珺, 徐东光, 韩志熔, 李斯, 黄雄. 结冰风洞过冷大水滴试验中混合翼设计[J]. 航空学报, 2023, 44(1): 627182-627182.
[7]	赵瑾, 孙向春, 张俊, 唐志共, 文东升. 热防护材料气固界面传热传质问题研究进展[J]. 航空学报, 2022, 43(10): 527577-527577.
[8]	郭向东, 柳庆林, 刘森云, 王梓旭, 李明. 结冰风洞中过冷大水滴云雾演化特性数值研究[J]. 航空学报, 2020, 41(8): 123655-123655.
[9]	郭向东, 张平涛, 赵照, 赖庆仁, 郭龙. 大型结冰风洞云雾场适航应用符合性验证[J]. 航空学报, 2020, 41(10): 123879-123879.
[10]	郭向东, 王梓旭, 李明, 刘蓓. 结冰风洞中液滴相变效应数值模拟[J]. 航空学报, 2018, 39(3): 121586-121586.
[11]	李伟斌, 魏东, 杜雁霞, 易贤, 杨肖峰. 动态结冰微观孔隙结构定量分析[J]. 航空学报, 2018, 39(2): 121536-121536.
[12]	雷桂林, 郑梅, 董威, 周志翔, 董奇. 航空发动机进气支板电热防冰试验[J]. 航空学报, 2017, 38(8): 121066-121066.
[13]	李岩, 王绍龙, 易贤, 周志宏, 郭龙. 绕轴旋转圆柱结冰特性结冰风洞试验[J]. 航空学报, 2017, 38(2): 520693-520703.
[14]	易贤, 王斌, 李伟斌, 郭龙. 飞机结冰冰形测量方法研究进展[J]. 航空学报, 2017, 38(2): 520700-520711.
[15]	郭向东, 王梓旭, 李明, 肖春华. 结冰风洞中液滴过冷特性数值研究[J]. 航空学报, 2017, 38(10): 121254-121254.