مطالعه پارامتري اثر تجمع يخ بر عملكرد آئروديناميك ايرفويل NACA0012 با استفاده از FENSAP-ICE

Parametric Study of the Effects of Ice Accretion on the Aerodynamic Performance of NACA0012 Airfoil Using FENSAP-ICE

در اين مقاله به مطالعه پارامتري اثر تجمع يخ بر عملكرد آئروديناميك ايرفويل NACA0012، با استفاده از نرم افزار FENSAP-ICE پرداخته شده است. نتايج اين تحقيق نشان مي دهد كه پروفيل يخ توليدشده روي لبه حمله ايرفويل با زاويه صفر درجه، تقريبا به صورت متقارن، داراي دو ناحيه فرورفتگي با انحناي مشخص است. فرورفتگي هاي بالا و پايين توده يخ، باعث ايجاد جريان چرخشي در اين نواحي شده كه اين امر سبب افزايش درگ فشاري و كاهش ضريب ليفت، نسبت به حالت بدون يخ مي شود. به طوري كه در زاويه حمله 12درجه، ضريب ليفت به ميزان 20/58% كاهش و ضريب درگ، 15/92% افزايش يافته است. همچنين تاثير دما و سرعت جريان هوا بر فرآيند تجمع يخ روي ايرفويل NACA0012، در دو رژيم مات و شفاف با جريان برگشتي، مورد مطالعه قرار گرفته است. با كاهش دما از صفر تا °C14- رژيم يخ شفاف با لبه هاي تيز ايجاد شده و ضخامت يخ افزايش مي يابد. در محدوده دمايي 14- تا°C 16-، فرآيند گذار از رژيم يخ شفاف به مات اتفاق افتاده و در دماهاي كمتر از °C16- يخ مات ايجاد مي شود. در پايان نيز يك سيستم يخ زداي حرارتي به منظور رفع كامل توده يخ، شبيه سازي شده است. نتايج نشان مي دهند كه با اعمال توان حرارتي به ميزان 30وات، توسط سيستم يخ زدا، فرآيند ذوب توده يخ شاخي 21/41گرمي آغاز شده و ضمن تشكيل جريان آب بازگشتي، كل يخ روي سطح ايرفويل به صورت كامل ذوب مي شود.

In the present study, a parametric study has been carried out to investigate the influence of ice accretion on the aerodynamic performance of NACA0012 airfoil through numerical simulations using FENSAP-ICE. The results reveal that at zero angle of attack the ice profile created on the leading edge of the airfoil is symmetric. The most dominant feature in the flowfield of an iced airfoil is a recirculation zone that forms due to concavity regions created on both upper and lower surfaces of the airfoil. The numerical simulations show that the appearance of the recirculation zone alters significantly the aerodynamic coefficients. At the angle of attack 12°, lift coefficient decreases by %20.58 and the drag coefficient increases by %15.92 in comparison with the clean airfoil. The effects of temperature and air flow velocity on the ice accretion created on the NACA0012 were investigated for glaze ice and rime ice. The thickness of ice increases with decreasing temperature, and glaze ice with the sharp horn is created at the temperatures ranging from 0°C to -14°C. Making the transition from glaze ice to rime ice occurs at temperatures varied from -14°C to -16°C and at temperatures below -16°C rime ice is created. In order to eliminate the ice accretion, a thermal de-icing system is simulated. By applying a heat power of 30 watts, the melting of 21.41 gr horn ice starts and the created ice on the airfoil surface is completely melted. It should be noted that with the introduction of thermal de-icing system the runback water flow on the airfoil’s surface occurs.