Ice Accretin Predictin Cde Vladimír Hrák University f Defence Institute f Hydrdynamics, Academy f Sciences Czech Republic 1
Dedicated t the memry f Dr. Bhumír Hření 2
Mtivatin The in-flight icing may affect all types f aircraft. Presence f ice n an aircraft surface can lead t a number f perfrmance degradatins: changes in pressure distributin decreased maximum lift and increased drag stall ccurring at lwer angles f attack and increased stall speed reduced cntrllability It is imprtant t t understand d hw the different ice shapes affect aircraft aerdynamics. It can be studied by flight tests, wind tunnel measurements, and cmputatinal simulatins. 3
Cmputatinal simulatin f ice accretin is an essential tl in design, develpment and certificatin f aircraft fr flight int icing cnditins. Currently, there exist several apprved ice accretin cdes: LEWICE (LEWis ICE accretin prgram) is sftware develped by the Icing Branch at NASA Glenn Research Center CANICE cde develped at the Ecle Plytechniue de Mntreal ONERA (Office Natinal d'etudes et de Recherches Aérspatiales) cde in France TRAJICE cde which was develped by DERA (Defence Evaluatin and Research Agency) in United Kingdm CIRA cde frm Italian Aerspace Research Center. 4
Ice Accretin Predictin Cde Cde was develped as a tl fr simulating flight int icing cnditins Presented sftware was subseuently develped and imprved. There are three main cde versins: R-Ice 1.1 Rime ice accretin predictin Ice 3.1 Glaze ice accretin predictin Ice 4.1 Multi-element airfils icing 5
Trajectry f water drplets The ptential flw field is calculated using 2-D panel methd. The relatin fr any pint inside the cntrl area is in frm Ptential flw field is then used t determine the trajectries f water drplets and the impingement pints n the bdy. Drplets passing thrugh the atmsphere are cnsidered as spherical elements n that the surrunding fluid frces and gravitatin act. Small water drplets have trajectries similar t streamlines, vice versa large water drplets trajectries are affected by the airfil inherency nly slightly. 6
Versin 1.1 Airfil rime ice accretin predictin Impinging super-cled water drplets freeze immediately upn impact. Cde applies a time-stepping prcedure t calculate the shape f an ice accretin. 0,050 0,050 y/c [1] y/c [1] 0,025 0,025 0,000 Chrd=0.9144 m; Vext=92.54 m.s -1 ; FluidP=100000 Pa; FluidT=257.6 K; Alpha=0 ; PartD=0.00002 m; PartCntent=0.00033 kg.m -3 ; RimeIceRh=900 kg.m -3 ; Steps=5; TimeStep=64.8 s; Time=324 s. 0,000 Chrd=0.9144 m; Vext=92.54 m.s -1 ; FluidP=100000 Pa; FluidT=257.6 K; Alpha=0 ; PartD=0.00002 m; PartCntent=0.00033 kg.m -3, RimeIceRh=900kg.m -3 ; Steps=5; TimeStep=244.8 s; Time=1224 s. -0,025-0,025-0,050-0,025 0,000 0,025 0,050 x/c [1] -0,050-0,050-0,025 0,000 0,025 0,050 x/c [1] 7
Quantitative cmparisn f using the current cmputatinal iceaccretin simulatin methds. 0,050 y/c [1] 0,025 0,000-0,025 Clean Airfil Experimental ONERA1990 Simn CANICE Paraschiviu TRAJICE ADSE ONERA2000 Duprat NASA R-ICE Generally, current ice accretin cdes give satisfied results f fthe rime ice simulatin. -0,050-0,025 0,000 0,025 0,050 x/c [1] 0,075 0,050 y/c [1] 0,025 0,000-0,025 Clean Airfil Experimental ONERA1990 Simn CANICE Paraschiviu TRAJICE ADSE ONERA2000 Duprat NASA R-ICE -0,050-0,050-0,025 0,000 0,025 x/c [1] 0,050 8
Versin 3.1 Airfil glaze ice accretin predictin Glaze ice creates at cmbinatins f temperature clse t freezing. Evapratin In that case, nt all f the impinging ev p water freezes n impact. Liuid water Thin layer f water is flwing very fr slwly alng the surface and freeze at ther lcatins. Ice accretin Bundary layer There are used theretical appraches are generally called as a shallw water thery. The cnservative euatins using fr the slutin f water flw in pen channels are frmally arranged. The flux terms are evaluated using a discntinuus Galerkin methd based n finite-vlume frmulatin. 9
Mass, mtin, and energy cnservative euatins culd be written in the general frm x t S S F Q general frm Q A Vectrs f variables Q, flw F and surces S, S are given by relatins: flw crss-sectin mass flux A EQ I g A Q Q E Q A n 1 2, F Q flw crss-sectin flw vlume thermal energy mass flux mmentum flux flux f energy A EQ E A dx dp A g I g t 2 0 S Internal T T T T A dx dp A g I g e e e w w w e e w w t n 2 S surces: f ev w fr w p e v v v S External surces: 10 ev ev w fr fr w p p e ev w fr w px p e L L ct ct ct v v v S surces:
Example f glaze ice-accretin simulatin. Airfil NAC 0012, chrd 0.45 m, angle f attack α = 0, velcity v =772ms -1 77.2, MVD = 18 μm, LWC = 0.32 g m -3, air temperature T = 270.5 K, icing expsitin 300 s. 006 0,06 y/b 0,04 0,02 0,00-0,02-0,04-0,06-0,02 0,00 0,02 0,04 0,06 0,08 x/b 0,10 Clean Airfil Experimental Paraschiviu CANICE ADSE TRAJICE Duprat ONERA NASA LEWICE Cmparisn f current ice-accretin cdes shws there is still rm fr imprvement in the uality f predictins. 11
Air temperature influence n iced airfil shapes Stream-wise shape (b), (c) Duble-hrn shape (d), (e) Span-wise ridge shape (f) (a) Rime ice (b) T = 269.6565 K (c) T = 270.15 K (d) T = 270.65 K (e) T = 271.65 K (f) T = 272.65 K Airfil NFL0414, chrd 0.45 m, angle f attack α = 0, velcity v = 77.2 m s -1, MVD = 18 μm, LWC = 0.32 g m -3, icing expsitin 900 s. 12
Versin 4.1 Multi-element airfils icing The latest cde versin enables slutin f multi-element airfils up t eight separate parts. Mutual flw verlap f circumfluent bdies ccurs. Drplet trajectries near an airfil with a sltted flap Drplet trajectries near an airfil with a sltted flap in landing psitin 13
Example f flapped airfil icing Ice accretin n the flap causes the reductin f the gap size between main element and flap. Cnseuently, it can have a large impact n the perfrmance degradatin f iced multi-element airfils. There is a ptential mechanical prblem in the elevatr mechanism itself. 14
Input file Cnfiguratin file: _ice.cfg Icing parameters # cnfiguratin file # ------------------ PrfDta=NACA0018 # NACA0018 Chrd=1 # 1.0;<0.1,10> Vext=50 # 50.0;<10,200> Alpha=5 # 5.0;<-20,20> FluidP=100000 # 1.e5;<0.1e5,1.2e5> 1e5 1 FluidT=263 # 263;<200,300> PartCntent=0.001 # 1.e-3;<0.05e-3,5.e-3> PartD=0.0001 # 100e-6;<5.e-6,5000.e-6> PartT=274.15 # 274.15;<200,300> BdyT=263.15 # 263.15;<200,300> RimeIceRh=900 # 900;<700,1000> GlaseIceRh=917 # 917;<700,1000> # ------------------ TimeStep=120 # 120;<1,600> Steps=5 # 5;<1,1000> MinBdyPints=50 # 50;<20,500> 15
Graphical utput files Airfil gemetry gem004.eps Successive ice accretin by steps mdif004.eps Drplet trajectries near iced airfil part004.eps p 16
beta004.eps Water drplets lcal cllectin efficiency β, nrmal v n and tangential v t velcity with regard t the airfil surface at impact lcatin 17
blay004.eps Parameters f bundary layer: Lcal frictin cefficient c f, displacement bundary layer thickness δ 1 (mass) and δ 2 (impulse), and Nusselt number Nu 18
flw004.eps Distributin f pressure cefficient c p and relative velcity v/v 19
flw005.dta Text utput file 1.000000 Chrd [m] 2.087219 Length [m] 5.000000 Alpha [st] 50.000000 Vext [m/s] 1.000000e+005 FluidP [Pa] 263.000000 FluidT [K] 1324643 1.324643 FluidRh [kg/m^3] 1.246121e-005 FluidNu [m^2/s] 1.000000e-004 PartD [m] 274.150000 PartT [K] 1.000000e-003 PartCntent [kg/m^3] 1000.000000 PartRh [kg/m^3] 263.150000 BdyT [K] # s[m] x[m] y[m] v[m/s] beta[-] vn[m/s] vt[m/s] cp[-] cf[-] Nu[-] -1.0395e+000 +9.9210e-001-8.7863e-002-3.2739e+001 0.0000e+000 +0.0000e+000 +0.0000e+000 +5.7127e-001 +0.0000e+000 6.7548e+002-1.0353e+000 +9.8800e-001-8.8559e-002-3.2634e+001 0.0000e+000 +0.0000e+000 +0.0000e+000 +5.7400e-001 +0.0000e+000 6.7548e+002-1.0311e+000 +9.8390e-001-8.9234e-002-3.4098e+001 0.0000e+000 +0.0000e+000 +0.0000e+000 +5.3493e-001 +0.0000e+000 6.7548e+002-1.0270e+000 +9.7979e-001-8.9889e-002-3.6759e+001 0.0000e+000 +0.0000e+000 +0.0000e+000 +4.5952e-001 +0.0000e+000 6.7548e+002-1.0228e+000 +9.7568e-001-9.0526e-002-3.8835e+001 0.0000e+000 +0.0000e+000 +0.0000e+000 +3.9672e-001 +0.0000e+000 6.7548e+002-1.0187e+000 +9.7157e-001-9.1144e-002-4.0328e+001 0.0000e+000 +0.0000e+000 +0.0000e+000 +3.4946e-001 +0.0000e+000 6.7548e+002...... +1.0187e+000 +000 +9.7264e-001-7.8940e-002 78940 002 +4.1296e+001 +001 0.0000e+0000000 +0.0000e+000 0000 +000 +0.0000e+000 0000 +000+31785 +3.1785e-001 +0.0000e+000 0000 +000 3.2149e+003 +1.0228e+000 +9.7658e-001-8.0263e-002 +3.9687e+001 0.0000e+000 +0.0000e+000 +0.0000e+000 +3.6997e-001 +0.0000e+000 3.2149e+003 +1.0270e+000 +9.8052e-001-8.1603e-002 +3.7486e+001 0.0000e+000 +0.0000e+000 +0.0000e+000 +4.3793e-001 +0.0000e+000 3.2149e+003 +1.0311e+000 +9.8445e-001-8.2961e-002 +3.4692e+001 0.0000e+000 +0.0000e+000 +0.0000e+000 +5.1859e-001 +0.0000e+000 3.2149e+003 +1.0353e+000 +9.8837e-001-8.4338e-002 +3.3053e+001 0.0000e+000 +0.0000e+000 +0.0000e+000 +5.6299e-001 +0.0000e+000 3.2149e+003 +1.0395e+000 +9.9229e-001 9229e-001-8.5735e-002 +3.2948e+001 0.0000e+000 0000e+000 +0.0000e+000 0000e+000 +0.0000e+000 0000e+000 +5.6577e-001 +0.0000e+000 0000e+000 3.2149e+003 #end 20
Clsing remarks ICE cde enables cmputatinal rime ice and glaze ice accretin predictin n single and multi-element airfils in acceptable time f slutin. Mathematical ti mdel has recently been mdified d fr variable wall temperature alng the airfil surface. The cde was als imprved fr the better apprximatin f transitin bundary layer lcatin. Presented cde culd be cnsidered at least as a fully cmparable with the current ice accretin predictin cdes. 21
Thank yu fr yur attentin 22