A numerical investigation was performed to determine the effect of airfoil on the optimum flap height using NACA0015 Wells turbine. The five double flaps which have 0.5% difference were selected. A Navier-Stokes code, CFX-TASCflow, was used to calculate the flow field of the Wells turbine. The basic feature of the Wells turbine is that even though the cyclic airflow produces oscillating axial forces on the airfoil blades, the tangential force on the rotor is always in the same direction. Geometry used to define the three dimension numerical grid is based upon that of an experimental test rig. This paper tries to disign the double flap of Wells turbine with the numerical analysis.