Abstract. Maximising the power production of large wind farms is key to the transition towards net zero. The overarching goal of this paper is to propose a computational method to maximise the power production of wind farms with two practical design strategies. First, we propose a gradient-free method to optimise the wind farm power production with high-fidelity surrogate models based on large-eddy simulations and a Bayesian framework. Second, we apply the proposed method to maximise wind farm power production by both micro-siting (layout optimisation) and wake steering (yaw angle optimisation). Third, we compare the optimisation results with the optimisation achieved with low-fidelity wake models. Finally, we propose a simple multi-fidelity strategy by combining the inexpensive wake models with the high-fidelity framework. The proposed gradient-free method can effectively maximise wind farm power production. Performance improvements relative to wake-model optimisation strategies can be attained, particularly in scenarios of increased flow complexity, such as in the wake steering problem, in which some of the assumptions in the simplified flow models become less accurate. The optimisation with high-fidelity methods takes into account nonlinear and unsteady fluid mechanical phenomena, which are leveraged by the proposed framework to increase the farm output. This paper opens up opportunities for wind farm optimisation with high-fidelity methods and without adjoint solvers.