• Wind and Structures
• /
• v.32 no.3
• /
• pp.205-217
• /
• 2021

Cold regions with high air density and wind speed attract wind energy producers across the globe exhibiting its potential for wind exploitation. However, exposure of wind turbine blades to such cold conditions bring about devastating impacts like aerodynamic degradation, production loss and blade failures etc. A series of wind tunnel tests were performed to investigate the effect of icing on the aerodynamic properties of wind turbine blades. A baseline clean wing configuration along with four different ice accretion geometries were considered in this study. Aerodynamic force coefficients were obtained from the surface pressure measurements made over the test model using MPS4264 Simultaneous pressure scanner. 3D printed Ice templates featuring different ice geometries based on Icing Research Tunnel data is utilized. Aerodynamic characteristics of both the clean wing configuration and Ice accreted geometries were analysed over a wide range of angles of attack (α) ranging from 0° to 24° with an increment of 3° for three different Reynolds number in the order of 105. Results show a decrease in aerodynamic characteristics of the iced aerofoil when compared against the baseline clean wing configuration. The key flow field features such as point of separation, reattachment and formation of Laminar Separation Bubble (LSB) for different icing geometries and its influence on the aerodynamic characteristics are addressed. Additionally, attempts were made to understand the influence of Reynolds number on the iced-aerofoil aerodynamics.

• Journal of computational fluids engineering
• /
• v.18 no.2
• /
• pp.41-48
• /
• 2013

Ice accretion on the solid surface is an importance factor in assessing the performance of aircraft and wind turbine blade. Changes in the external shape due to ice accretion can greatly deteriorate the aerodynamic performance. In this study, a three-dimensional upwind-type second-order positivity-preserving finite volume CFD scheme based on the unstructured mesh topology is developed to simulate two-phase flow in atmospheric icing condition. The code is then validated by comparing with NASA IRT experimental data on the sphere. The present results of the collection efficiency are found to be in close agreement with experimental data and show improvement near the stagnation region.

• Journal of computational fluids engineering
• /
• v.20 no.2
• /
• pp.23-36
• /
• 2015

The $2^{nd}$ generation ice accretion analysis program has been developed and validated for various icing conditions. The essential feature of the $2^{nd}$ generation code lies in its capability of handling general 3-D geometry and improved accuracy. The entire velocity fields are obtained based on Navier-Stokes equations in order to take the massively separated flow field into account. Unlike $1^{st}$ generation code, the droplet trajectories are calculated using Eulerian approach, which is adopted to yield appropriate collection efficiency even in the shadow region. For improved thermodynamic analysis on the surfaces, water film model and modified Messinger model are newly included in the present analysis. The ice shape for a given time step is obtained by considering the exact amount of ice accreted on the surface. Each module of the icing analysis code has been seamlessly integrated on the OpenFOAM platform. The developed code was validated against available experimental data for 2D airfoils and 3D DLR-F4. Due to the lack of experimental data, the computed results of DLR-F4 were compared with those obtained from FENSAP-ICE, which is state-of-the-art 3D icing analysis code. It was clearly shown that the present code produces comparable results to those of FENSAP-ICE, in terms of prediction accuracy and the capability of handling general 3-D geometries.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.42 no.2
• /
• pp.134-143
• /
• 2014

A significant degradation in the aerodynamic performance of wind turbine system can occur by ice accretion on the surface of blades operated in cold climate. The ice accretion can result in performance loss, overloading due to delayed stall, excessive vibration associated with mass imbalance, ice shedding, instrumental measurement errors, and, in worst case, wind turbine system shutdown. In this study, the effects of ice accretions on the aerodynamic characteristics of wind turbine blade sections are investigated on the basis of modern CFD method. In addition, the computational results are used to predict the performance of three-dimensional wind turbine blade system through the blade element momentum method. It is shown that the thickness of ice accretion increases from the root to the tip and the effects of icing conditions such as relative wind velocity play significant role in the shape of ice accretion.

• International Journal of Aeronautical and Space Sciences
• /
• v.18 no.2
• /
• pp.352-364
• /
• 2017

