航空发动机旋转帽罩结冰表面换热系数研究

推进技术 ›› 2017, Vol. 38 ›› Issue (4): 853-859.

航空发动机旋转帽罩结冰表面换热系数研究

张丽芬,刘振侠,张斐,张美华

西北工业大学动力与能源学院，陕西西安 710072,西北工业大学动力与能源学院，陕西西安 710072,西北工业大学动力与能源学院，陕西西安 710072,西北工业大学动力与能源学院，陕西西安 710072

发布日期:2021-08-15
作者简介:张丽芬，女，博士，讲师，研究领域为航空发动机结冰、防冰，以及冰颗粒摄入的动力学和热力学。
基金资助:
国家自然基金青年基金(51406167)；中央高校基本科研业务费专项资金(3102014JCQ01051)。

Numerical Simulation of Heat Transfer Coefficient on Icing Surface of Rotating Cone in Aero-Engine

College of Power and Energy，Northwestern Polytechnical University，Xi’an 710072，China,College of Power and Energy，Northwestern Polytechnical University，Xi’an 710072，China,College of Power and Energy，Northwestern Polytechnical University，Xi’an 710072，China and College of Power and Energy，Northwestern Polytechnical University，Xi’an 710072，China

Published:2021-08-15

摘要/Abstract

摘要： 旋转帽罩表面的对流换热系数是结冰、防冰研究的重要参数，旋转帽罩结冰模拟中对流换热量和蒸发量都需要已知对流换热系数。为了获取旋转帽罩表面的对流换热系数，采用数值模拟对旋转帽罩结冰表面的对流换热系数进行了研究。首先验证了网格与计算方法的合理性和可靠性，在此基础上，对影响表面换热的因素进行了分析，计算了不同转速、来流速度、来流温度对锥角40°，锥高176mm的旋转帽罩表面换热系数的影响。结果表明:锥尖区域，来流速度的影响与转速影响相比占主导；除锥尖外的其他区域，转速对换热系数的影响占据了主导地位；来流温度对整个表面的换热系数均有影响。采用不同工况下的数值模拟结果，在锥尖区域建立了来流雷诺数与努赛尔数之间的关联式；在除锥尖外的其他区域建立了旋转雷诺数与努赛尔数之间的关联式。

关键词: 航空发动机；旋转帽罩；对流换热系数；关联式

Abstract: Heat transfer coefficient is an important parameter for icing and anti-icing of rotating cone. It needs to be known for computing heat transfer and evaporation in numerical simulation of ice accretion. In order to obtain this parameter，heat transfer coefficient on icing surface of a rotating cone has been simulated by numerical method. Firstly，grid independence and numerical approach were confirmed. Then，factors that affect heat transfer were analyzed，including revolving velocity，inflow velocity and inflow temperature. Heat transfer coefficients influenced by these factors on icing surface of a rotating cone with half angle of 40°and height of 176 mm were predicted. The results show that in the cone tip area inflow velocity is more important than revolving velocity，and in other area the revolving velocity is dominated，compared with inflow velocity，and the influence of inflow temperature is also apparent for all icing surface. Finally，heat transfer coefficient correlations are obtained on the base of numerical results under different conditions，in the tip area heat transfer coefficient correlated formula is established using Nu number and inflow Re number，and in other area the correlation is based on Nu number and involving Re number.

Key words: Aero-engine；Rotating cone；Heat transfer coefficient；Correlation

张丽芬,刘振侠,张斐,张美华. 航空发动机旋转帽罩结冰表面换热系数研究[J]. 推进技术, 2017, 38(4): 853-859.

ZHANG Li-fen,LIU Zhen-xia,ZHANG Fei,ZHANG Mei-hua. Numerical Simulation of Heat Transfer Coefficient on Icing Surface of Rotating Cone in Aero-Engine[J]. Journal of Propulsion Technology, 2017, 38(4): 853-859.

参考文献

