Numerical simulation of ice shedding motion characteristic on airfoil surface

In order to accurately predict the motion trajectory of ice shedding and ensure the safe flight of aircraft, the motion characteristics of ice shedding under two-dimensional (2D) and three-dimensional (3D) conditions are simulated and analyzed. Considering the influence of any possible shape of shedding ice and its rotation and the magnitude and direction of acceleration with time under aerodynamic force, a six degree-of-freedom analysis method is introduced in this paper. This paper proposes a theoretical model, which can be used to calculate the 3D trajectory of ice shedding with arbitrary shape. The dynamic analysis of real 3D shedding ice is carried out, and the motion behavior of shedding ice with different positions and shapes is calculated. The results show that the movement mode of the shedding ice after leaving the aircraft is translation and rotation. The shape of the low-speed region on the leeward side of the shedding ice will first increase, then decrease, and then increase with the rotation of the ice body. The influence of ice shape on ice shedding trajectory is mainly that the shedding ice continues to flip during the downstream movement of the flow field, and the projected area of the effective windward area in the lift and the drag direction changes with time. The average deviation of the shedding ice at position 5 along the spanwise is only 22.9% of that at position 1. Finally, the closer the initial position of ice shedding is to the airfoil root, the greater the probability of ice shedding hitting the aircraft fuselage. In this paper, the probability of ice shedding hitting the aircraft fuselage is 8%, which all occurred in the case with position 1 as the initial position of ice shedding.