-
摘要:
相似性是风洞实验的一项基本要求，结冰风洞实验也不例外。为了系统性地研究结冰风洞实验中的相似性问题，首先，通过对飞行结冰问题进行分析和总结，找出了其中涉及的各种变量，并分析了其物理意义和量纲。其次，采用相似理论分析方法对结冰问题进行分析，得到了若干无量纲参数，并分析了这些无量纲参数的意义。再次，通过忽略一些不重要的影响因素，减少无量纲参数，并对飞行结冰中涉及的影响因素进行了合理地简化，最终得出了结冰实验的相似准则。最后，结合结冰风洞的运行参数，按照相似准则的要求得到了缩比模型结冰实验运行参数选取方法，并采用CFD方法进行了验证。结果表明：得到的缩比模型结冰实验运行参数选取方法是可行的。
-
关键词:
- 结冰 /
- 风洞实验 /
- 相似理论 /
- 缩比模型 /
- 计算流体力学(CFD)
Abstract:Similarity is the basic requirements of all kinds of wind tunnel tests including icing wind tunnel tests. To systematically study the similarity problem of icing wind tunnel test, the in-flight icing problem is analyzed and summed up, and then the variables involved are summarized, including both their physical meanings and dimensions. Then some dimensionless variables are obtained by applying the similarity theory analysis method to the in-flight icing problem and the physical meanings of the dimensionless variables are analyzed. Some simplifications to the dimensionless variables of the in-flight icing problem are conducted by ignoring the unimportant factors, and then the dominant similarity parameters of icing wind tunnel tests are obtained. Finally, the method of selecting the operating parameters of the icing wind tunnel test on scale models is achieved by applying the similarity laws, it is verified by CFD and its feasibility is confirmed.
-
Key words:
- icing /
- wind tunnel test /
- similarity theory /
- scale model /
- computational fluid dynamics (CFD)
-
图 5 霜冰冰形计算结果与实验结果
Figure 5. Computed and experimental icing results and experiment results under rime ice situation
图 6 明冰冰形计算结果与实验结果
Figure 6. Computed and experimental icing results and experiment results under glaze ice situation
图 7 霜冰冰形相似验证的计算结果
Figure 7. Similarity verification calculation results of icing under rime ice situation
图 8 明冰冰形相似验证的计算结果
Figure 8. Similarity verification calculation results of icingunder glaze ice situation
表 1 结冰问题中相关物理量及其量纲
Table 1. Relevant physical variables and their dimensions in icing problem
物理量性质 参数 意义 量纲 空间 δi(s, t) 积冰厚度 L δ 过冷水滴直径 L d 过冷水滴平均间距 L c 飞行器特征尺度 L s 模型表面曲线坐标 L hfilm 水膜厚度 L X(xi), i=1, 2, 3 水滴的空间位置矢量 L 时间 ti 积冰时间 T 物性参数及系统参数 νa 空气的运动黏性系数 L2T-1 νw 水的运动黏性系数 L2T-1 ρa 空气的密度 ML-3 ρw 水的密度 ML-3 ρi 冰的密度 ML-3 ka 空气的热传导系数 MLT-3Θ-1 kw 水的热传导系数 MLT-3Θ-1 ki 冰的热传导系数 MLT-3Θ-1 cp, a 空气的比热容 L2T-2Θ-1 cp, w 水的比热容 L2T-2Θ-1 cp, i 冰的比热容 L2T-2Θ-1 hf 水的凝固潜热 L2T-2 h 对流换热系数 MT-3Θ-1 流场参数 U∞ 来流速度(飞行速度) LT-1 ua(ua, i), i=1, 2, 3 空气气流速度 LT-1 ud(ud, i), i=1, 2, 3 过冷水滴速度 LT-1 uf(uf, i), i=1, 2, 3 水膜速度 LT-1 p 空气气流静压 ML-1T-2 Ta 空气流场的温度 Θ Td 水滴的温度 Θ 结冰表面参数 Ts 结冰表面温度 Θ σw/a 水-气界面张力系数 MT-2 σw/s 水-固表面张力系数 MT-2 其他变量 g 重力加速度 LT-2 Tfre 冰点温度 Θ 
下载: 导出CSV
表 2 结冰问题中的无量纲参数
Table 2. Dimensionless parameters in icing problem
参数 意义 
无量纲结冰厚度 
过冷水滴尺寸与模型尺寸之比 
液态水含量与空气密度的比值 
模型几何相似 
无量纲水膜厚度 
水滴空间位置相似 
无量纲积冰时间 
物性参数的比值 
普朗特数 
马赫数 
雅克伯数 
努赛尔数 
来流雷诺数 
速度场相似 
欧拉数 
温度场相似 
韦伯数 
弗劳德数
下载: 导出CSV
表 3 传统相似准则要求与本文相似准则的关系
Table 3. Relationship between traditional similarity laws and proposed similarity law
传统相似准则的要求 相关的无量纲参数 应用中满足相似性要求的条件 几何相似 
采用缩比模型，自然满足 流场相似 Π14=Ma 选取合理的速度和温度，容易满足 Π19=Re∞ 往往不容易满足，尽量配合其他条件，使其他条件优先满足 水滴运动相似 
令K保持不变 Π19=Re∞ 水滴收集相似 
选取合理液态水含量，容易满足 
取合理的积冰时间，容易满足 β 水滴运动相似的前提下即可满足 η 对流换热相似，依赖于Π18 水滴-结冰表面动力学相似 
明冰条件下非常重要 结冰表面热力学相似 Π18=Nu 很难满足，取决于Π4、Π9、Π11 
采用缩比模型，自然满足 Π19=Re∞ 往往不容易满足，尽量配合其他条件，使其他条件优先满足 Π11=Pr 自然满足
下载: 导出CSV
表 4 积冰冰形预测的参数
Table 4. Parameters in icing prediction
变量 霜冰情形 明冰情形 平均气动弦长c/m 0.533 4 0.533 4 来流速度U∞/(m·s-1) 102.8 67.1 来流静温T∞/K 256.49 265.07 来流静压p∞/Pa 90 760 90 760 迎角α/(°) 3.5 3.5 液态水含量LWC/(g·m-3) 0.55 1.30 平均水滴直径MVD/μm 20 30 积冰时间ti/s 420 360 积冰时间步数n 3 2 积冰时间步长Δt/s 140 180
下载: 导出CSV
表 5 霜冰情形相似性验证计算参数选取
Table 5. Calculation parameter selection in verification of similarity under rime ice situation
变量 全尺寸模型 缩比模型 缩比比例k0 1 1/3 平均气动弦长c/m 0.533 4 0.177 8 来流速度U∞/(m·s-1) 102.8 102.8 来流静温T∞/K 256.49 256.49 来流静压p∞/Pa 90 760 90 760 迎角α/(°) 3.5 3.5 液态水含量LWC/(g·m-3) 0.55 0.55 平均水滴直径MVD/μm 20 11.547 积冰时间ti/s 420 140 积冰时间步数n 3 3 积冰时间步长Δt/s 140 46.667
下载: 导出CSV
表 6 明冰情形相似性验证计算参数选取
Table 6. Calculation parameter selection in verification of similarity under glaze ice situation
变量 全尺寸模型 缩比模型 缩比比例k0 1 1/3 平均气动弦长c/m 0.533 4 0.177 8 来流速度U∞/(m·s-1) 67.1 96.775 来流静温T∞/K 265.07 265.07 来流静压p∞/Pa 90 760 90 760 迎角α/(°) 3.5 3.5 液态水含量LWC/(g·m-3) 1.30 1.30 平均水滴直径MVD/μm 30 14.423 积冰时间ti/s 360 83.2 积冰时间步数n 2 2 积冰时间步长Δt/s 180 41.6
下载: 导出CSV
-
[1] THOMAS S K, CASSONI R P, MACARTHUR C D.Aircraft anti-icing and de-icing techniques and modeling[J].Journal of Aircraft, 1996, 33(5):841-854. doi: 10.2514/3.47027 [2] 裘燮刚, 韩凤华.飞机防冰系统[M].北京:国防工业出版社, 2004.QIU X G, HAN F H.Anti-icing systems of aircraft[M].Beijing:National Defense Industry Press, 2004(in Chinese). [3] 朱春玲, 朱程香.飞机结冰及其防护[M].北京:科学出版社, 2016.ZHU C L, ZHU C X.Aircraft icing and its protection[M].Beijing:Science Press, 2016(in Chinese). [4] 林贵平, 卜雪琴, 申晓斌, 等.飞机结冰与防冰技术[M].北京:北京航空航天大学出版社, 2016.LIN G P, BU X Q, SHEN X B, et al.Aircraft icing and anti-icing/deicing technologies[M].Beijing:Beihang University Press, 2016(in Chinese). [5] 任鹏飞, 徐宇, 宋娟娟, 等.霜冰条件对翼型气动性能影响数值研究[J].工程热物理学报, 2014, 35(4):663-668. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gcrwlxb201404010REN P F, XU Y, SONG J J, et al.Numerical study about the influence of rime ice conditions on airfoil[J].Journal of Engineering Thermophysics, 2014, 35(4):663-668(in Chinese). http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gcrwlxb201404010 [6] CAO Y, HUANG J, XU Z, et al.Insight into rime ice accretion on an aircraft wing and corresponding effects on aerodynamic performance[J].The Aeronautical Journal, 2016, 120(1229):1101-1122. doi: 10.1017/aer.2016.49 [7] BRAGG M B, BROEREN A P, BLUMENTHAL L A.Iced-airfoil aerodynamics[J].Progress in Aerospace Sciences, 2005, 41(5):323-362. doi: 10.1016/j.paerosci.2005.07.001 [8] POURYOUSSEFI S G, MIRZAEI M, NAZEMI M M, et al.Experimental study of ice accretion effects on aerodynamic performance of an NACA 23012 airfoil[J].Chinese Journal of Aeronautics, 2016, 29(3):585-595. doi: 10.1016/j.cja.2016.03.002 [9] 李国知, 曹义华.旋翼结冰对直升机飞行动力学特性的影响[J].航空学报, 2011, 32(2):187-194. http://d.old.wanfangdata.com.cn/Periodical/hkxb201102001LI G Z, CAO Y H.Effects of rotor icing on helicopterflight dynamic characteristics[J].Acta Aeronautica et Astronautica Sinica, 2011, 32(2):187-194(in Chinese). http://d.old.wanfangdata.com.cn/Periodical/hkxb201102001 [10] 张恒, 李杰, 龚志斌.多段翼型缝翼前缘结冰大迎角分离流动数值模拟[J].航空学报, 2017, 38(2):520733. http://d.old.wanfangdata.com.cn/Periodical/hkxb201702005ZHANG H, LI J, GONG Z B.Numerical simulation of separated flow around a multi-element airfoil at high angle of attack with iced slat[J].Acta Aeronautica et Astronautica Sinica, 2017, 38(2):520733(in Chinese). http://d.old.wanfangdata.com.cn/Periodical/hkxb201702005 [11] 周莉, 徐浩军, 龚胜科, 等.飞机结冰特性及防除冰技术研究[J].中国安全科学学报, 2010, 20(6):105-110. doi: 10.3969/j.issn.1003-3033.2010.06.018ZHOU L, XU H J, GONG S K, et al.Research of aircraft icing characteristics and anti-icng and deicing technology[J].China Safty Science Journal, 2010, 20(6):105-110(in Chinese). doi: 10.3969/j.issn.1003-3033.2010.06.018 [12] BEAUGENDRE H, MORENCY F, HABASHI W G.FENSAP-ICE's three-dimensional in-flight ice accretion module:ICE3D[J].Journal of Aircraft, 2003, 40(2):239-247. doi: 10.2514/2.3113 [13] IULIANO E, MINGIONE G, PETROSINO F, et al.Eulerian modeling of large droplet physics toward realistic aircraft icing simulation[J].Journal of Aircraft, 2011, 48(5):1621-1632. doi: 10.2514/1.C031326 [14] JUNG S K, SHIN S, MYONG R S, et al.An efficient CFD-based method for aircraft icing simulation using a reduced order model[J].Journal of Mechanical Science and Technology, 2011, 25(3):703-711. doi: 10.1007/s12206-011-0118-4 [15] PETROSINO F, MINGIONE G, CAROZZA A, et al.Ice accretion model on multi-element airfoil[J].Journal of Aircraft, 2011, 48(6):1913-1920. doi: 10.2514/1.C031346 [16] 易贤.飞机积冰的数值计算与积冰实验相似准则研究[D].绵阳: 中国空气动力研究与发展中心, 2007. http://cdmd.cnki.com.cn/Article/CDMD-90113-2007189443.htmYI X.Numerical computation of aircraft icing and study on icing test scaling law[D].Mianyang: Chinese Aerodynamic Research and Development Center, 2007(in Chinese). http://cdmd.cnki.com.cn/Article/CDMD-90113-2007189443.htm [17] 雷梦龙, 常士楠, 杨波.基于Myers模型的三维结冰数值仿真[J].航空学报, 2018, 39(9):21952. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=hkxb201809003LEI M L, CHANG S N, YANG B.Three-dimensional numerical simulation of icing using Myers model[J].Acta Aeronautica et Astronautica Sinica, 2018, 39(9):21952(in Chinese). http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=hkxb201809003 [18] 常士楠, 苏新明, 邱义芬.三维机翼结冰模拟[J].航空学报, 2011, 32(2):212-222. http://d.old.wanfangdata.com.cn/Periodical/hkxb201102004CHANG S N, SU X M, QIU Y F.Ice accretion simulation on three dimensional wings[J].Acta Aeronautica et Astronautica Sinica, 2011, 32(2):212-222(in Chinese). http://d.old.wanfangdata.com.cn/Periodical/hkxb201102004 [19] HUANG X, MYERS B, D'AVIRRO J, et al.Icing wind-tunnel icing test on a contaminated full-scale Wing-Model at takeoff conditions: AIAA-2008-6417[R].Reston: AIAA, 2008. [20] 战培国.结冰风洞研究综述[J].实验流体力学, 2007, 21(3):92-96. doi: 10.3969/j.issn.1672-9897.2007.03.019ZHAN P G.A review of icing wind tunnel[J].Journal of Experiments in Fluid Mechanics, 2007, 21(3):92-96(in Chinese). doi: 10.3969/j.issn.1672-9897.2007.03.019 [21] 李周复.风洞特种实验技术[M].北京:航空工业出版社, 2010.LI Z F.Techniques in special wind tunnel tests[M].Beijing:Aviation Industry Press, 2010. [22] 李周复.风洞试验手册[M].北京:航空工业出版社, 2015.LI Z F.Handbook of wind tunnel test[M].Beijing:Aviation Industry Press, 2015(in Chinese). [23] 赵凯华.定性与半定量物理学[M].北京:高等教育出版社, 2008.ZHAO K H.Qualitative and semi-quantitative physics[M].Beijing:Higher Education Press, 2008(in Chinese). [24] 李之光.相似与模化[M].北京:国防工业出版社, 1982.LI Z G.Similarity and modeling[M].Beijing:National Defense Industry Press, 1982(in Chinese). [25] ANDERSON D N.Manual of scaling method: NASA/CR-2004-212875[R].Washington, D.C.: NASA, 2004. [26] 易贤, 周志宏, 杜雁霞, 等.考虑相变时间效应的结冰实验相似参数[J].实验流体力学, 2016, 30(2):14-19. http://www.cnki.com.cn/Article/CJFDTotal-LTLC201602003.htmYI X, ZHOU Z H, DU Y X, et al.An icing scaling parameter with the effects of phase change time[J].Journal of Experiments in Fluid Mechanics, 2016, 30(2):14-19(in Chinese). http://www.cnki.com.cn/Article/CJFDTotal-LTLC201602003.htm [27] 周志宏, 易贤, 桂业伟, 等.考虑水滴动力学效应的结冰实验相似准则[J].实验流体力学, 2016, 30(2):20-25. http://www.cnki.com.cn/Article/CJFDTotal-LTLC201602004.htmZHOU Z H, YI X, GUI Y W, et al.Icing scaling law with the dynamics effects of water droplets[J]Journal of Experiments in Fluid Mechanics, 2016, 30(2):20-25(in Chinese). http://www.cnki.com.cn/Article/CJFDTotal-LTLC201602004.htm [28] 李静.发动机进口静止/旋转部件结冰及结冰实验的相似研究[D].西安: 西北工业大学, 2015. http://cdmd.cnki.com.cn/Article/CDMD-10699-1016909385.htmLI J.Numerical simulation of ice accretion and study for icing scaling law on aero-engine entry components[D].Xi'an: Northwestern Polytechnical University, 2015(in Chinese). http://cdmd.cnki.com.cn/Article/CDMD-10699-1016909385.htm [29] 张丽芬, 张美华, 吴丁毅, 等.旋转帽罩结冰相似准则的研究[J].推进技术, 2015, 36(8):1164-1169. http://d.old.wanfangdata.com.cn/Periodical/tjjs201508007ZHANG L F, ZHANG M H, WU D Y, et al.Research on icing scaling law for rotating cone[J].Journal of Propulsion Technology, 2015, 36(8):1164-1169(in Chinese). http://d.old.wanfangdata.com.cn/Periodical/tjjs201508007 [30] BILANIN, A J.Proposed modifications to ice accretion/icing scaling theory[J].Journal of Aircraft, 1991, 28(6):353-359. doi: 10.2514/3.46034 [31] 杨世铭, 陶文铨.传热学[M].北京:高等教育出版社, 2006.YANG S M, TAO W Q.Heat transfer[M].Beijing:Higher Education Press, 2006(in Chinese). [32] 孔维梁.飞机异常结冰的过冷凝固机理及理论研究[D].上海: 上海交通大学, 2015. http://cdmd.cnki.com.cn/Article/CDMD-10248-1016787495.htmKONG W L.Mechanism and fundamental theory of supercooled water solidification on the abnormal aircraft icing[D].Shanghai: Shanghai Jiao Tong University, 2015(in Chinese). http://cdmd.cnki.com.cn/Article/CDMD-10248-1016787495.htm [33] WRIGHT W B, RUTKOWSKI A.Validation results for LEWICE 2.0: NASA/CR-1999-208690[R].Washington, D.C.: NASA, 1999. -
下载:
点击查看大图
计量
- 文章访问数: 786
- HTML全文浏览量: 75
- PDF下载量: 152
- 被引次数: 0
