基于欧拉壁面液膜模型的三维热气防冰腔数值仿真

doi:10.13700/j.bh.1001-5965.2017.0362

北京航空航天大学学报 ›› 2018, Vol. 44 ›› Issue (5): 959-966.doi: 10.13700/j.bh.1001-5965.2017.0362

基于欧拉壁面液膜模型的三维热气防冰腔数值仿真

李延, 郭涛, 常红亮

中国航空工业集团公司第一飞机设计研究院, 西安 710089

收稿日期:2017-05-31 出版日期:2018-05-20 发布日期:2018-05-29
通讯作者: 李延 E-mail:276637284@qq.com
作者简介:李延,男,硕士,工程师。主要研究方向:飞机防除冰系统设计与仿真计算;郭涛,男,博士,高级工程师。主要研究方向:飞机防除冰系统设计;常红亮,男,学士,研究员。主要研究方向:飞机防除冰系统设计。

Numerical simulation of 3D hot-air anti-icing chamber based on Eulerian wall film model

LI Yan, GUO Tao, CHANG Hongliang

The First Aircraft Institute, Aviation Industry Corporation of China, Xi'an 710089, China

Received:2017-05-31 Online:2018-05-20 Published:2018-05-29

摘要/Abstract

摘要： 提出了一种基于欧拉壁面液膜（EWF）模型的热气防冰腔性能仿真计算的新方法。通过FLUENT软件用户自定义标量（UDS）框架求解水滴控制方程获取三维表面水滴撞击特性。通过对各微元的水收集率、水膜蒸发率等进行质量平衡分析得到了通过该微元的质量流量，并以此作为EWF模型质量流量边界条件进行空气驱动下三维水膜厚度分布的计算，进而建立了防冰表面水膜流动动态模型。在此基础上建立了适用于三维防冰表面的耦合换热模型，通过引入亚松弛因子实现了内外流场、水膜流动及蒙皮导热的松散耦合求解。通过对某发动机短舱模型三维算例计算结果的分析和对比，结果表明所采用的计算方法是合理可信的，可以用于三维防冰腔性能的计算。

关键词: 欧拉壁面液膜(EWF)模型, 热气防冰, 松散耦合, 防冰热载荷, 水膜, 数值仿真, 三维

Abstract: A new computation method of hot-air anti-icing chamber performance based on Eulerian wall film (EWF) model was presented in this paper. User defined scalar (UDS) of FLUENT software was used to calculate water droplet impingement efficiency by solving droplet governing equation. Mass balance analysis of water collection rate and film evaporation rate was performed to get the mass flow rate of each micro unit, which was used as boundary condition in calculating 3D film thickness distribution driven by air, and then the dynamic model of film flow on anti-icing surface was set up. Based on the above work, conjugate heat transfer model for 3D anti-icing surface was built, and under-relaxation factor was used in solving the loose coupling of inner/outer flow field, water film flow and skin heat conduction. The method of this paper was used in the calculation of a nacelle anti-icing chamber performance, and the results show good compliance with the physical phenomenon. The method in this paper can be used in 3D hot-air anti-icing chamber performance calculation.

Key words: Eulerian wall film (EWF) model, hot-air anti-icing, loose couple, anti-icing heat load, water film, numerical simulation, three-dimensional

中图分类号:

V244.1⁺.5

李延, 郭涛, 常红亮. 基于欧拉壁面液膜模型的三维热气防冰腔数值仿真[J]. 北京航空航天大学学报, 2018, 44(5): 959-966.

LI Yan, GUO Tao, CHANG Hongliang. Numerical simulation of 3D hot-air anti-icing chamber based on Eulerian wall film model[J]. JOURNAL OF BEIJING UNIVERSITY OF AERONAUTICS AND A, 2018, 44(5): 959-966.

参考文献

