The response of turbulent flow to time-modulated forcing is studied by direct numerical simulation of the Navier–Stokes equations. The forcing is modulated via periodic energy-input variations at a frequency x. Harmonically modulated forcing of the large scales is shown to yield a response maximum at frequencies in the range of the inverse of the large-eddy turnover time, as well as a characteristic rapid change of the phase-angle between forcing and response. Harmonically modulated broadband forcing is also studied in case a wide spectrum of length-scales is forced simultaneously. If smaller length-scales are also explicitly agitated by the forcing, the response maximum is found to occur at higher frequencies and to become less pronounced....