Advanced numerical methods for mantle convection models