I am trying to get the gradient of function f using the following coordinate-free del operator.
$\nabla = \lim_{\Delta v \to 0} \frac{\iint_S dS}{\Delta v} $
For infinitesimal volume from $ (r, \phi , z) $ to $(r+ \Delta, \phi + \Delta \phi, z+\Delta z) $
here,
$$ \Delta v = r \Delta r \Delta \phi \Delta z $$
Through r-direction we can calculate the following
$$ \iint_{S_r} f dS = f(r+\Delta r, \phi, z)(r+\Delta r)\Delta \phi \Delta z - f(r,\phi, z)r\Delta \phi \Delta z $$ Then, (when $ \Delta v \to 0 $) $$ \frac{\iint_{S_r} f dS}{\Delta v} = \frac{\partial f}{\partial r} + \frac{1}{r}f = \frac{1}{r}\frac{\partial}{\partial r} (rf)\ \cdots \ (1) $$
But, as you know, the r component of the gradient is $ \frac{\partial f}{\partial r} $
Surprisingly, the calculation is correct in getting divergence!
Can somebody tell me what's wrong with the calculation in gradient?