Numerical optimization of rotor blade airfoils is performed with a response surface method for helicopter rotor. For the baseline airfoils, OA 312, OA 309, and OA 407 airfoils are selected and optimized to improve aerodynamic performance. Aerodynamic coefficients required for the response surface method are obtained by using Navier-Stokes solver with k-$omega$ Shear Stress Transport turbulence model. An optimized airfoil has increased drag divergence Mach number. The present design optimization method can generate an optimized airfoil with multiple design constraints, whenever it is designed from different baseline airfoils at the same design condition.