The uncertainties such as manufacturing errors and environmental variations are inevitably encountered in engineering design, therefore in order to find robust solutions which keep high performance over a wide range, a meanvariance multi-objective robust optimization (RO) framework was proposed. Specifically, a kriging-based uncertainty quantification (UQ) formulation was proposed, which formulates uncertainty parameters and optimization design variables in product form, thereby the optimization search and UQ are conducted simultaneously in RO. Before RO, the sequential sampling strategies were employed to build high-accurate surrogate, which was shown effective by our five-variable engineering problem. Through aero-thermal design of vortex generators (VG) in U-bend channel of hightemperature blades, some issues such as constraints selection and knowledge mining in RO space were addressed. The RO of VG was carried out with one Gaussian-distributed and one uniformly-distributed uncertainty parameter, and three optimization design variables. One deterministic optimization (DO) process at the nominal condition and several RO process with different constraints were conducted. It was shown that the proposed RO was able to find robust solutions that have high-performance and are not sensitive to uncertainty fluctuations; which were validated by CFD. Meanwhile, the solution of DO was found helpful to exclude variance-dominated solutions which have much worse performance. But some constraints based on solution of DO would make the RO solution set to be null, therefore the RO can be an iterative process, in this regard, it is attractive to build high-accurate surrogates before RO. From another perspective, the several RO with different constraints are useful to group the solutions and get a insight into black-box design space like VG, which is valuable to select the final solution from many Pareto solutions.