关于几个公式的疑问
-
第一个公式是An Introduction to Computational Fluid Dynamics The Finite Volume Method 2nd Edition中的13页式2-10
https://pan.baidu.com/s/1nvqBgtb
的最后一项是不是应该为$D(\phi\rho)/Dt$第二个公式是
http://openfoam.org/release/2-2-0/numerics-boundedness/
中的第二个公式
左边项应该为$\rho Dh/Dt$,右边第二项应该是div(\rho U h) -
方程实际要求的是:
\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)
包含进来的原因为这有利于使求解变量有界以更好地收敛。