Ice Accretion Prediction for Single Airfoil and Multi-Element Airfoil

2011 ◽  
Vol 138-139 ◽  
pp. 140-145
Author(s):  
Zhi Guo Sun ◽  
Cheng Xiang Zhu ◽  
Chun Ling Zhu

Ice accretion on aircraft components is an enormous threat to flight safety. In this paper, ice accretions on the leading edge of the NACA 0012 airfoil and the NLR 7301 multi-element airfoil with flap are predicted using the icing code developed by us. This code mainly contains five modules which are grid module, airflow module, droplet module, heat module, and boundary reconstruction module. The effectiveness and robustness of this code are tested by executing the five modules orderly and repeatedly. The Spalart-Allmaras one-equation turbulence model is adopt to calculate the viscous airflow field and the four-order Runge-Kutta method is used to solve the droplet trajectory equations. In order to enhance the efficiency of the icing calculations, the multi-block grid technique is integrated into the grid module. Based on the above methods, numerical results in both two cases are presented and the necessary comparisons with the experimental data are given in corresponding chapters. The computational results show that performance of the icing code is very good for the wide range of icing conditions.

2012 ◽  
Vol 249-250 ◽  
pp. 40-45
Author(s):  
Bo An ◽  
Wei Min Sang

Aircraft icing cause significant degradation in aerodynamics performance and flight safety. Numerical methods are developed and presented to simulate two icing-related problems for airfoils, namely ice accretion and icing effects. Ice accretion on the leading edge of the NACA 0012 airfoil is predicted using CFD method based on spring analogy. A four-order Runge-Kutta method is used to solve the droplet trajectory equation. Besides, we use the integral form of Navier-Stokes equations and the Spalart-Allmaras turbulence model to study the icing effects. Designing three different icing models, the flow fields are analyzed. The results are in good agreement with the experimental data and show preliminarily that numerical method is feasible and effective.


Aerospace ◽  
2020 ◽  
Vol 7 (9) ◽  
pp. 123
Author(s):  
Sho Uranai ◽  
Koji Fukudome ◽  
Hiroya Mamori ◽  
Naoya Fukushima ◽  
Makoto Yamamoto

Ice accretion is a phenomenon whereby super-cooled water droplets impinge and accrete on wall surfaces. It is well known that the icing may cause severe accidents via the deformation of airfoil shape and the shedding of the growing adhered ice. To prevent ice accretion, electro-thermal heaters have recently been implemented as a de- and anti-icing device for aircraft wings. In this study, an icing simulation method for a two-dimensional airfoil with a heating surface was developed by modifying the extended Messinger model. The main modification is the computation of heat transfer from the airfoil wall and the run-back water temperature achieved by the heater. A numerical simulation is conducted based on an Euler–Lagrange method: a flow field around the airfoil is computed by an Eulerian method and droplet trajectories are computed by a Lagrangian method. The wall temperature distribution was validated by experiment. The results of the numerical and practical experiments were in reasonable agreement. The ice shape and aerodynamic performance of a NACA 0012 airfoil with a heater on the leading-edge surface were computed. The heating area changed from 1% to 10% of the chord length with a four-degree angle of attack. The simulation results reveal that the lift coefficient varies significantly with the heating area: when the heating area was 1.0% of the chord length, the lift coefficient was improved by up to 15%, owing to the flow separation instigated by the ice edge; increasing the heating area, the lift coefficient deteriorated, because the suction peak on the suction surface was attenuated by the ice formed. When the heating area exceeded 4.0% of the chord length, the lift coefficient recovered by up to 4%, because the large ice near the heater vanished. In contrast, the drag coefficient gradually decreased as the heating area increased. The present simulation method using the modified extended Messinger model is more suitable for de-icing simulations of both rime and glaze ice conditions, because it reproduces the thin ice layer formed behind the heater due to the runback phenomenon.


2017 ◽  
Vol 863 ◽  
pp. 229-234
Author(s):  
Muhammad S. Virk

