方程实际要求的是：
\begin{equation}
\frac{\partial \rho\mathbf{U}}{\partial t}+\nabla\cdot(\rho\mathbf{U}\mathbf{U})=0
\end{equation}
对应：
连续性方程：
\begin{equation}
\frac{\partial \rho}{\partial t}+\nabla\cdot(\rho\mathbf{U})=0
\end{equation}
对应：
因此实际求解的就是：
fvVectorMatrix UEqn ( fvm::ddt(rho, U) + fvm::div(rhoPhi, U) + turbulence->divDevRhoReff(rho, U) );把- fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)包含进来的原因为这有利于使求解变量有界以更好地收敛。
