Multimode modelling of the Rayleigh-Taylor instability