The selection of the optimal position of the flap was performed in order to improve the aerodynamic performance during the take-off and landing processes of aircraft. For this, the existing airfoils of the main wing and flap are selected as the baseline model and the lift coefficients (cl) according to angle of attacks (AOA) were calculated with the change of the position of flap airfoil. The objective function was defined as the consideration of the maximum cl, lift to drag ratio and cl at certain AOA. Then, at 121 experimental points within $20mm{ imes}20mm$ domain, two dimensional flow simulations with Spalart-Allmaras turbulence model were performed concerning the AOA from 0 to 15 degree. If the optimal position was located at the domain boundary, the domain moved to the optimal position. These processes were iterated until the position was included in the inside of the domain. From these processes, the flow separation at low AOA was removed and cl increased linearly comparing with that of the baseline model.