(((∂ρ)/(∂t))+∇·ρ u=0(16)(∂/(∂t))(ρ u)+∇·(ρ uu)=∇·σ+ρ g(17)σ=-PI+τ
τ=2η D(18)ρ (Subscript[C, P])(((∂T)/(∂t))+v·∇T)=∇·(k∇T)+η((OverDot[γ])^2)(19))