[1] 林贵平,卜雪琴,申晓斌,等.飞机结冰与防冰技术[M].北京:北京航空航天大学出版社,2016:30-46.LIN G P,BU X Q,SHEN X B,et al.Aircraft icing and anti-icing technology[M].Beijing:Beihang University Press,2016:30-46(in Chinese).
[2] DONG W,ZHU J,MIN X H.Calculation of the heat transfer and temperature on the aircraft anti-icing surface[C]//27th International Congress of the Aeronautical Sciences,2010:1-9.
[3] DOMINGOS R H,PAPADAKIS M,ZAMORA A O.Computational methodology for bleed air ice protection system parametric analysis:AIAA-2010-7834[R].Reston:AIAA,2010.
[4] MORENCY F,BRAHIMI M T,TEZOK F,et al.Anti-icing system simulation using CANICE[J].Journal of Aircraft,1999,36(6):999-1006.
[5] PELLISSIER M,HABASHI W G,PUEYO A.Design optimization of hot-air anti-icing systems by FENSAP-ICE:AIAA-2010-1238[R].Reston:AIAA,2010.
[6] WANG H Z,TRAN P,HABASHI W G.Anti-icing simulation in wet air of a piccolo system using FE NSAP-ICE:SAE 2007-01-3357[R].Warrendale:SAE International,2007.
[7] 李延,郭涛.基于松散耦合的三维热气防冰腔数值仿真[C]//第六届中国航空学会青年科技论坛文集.北京:航空工业出版社,2014:1245-1250.LI Y,GUO T.Numerical simulation of three dimensional thermal anti-icing chamber by loose couple method[C]//CSAA the sixth Youth Science and Technology Forum.Beijing:Aviation Industry Press,2014:1245-1250(in Chinese).
[8] 卜雪琴,林贵平,郁嘉.三维内外热耦合计算热气防冰系统表面温度[J].航空动力学报,2009,24(11):2495-2500.BU X Q,LIN G P,YU J.Three dimensional conjugate heat transfer for the surface temperature of wing hot-air anti-icing system[J].Journal of Aerospace Power,2009,24(11):2495-2500(in Chinese).
[9] 侯盼雪,林贵平,卜雪琴,等.后掠翼热气防冰系统数值仿真[J].航空学报,2012,33(5):809-817.HOU P X,LIN G P,BU X Q,et al.Numerical simulation of a swept wing hot-air anti-icing system[J].Acta Aeronautica et Astronautica Sinica,2012,33(5):809-817(in Chinese).
[10] 郁嘉,赵柏阳,卜雪琴,等.某型飞机发动机短舱热气防冰系统性能数值模拟[J].空气动力学学报,2016,34(3):302-307.YU J,ZHAO B Y,BU X Q,et al.Numerical simulation of the performance of an engine nacelle hot-air anti-icing system[J].Acta Aerodynamica Sinica,2016,34(3):302-307(in Chinese).
[11] Ansys Inc.ANSYS FLUENT 14 User Guide[Z].New Hampshire:ANSYS Inc.,2011.
[12] 毛晓东,卜雪琴,赵国昌,等.基于UDS框架的水撞击特性数值计算方法[J].沈阳航空航天大学学报,2016,33(1):8-12.MAO X D,BU X Q,ZHAO G C,et al.Calculation of water droplet impingement characteristics based on UDS frame[J].Journal of Shenyang Aerospace University,2016,33(1):8-12(in Chinese).
[13] WIROGO S,SRIRAMBHATLA S.An eulerian method to calculate the collect efficiency on two and three dimensional bodies:AIAA-2003-1073[R].Reston:AIAA,2003.
[14] TONG X L,LUKE E A.Eulerian simulations of icing collection efficiency using a singularity diffusion model:AIAA-2005-1246[R].Reston:AIAA,2005.
[15] 杨胜华,林贵平,申晓斌.三维复杂表面水滴撞击特性计算[J],航空动力学报,2010,25(2):284-290.YANG S H,LIN G P,SHEN X B.Water droplet impingement prediction for three dimensional complex surfaces[J].Journal of Aerospace Power,2010,25(2):284-290(in Chinese).
[16] 王超,常士楠,吴梦龙,等.过冷大水滴飞溅特性数值分析[J].航空学报,2014,35(4):1004-1011.WANG C,CHANG S N,WU M L,et al.Numerical investigation of splashing characteristics in super-cooled large droplet regime[J].Acta Aeronautica et Astronautica Sinica,2014,35(4):1004-1011(in Chinese).
[17] MESSINGER B L.Equilibrium temperature of an unheated icing surface as a function of air speed[J].Journal of the Aeronautical Sciences,1953,20(1):29-42.
[18] VELAZQUEZ M T,HANSMAN R J.Implementation of combined feather and surface-normal ice growth models in LEWICE/X:AIAA-1995-0735[R].Reston:AIAA,1995.
[19] 陈鹏,邓飞,刘思亭.三维空间属性插值方法的研究[J].电脑知识与技术,2015,11(7):235-239.CHEN P,DENG F,LIU S T.Research of three dimensional space attribute interpolation methods[J].Computer Knowledge and Technology,2015,11(7):235-239(in Chinese).

