With [tex]\mathbf f=\mathbf r[/tex] (?) you have divergence
[tex]\nabla\cdot\mathbf f=\dfrac{\partial x}{\partial x}+\dfrac{\partial y}{\partial y}+\dfrac{\partial z}{\partial z}=3[/tex]
Then by the divergence theorem, the flux is
[tex]\displaystyle\iint_S\mathbf f\cdot\mathrm d\mathbf S=3\iiint_R\mathrm dV[/tex]
where [tex]R[/tex] denotes the interior of the region with boundary [tex]S[/tex].
Convert to spherical coordinates to set up the volume integral as
[tex]\displaystyle3\int_{\varphi=0}^{\varphi=\pi/2}\int_{\theta=0}^{\theta=2\pi}\int_{\rho=0}^{\rho=4}\rho^2\sin\varphi\,\mathrm d\rho\,\mathrm d\theta\,\mathrm d\varphi[/tex]
[tex]=\displaystyle6\pi\left(\int_0^4\rho^2\,\mathrm d\rho\right)\left(\int_0^{\pi/2}\sin\varphi\,\mathrm d\varphi\right)[/tex]
[tex]=128\pi[/tex]