Natural ice-accretion poses threat to aviation flight safety. When an aircraft passes through the clouds containing supercooled water droplet, ice could be accumulated according to certain physical laws. This paper introduces a three-dimensional icing simulation software AERO-ICE. The software is composed of four modules: automatic grid generation, RANS calculation of air flow field, Euler calculation of water drop field and thermodynamic analysis of icing. In the calculation of the air flow field, the SPF k-v2-ω turbulence model is adopted, which introduces the fixing term of non-equilibrium characteristics of turbulence. The model is developed to improve the prediction of non-streamlined large separation flow. Because of the large source term, the solving process of the Euler equation of the water field can be hard to converge. AERO-ICE software uses flow field smoothing, second-order MUSCL space discretization and LU-SGS implicit time marching method to improve numerical stability. In the icing thermodynamic analysis module, the AERO-ICE software has both Messinger and Myers models. The wall temperature predicted by the Messinger model is used as the boundary condition of the Myers model. This method solves the empirical problem of temperature setting in the Myers model. AERO-ICE software supports multi-block gird, multi-grid acceleration technology and large-scale parallel computing. The calculation of the ice-accretion has been preliminarily verified.
LI Haoran
,
DUAN Yuyu
,
ZHANG Yufei
,
CHEN Haixin
. Numerical method of ice-accretion software AERO-ICE[J]. ACTA AERONAUTICAET ASTRONAUTICA SINICA, 2021
, 42(S1)
: 726371
-726371
.
DOI: 10.7527/S1000-6893.2021.26371
[1] 中国民用航空规章 第25部 运输类飞机适航标准CCAR-25-R4, 2011[EB/OL].[2011-12-07]. https://www.docin.com/p-1963805781.html. CIVIL AVIATION ADMINISTRATION OF CHINA 25 AIRWORTHINESS STANDARDS FOR TRANSPORT AIRCRAFT CCAR-25-R4, 2011[EB/OL].[2011-12-07]. https://www.docin.com/p-1963805781.html(in Chinese).
[2] LI H R, ZHANG Y F, CHEN H X. Optimization design of airfoils under atmospheric icing conditions for UAV[J]. Chinese Journal of Aeronautics.doi: 10.1016/j.cja.2021.04.031.
[3] 陈维建, 张大林. 飞机机翼结冰过程的数值模拟[J]. 航空动力学报, 2005, 20(6): 1010-1017. CHEN W J, ZHANG D L. Numerical simulation of ice accretion on airfoils[J]. Journal of Aerospace Power, 2005, 20(6): 1010-1017 (in Chinese).
[4] 易贤. 飞机积冰的数值计算与积冰试验相似准则研究[D].绵阳: 中国空气动力研究与发展中心, 2007. YI X. Numerical computation of aircraft icing and study on icing test scaling law[D].Mianyang: China Aerodynamics Research and Development Center, 2007 (in Chinese).
[5] 易贤, 朱国林. 考虑传质传热效应的翼型积冰计算[J]. 空气动力学学报, 2004, 22(4): 490-493. YI X, ZHU G L. Computation of glaze ice accretion on airfoil[J]. Acta Aerodynamica Sinica, 2004, 22(4): 490-493 (in Chinese).
[6] 桑为民, 李凤蔚, 施永毅. 结冰对翼型和多段翼型绕流及气动特性影响研究[J]. 西北工业大学学报, 2005, 23(6): 729-732. SANG W M, LI F W, SHI Y Y. Icing research on airfoil and multi-element airfoil based on flow field and aerodynamic performance[J]. Journal of Northwestern Polytechnical University, 2005, 23(6): 729-732 (in Chinese).
[7] 卜雪琴, 林贵平. 基于CFD的水收集系数及防冰表面温度预测[J]. 北京航空航天大学学报, 2007, 33(10): 1182-1185. BU X Q, LIN G P. Predictions of collection efficiency and anti-icing surface temperature[J]. Journal of Beijing University of Aeronautics and Astronautics, 2007, 33(10): 1182-1185 (in Chinese).
[8] 杨倩, 常士楠, 袁修干. 水滴撞击特性的数值计算方法研究[J]. 航空学报, 2002, 23(2): 173-176. YANG Q, CHANG S N, YUAN X G. Study on numerical method for determining the droplet trajectories[J]. Acta Aeronautica et Astronautica Sinica, 2002, 23(2): 173-176 (in Chinese).
[9] 孙志国. 飞机结冰数值计算与冰风洞部件设计研究[D].南京: 南京航空航天大学, 2012. SUN Z G. Research on numerical simulation of ice accertion and design for icing research tunnel parts[D].Nanjing: Nanjing University of Aeronautics and Astronautics, 2012 (in Chinese).
[10] 陶文铨. 数值传热学[M].2版. 西安: 西安交通大学出版社, 2001: 442. TAO W Q. Numerical heat transfer[M].2nd edition. Xi'an: Xi'an Jiaotong University Press, 2001: 442 (in Chinese).
[11] BRAGG M B, BROEREN A P, BLUMENTHAL L A. Iced-airfoil aerodynamics[J]. Progress in Aerospace Sciences, 2005, 41(5): 323-362.
[12] LOPEZ M, WALTERS D K. Prediction of transitional and fully turbulent flow using an alternative to the laminar kinetic energy approach[J]. Journal of Turbulence, 2016, 17(3): 253-273.
[13] RUMSEY C. Exploring a method for improving turbulent separated-flow predictions with k-ω models[R]. Washington, D.C.: NASA, 2009.
[14] LI H R, ZHANG Y F, CHEN H X. Aerodynamic prediction of iced airfoils based on modified three-equation turbulence model[J]. AIAA Journal, 2020, 58(9): 3863-3876.
[15] LI H R, ZHANG Y F, CHEN H X. Numerical simulation of iced wing using separating shear layer fixed turbulence models[J]. AIAA Journal, 2021, 59(9): 3667-3681.
[16] BRADSHAW P, FERRISS D H, ATWELL N P. Calculation of boundary-layer development using the turbulent energy equation[J]. Journal of Fluid Mechanics, 1967, 28(3): 593.
[17] MENTER F R. Two-equation eddy-viscosity turbulence models for engineering applications[J]. AIAA Journal, 1994, 32(8): 1598-1605.
[18] LI H R, ZHANG Y F, CHEN H X. Optimization of supercritical airfoil considering the ice-accretion effects[J]. AIAA Journal, 2019, 57(11): 4650-4669.
[19] MELE B, TOGNACCINI R, CATALANO P. Performance assessment of a transonic wing-body configuration with riblets installed[J]. Journal of Aircraft, 2015, 53(1): 129-140.
[20] WRIGHT W B. User’s manual for LEWICE 3.2: NACA CR-2008-214255[R].Boston:NACA, 2008.
[21] 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.
[22] MYERS T G, CHARPIN J P F. A mathematical model for atmospheric ice accretion and water flow on a cold surface[J]. International Journal of Heat and Mass Transfer, 2004, 47(25): 5483-5500.
[23] MYERS T G. Extension to the messinger model for aircraft icing[J]. AIAA Journal, 2001, 39(2): 211-218.
[24] CAO Y H, HOU S. Extension to the myers model for calculation of three-dimensional glaze icing[J]. Journal of Aircraft, 2015, 53(1): 106-116.
[25] WRIGHT W. A summary of validation results for LEWICE 2.0[C]//37th Aerospace Sciences Meeting and Exhibit. Reston: AIAA, 1999.
CJA