Numerical methods for coupled thermal and electric current fields in axisymmetric devices