I have recently come to realise that many of the most fundamental theorems can be reduced to a continuity equation. Doing some research on the topic of said equations, I have found out they have stronger implications than just the conservation of a given magnitude, since they also impose locality. I am specially interested in Poynting's theorem:
$$-\frac{\partial U}{\partial t}=\int_V\vec{J}\cdot \vec{E}dV + \oint \vec{S}\cdot d\vec{A}$$
Using Gauss' theorem to transform the surface integral into a volume integral and taking a region of space in which there is no free charge, we get:
$$\frac{du}{dt}+\vec{\nabla}\cdot \vec{S}=0$$
where $u$ is the electromagnetic energy density, which fullfils: $U=\int_V u\ dV$
What are the implications of this? My best guess is the electromagnetic energy density in a given region diminishes when EM energy escapes said region in the form of radiation, represented by the divergence of Poynting's vector. Is there a general (physical) interpretation for what happens to the system for two general magnitudes $f,\vec{G}$ such that they fulfill the continuity equation?
$$\frac{\partial f}{\partial t}+\vec{\nabla}\cdot\vec{G}=0$$
I have also found that, regardless of which theorem I choose to reduce to this form, I always end up with a density on the time derivative and a flux on the gradient. How so?