Abstract. Given a wind farm with known dimensions and number of wind turbines, we try to find the optimum positioning of wind turbines that maximises wind-farm energy production. In practice, given that optimisation has to be performed for many wind directions, and taking into account the yearly wind distribution, such an optimisation is computationally only feasible using fast engineering wake models such as the Jensen model. These models are known to have accuracy issues, in particular since their representation of wake interaction is very simple. In the present work, we propose an optimisation approach that is based on a hybrid combination of large-eddy simulation (LES) and the Jensen model; in this approach, optimisation is mainly performed using the Jensen model, and LES is used at a few points only during optimisation for online tuning of the wake-expansion coefficient in the Jensen model, as well as for validation of the results. An optimisation case study is considered, in which the placement of 30 turbines in a 4 km by 3 km rectangular domain is optimised in a neutral atmospheric boundary layer. Optimisation for both a single wind direction and multiple wind directions is discussed.