Considering the mass and energy sources carried by the accumulated ice layer, an unsteady heat and mass transfer model of the runback water film on the deicing surface is established to simulate aircraft electro-thermal deicing process. With the extension of the freezing coefficient to the transient calculation, the coupled heat transfer of the runback water and the solid skin is solved at each time step by a temperature-based method. Unsteady numerical simulation is carried out for the electro-thermal deicing system of a NACA 0012 airfoil. The temperature variations with time are in acceptable agreement with the literature data, and the unsteady temperature-based deicing model is verified. The calculation results of temperature, runback water flux and ice thickness on the deicing surface are analyzed at different time points, and it is shown that the unsteady electro-thermal deicing model can capture the main features of the icing, ice melting and re-freezing processes in the transient deicing simulations.
Heinrich A, Ross R and Zumwalt G. Aircraft icing handbook. Vol. 1. Egg Harbor Township: Atlantic City International Airport: Department of Transportation, Federal Aviation Administration Technical Center, DOT/FAA/CT-88/8-1.
2.
CaoYTanWWuZ. Aircraft icing: an ongoing threat to aviation safety. Aerosp Sci Technol2018; 75: 353–385.
3.
ShinkafiALawsonC. Enhanced method of conceptual sizing of aircraft electro-thermal de-icing system. Int J Aerosp Mech Eng2014; 8: 1073–1080.
4.
BuXQLinGPYuJ, et al.Numerical analysis of a swept wing hot air ice protection system. Proc IMechE, Part G: J Aerospace Engineering2014; 228: 1507–1518.
5.
SilvaGALSilvaresOMZerbiniEJGJ. Numerical simulation of airfoil thermal anti-ice operation, part 2: implementation and results. J Aircr2007; 44: 634–641.
6.
BennaniLVilledieuPSalaunM, et al.Numerical simulation and modeling of ice shedding: process initiation. Comput Struct2014; 142: 15–27.
7.
PourbagianMHabashiWG. Aero-thermal optimization of in-flight electro-thermal ice protection systems in transient de-icing mode. Int J Heat Fluid Flow2015; 54: 167–182.
8.
Aircraft Engineering and Aerospace Technology. Electro-thermal ice protection system for the B-787, Bingley: Emerald Group Publishing Limited, 2007.
9.
Reid T, Baruzzi GS, Aliaga C, et al. FENSAP-ICE: application of unsteady CHT to de-icing simulations on a wing with inter-cycle ice formation. AIAA paper, no. 2010-7835, 2010.
10.
ThomasSKCassoniRPMacArthurCD. Aircraft anti-icing and de-icing techniques and modeling. J Aircr1996; 33: 841–854.
11.
Stallabrass JR. Thermal aspects of de-icer design. In: The international helicopter icing conference, Ottawa, Canada, 23–26 May 1972.
12.
RoelkeRJKeithTGDewittKJ, et al.Efficient numerical simulation of a one-dimensional electrothermal deicer pad. J Aircr1988; 25: 1097–1101.
13.
Yaslik AD, DeWitt KJ and Keith TG. Further developments in three-dimensional numerical simulation of electrothermal deicing systems. AIAA paper, no. 92-0528, 1992.
14.
HuangJRKeithTGDewittKJ. Effect of curvature in the numerical simulation of an electrothermal de-icer pad. J Aircr1995; 32: 84–91.
15.
LeiGLDongWZhengM, et al.Numerical investigation on heat transfer and melting process of ice with different porosities. Int J Heat Mass Transf2017; 107: 934–944.
16.
Wright WB, Keith TGJ and Dewitt KJ. Numerical analysis of a thermal deicer. AIAA paper, no. 92-0527, 1992.
17.
Wright WB, Al-Khalil KM and Miller DR. Validation of NASA thermal ice protection computer codes: part 2: the validation of LEWICE/Thermal. AIAA paper, no. 97-0050, 1997.
18.
MyersTGCharpinJPF. A mathematical model for atmospheric ice accretion and water flow on a cold surface. Int J Heat Mass Transf2004; 47: 5483–5500.
19.
HarirecheOVerdinPThompsonCP, et al.Explicit finite volume modeling of aircraft anti-icing and de-icing. J Aircr2008; 45: 1924–1936.
20.
ReidTBaruzziGSHabashiWG. FENSAP-ICE: unsteady conjugate heat transfer simulation of electrothermal de-icing. J Aircr2012; 49: 1101–1109.
21.
DingLChangSYangS, et al.Study on heating strategies of nose cone electrothermal anti-icing. J Aircr2018; 55: 1310–1314.
22.
MuZDLinGPShenXB, et al.Numerical simulation of unsteady conjugate heat transfer of electrothermal deicing process. Int J Aerosp Eng2018; 2018: 1–16.
Shen X, Lin G, Yang S, et al. Eulerian method for droplets impingement calculation. In: Proceedings of 48th SAFE symposium, San Diego, California, United States, 8–10 November 2011.
25.
SilvaGALSilvaresOMZerbiniEJGJ. Numerical simulation of airfoil thermal anti-ice operation, part 1: mathematical modelling. J Aircr2007; 44: 627–633.
26.
MessingerBL. Equilibrium temperature of an unheated icing surface as a function of airspeed. J Aeronaut Sci2012; 20: 29–42.
27.
BuXQLinGPYuG, et al.Numerical simulation of an airfoil electrothermal anti-icing system. Proc IMechE, Part G: J Aerospace Engineering2013; 227: 1608–1622.
28.
Domingos RH, Papadakis M and Zamora AO. Computational methodology for bleed air ice protection system parametric analysis. AIAA paper, no. 2010-7834, 2010.
29.
Wright WB. User manual for the NASA Glenn ice accretion code LEWICE version 2.2.2. NASA-CR-2002-211793, 2002. Washington, DC: NASA.
30.
Shen X, Liu X, Lin G, et al. Effects of anisotropic composite skin on electro-thermal anti-icing system. Proc IMechE, Part G: J Aerospace Engineering 2019; 0: 1–11.
