Aerodynamic design optimization of rotor airfoil has been performed with advanced design method for improved aerodynamic characteristics of ONERA airfoils. A multiple response surface method is used to consider various requirements in rotor airfoil design. Shape functions for mean camber line are proposed to extend possible design domain. Numerical simulations are performed using KFLOW, a Navier-Stokes solver with shear stress transport turbulence model. The present design method provides favorable configurations for the high performance rotor airfoil. Resulting optimized airfoils give better aerodynamic performance than the baseline airfoils.