﻿ 基于欧拉法模拟旋转帽罩水滴撞击特性
 文章快速检索 高级检索

Simulation of droplet impingement characteristics of spinner based on Eulerian method
Wu Menglong, Chang Shinan, Leng Mengyao, Wang Chao
School of Aeronautic Science and Engineering, Beijing University of Aeronautics and Astronautics, Beijing 100191, China
Abstract:Based on Eulerian theory, the droplet impingement characteristics of the spinner were simulated. The droplet impingement characteristics of the spinner were achieved by calculating the droplets' flow field with numerical diffusion. The rotate speed impact on the droplet impingement characteristics of spinner was analyzed. Results show that the greater the rotate speed, the smaller the droplet impingement limit. In typical flight conditions and weather conditions, the centrifugal force on the trajectory of droplets is much smaller than the inertial force. The rotate speed has little influence on the droplet impingement characteristics of the spinner.
Key words: aircraft icing     multiphase flow     Eulerian method     droplet impingement characteristics     spinner

1) 把气体当作连续介质,而将颗粒视为离散体系,即拉格朗日法;

2) 把气体与颗粒都看成共同存在且相互渗透的连续介质,把颗粒作为拟流体,即欧拉法.

1.2 边界条件

 图 1 水滴相壁面边界条件 Fig. 1 Wall boundary condition of droplet flowfield
1.3 耗散函数

 图 2 旋转帽罩计算域网格 Fig. 2 Computational mesh of the spinner
2.2 计算条件

 状态 旋转速度(r/min) 压强/Pa 飞行速度/(m/s) 温度/K MVD/μm LWC/(g/m3) 状态1 1000 89867 78 260 20 1 状态2 3000 89867 78 260 20 1 状态3 5000 89867 78 260 20 1 状态4 5000 89867 170 260 20 1 状态5 5000 89867 50 260 20 1 状态6 5000 89867 78 260 10 1 状态7 5000 89867 78 260 40 1 状态8 0 89867 50 260 10 1 状态9 5000 89867 50 260 10 1 注：MVD—水滴平均容积直径；LWC—液态水含量.
3 计算结果与分析

 图 3 状态3，y=0截面局部水滴收集系数分布曲线 Fig. 3 Case 3, cross-section y=0,local water droplets collect coefficient distribution curves

 图 4 状态1~3,y=0截面局部水滴收集系数分布曲线 Fig. 4 Case 1~3,cross-section y=0,local water droplets collect coefficient distribution curves

 图 5 状态3,速度78m/s,转速5000r/min水滴轨迹图 Fig. 5 Case 3,velocity=78m/s,rotate speed= 5000r/min,the scheme of droplet trajectory
3.2 来流速度对水滴撞击特性影响

 图 6 状态3~5,y=0截面局部水滴收集系数分布曲线 Fig. 6 Case 3~5,cross-section y=0,local water droplets collect coefficient distribution curves
3.3 水滴平均容积直径对水滴撞击特性影响

 图 7 状态3,6,7,y=0截面局部水滴收集系数分布曲线 Fig. 7 Case 3,6,7,cross-section y=0,local water droplets collect coefficient distribution curves

 图 8 状态8和9,y=0截面局部水滴收集系数分布曲线 Fig. 8 Case 8 and 9,cross-section y=0,local water droplets collect coefficient distribution curves

1) 欧拉法求解旋转帽罩的空气-水滴两相流流场,继而得到旋转帽罩的水滴撞击特性是可行的.

2) 在飞行条件与气象条件保持不变的情况下,随着旋转帽罩转速增加,水滴撞击极限稍有减小.

