The mathematical models and a numerical code for numerical simulation of a thermal anti-icing system are presented in this article. Mass conservation is applied to the runback water flow. An energy balance is imposed on the airfoil skin including the water flow. The heat transfer coefficient distributions are obtained using the boundary layer integral method. The external flowfield and the local water collection efficiency data are predicted using an Eulerian method, based on a computation fluid dynamic code and its user-defined functions. Given the input of the electrical power density distribution, the numerical code is able to calculate airfoil equilibrium surface temperature, mass flux of runback water, runback ice mass flux, and range if happens. A user interface is developed to integrate the computation fluid dynamic code to achieve a method for the analysis of a thermal anti-icing system. All the numerical results are compared with both experimental data and other numerical results presented in the literature.
GentRWDartNPCansdaleJT. Aircraft icing. Philos Trans R Soc London Ser A2000; 358: 2873–2911.
2.
SilvaGALSilvaresOMZerbiniEJGJ. Numerical simulation of airfoil thermal anti-ice operation part 1: mathematical modeling. J Aircr2007; 44(2): 627–634.
3.
WrightWB. User manual for the NASA Glenn ice accretion code LEWICE. 2.2.2 Final report. NASA/CR-2002-211793, 2002. Washington, DC: NASA.
4.
Gent RW. TRAJICE2-A combined water droplet trajectory and ice accretion prediction program for aerofoils. TR90054, 1990. Farnborough, UK: Royal Aircraft Establishment.
MorencyFBrahimiMTTezokF. Hot air anti-icing system modelization in the ice prediction code CANICE. AIAA paper no. 98-0192, 1998. Reston, VA: AIAA.
7.
BeaugendreHMorencyFHabashiWG. Development of a second generation in-flight icing simulation code. J Fluids Eng2006; 128(2): 378–387.
8.
MessingerBL. Equilibrium temperature of an unheated icing surface as a function of air speed. J Aeronaut Sci1953; 20(1): 29–42.
9.
WrightWBGentRWGuffondD. DRA/NASA/ONERA collaboration on icing research part 2: prediction of airfoil ice accretion. NASA CR 202349, 1997. Washington, DC: NASA.
10.
OlsenWAJShawRJNewtonJ. Ice shapes and the resulting drag increase for a NACA 0012 airfoil. NASA TM 83556, 1984. Washington, DC: NASA.
11.
Morency F, Tezok F and Paraschivoiu I. Heat and mass transfer in the case of an anti-icing system modelisation. AIAA paper no. 99-0623, 1999. Reston, VA: AIAA.
12.
Silva GAL, Silvares OM and Zerbini EJGJ. Boundary-layers integral analysis-heated airfoils in icing conditions. AIAA paper no. 2008-475, 2008. Reston, VA: AIAA.
13.
Wright WB. An evaluation of jet impingement heat transfer correlations for piccolo tube application. AIAA paper no. 2004-62, 2004. Reston, VA: AIAA.
14.
Saeed F and Al-Garni AZ. Numerical simulation of surface heat transfer from an array of hot air jets. AIAA paper no. 2007-4287, 2007. Reston, VA: AIAA.
BuXQLinGPPengYX. Method for calculation of anti-icing heat loads. Acta Aeronaut Astronaut Sin2006; 27(2): 208–212.
17.
Al-Khalil KM. Numerical simulation of an aircraft anti-icing system incorporating a rivulet model for the runback water. PhD Thesis, The University of Toledo, Toledo, OH, 1991.
18.
SilvaGALSilvaresOMZerbiniEJGJ. Numerical simulation of airfoil thermal anti-ice operation part 2: implementation and results. J Aircr2007; 44(2): 635–641.
19.
BuXQLinGP. Predictions of collection efficiency and anti-icing surface temperature. J Beijing Univ Aeronaut Astronaut2007; 33(10): 1182–1185.
20.
Silva GAL, Silvares OM and Zerbini EJGJ. Water film breakdown and rivulets formation effects on thermal anti-ice operation simulation. AIAA paper no. 2006-3785, 2006. Reston, VA: AIAA.
21.
FortinGLaforteJ-LIlincaA. Heat and mass transfer during ice accretion on aircraft wings with an improved roughness model. Int J Therm Sci2006; 45(6): 595–606.
22.
Wirogo S and Srirambhatla S. An Eulerian method to calculate the collection efficiency on two and three dimensional bodies. AIAA paper no. 2003-1073, 2003. Reston, VA: AIAA.
23.
Yang SH, Lin GP and Shen XB. An Eulerian method for water droplet impingement prediction and its implementations. In: Proceedings of the 2009 international symposium on aircraft airworthiness, Beijing, China, 2–4 November 2009, pp.72–81.
24.
IncroperaFPDeWittDP. Fundamentals of heat and mass transfer, 5th ed. New York: John Wiley & Sons, 2002.
25.
SmithAGSpaldingDB. Heat transfer in a laminar boundary layer with constant fluid properties and constant wall temperature. J Royal Aeronaut Soc1958; 62(1): 60–64.
26.
AmbrokGS. Approximate solution of equations for the thermal boundary layer with variations in boundary layer structure. Sov Phys – Tech Phys1957; 2: 1979–1986.
27.
Abu-GhannamBShawR. Natural transition of boundary layers – the effects of turbulence, pressure gradient and flow history. J Mech Eng Sci1980; 22(5): 213–228.
28.
KaysWMCrawfordME. Convective heat and mass transfer, 4th ed. New York: McGraw-Hill, 2005.
29.
Al-Khalil KM, Horvath C, Miller DR, et al. Validation of NASA thermal ice protection computer codes part3: the validation of ANTICE. NASA TM-2001-210907, Cleveland, OH, 2001. Washington, DC: NASA.
30.
Silva GAL, Silvares OM and Zerbini EJGJ, et al. Differential boundary-layer analysis and runback water flow model applied to flow around airfoils with thermal anti-ice. AIAA paper no. 2009-3967, 2009. Reston, VA: AIAA.
