ρ(∂v/∂t+v⋅∇v)=−∇P+μ∇^2v+ρf