(((∂ρ)/(∂t))+∇·ρ u=0(1)(∂/(∂t))(ρ u)+∇·(ρ uu)=-∇P+∇·τ+ρ g(2)ρ (Subscript[C, P])(((∂T)/(∂t))+u·∇T)=∇·(k∇T)+τ:D(3))