Surrogate-based acceleration of the coupling of black-box solvers