Numerical Simulation of Ice Crystal Trajectory and Its Influencing Factors Based on Lagrangian Method

Ice crystals are one of the important factors causing aircraft icing, and icing prediction is a major problem in the aviation field. The calculation of ice crystal trajectory is a precondition of the accurate simulation of ice accretion. In this paper, an ice crystal trajectory simulation method based on the Lagrangian method is proposed, and the effects of various parameters on the calculation results are discussed. Control equations of the ice crystal motion trajectory were established, and the solution methods of the equations are given. The method was used to calculate the ice crystal motion characteristics under various conditions, and the calculation results were compared with experimental results and reference results, to verify the correctness of the method. The effects of different particle densities, drag coefficient models, ice crystal diameters, inflow velocities, and angles of attack on the ice crystal collection coefficient were studied. The results showed that the smaller the particle density, the smaller the collection coefficient; the collection coefficient increases with the increase of ice crystal diameter and inflow velocity.