3) 在典型飞行条件与气象条件下,离心力对水滴运动轨迹的影响远小于惯性力对其的影响.旋转帽罩转速对水滴的运动轨迹影响较小,水滴撞击特性变化不大.

 [1] 裘燮纲,韩凤华.飞机防冰系统[M].北京:航空专业教材编审组,1985:50-53Qiu Xiegang,Han Fenghua.Aircraft anti-icing system[M].Beijing:Compilation and Examination Group of Aero Specialized Teaching Materials,1985:50-53(in Chinese) [2] 周力行.湍流气粒两相流动和燃烧的理论与数值模拟[M].北京:科学出版社,1994:293-298Zhou Lixing.Theory and numerical modeling of turbulent gas-particle flows and combustion[M].Beijing:Science Press,1994:293-298(in Chinese) [3] 杨倩,常士楠,袁修干.水滴撞击特性的数值计算方法研究[J].航空学报,2002,23(2):173-176 Yang Qian,Chang Shinan,Yuan Xiugan.Study on numerical method for determining the droplet trajectories[J].Acta Aeronuatica et Astronautica Sinca,2002,23(2):173-176(in Chinese) Cited By in Cnki (69) [4] 常士楠,苏新明,邱义芬.三维机翼结冰模拟[J].航空学报,2011,32(2):212-222 Chang Shinan,Su Xinming,Qiu Yifen.Ice accretion simulation on three dimensional wings[J].Acta Aeronautica et Astronautica Sinca,2011,32(2):212-222(in Chinese) Cited By in Cnki (9) [5] 申晓斌,林贵平,杨胜华.三维发动机进气道水滴撞击特性分析[J].北京航空航天大学学报,2011,37(1):1-5 Shen Xiaobin,Lin Guiping,Yang Shenghua.Analysis on three dimensional water droplets impingement characteristics of engine inlet[J].Journal of Beijing University of Aeronautics and Astronautics,2011,37(1):1-5(in Chinese) Cited By in Cnki (3) [6] Bourgault Y,Boutanios Z,Habashi W G.Three-dimensional eulerian approach to droplet impingement simulation using FENSAP-ICE,part 1:model, algorithm, and validation[J].Journal of Aircraft,2000,37(1):95-103 Click to display the text [7] 赵秋月,董威,朱剑鋆.发动机旋转整流帽罩的水滴撞击特性分析[J].燃气涡轮试验与研究,2011,24(4):32-35 Zhao Qiuyue,Dong Wei,Zhu Jianjun.Droplets impinging characteristic analysis of the rotating fairing of aero-engine[J].Gas Turbine Experiment and Research,2011,24(4):32-35(in Chinese) Cited By in Cnki (2) [8] Das K,Hamed A,Basu D.Ice shape prediction for turbofan rotating blades[R].AIAA 2006-0209,2006 [9] 杨胜华,林贵平,申晓斌.三维复杂表面水滴撞击特性计算[J].航空动力学报,2010,25(2):284-290 Yang Shenghua,Lin Guiping,Shen Xiaobin.Water droplet impingement prediction for three-dimensional complex surfaces[J].Journal of Aerospace Power,2010,25(2):284-290(in Chinese) Cited By in Cnki (10) [10] Slater S A,Young J B.The calculation of inertial particle transport in dilute gas-particle flows[J].International Journal of Multiphase Flow,2001,27(1):61-87 Click to display the text [11] Wirogo S,Srirambhatla S.An Eulerian method to calculate the collection efficiency on two and three dimensional bodies[R].AIAA 2003-1073,2003 [12] Tong X,Luke E A.Eulerian simulation of icing collection efficiency using a singularity diffusion model[R].AIAA 2005-1246,2005 [13] Cao Y,Ma C,Zhang Q,et al.Numerical simulation of ice accretions on an aircraft wing[J].Aerospace Science and Technology,2012,23(1):296-304 Click to display the text [14] Linke D A.Systems of commercial turbofan engines[M].Berlin:Springer,2008:179-183 [15] 李西园.旋转部件结冰缩比试验相似性研究[D].北京:北京航空航天大学,2010 Li Xiyuan.Similarity of icing scaling test for aircraft rotating components[D].Beijing:Beijing University of Aeronautics and Astronautics,2010(in Chinese)

#### 文章信息

Wu Menglong, Chang Shinan, Leng Mengyao, Wang Chao

Simulation of droplet impingement characteristics of spinner based on Eulerian method

Journal of Beijing University of Aeronautics and Astronsutics, 2014, 40(9): 1263-1267.
http://dx.doi.org/10.13700/j.bh.1001-5965.2013.0559