Evaluate the flux of the vector field $\mathbf{F} = \langle x, y, z \rangle$ through the boundary of the region defined by $x^2 + y^2 + z^2 \le a^2$ using the Divergence Theorem.