We aim to develop fast and reliable surrogate models for yaw-based wind farm control. The surrogates, based on polynomial chaos expansion, are built using high-fidelity flow simulations combined with aeroelastic simulations of the turbine performance and loads. Optimization results performed using two Vestas V27 turbines in a row for a specific atmospheric condition suggest that a power gain of almost 3 % ± 1 % can be achieved at close spacing by yawing the upstream turbine more than 15°.