In this paper the incompressible subsonic flow distribution and aerodynamic characteristics of NACA 0015 airfoil in Reynolds number 2e6 is verified. Contours of pressure distribution, velocity and streamlines are shown around the airfoil. Pressure coefficients, lift coefficients and lift to drag coefficient are obtained in different angles of attack. The simulation is done by Abaqus finite element software and XFLR5. The domain of analysis includes 24620 elements and 27775 nodes. The results were repeated as the number of elements increased. It means Increasing the number of elements does not affect the accuracy of the results and only increases the analysis time. However, In order to increase the accuracy and more detailed study of the boundary layer, meshing is done more carefully near the surface. The results show that by increasing the angle of attack, pressure coefficient increases. The stall angle is 18.3 degrees in Reynolds number 2e6.