We report on molecular dynamics studies of heat flow in superlattices. The computer simulations are performed using classical mechanics with periodic boundary conditions. The heat flow is in the direction normal to the layers. We have studied the variation of the conductivity with the repeat distance and the effect of interfacial roughness. We discuss the relation of these results to experimental data in the literature.