A novel 3D mathematical model for water film runback and icing on a rotating surface is established in this work, where both inertial forces caused by the rotation and shear forces due to the air flow are taken into account. The mathematical model of the water film runback and energy conservation of phase transition process is established, with a cyclical average method applied to simulate the unsteady parameters variation at angles of attack. Ice accretion on a conical spinner surface is simulated and the results are compared with the experimental data to validate the presented model. Then Ice accretion on a cowling surface is numerically investigated. Results show that a higher temperature would correspond to a larger runback ice area and thinner ice layer for glaze ice. Rotation would enhance the icing process, while it would not significantly affect the droplet collection efficiency for an axi-symmetric surface. In the case at angle of attack, the effect of rotation on ice shape is appreciable, ice would present a symmetric shape, while in a stationary case the shape is asymmetric.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.8
• /
• pp.747-757
• /
• 2010

The ice accreted on the airfoil is one of the critical drivers that causes the degradation of aerodynamic performance as well as aircraft accidents. Hence, an efficient numerical code to predict the accreted ice shape is crucial for the successful design of de-icing and anti-icing devices. To this end, a numerical code has been developed for the prediction of glaze ice accretion shape on 2D airfoil. Constant Source-Doublet method is used for the purpose of computational efficiency and heat transfer in the icing process is accounted for by Messinger model. The computational results are thoroughly compared against available experiments and other computation codes such as LEWICE and TRAJICE. The direction and thickness of ice horn are shown to yield similar results compared to the experiments and other codes. In addition, the effects of various parameters - temperature, free-stream velocity, liquid water contents, and droplet diameter - on the ice shape are systematically analyzed through parametric studies.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.11
• /
• pp.1105-1114
• /
• 2010

Ice shapes accumulated on the aircraft surfaces are categorized into rime and glaze ice, which are highly dependent on various parameters such as ambient temperature, liquid water contents (LWC), mean volumetric droplet diameter and freestream velocity. In this study, quantitative analyses on the ice accretion have been attempted in a systematical manner and the key findings are as follows. First, the increase of freestream velocity can cause tremendous change in the ice accumulation such as the growth of ice accretion area, ice heading direction and maximum thickness of ice horn. Second, LWC is found to be linearly proportional to the ice accretion area. Third, the effects of ambient temperature on incoming water mass seem to be relatively small in comparison with LWC and freestream velocity. Finally, it was shown that MVD has only a little influence on ice shapes. However, it may increase the ice accretion area by increasing the droplet impacting range.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.11
• /
• pp.1094-1104
• /
• 2010

In the previous studies, the validation of numerical codes has been conducted based on the qualitative comparison of predicted ice shapes with experiments, which poses a significant limit on the systematic analysis of ice shapes due to the variation of meteorological conditions. In response to this, the numerical code has been quantitatively validated against available experiment for the ice accretion on cylinders and airfoils in the present study. Ice shapes accumulated on the bodies are systematically investigated with respect to various icing parameters. To this end, maximum thickness, heading direction and ice thickness are quantified and expressed in the polar coordinate system for the comparison with other numerical results. By applying the quantitative analysis, similar shapes are intuitively distinguished. The developed numerical code underestimates the ice accretion area and the ice thickness of lower surface. In order to improve the accuracy, further accurate aerodynamic solver is required for the water droplet trajectories.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.6
• /
• pp.530-536
• /
• 2010

In-flight icing is a critical technical issue for aircraft safety and, in particular, ice accretions on aircraft surfaces can drastically impair aerodynamic performances and control authority. In order to investigate icing effects on the aerodynamic characteristics of KC-100 aircraft, a state-of-the-art CFD code, FENSAP-ICE, was used. A main wing section and full configuration of KC-100 aircraft were considered for the icing analysis. Also, shapes of iced area were calculated for the design of anti-/de-icing devices. The iced areas around leading edge of main wing and horizontal tail wing were observed maximum 7.07% and 11.2% of the chord length of wing section, respectively. In case of wind shield, 16.7% of its area turned out to be covered by ice. The lift of KC-100 aircraft were decreased to 64.3%, while the drag was increased to 55.2%.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.27-30
• /
• 2009

The aerodynamic performance of aircraft in icing condition can deteriorate considerably by contamination of aerodynamic and propulsive systems due to icing accretions on aircraft surfaces. A computational analysis based on the Eulerian description was performed on an airfoil to investigate effects of ice accretions on airfoil aerodynamics. A water droplet with liquid water concentration ($0.00075kg/m^3$) and mean volume diameter ($20{\mu}m$) was considered and applied to various angles of attack to investigate the stall angle decrease and the drag increment.