The current study deals with numerical simulations of a National Advisory Committee of Aeronautics (NACA) airfoil, NACA 0012, using a commercial Computational Fluid Dynamics (CFD) Code, in air phase flow as well as in two-phase flow of air and sand particles, which consisted of 1 percent and 5 percent sand particles in air. The simulations where accomplished for Reynolds number of Re=1.76×10 6 , at various angles of attack, on Realizable k-ε turbulence model and the injection of the particles was succeeded using the Discrete Phase Model (DPM). The validation of the obtained numerical results was achieved by comparing them with reliable experimental data from other researchers and it was shown that the existence of sand particles in the air influences the aerodynamic performance of the airfoil. In particular, the predicted lift coefficient was decreased and at the same time the drag coefficient was increased. By the help of contours of static pressure and DPM concentration, it was possible to study the region around NACA 0012 airfoil.