The motivation for the divergence being interpreted as the flux of stuff used the following:
$$\text{div} F(a) = \lim_{r\to0}\frac{3}{4\pi r^3}\int_{|x-a|=r} F\cdot n dA$$
Without the $r^3$ in the denominator the right hand side would converge to $0$. But in the context of the divergence theorem what happens if it is $r^2$ instead, can the limit be nonzero?
$\mathbf F=1 \mathbf e_r$. Then the integral is just $4\pi r^2$ and the limit is 3.
I'm not clear on what you mean by "in the context of the divergence theorem", though. In three dimensions, the size of the sphere of radius $r$ is always proportional to $r^3$, and there is no way you would get the above formula with $r^2$ instead.
Of course, in two dimensions, this is exactly what one does find.
$\nabla\cdot \mathbf F=\lim \frac 1 {\pi r^2} \int_{\text{circle}} \mathbf {F\cdot dn}$
Try $\mathbf F = r \mathbf e_r$ for a simple example. Both sides should give 2.