As a step toward accurate prediction of droplet impingement and ice accretion on aircraft, an Eulerian-based droplet impingement and ice accretion code for air flows around an airfoil containing water droplets is developed. A CFD solver based on the finite volume method was also developed to solve the clean airflow. The finite-volume-based approach for simulating droplet impingement on an airfoil was employed owing to its compatibility with the CFD solver and robustness. For ice accretion module, a simple model based on the control volume is combined with the droplet impingement module that provides the collection efficiency. To validate the present code, it is compared with NASA Glenn IRT (Icing Research Tunnel) experimental data and other well-known icing codes such as LEWICE and FENSAP-ICE. It is shown that the collection efficiency and shape of ice accretion are in good agreement with previous experimental and simulation results.