\rho\left( \partial_t \vec{u} + \vec{u} \cdot \nabla\vec{u} \right) = - \nabla p + \mu \nabla^2 \vec{u}