Turbulent flow over a series of increasingly high, two-dimensional bumps is studied by well-resolved large-eddy simulation. The mean flow and Reynolds stresses for the lowest bump are in good agreement with experimental data. The flow encounters a favourable pressure gradient over the windward side of the bump, but does not relaminarize, as is evident from near-wall fluctuations. A patch of high turbulent kinetic energy forms in the lee of the bump and extends into the wake. It originates near the surface, before flow separation, and has a significant influence on flow development. The highest bumps create a small separation bubble, whereas flow over the lowest bump does not separate. The log law is absent over the entire bump, evidencing strong disequilibrium. This dataset was created for data-driven modelling. An optimization method is used to extract fields of variables that are used in turbulence closure models. From this, it is shown how these models fail to correctly predict the behaviour of these variables near to the surface. The discrepancies extend further away from the wall in the adverse pressure gradient and recovery regions than in the favourable pressure gradient region.