This paper describes the numerical optimization of an axial fan focused on the blade and guide vane (GV). For numerical analysis, three-dimensional (3D) steady-state Reynolds-averaged Navier-Stokes (RANS) equations with the shear stress transport (SST) turbulence model are discretized by the finite volume method (FVM). The objective function is enhancement of aerodynamic performance with specified total pressure. To select the design variables which have main effect to the objective function, 2k factorial design is employed as a method for design of experiment (DOE). In addition, response surface method (RSM) based on the central composite design applied to carry out the single-objective optimization. Effects on the components such as bell mouth and hub cap are considered with previous analysis. The internal flow characteristics between base and optimized model are analyzed and discussed.