结冰风洞中过冷大水滴云雾演化特性数值研究

doi:10.7527/S1000-6893.2019.23655

流体力学与飞行力学

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

前一篇 | 后一篇

结冰风洞中过冷大水滴云雾演化特性数值研究

郭向东^1,2, 柳庆林², 刘森云², 王梓旭², 李明²

1. 中国空气动力研究与发展中心空气动力学国家重点实验室, 绵阳 621000;
2. 中国空气动力研究与发展中心飞行器结冰与防/除冰重点实验室, 绵阳 621000

收稿日期:2019-11-13 修回日期:2019-12-27 出版日期:2020-08-15 发布日期:2020-01-02
通讯作者: 郭向东 E-mail:easterkuo@163.com
基金资助:
国家"973"计划（2015CB755800）

Numerical study of supercooled large droplet cloud evolution characteristics in icing wind tunnel

GUO Xiangdong^1,2, LIU Qinglin², LIU Senyun², WANG Zixu², LI Ming²

1. State Key Laboratory of Aerodynamics, China Aerodynamics Research and Development Center, Mianyang 621000, China;
2. Key Laboratory of Aircraft Icing and Anti/De-Icing, China Aerodynamics Research and Development Center, Mianyang 621000, China

Received:2019-11-13 Revised:2019-12-27 Online:2020-08-15 Published:2020-01-02
Supported by:
National Basic Research Program of China (2015CB755800)

摘要/Abstract

摘要： 为明晰结冰风洞中过冷大水滴（SLD）云雾演化特性，发展了基于欧拉法的SLD液滴运动、传热和传质耦合计算方法，针对3 m×2 m结冰风洞主试验段水平收缩构型，分析了SLD云雾沉降收缩特性、动量平衡特性和热平衡特性，探索了液滴变形破碎的影响，评估了构型出口处SLD液滴动量平衡和热平衡状态。研究结果表明：直径超过250 μm的SLD液滴在构型内会发生显著形变，液滴尺寸越大则变形程度越强，尤其在160 m/s工况下，当液滴直径超过750 μm后，SLD液滴甚至会发生破碎；液滴变形破碎效应会增大液滴加速度和液滴温度下降率，促使SLD液滴趋近动量平衡和热平衡状态；SLD云雾（最大液滴直径小于1 000 μm）在构型出口处会出现显著的粒径浓度分层、动量分层和热分层现象，其中直径小于100 μm的小尺寸液滴速度快、温度低且不断凝结，趋于平衡态，但直径超过500 μm的大尺寸SLD液滴速度慢、温度高且不断蒸发，则显著偏离平衡态；增大试验段气流速度尽管会减弱SLD云雾粒径浓度分层程度，但会增强动量分层和热分层程度，尤其在160 m/s工况下，SLD云雾会均匀分布在构型出口中心区域内（-0.75 m < Y < 0.75 m且-0.5 m < Z < 0.5 m），与其平衡态间的最大速度差和温度差将分别超过18 m/s和20℃。

关键词: 飞机结冰, 结冰风洞, 过冷大水滴, 云雾演化, 欧拉法

Abstract: In order to understand the Supercooled Large Droplet (SLD) cloud evolution characteristics in icing wind tunnels, a method based on Eulerian theory is developed to simulate SLD cloud flows coupling of momentum, mass and heat transfer. Using this method, the evolution of SLD cloud is investigated for the horizonal contraction configuration of the main section in a 3 m×2 m icing wind tunnel. The SLD cloud evolution characteristics, including cloud sinking and contraction, cloud momentum equilibrium and cloud thermal equilibrium, are analyzed. The effects of droplet deformation and breakup on the SLD mechanical equilibrium characteristic are explored. The states of cloud momentum equilibrium and thermal equilibrium in the main test section are evaluated. Results show that SLD droplets with diameter larger than 250 μm would experience significant deformation in the configuration. The increased droplet size could enhance the degree of deformation. Particularly at the test section velocity of 160 m/s, SLD droplets with the diameter larger than 750 μm would break up. Then the effects of droplet deformation and breakup could increase the rates of droplet acceleration and temperature reduction, so that the SLD droplets would approach the momentum and thermal equilibrium states. The SLD cloud with the maximum diameter of 1000 μm shows size and concentration stratification, momentum stratification and thermal stratification at the configuration outlet. In the SLD cloud, the small droplets with diameter smaller than 100 μm would reach the equilibrium states with higher speed, lower temperature and constant condensation. However, large SLD droplets with diameter larger than 500 μm would deviate from equilibrium states significantly with lower speed, higher temperature and constant evaporation. Increased test section velocity could reduce the degree of SLD cloud concentration stratification, but enhance the degree of momentum and thermal stratification. In particular, in the condition with the test section of 160 m/s, the SLD cloud will be uniform in the central area of the outlet (almost -0.75 m < Y < 0.75 m and -0.5 m < Z < 0.5 m), but the maximum velocity difference and temperature difference compared with the equilibrium states are higher than 18 m/s and 20 ℃, respectively.

