PurposeThis paper aims to present a multi-fidelity surrogate-based optimization (MFSBO) method for computationally accurate and efficient design of microfluidic concentration gradient generators (µCGGs).Design/methodology/approachCokriging-based multi-fidelity surrogate model (MFSM) is constructed to combine data with varying fidelities and computational costs to accelerate the optimization process and improve design accuracy. An adaptive sampling approach based on parallel infill of multiple low-fidelity (LF) samples without notably adding computation burden is developed. The proposed optimization framework is compared with a surrogate-based optimization (SBO) method that relies on data from a single source, and a conventional multi-fidelity adaptive sampling and optimization method in terms of the convergence rate and design accuracy.FindingsThe results demonstrate that proposed MFSBO method allows faster convergence and better designs than SBO for all case studies with 49% more reduction in the objective function value on average. It is also found that parallel infill (MFSBO-4) with four LF samples, enables more robust, efficient and accurate designs than conventional multi-fidelity infill (MFSBO-1) that only adopts one LF sample during each iteration for more complex optimization problems.Originality/valueA MFSM based on cokriging method is constructed to utilize data with varying fidelities, accuracies and computational costs for µCGG design. A parallel infill strategy based on multiple infill criteria is developed to accelerate the convergence and improve the design accuracy of optimization. The proposed methodology is proved to be a feasible method for µCGG design and its computational efficiency is verified.