• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.39 no.8
• /
• pp.689-701
• /
• 2011

The relations between ambient conditions and ice accretion shapes are quantitatively analyzed by employing self-organization maps and analysis of variance. Liquid water contents(LWC), mean volumetric droplet diameter(MVD), ambient temperature and free-stream velocity are chosen as ambient conditions which change ice accretion shapes. The parameters of ice accretion shape are selected as maximum thickness, icing limits, ice heading, and ice accretion area. Qualitative analysis was conducted by employing self-organization maps which show the qualitative relations between ice shapes and ambient conditions. The quantitative results of analysis of variance yield intensity of ambient conditions to the parameters of ice accretion shapes.

• 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.

• 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.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%.

• Journal of Mechanical Science and Technology
• /
• v.19 no.12
• /
• pp.2312-2320
• /
• 2005

The aerodynamic effects of leading-edge accretion can raise important safety concerns since the formulation of ice causes severe degradation in aerodynamic performance as compared with the clean airfoil. The objective of this study is to develop a numerical simulation strategy for predicting the particle trajectory around an MS-0317 airfoil in the test section of the NASA Glenn Icing Research Tunnel and to investigate the impingement characteristics of droplets on the airfoil surface. In particular, predictions of the mean velocity and turbulence diffusion using turbulent flow solver and Continuous Random Walk method were desired throughout this flow domain in order to investigate droplet dispersion. The collection efficiency distributions over the airfoil surface in simulations with different numbers of droplets, various integration time-steps and particle sizes were compared with experimental data. The large droplet impingement data indicated the trends in impingement characteristics with respect to particle size ; the maximum collection efficiency located at the upper surface near the leading edge, and the maximum value and total collection efficiency were increased as the particle size was increased. The extent of the area impinged on by particles also increased with the increment of the particle size, which is similar as compared with experimental data.