Mantle models with high viscosity plates and a low viscosity asthenosphere are investigated.
Free-slip and force-balance surface boundary conditions are used.
Long wavelength flows are found with changing heating rate and Rayleigh number.
Mean heat flux decreases with aspect ratio in force-balance models.