旋转整流罩表面球形脱落冰运动轨迹的数值模拟

推进技术 ›› 2016, Vol. 37 ›› Issue (9): 1609-1616.

旋转整流罩表面球形脱落冰运动轨迹的数值模拟

胡娅萍,吉洪湖,吴铁鹰,王健

南京航空航天大学能源与动力学院，江苏南京 210016,南京航空航天大学能源与动力学院，江苏南京 210016,南京航空航天大学能源与动力学院，江苏南京 210016,南京航空航天大学能源与动力学院，江苏南京 210016

发布日期:2021-08-15
作者简介:胡娅萍，女，博士，副教授，研究领域为航空发动机积冰、防冰以及密封技术。

Numerical Simulation of Spherical Shed Ice

College of Energy and Power Engineering，Nanjing University of Aeronautics and Astronautics， Nanjing 210016，China,College of Energy and Power Engineering，Nanjing University of Aeronautics and Astronautics， Nanjing 210016，China,College of Energy and Power Engineering，Nanjing University of Aeronautics and Astronautics， Nanjing 210016，China and College of Energy and Power Engineering，Nanjing University of Aeronautics and Astronautics， Nanjing 210016，China

Published:2021-08-15

摘要/Abstract

摘要： 为了合理地设计旋转整流罩的防冰系统，以旋转坐标系为参考系，建立了旋转整流罩表面球形脱落冰在涡扇发动机进口流场中运动的数学模型，考虑了脱落冰所受的气动曳力、压差阻力、离心力以及哥氏力，并据此开发了球形脱落冰运动轨迹的计算程序。采用该程序与商业软件NUMECA计算得到的发动机进口空气流场数据对球形脱落冰的运动过程进行了数值模拟。计算结果表明，球形冰的运动轨迹线在靠近整流罩处接近于沿径向的直线，此后球形冰在远离旋转轴的同时也向下游弯曲；直径较小和脱落位置靠近下游的球形冰更易进入发动机内，并且其跟随旋转轴转动的趋势越显著。

关键词: 防冰；旋转整流罩；球形脱落冰；运动轨迹；数值模拟

Abstract: For designing the anti-icing system of rotating spinner，the mathematical model of motion of a spherical shed ice in the flowfield of turbofan engine entrance was established under the rotating coordinates system，which took the aerodynamic drag force，the pressure drag force，the centrifugal force and Coriolis force into account. Subsequently computation code of motion trajectory is developed based on the model. Using the code and the flowfield data of some turbofan engine entrance computed by NUMECA，the motion trajectory of spherical shed ice was simulated. The results show that the motion trajectory is close to radial line near the spinner. And then the trajectory is bent towards downstream along the axis direction as well as the shed ice is far away from the axis. The shed ice with smaller diameter and more downstream shedding position would tend to move around the axis and enter into the engine.

Key words: Anti-icing；Rotating spinner；Spherical shed ice；Motion trajectory；Numerical simulation

胡娅萍,吉洪湖,吴铁鹰,王健. 旋转整流罩表面球形脱落冰运动轨迹的数值模拟[J]. 推进技术, 2016, 37(9): 1609-1616.

HU Ya-ping,JI Hong-hu,WU Tie-ying,WANG Jian. Numerical Simulation of Spherical Shed Ice[J]. Journal of Propulsion Technology, 2016, 37(9): 1609-1616.

参考文献

