Ice accretion on aircraft surface plays a critical role in the performance and safety of aircraft. Especially, changes in the external shape of the aircraft due to ice accretion can greatly deteriorate the aerodynamic performance of aircraft. In the study, air flow and droplet impingement codes for prediction of ice accretion on aircraft are developed. A finite volume method CFD solver based on unstructured grids is developed to solve the dean air flow and Eulerian-based droplet field. The codes are then compared with NASA IRT experimental data and computational results obtained by a state-of-the-art icing code. The present prediction is shown to be in close agreement with IRT data.