In this study, an unsteady numerical simulation method based on dynamic grids was established to solve unsteady aircraft icing processes. The dual-time step method was used to achieve the time advancement of the transient airflow and droplets fields. The droplets trajectory was calculated using the Euler method, which incorporates a supercooled large droplets model. The transient evolution of the water film and ice layer on the substrate surface was considered in the icing thermodynamic model. The accuracy of the unsteady method for ice-shape prediction was verified via comparison and analysis of numerical examples. Additionally, the rime and glaze ice formation processes were explored, and the coupling effects of ice horn growth and the air convective heat transfer coefficient, water collection coefficient, and water film flow on the substrate surface were studied. The ice shapes calculated by unsteady and quasi-steady methods were compared. The results indicated that the latter predicts a larger icing range than the former. Moreover, there may be more obvious distinctions between these two methods in predicting glaze ice.
ColeJSandW. Statistical study of aircraft icing accidents. In: Proceedings of the 29th Aerospace Sciences Meeting, Reno,NV,U.S.A, 07 January 1991 - 10 January 1991.
2.
AshendenRMarwitzJ. Characterizing the supercooled large droplet environment with corresponding turboprop aircraft response. J Aircraft1998; 35: 912–920. DOI: 10.2514/2.2386.
3.
KimHBraggM. Effects of leading-edge ice accretion geometry on airfoil performance. In: Proceedings of the 17th Applied Aerodynamics Conference, Norfolk,VA,U.S.A, 28 June 1999 - 01 July 1999.
4.
KreegerREBroerenAP. Collaborative experiments and computations in aircraft icing. In: Proceedings of the 2018 Applied Aerodynamics Conference, Atlanta, Georgia, June 25-29, 2018.
5.
WrightW. A summary of validation results for LEWICE 2.0. In: Proceedings of the 37th Aerospace Sciences Meeting and Exhibit, Reno,NV,U.S.A, 11 January 1999 - 14 January 1999.
6.
GhenaiCLinCX. Verification and validation of NASA LEWICE 2.2 icing software code. J Aircraft2006; 43: 1253–1258. DOI: 10.2514/1.18977.
BaumertABansmerSTrontinP, et al.Experimental and numerical investigations on aircraft icing at mixed phase conditions. Int J Heat Mass Tran2018; 123: 957–978. DOI: 10.1016/j.ijheatmasstransfer.2018.02.008.
10.
PapadakisMRachmanAWongS-C, et al.Water impingement experiments on a NACA 23012 airfoil with simulated glaze ice shapes. In: Proceedings of the 42nd AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada, 05 January 2004 - 08 January 2004.
11.
BourgaultYHabashiWGDompierreJ, et al.A finite element method study of Eulerian droplets impingement models. Int J Numer Methods Fluid1999; 29: 429–449.
12.
MessingerBL. Equilibrium temperature of an unheated icing surface as a function of air speed. J Aeronaut Sci1953; 20: 29–42. DOI: 10.2514/8.2520.
13.
WangZZhuC. Numerical simulation for in-cloud icing of three-dimensional wind turbine blades. Simulation2017; 94: 31–41. DOI: 10.1177/0037549717712039.
14.
CaoYMaCZhangQ, et al.Numerical simulation of ice accretions on an aircraft wing. Aero Sci Technol2012; 23: 296–304. DOI: 10.1016/j.ast.2011.08.004.
15.
ChenXZhaoQ-J. Computational investigations of water collection efficiency on blades in unsteady vortex flowfield of rotor. Aero Sci Technol2018; 79: 482–491. DOI: 10.1016/j.ast.2018.06.016.
16.
BourgaultYBeaugendreHHabashiWG. Development of a shallow-water icing model in FENSAP-ICE. J Aircraft2000; 37: 640–646. DOI: 10.2514/2.2646.
17.
MyersTG. Extension to the messinger model for aircraft icing. AIAA J2001; 39: 211–218. DOI: 10.2514/2.1312.
18.
MyersTGCharpinJPF. A mathematical model for atmospheric ice accretion and water flow on a cold surface. Int J Heat Mass Tran2004; 47: 5483–5500. DOI: 10.1016/j.ijheatmasstransfer.2004.06.037.
19.
YanxiaDYeweiGChunhuaX, et al.Investigation on heat transfer characteristics of aircraft icing including runback water. Int J Heat Mass Tran2010; 53: 3702–3707. DOI: 10.1016/j.ijheatmasstransfer.2010.04.021.
20.
CaoYHouS. Extension to the Myers model for calculation of three-dimensional glaze icing. J Aircraft2016; 53: 106–116. DOI: 10.2514/1.C033212.
ShinJBondT. Results of an icing test on a NACA 0012 airfoil in the NASA Lewis icing research Tunnel. In: Proceedings of the 30th Aerospace Sciences Meeting and Exhibit, Reno,NV,U.S.A., 06 January 1992 - 09 January 1992.
23.
LiLLiuYZhangZ, et al.Effects of thermal conductivity of airframe substrate on the dynamic ice accretion process pertinent to UAS inflight icing phenomena. Int J Heat Mass Tran2019; 131: 1184–1195. DOI: 10.1016/j.ijheatmasstransfer.2018.11.132.
24.
GoriGParmaGZoccaM, et al.Local solution to the unsteady stefan problem for in-flight ice accretion modeling. J Aircraft2018; 55: 251–262. DOI: 10.2514/1.C034412.
25.
LiuTQuKCaiJ, et al.A three-dimensional aircraft ice accretion model based on the numerical solution of the unsteady Stefan problem. Aero Sci Technol2019; 93: 105328. DOI: 10.1016/j.ast.2019.105328.
26.
XiaoMZhangYZhouF. Numerical investigation of the unsteady flow past an iced multi-element airfoil. AIAA J2020; 58: 3848–3862. DOI: 10.2514/1.J059114.
27.
VerdinPGCharpinJPF. Multi-stepping ice prediction on cylinders and other relevant geometries. J Aircraft2013; 50: 871–878. DOI: 10.2514/1.C032033.
DonizettiABellostaTRausaA, et al.Level-set mass-conservative front-tracking technique for multistep simulations of in-flight ice accretion. J Aircraft2023; 60: 1050–1060. DOI: 10.2514/1.C037027.
30.
ZhuCZhaoHZhaoN, et al.An adaptive Cartesian method for prediction of the unsteady process of aircraft ice accretion. Nanjing Univ Aeronaut & Astronaut, Key Lab Fundamental Sci Natl Def Adv Design Techn; 30: 515–535. DOI: 10.4208/cicp.OA-2018-0228.
31.
WangYQinNZhaoN. Delaunay graph and radial basis function for fast quality mesh deformation. J Comput Phys2015; 294: 149–172. DOI: 10.1016/j.jcp.2015.03.046.
32.
KnoppTEisfeldBCalvoJB. A new extension for k– turbulence models to account for wall roughness. Int J Heat Fluid Flow2009; 30: 54–65. DOI: 10.1016/j.ijheatfluidflow.2008.09.009.
33.
ShinJBerkowitzBChenH, et al.Prediction of ice shapes and their effect on airfoil performance. In: Proceedings of the 29th Aerospace Sciences Meeting, Reno,NV,U.S.A., 07 January 1991 - 10 January 1991.
34.
BidwellCSPotapzcukMG. Users manual for the NASA Lewis three-dimensional ice accretion code (LEWICE 3D). Washington, D.C.: Icing, 1993.
35.
FornasierS. Ice accretion simulation over 2D airfoils by means of a methodology embedded in a CFD solver. M.Sc. thesis, Università degli Studi di Padova. 2015.
36.
ToroEFSpruceMSpearesW. Restoration of the contact surface in the HLL-Riemann solver. Shock Waves1994; 4: 25–34.
37.
BattenPClarkeNLambertC, et al.On the choice of wavespeeds for the HLLC riemann solver. SIAM J Sci Comput1997; 18: 1553–1570.
38.
SeokkwanYAntonyJ. Lower-upper symmetric-gauss-seidel method for the euler and Navier-Stokes equations. AIAA J1988; 26(9): 10007.
39.
IulianoEMingioneGPetrosinoF, et al.Eulerian modeling of SLD physics towards more realistic aircraft icing simulation. In: AIAA Atmospheric and Space Environments Conference, Toronto, Ontario, Canada, 02 August 2010 - 05 August 2010.
WrightWBRutkowskiA. Validation results for LEWICE 2.0. NASA CR-208690. 1999. DOI: 10.2514/6.1999-249.
42.
WrightWPotapczukM. Semi-empirical modelling of SLD physics. In: Proceedings of the 42nd AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada, 05 January 2004 - 08 January 2004.
43.
HanYPalaciosJ. Validation of a LEWICE-based icing code with coupled heat transfer prediction and aerodynamics performance determination. In: Proceedings of the 9th AIAA Atmospheric and Space Environments Conference, Denver, Colorado, 5-9 June 2017.
WirogoSSrirambhatlaS. An eulerian method to calculate the collection efficiency on two and three dimensional bodies. In: Proceedings of the 41st Aerospace Sciences Meeting and Exhibit, Reno, Nevada, 06 January 2003 - 09 January 2003.
