A computational method was developed for investigating boundary layer control. Solutions of the Reynolds-averaged Navier-Stokes equations were obtained using the two-equation k-∈ turbulence model which includes the low-Reynolds-number effect in the near-wall region. Stream function and vorticity together with the turbulent kinetic energy and its dissipation rate were calculated for the flowfield in a body-fitted coordinate system. By increasing the amount of suction on the upper surface, flow separation could be totally eliminated. Transition from laminar to turbulent flow was delayed. Aerodynamic performance was substantially improved.