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