I went back and checked everything running models using all the default parameters for the visco plastic material model. For the boundary conditions, I prescribe the outer spherical boundary at 273 K and the inner at 2973 K with no free surface. I looked and the residuals are small and the compositional fields are similar between the models.