where ∂/∂t is the Eulerian time derivative; ρ is the local density, which characterises
the amount of mass per unit volume (kg/m3); v is local velocity (m/s) and div() or
∇·denotes the divergence operator. The divergence is a scalar function of a vector