Key words: aircraft icing, icing wind tunnel, supercooled large droplet, cloud evolution, Eulerian method

中图分类号:

V221.3

郭向东, 柳庆林, 刘森云, 王梓旭, 李明. 结冰风洞中过冷大水滴云雾演化特性数值研究[J]. 航空学报, 2020, 41(8): 123655-123655.

GUO Xiangdong, LIU Qinglin, LIU Senyun, WANG Zixu, LI Ming. Numerical study of supercooled large droplet cloud evolution characteristics in icing wind tunnel[J]. ACTA AERONAUTICAET ASTRONAUTICA SINICA, 2020, 41(8): 123655-123655.

参考文献 33

[1]	杜雁霞,李明,桂业伟,等. 飞机结冰热力学行为研究综述[J]. 航空学报,2017,38(2):520706. DU Y X, LI M, GUI Y W, et al. Review of thermodynamic behaviors in aircraft icing process[J]. Acta Aeronautica et Astronautica Sinica, 2017, 38(2):520706(in Chinese).
[2]	桂业伟,周志宏,李颖晖,等. 关于飞机结冰的多重安全边界问题[J]. 航空学报,2017,38(2):520723. GUI Y W, ZHOU Z H, LI Y H, et al. Multiple safety boundaries protection on aircraft icing[J]. Acta Aeronautica et Astronautica Sinica, 2017, 38(2):520723(in Chinese).
[3]	COBER S, BERNSTEIN B, JECK R, et al. Data and analysis for the development of an engineering standard for supercooled large drop conditions:DOT/FAA/AR-09/10[R]. Washington, D. C.:FAA, 2009.
[4]	Federal Aviation Administration. Airplane and engine certification requirements in supercooled large drop, mixed phase, and ice crystal icing conditions[S]. Washington, D.C.:FAA,2014.
[5]	STEEN L E, IDE R F, VAN ZANTE J F, et al. NASA Glenn icing research tunnel:2014 and 2015 cloud calibration procedures and results:NASA/TM-2015-218758[R]. Washington, D. C.:NASA, 2015.
[6]	BELLUCCI M, ESPOSITO B M, MARRAZZO M, et al. Calibration of the CIRA IWT in the low speed configuration[c]//45th AIAA Aerospace Science Meeting and Exhibit. Reston:AIAA, 2007.
[7]	TAN S C, PAPADAKIS M. General effects of large droplet dynamics on ice accretion modeling[C]//41st AIAA Aerospace Science Meeting and Exhibit. Reston:AIAA, 2003.
[8]	WILLBANKS C E, SCHULZT R J. Analytical study of icing simulation for turbine engines in altitude test cells[J]. Journal of Aircraft, 1975, 12(12):960-967.
[9]	GATES E M, LAM W, LOZOWSKI E P. Spray evolution in icing wind tunnels[J]. Cold Regions Science and Technology, 1988, 15:65-74.
[10]	MILLER D R, ADDY H E, IDE R F. A study of large droplet ice accretions in the NASA GLENN IRT at near-freezing conditions:NASA TM-1996-107142-REV1[R]. Washington, D. C.:NASA, 1996.
[11]	IDE R F, OLDENBURG J R. Icing cloud calibration of the NASA Glenn icing research tunnel:NASA/TM-2001-210689[R]. Washington, D. C.:NASA, 2001.
[12]	VAN ZANTE J F, IDE R F, STEEN L E. NASA Glenn icing research tunnel:2012 cloud calibration procedure and results[C]//4th AIAA Atmospheric and Space Environments Conference. Reston:AIAA, 2012.
[13]	KING-STEEN L E, IDE R F. Creating a bimodal drop-size distribution in the NASA Glenn icing research tunnel[C]//9th AIAA Atmospheric and Space Environments Conference. Reston:AIAA, 2017.
[14]	ORCHARD D M, SZILDER K, DAVISON C R. Design of an icing wind tunnel contraction for supercooled large drop conditions[C]//2018 AIAA Atmospheric and Space Environments Conference. Reston:AIAA, 2018.
[15]	易贤, 马洪林, 王开春, 等. 结冰风洞液滴运动及传质传热特性分析[J]. 四川大学学报, 2012, 44(2):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]. Journal of Sichuan University, 2012, 44(2):132-135(in Chinese).
[16]	DU Q, WANG Z Z, ZHU C L. Analysis and calculation of droplet-air mixed phase flow model in icing wind tunnel[C]//32nd AIAA Aerodynamic Measurement Technology and Ground Testing Conference. Reston:AIAA, 2016.
[17]	郭向东,王梓旭,李明,等.结冰风洞中液滴过冷特性数值研究[J].航空学报, 2017, 38(10):121254. GUO X D, WANG Z X, LI M, et al. Numerical study of droplet supercooling in an icing wind tunnel[J]. Acta Aeronautica et Astronautica Sinica, 2017, 38(10):121254(in Chinese).
[18]	郭向东,王梓旭,李明,等.结冰风洞中液滴相变效应数值研究[J].航空学报, 2018, 39(1):121586. GUO X D, WANG Z X, LI M, et al. Numerical investigation of phase transition effects of droplet in icing wind tunnel[J]. Acta Aeronautica et Astronautica Sinica, 2018, 39(1):121586(in Chinese).
[19]	LI H X, CHEN F, HU H. Simultaneous measurements of droplet size, flying velocity and transient temperature of in flight droplets by using a molecular tagging technique[J]. Experiments in Fluids, 2015, 56(10):194.
[20]	NORDE E. Eulerian method for ice crystal icing in turbofan engines[D]. The Netherlands:University of Twente, 2017.
[21]	NORDE E, VAN DER WEIDE E T A, HOEIJMAKERS H W M. Eulerian method for ice crystal icing[J]. AIAA Journal, 2018, 56(1):222-234.
[22]	LUXFORD G. Experimental and modelling investigation of the deformation, drag and break-up of drizzle droplets subjected to strong aerodynamic forces in relation to SLD aircraft icing[D]. Cranfield:Cranfield University, 2005.
[23]	LUXFORD G, HAMMOND D, IVEY P. Modelling, imaging and measurement of distortion, drag and break-up of aircraft-icing droplets[C]//43rd AIAA Aerospace Science Meeting and Exhibit. Reston:AIAA, 2005.
[24]	CLIFT R, GRACE J R, WEBER M E. Bubbles, drops and particles[M]. New York:Academic Press, 1978.
[25]	KIM I, BACHCHAN N, PEROOMIAN O. Supercooled large droplet modeling for aircraft icing using an Eulerian-Eulerian approach[J]. Journal of Aircraft, 2016, 53(2):487-500.
[26]	HONSEK R, HABASHI W G. Eulerian modeling of in-flight icing due to supercooled large droplets[J]. Journal of Aircraft, 2008, 45(4):1290-1296.
[27]	PILCH M, ERDMAN C A. Use of breakup time data and velocity history data to predict the maximum size of stable fragments for acceleration-induced breakup of a liquid drop[J]. International Journal of Multiphase Flow, 1987, 13(6):741-757.
[28]	CENGEL Y A, GHAJAR A J. Heat and mass transfer:fundamentals & applications[M]. 5th ed. New York:McGraw-Hill Education, 2015.
[29]	VILLEDIEU P, TRONTIN P, CHAUVIN R. Glaciated and mixed-phase ice accretion modeling using ONERA 2D icing suite[C]//6th AIAA Atmospheric and Space Environments Conference. Reston:AIAA, 2014.
[30]	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.
[31]	MILLER R S, HARSTAD K, BELLAN J. Evaluation of equilibrium and non-equilibrium evaporation models for many-droplet gas-liquid flow simulations[J]. International Journal of Multiphase Flow, 1998, 24:1025-1055.
[32]	HINDMARSH J P, RUSSEL A B, CHEN X D. Experimental and numerical analysis of the temperature transition of a suspended freezing water droplet[J]. International Journal of Heat and Mass Transfer, 2003, 46:1199-1213.
[33]	KING M C, BACHALO W D, BELL D, et al. Weber number tests in the NASA icing research tunnel[C]//2018 AIAA Atmospheric and Space Environments Conference. Reston:AIAA, 2018.

