Supercooled large droplet issues in aircraft icing have been continually reported due to the important safety considerations. In order to simulate the impingement behavior of large droplets, a two-dimensional and compressible Navier-Stokes code was developed to determine the flow field around the test model. Also, the Eulerian-based droplet impingement model including a semi-empirical approach for the droplet-wall interaction process and droplet break-up was developed. In particular, the droplet-wall interactions were considered as numerical boundary conditions for the droplet impingement simulation in the supercooled large droplet conditions. Finally, the present results were compared with the experimental test data and the LEWICE results. The droplet impingement area and maximum collection efficiency values between present results and wind tunnel data were in good agreements. Otherwise, the inclination of collection efficiency of the present result is over-predicted than the wind tunnel data around a lower surface of the NACA 23012 airfoil.