Each component gives, e.g., \(\partial_x\big((x^2+y^2+z^2)x\big)=3x^2+y^2+z^2\). Adding the three yields \(5(x^2+y^2+z^2)\). By the product rule: \(\nabla f\cdot\mathbf r+f\,\operatorname{div}\mathbf r=2r^2+3r^2\).
The divergence of a vector field \(\mathbf F=\langle P,Q,R\rangle\) is the scalar field that measures, at each point, how much the flow of the field expands or compresses.