Investigating Bayesian optimization for expensive-to-evaluate black box functions: Application in fluid dynamics