[1] Aliaga C N, Aube M S, Baruzzi G S, et al. Fensap-Ice-Unsteady: Unified Inflight Icing Simulation Methodology for Aircraft, Rotorcraft, and Jet Engines[J]. Journal of Aircraft, 2011, 48(1): 119-126.
[2] Verdin P G, Charpin J P F. Multi-Stepping Ice Prediction on Cylinders and other Relevant Geometries[J]. Journal of Aircraft, 2013, 50(3): 871-878.
[3] 易贤, 桂业伟, 朱国林. 飞机三维结冰模型及其数值求解方法[J]. 航空学报, 2010, 31(11): 2153-2158.
[4] 常士楠, 苏新明, 邱义芬. 三维机翼结冰模拟[J]. 航空学报, 2011, 32(2): 212-222.
[5] Cao Y H, Ma C, Zhang Q, et al. Numerical Simulatin of Ice Accretions on an Aircraft Wing[J]. Aerospace Science and Technology, 2012, 23(1): 296-304.
[6] 宋馨, 林贵平, 卜雪琴. 防冰表面的对流换热计算分析[J]. 北京航空航天大学学报, 2011, 37(1): 49-53.
[7] 卜雪琴, 林贵平, 郁嘉. 三维内外热耦合计算热气防冰系统表面温度[J]. 航空动力学报, 2009, 24(11): 2495-2500.
[8] 申晓斌, 林贵平, 卜雪琴, 等. 发动机进气道短舱前缘结冰三维模拟研究[J]. 航空学报, 2013, 34(3):517-524.
[9] 卜雪琴, 林贵平, 彭又新, 等. 防冰热载荷计算的一种新方法[J]. 航空学报, 2006, 27(2): 208-212.
[10] 卜雪琴, 林贵平. 基于CFD的水收集系数及防冰表面温度预测[J]. 北京航空航天大学学报, 2007, 33(10): 1182-1185.
[11] 赵曙, 朱惠人, 郭涛, 等. 旋转带肋回转通道流动换热数值模拟[J]. 西安交通大学学报, 2014, 48(2): 125-130.
[12] 徐磊, 常海萍, 潘金栋. 旋转条件下带出流孔的受限空间内冲击换热[J]. 推进技术, 2008, 29(2): 149-152. (XU Lei, CHANG Hai-ping, PAN Jin-dong. Heat Transfer of Impingement in Limited Room Configuration with Effusion in Rotating State[J]. Journal of Propulsion Technology, 2008, 29(2): 149-152.)
[13] 陈豪, 邓宏武, 李洋, 等. 高旋转数下45°斜肋回转通道平均换热热性研究[J]. 推进技术, 2015, 36(9): 1339-1346. (CHEN Hao, DENG Hong-wu, LI Yang, et al. Heat Transfer in Rotating Two-Pass Square Duct with 45 Ribbed Walls under High Rotation Numebers[J]. Journal of Propulsion Technology, 2015, 36(9): 1339-1346.)
[14] Bourgault Y, Habashi W G, Beaugendre. Development of a Shallow Water Icing Model in Fensap-Ice[R]. AIAA 99-0246.
[15] Hansman R H, Turnock S T. Investigation of Surface Water Behavior During Glaze Ice Accretion[J]. Journal of Aircraft, 1989, 26(2): 140-147.
[16] 陶文铨. 数值传热学(2版)[M]. 西安:西安交通大学出版社, 2001.
[17] 裘燮纲, 韩凤华. 飞机防冰系统[M]. 北京:航空专业教材编审组出版, 1985.
[18] 杨世铭, 陶文铨. 传热学(第三版)[M]. 北京:高等教育出版社, 1998.　　　　　　　　　　*　收稿日期:2015-04-26；修订日期:2016-02-28。基金项目:国家自然基金青年基金(51406167)；中央高校基本科研业务费专项资金(3102014JCQ01051)。作者简介:张丽芬,女,博士,讲师,研究领域为航空发动机结冰、防冰,以及冰颗粒摄入的动力学和热力学。 E-mail: zhanglifen@nwpu.edu.cn(编辑:梅瑛)