Free surface flow of a layered medium of simple fluids driven by a torsionally oscillating, cylindrical rod is investigated. The stress, expressed as a series of multiple integrals of polynomials in the strain histories, is expanded into a Fréchet series in terms of the amplitude of the oscillation of the rod. A domain perturbation method is used to perturb the rest state. The mean shapes and deviations from flatness at the rod of the interfaces between layers are determined together with the flow field up to and including second order in the amplitude of the oscillation. © 1986.