基于欧拉壁面液膜模型的三维热气防冰腔数值仿真

Numerical simulation of 3D hot-air anti-icing chamber based on Eulerian wall film model

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

本文评价

[1]	赵蓓蕾, 赵继广, 崔村燕, 王岩, 孙武华. 贮箱爆炸碎片初始速度及影响因素[J]. 北京航空航天大学学报, 2019, 45(9): 1797-1804.
[2]	史贤俊, 王康, 肖支才, 龙玉峰, 陈垚. 面向复杂系统的三维Bayes网络测试性验证模型[J]. 北京航空航天大学学报, 2019, 45(7): 1303-1313.
[3]	史绍琨, 赵久奋, 崇阳, 杨奇松, 尤浩. 带落角约束的新型二阶滑模三维制导律[J]. 北京航空航天大学学报, 2019, 45(3): 614-623.
[4]	杨彬, 李和平, 曾慧. 基于视频的三维人体姿态估计[J]. 北京航空航天大学学报, 2019, 45(12): 2463-2469.
[5]	张浩悦, 程晓琦, 刘畅, 孙军华. 基于全局稀疏地图的AGV视觉定位技术[J]. 北京航空航天大学学报, 2019, 45(1): 218-226.
[6]	邱丽芳, 陈海翔, 吴友炜. 新型单轴柔性铰链拓扑结构设计与柔度分析[J]. 北京航空航天大学学报, 2018, 44(6): 1133-1140.
[7]	翟敏, 刘华平, 张天昱, 卢山, 许静文. 基于智能手机APP的非合作目标三维重建与交互[J]. 北京航空航天大学学报, 2018, 44(12): 2637-2643.
[8]	张梦颖, 王华, 葛霖, 程浩. 大型铁磁性物体近场磁异常场数值仿真[J]. 北京航空航天大学学报, 2018, 44(11): 2454-2462.
[9]	蔡军, 赵原, 李宇豪, 解杨敏. 一种三维激光扫描系统的设计及参数标定[J]. 北京航空航天大学学报, 2018, 44(10): 2208-2216.
[10]	闫攀运, 梁国柱, 吕永志, 齐志会, 王玺. 车载冷发射系统多刚体动力学快速仿真研究[J]. 北京航空航天大学学报, 2017, 43(1): 71-78.
[11]	蒙志君, 平学寿, 陈旭智. 无人直升机三维避障方法及仿真[J]. 北京航空航天大学学报, 2016, 42(8): 1619-1626.
[12]	毕严先, 魏少明, 王俊, 张耀天, 孙忠胜, 袁常顺. 基于多假设跟踪的散射点关联和三维重构方法[J]. 北京航空航天大学学报, 2016, 42(6): 1219-1227.
[13]	毛晓东, 林贵平, 郁嘉. 弹射座椅不利姿态控制规律设计[J]. 北京航空航天大学学报, 2016, 42(3): 426-434.
[14]	刘振国, 亚纪轩, 林强, 王一博, 胡龙. 三维全五向编织复合材料细观结构实验分析[J]. 北京航空航天大学学报, 2016, 42(2): 211-217.
[15]	高霄鹏, 孙培成, 董祖舜, 吴彬, 魏可可. 支柱式沉浸水橇降载性能相关影响因素[J]. 北京航空航天大学学报, 2016, 42(2): 236-242.