Abstract: We present systematic study on droplet trajectory calculation and impingement on wings.Fan-shaped partition and ceil topology methods ate applied to water droplet trajectory calculation.The droplet momentum equation is solved with fourth Runge-Kutta method.For a three-dimensional wing,position of droplet and velocity interpolation of collection efficiency on wing surface are discussed.Three discretization methods for control equation are used and compared.They are Euler method,predictor-corroctor,and Runge-Kutta method.Impingement of droplet is analyzed and droplet collection efficiency on a NACA0012 wing is presented.Effects of droplet diameter on collection efficiency are considered.A particle trajectory module DROP3D is integrated to ice accretion software. which gives droplet impingements and ice accretion on wing surface in different icing conditions.Accuracy and reliability of the code is verified.

Key words: aircraft icing, water droplet trajectory, collection efficiency, discretization, grid partition, wing