[1] 裘燮纲. 飞机防冰系统[M]. 南京:航空专业教材编审组, 1996.
[2] Britton R K, Bond T H. An Overview of Shed Ice Impact Studies in the NASA Lewis Icing Research Tunnel[R]. AIAA 93-0301.
[3] Howard B, Jay O, James J S, et al. Determination of Shed Ice Particle Size Using High Speed Digital Imaging[R]. NASA TM-07406, 1996.
[4] Jonathan W N, Ethiraj V, Mark P L. Particle Impact Risk Assessment for Ablative Thermal Protection Systems[R]. AIAA 98-0166.
[5] Wright W B, Keith T G, DeWitt K J. Numerical Simulation of Icing, Deicing and Shedding[R]. AIAA 91-0665.
[6] Kohlman D L, Winn R C. Analytical Prediction of Trajectories of Ice Pieces after Release in an Airstream[R].AIAA 2001-0680.
[7] Chandrasekharan R, Hinson M. Trajectory Simulation of Ice Shed from a Business Jet[R]. SAE-2003-01-3023.
[8] Papadakis M, Yeong H W, Suares I G. Experimental and Computational Investigation of Ice Shedding from Aircraft Surfaces[R]. AIAA 2006-1010.
[9] Papadakis M, Yeong H W, Suares I G. Simulation of Ice Shedding from a Business Jet Aircraft[R]. AIAA 2007-506.
[10] Papadakis M, Yeong H W, Suares I G. Parametric Investigation of Ice Shedding from a Business Jet Aircraft[R]. SAE-2007-01-3359.
[11] Papadakis M, Yeong H W, Shimoi K, et al. Ice Shedding Experiments with Simulated Ice Shapes[R]. AIAA 2009-3972.
[12] Yeong H W, Papadakis M, Shimoi K. Ice Trajectory and Monte Carlo Analyses for a Business Jet[R]. AIAA 2009-3973.
[13] Baruzzi G S, Lagacé P, Aubé M S, et al. Development of a Shed-Ice Trajectory Simulation in FENSAP-ICE[R]. SAE-2007-01-3360.
[14] Widhalm M. Lagrangian Trajectory Simulation of Rotating Regular Shaped Ice Particles[R]. SAE-2015-01-2141.
[15] 肖春华, 林贵平, 桂业伟, 等. 冰脱落对电热除冰传热特性的影响研究[J]. 空气动力学学报, 2012, 30(4): 551-556.
[16] 王建涛, 易贤, 肖中云, 等. ARJ21-700飞机冰脱落数值模拟[J]. 空气动力学学报, 2013, 31(4): 430-436.
[17] 胡娅萍, 吉洪湖, 王健, 等. 锥角对旋转整流罩积冰影响的模拟实验[J]. 航空动力学报, 2014, 29(3): 495-503.
[18] 王健, 胡娅萍, 吉洪湖, 等. 旋转整流罩积冰生长与脱落过程的实验[J]. 航空动力学报, 2014, 29(6):1352-1357.
[19] 张丽芬, 张美华, 吴丁毅, 等. 旋转帽罩结冰相似准则的研究[J]. 推进技术, 2015, 36(8): 1164-1169. (ZHANG Li-fen, ZHANG Mei-hua, WU Ding-yi, et al. Research on Icing Scaling Law for Rotating Cone[J].Journal of Propulsion Technology, 2015, 36(8): 1164-1169.)
[20] 王云. 航空发动机原理[M]. 北京:北京航空航天大学出版社, 2009.
[21] 李晓娟, 桂幸民. 风扇∕增压级设计与非设计性能数值模拟[J]. 推进技术, 2005, 26(6): 522-525. (LI Xiao-juan, GUI Xing-min. Numerical Simulation of Fan/Compressor Design and Off Design Performance[J]. Journal of Propulsion Technology, 2005, 26(6): 522-525.)
[22] 金东海, 陈佳, 桂幸民. 压气机叶栅多点气动优化设计[J]. 推进技术, 2007, 28(4): 367-372. (JIN Dong-hai, CHEN Jia, GUI Xing-min. Multi-Point Aerodynamic Design Optimization for Compressor Cascade Airfoils[J]. Journal of Propulsion Technology, 2007, 28(4): 367-372.)(编辑:朱立影)　　　　　　　　　　　　　*　收稿日期:2015-01-25；修订日期:2015-08-15。作者简介:胡娅萍,女,博士,副教授,研究领域为航空发动机积冰、防冰以及密封技术。E-mail: hyp@nuaa.edu.cn