Bimodal SLD Ice Accretion on Swept NACA 0012 Airfoil Models

2020 ◽  
Author(s):  
Mark G. Potapczuk ◽  
Jen-Ching Tsao
Keyword(s):  
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.


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.


Author(s):  
Z. Chara ◽  
V. Horak ◽  
D. Rozehnal

The phenomenon of in-flight icing may affect all types of aircraft. Presence of ice on wings can lead to a number of aerodynamic degradation problems. Thus, it is important to understand the different ice shapes that can form on the wings and how they affect aerodynamics. When compared to wings without ice, wings with ice indicate decreased maximum lift, increased drag, changes in pressure distribution, stall occurring at much lower angles of attack, increased stall speed, and reduced controllability. The in-house ice accretion prediction code R-ICE using 2-D panel method was developed. The CFD simulation with the software ANSYS CFX 11.0 was used to simulate flow around iced airfoils NACA 0012. These airfoils were experimentally investigated in a wind tunnel. The paper presents a comparison of lift and drag coefficients experimentally observed and numerically simulated.


Author(s):  
Mark G. Potapczuk ◽  
Jen-Ching Tsao ◽  
Laura-Cheri King-Steen
Keyword(s):  

2020 ◽  
Vol 34 (14n16) ◽  
pp. 2040068
Author(s):  
Tong Liu ◽  
Jin-Sheng Cai ◽  
Kun Qu ◽  
Shu-Cheng Pan

This paper presents a comprehensive aircraft icing simulation tool implemented in an in-house Navier–Stokes parallel multi-block solver. In detail, the droplet flow field is solved by Eulerian approach, and a Partial Differential Equation (PDE)-based ice accretion model is adopted to determine the runback water flow and icing rate. Numerical validations are performed on the two-dimensional (2D) NACA 0012 airfoil, where good agreements with the literature are observed. Additionally, the paper investigates the influence of droplet size on the final ice shape. Results show that droplets with greater Median Volume Diameter (MVD) are more likely to impact on the wall, which results in larger droplet impingement limit and icing limit.


Sign in / Sign up

Export Citation Format

Share Document