Collection efficiency and ice accretion calculations for a sphere, a swept MS(1)-317 wing, a swept NACA-0012 wing tip, an axisymmetric inlet, and a Boeing 737-300 inlet

Author(s):  
Colin Bidwell ◽  
Stanley Mohler, y R, Jr
1983 ◽  
Vol 4 ◽  
pp. 174-179
Author(s):  
P. McComber ◽  
J.-L. Laforte ◽  
D. Bouchard ◽  
D. D. Nguyen

There is at present a need to develop a better technique for measuring the rate of icing on structures such as, for example, overhead transmission lines. For aircraft and helicopter icing, the most widely used method of measurement is the rotating cylinder. However, for measuring the icing of structures, this method is difficult to apply and also less accurate due to lower wind velocities. Different approaches are now being developed using fixed cylinders.Icing tests were conducted with fixed and rotating cylinders in a wind tunnel. The rate of icing was obtained through measurements of volume, accretion cross-section and time of deposition. Tests were made using five different liquid water contents and droplet diameter spectra, and four cylinder diameters, keeping the wind velocity and temperature constant. The rate of icing is presented as a function of the diameters of the fixed and rotating cylinders for each of the liquid water contents tested. Results indicate that at lower wind velocities the accretion rate is overestimated for the smaller rotating cylinders. This difference is probably due to the variation of the collection efficiency with diameter. From these results it is suggested that the rate of ice accretion on structures should be based on at least two fixed cylinders of different small sizes in order to take into account the effect of the collection efficiency.


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.


Author(s):  
Wei Dong ◽  
JianJun Zhu ◽  
Rui Wang ◽  
Yong Chen

The physical processes involved in ice accretion on the rotating blade are complex. It is important to develop high fidelity numerical method and simulate the icing process on the blade under icing conditions. This paper presents a numerical study on the icing process on the rotating blade. The flow field around the blade is obtained using ANSYS FLUENT. The trajectories of supercooled water droplets and the collection efficiency are calculated by Eulerian approach. Heat and mass balance on the rotating blade surface is taken into account in icing process simulations. The NASA Rotor 67 blade is chosen as the computational model. The collection efficiency on the blade surface is computed and the impingement characteristics are analyzed. The 3D icing accretion on Rotor 67 blade is predicted at design point. The ice shapes of accretion time of 5s, 10s and 15s are simulated and the ice shapes at different span positions of the rotating blade are compared.


2012 ◽  
Vol 512-515 ◽  
pp. 754-757
Author(s):  
Xian Yi ◽  
Kai Chun Wang ◽  
Hong Lin Ma

A three dimensional numerical method and its computer codes, which are suitable to predict the process of horizontal axis wind turbine icing, are presented. The method is composed of the Multiple Reference Frame (MRF) method to calculate flowfield of air, an Eulerian method to compute collection efficiency and a three dimensional icing model companying with an iterative arithmetic for solving the model. Ice accretion on a 1.5 MW horizontal axis wind turbine is then computed with the numerical method, and characteristics of droplet collection efficiency and ice shape/type are obtained. The results show that ice on the hub and blade root is slight and it can be neglected comparing with ice near blade tip. From blade tip to root, ice becomes thinner and glaze ice may changes into rime ice.


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):  
Hao Zhang ◽  
Chihyung Wen ◽  
Junwei Su

Droplet impingement is the basic module in both ice accretion and anti-icing numerical calculation. A three dimensional finite volume approach with the capacity of modeling the in-flight droplet impingement on a wide range of subsonic regime is therefore established in this research, using OpenFOAM®. The Eulerian model is applied to estimate the droplet flow field with the same computational grid sets as those of the air flow calculation. The roughness effect caused by ice accretion is considered in the wall function modeling. Thus, the collection efficiency could be investigated for further icing numerical simulations. This approach is validated on both cylinder and sphere benchmark cases. The results are compared with the corresponding experimental and LEWICE (LEWis ICE accretion program) simulation data.


Sign in / Sign up

Export Citation Format

Share Document