Ice accretion may occur when the sold surface passes through the clouds containing supercooled water droplets. In the case of aircraft, it can result in serious performance degradation and safety hazard. In this study, numerical analysis code has been developed to predict the rime ice shapes on a 2-D airfoil and the computation results are validated against experimental data of NASA and other computation results of well-known ice prediction code, LEWICE. In addition, the effects of various numerical parameters on the ice shape have been systematically investigated.