In this study, efficient global optimization (EGO) with a multi-fidelity hybrid surrogate model for multi-objective optimization is proposed to solve multi-objective real-world design problems. In the proposed approach, a design exploration is carried out assisted by surrogate models, which are constructed by adding a local deviation estimated by the kriging method and a global model approximated by a radial basis function. An expected hypervolume improvement is then computed on the basis of the model uncertainty to determine additional samples that could improve the model accuracy. In the investigation, the proposed approach is applied to two-objective and three-objective optimization test functions. Then, it is applied to aerodynamic airfoil design optimization with two objective functions, namely minimization of aerodynamic drag and maximization of airfoil thickness at the trailing edge. Finally, the proposed method is applied to aerodynamic airfoil design optimization with three objective functions, namely minimization of aerodynamic drag at cruising speed, maximization of airfoil thickness at the trialing edge and maximization of lift at low speed assuming a landing attitude. XFOILis used to investigate the low-fidelity aerodynamic force, and a Reynolds-averaged Navier-Stokes simulation is applied for high-fidelity aerodynamics in conjunction with a high-cost approach. For comparison, multi-objective optimization is carried out using a kriging model only with a high-fidelity solver (single fidelity). The design results indicate that the non-dominated solutions of the proposed method achieve greater data diversity than the optimal solutions of the kriging method. Moreover, the proposed method gives a smaller error than the kriging method.