Free surface computations in mantle convection models