E-mail：hkxb@buaa.edu.cn

关于我们

期刊社服务

专业学科

封面文章

友情链接

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

结冰风洞中过冷大水滴云雾演化特性数值研究

Numerical study of supercooled large droplet cloud evolution characteristics in icing wind tunnel

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

参考文献 33

相关文章 15

编辑推荐

Metrics

本文评价

[1]	熊华杰, 安怡竞, 吴主龙, 周志宏. 欧拉壁膜模型在结冰模拟中的扩展[J]. 航空学报, 2023, 44(S2): 729287-729287.
[2]	杨映麟, 常士楠, 石奇玉, 齐海峰. 气流剪切作用下的水膜实验[J]. 航空学报, 2023, 44(S2): 729290-729290.
[3]	石达志, 桑为民, 李世杰, 安博. 基于LBM的结冰表面水滴撞击特性数值模拟[J]. 航空学报, 2023, 44(S2): 729192-729192.
[4]	许岭松, 吴渊, 朱东宇, 张付昆, 刘昱. FL⁃61结冰风洞翼型俯仰振荡机构研制[J]. 航空学报, 2023, 44(S2): 729296-729296.
[5]	何磊, 钱炜祺, 董康生, 易贤, 柴聪聪. 基于卷积神经网络的结冰翼型气动特性建模[J]. 航空学报, 2023, 44(5): 126434-126434.
[6]	任靖豪, 王强, 李维浩, 刘宇, 易贤. 基于梯度下降的水滴收集率计算方法[J]. 航空学报, 2023, 44(4): 126381-126381.
[7]	陈勇, 孔维梁, 刘洪. 飞机过冷大水滴结冰气象条件运行设计挑战[J]. 航空学报, 2023, 44(1): 626973-626973.
[8]	束珺, 徐东光, 韩志熔, 李斯, 黄雄. 结冰风洞过冷大水滴试验中混合翼设计[J]. 航空学报, 2023, 44(1): 627182-627182.
[9]	牛俊杰, 桑为民, 李栋, 郝莲, 王泽林. 一种基于水滴收集量代理模型的结冰试飞空域确定方法[J]. 航空学报, 2023, 44(1): 627697-627697.
[10]	许文纲, 王志颖, 孙闯, 严如强, 陈雪峰. 转频调制下齿轮泵压力脉动机理[J]. 航空学报, 2022, 43(9): 625508-625508.
[11]	王昭力, 曾涛, 周志宏, 陈宇. TVD格式在水滴流场数值模拟中的应用[J]. 航空学报, 2022, 43(12): 627010-627010.
[12]	李浩然, 段玉宇, 张宇飞, 陈海昕. 结冰模拟软件AERO-ICE中的关键数值方法[J]. 航空学报, 2021, 42(S1): 726371-726371.
[13]	郭向东, 张平涛, 赵照, 赖庆仁, 郭龙. 大型结冰风洞云雾场适航应用符合性验证[J]. 航空学报, 2020, 41(10): 123879-123879.
[14]	魏扬, 徐浩军, 薛源, 郑无计, 李哲, 裴彬彬. 基于神经网络自适应动态逆的结冰飞机飞行安全边界保护方法[J]. 航空学报, 2019, 40(5): 122488-122488.
[15]	雷梦龙, 常士楠, 杨波. 基于Myers模型的三维结冰数值仿真[J]. 航空学报, 2018, 39(9): 121952-121962.