Process-constrained batch Bayesian approaches for yield optimization in multi-reactor systems