in Kugelkoordinaten schreibt sich das Feld wie folgt:
v=r1⎝⎛rsin(φ)sin(ϑ)−rcos(φ)sin(ϑ)1⎠⎞=⎝⎛sin(φ)sin(ϑ)−cos(φ)sin(ϑ)r1⎠⎞=−sin(ϑ)⎝⎛−sin(φ)cos(φ)0⎠⎞+r1⎝⎛001⎠⎞=−sin(ϑ)eφ+r1(cos(ϑ)er−sin(ϑ)eϑ)
Die Divergenz in Kugelkoordinaten findest du hier :
https://de.wikipedia.org/wiki/Divergenz_eines_Vektorfeldes#Zylinder-…
Damit ergibt sich:
div(v)=r21∂r∂(rcos(ϑ))+rsin(ϑ)1∂ϑ∂(r−sin2(ϑ))+rsin(ϑ)1∂φ∂sin(ϑ)=r2cos(ϑ)−2r2cos(ϑ)+0=−r2cos(ϑ)
Das ist hier also doch etwas umständlich, da ein Teil von v in Richtung z-Achse zeigt. In kartesischen Koordinaten hätte man auch
dxd(x2+y2+z2−y)+dyd(x2+y2+z2x)+dzd(x2+y2+z21)=−x2+y2+z23z=−cos(ϑ)/r2
ermitteln können.