A multiphase numerical study has been carried out to understand the effects of wind turbine blade profile (airfoil) symmetry on resultant ice accretion. Two symmetric (NACA 0006 & 0012) and two non-symmetric airfoils (NACA 23012 & N-22) were used for this preliminary study. Based upon the airflow field calculations and super cooled water droplets collision efficiency, the rate and shape of accreted ice was simulated for rime ice conditions. Analysis showed higher air velocity along top surface of the non-symmetric airfoils as compared to symmetrical airfoils that also effects the droplet behavior and resultant ice growth. Results show that change in blade profile symmetry effects the resultant ice accretion. For symmetric airfoils, more streamlines ice shapes were observed along leading edge as compared to non- symmetric airfoils.


Author(s):  
Ranjan Saha ◽  
Jens Fridh ◽  
Torsten Fransson ◽  
Boris I. Mamaev ◽  
Mats Annerfeldt

An experimental study of the hub leading edge contouring using fillets is performed in an annular sector cascade to observe the influence of secondary flows and aerodynamic losses. The investigated vane is a three dimensional gas turbine guide vane (geometrically similar) with a mid-span aspect ratio of 0.46. The measurements are carried out on the leading edge fillet and baseline cases using pneumatic probes. Significant precautions have been taken to increase the accuracy of the measurements. The investigations are performed for a wide range of operating exit Mach numbers from 0.5 to 0.9 at a design inlet flow angle of 90°. Data presented include the loading, fields of total pressures, exit flow angles, radial flow angles, as well as profile and secondary losses. The vane has a small profile loss of approximately 2.5% and secondary loss of about 1.1%. Contour plots of vorticity distributions and velocity vectors indicate there is a small influence of the vortex-structure in endwall regions when the leading edge fillet is used. Compared to the baseline case the loss for the filleted case is lower up to 13% of span and higher from 13% to 20% of the span for a reference condition with Mach no. of 0.9. For the filleted case, there is a small increase of turning up to 15% of the span and then a small decrease up to 35% of the span. Hence, there are no significant influences on the losses and turning for the filleted case. Results lead to the conclusion that one cannot expect a noticeable effect of leading edge contouring on the aerodynamic efficiency for the investigated 1st stage vane of a modern gas turbine.


2018 ◽  
Vol 72 ◽  
pp. 01007 ◽  
Author(s):  
Faizan Afzal ◽  
Muhammad S. Virk

This paper describes a brief overview of main issues related to atmospheric ice accretion on wind turbines installed in cold climate region. Icing has significant effects on wind turbine performance particularly from aerodynamic and structural integrity perspective, as ice accumulates mainly on the leading edge of the blades that change its aerodynamic profile shape and effects its structural dynamics due to added mass effects of ice. This research aims to provide an overview and develop further understanding of the effects of atmospheric ice accretion on wind turbine blades. One of the operational challenges of the wind turbine blade operation in icing condition is also to overcome the process of ice shedding, which may happen due to vibrations or bending of the blades. Ice shedding is dangerous phenomenon, hazardous for equipment and personnel in the immediate area.


Author(s):  
Ranjan Saha ◽  
Boris I. Mamaev ◽  
Jens Fridh ◽  
Björn Laumert ◽  
Torsten H. Fransson

Experiments are conducted to investigate the effect of the pre-history in the aerodynamic performance of a three-dimensional nozzle guide vane with a hub leading edge contouring. The performance is determined with two pneumatic probes (5 hole and 3 hole) concentrating mainly on the endwall. The investigated vane is a geometrically similar gas turbine vane for the first stage with a reference exit Mach number of 0.9. Results are compared for the baseline and filleted cases for a wide range of operating exit Mach numbers from 0.5 to 0.9. The presented data includes loading distributions, loss distributions, fields of exit flow angles, velocity vector and vorticity contour, as well as, mass-averaged loss coefficients. The results show an insignificant influence of the leading edge fillet on the performance of the vane. However, the pre-history (inlet condition) affects significantly in the secondary loss. Additionally, an oil visualization technique yields information about the streamlines on the solid vane surface which allows identifying the locations of secondary flow vortices, stagnation line and saddle point.


Sign in / Sign up

Export Citation Format

Share Document