流体力学与飞行力学

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

  • 齐海峰 ,
  • 常士楠 ,
  • 任战鹏 ,
  • 杨映麟
展开
  • 1.北京航空航天大学 航空科学与工程学院,北京 100191
    2.中国飞机强度研究所 强度与结构完整性全国重点实验室,西安 710065
.E-mail: sn_chang@buaa.edu.cn

收稿日期: 2025-05-13

  修回日期: 2025-06-08

  录用日期: 2025-07-11

  网络出版日期: 2025-07-15

基金资助

国家自然科学基金(12172031);强度与结构完整性国家重点实验室开放基金(LSSIKFJJ202405002)

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

  • Haifeng QI ,
  • Shinan CHANG ,
  • Zhanpeng REN ,
  • Yinglin YANG
Expand
  • 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 date: 2025-05-13

  Revised date: 2025-06-08

  Accepted date: 2025-07-11

  Online published: 2025-07-15

Supported by

National Natural Science Foundation of China(12172031);National Key Laboratory of Strength and Structural Integrity Open Fund(LSSIKFJJ202405002)

摘要

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

本文引用格式

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

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.

参考文献

[1] 钟富豪, 刘森云, 刘秀芳, 等. 机翼混合相冰晶结冰现象的数值研究[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).
[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] 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]. 航空发动机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] 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 Engineering2024238: 122044.
[6] 宋建宇, 张森, 涂杰, 等. 航空发动机冰晶结冰适航审定要点及符合性方法研究[J]. 燃气涡轮试验与研究202235(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 Research202235(3): 8-15 (in Chinese).
[7] 马乙楗, 柴得林, 易贤, 等. 考虑侵蚀效应的冰晶/混合相结冰计算方法[J]. 航空学报202344(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 Sinica202344(15): 528609 (in Chinese).
[8] 魏震, 刘秀芳, 钟富豪, 等. 微小冰晶粒子融化特性可视化实验[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(S2): 729301 (in Chinese).
[9] 丁军亮, 赵利利, 杨涛, 等. 自然结冰飞行试验技术综述[J]. 航空学报202344(17): 028270.
  DING J L, ZHAO L L, YANG T, et al. Flight test technology of natural icing[J]. Acta Aeronautica et Astronautica Sinica202344(17): 028270 (in Chinese).
[10] 胡应交, 徐峰, 杨志军. 航空发动机整机防结冰试验能力综述[J]. 航空学报202344(): 729449.
  HU Y J, XU F, YANG Z J. Overview of aero-engine ice testing capability[J]. Acta Aeronautica et Astronautica Sinica202344(S2): 729449 (in Chinese).
[11] 李斯, 束珺, 张志强, 等. 冰风洞过冷大水滴云雾水滴质量分布模拟[J]. 南京航空航天大学学报202355(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 & Astronautics202355(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 Science2024159: 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 Aeronautics202336(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 Engineering2024236: 121539.
[15] 蔡英磊, 董威. 冰风洞试验中水滴的传热传质计算研究[J]. 航空发动机201339(3): 36-40.
  CAI Y L, DONG W. Calculation and analysis of heat and mass transfer for water droplet in icing tunnel test[J]. Aeroengine201339(3): 36-40 (in Chinese).
[16] 柯鹏, 杨春信, 吴江浩, 等. 结冰实验段空气-水混合流动过程的数值模拟[J]. 航空动力学报200924(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 Power200924(7): 1449-1456 (in Chinese).
[17] 易贤, 马洪林, 王开春, 等. 结冰风洞液滴运动及传质传热特性分析[J]. 四川大学学报(工程科学版)201244(): 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 Sciences201244(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 Aircraft197512(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]. 计算机仿真201835(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 Simulation201835(10): 66-71 (in Chinese).
[21] 何鹏, 田骏丹, 田震, 等. 航空发动机试车台空气雾化加湿系统设计与研究[J]. 科学技术与工程202424(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 Engineering202424(11): 4495-4500 (in Chinese).
[22] 田小江, 孙玥, 吴锋, 等. 航空发动机结冰试验中水滴特性参数分析[J]. 燃气涡轮试验与研究202235(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 Research202235(1): 41-44, 57 (in Chinese).
[23] 孙澄川, 卢静, 陈东, 等. 悬浮液冷喷涂的液滴蒸发及加速过程研究[J]. 材料研究与应用202317(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 Application202317(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, 20131547(1): 659-670.
[25] WILLIS K D, ORME M E. Experiments on the dynamics of droplet collisions in a vacuum[J]. Experiments in Fluids200029(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 Transfer2020155: 119720.
[27] SAZHIN S S. Advanced models of fuel droplet heating and evaporation[J]. Progress in Energy and Combustion Science200632(2): 162-214.
[28] 潘艳秋, 杜宇杰, 俞路. 小型喷雾冷却装置两相区传热特性模拟及优化[J]. 天津大学学报(自然科学与工程技术版)202356(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)202356(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.
文章导航

/