Numerical study of ice crystal movement and melting in rotating blade channels

The motion and melting characteristics of ice crystals in the rotating blade channel are investigated. Firstly, the method of calculating local collection coefficient is proposed for rotating parts. Secondly, the numerical simulation of ice crystal movement and melting in the rotating blade channel is carried out to analyze the influence of ice crystal geometry parameters and working condition changes on the ice crystal impact location and ice crystal melting rate. The results show: ① the collection coefficient of ice crystal at the leading edge of the blade is the highest, the trailing edge of the pressure surface is also the area where the ice crystal may impact, while the root of the blade is less affected by centrifugal force; ② the larger the ice crystal content, the larger the collection rate of ice crystal at the same position on the blade surface; the larger the ice crystal diameter, the larger the collection rate of ice crystal at the suction surface, the smaller the collection rate of ice crystal at the pressure surface; the higher the non-spherical degree of ice crystal, the more likely the ice crystal will impact on the pressure surface; the higher the non-sphericity of ice crystals, the easier it is for the ice crystals to impact on the pressure surface; when the rotational speed decreases, the ice crystals are more likely to impact on the suction surface, and the impact area is also larger and closer to the blade root; ③ The melting ice crystals impacting on the pressure surface are more widely distributed, and the pressure surface is more prone to ice crystal adhesion and freezing than the suction surface; while, there are also melted ice crystals on the suction surface near the trailing edge, where ice crystals may also adhere; the larger the ice crystal content, the larger the diameter, the larger the non-sphericity, and when the rotational speed increases, the melting rate of ice crystals decreases.