ρ{∂v/∂t+v( ^)∇(^ )v}=-∇p+μ∇^2v+ρf