$\frac{d}{dt}\int \rho d\tau+\oint \rho V\cdot d\sigma$