@东岳 在 OpenFOAM-6的twoPhaseEulerFoam学习 中说：
是的，后面的kI放在了压力里面
最开始我以为是把kI放入了压力泊松方程中，但是我pEqn里一直都没有出现kI，后来我在xiaopingqiu的个人网站上看见了一种说法：动量方程中的压力是雷诺时均压力和雷诺应力的各向同性分量（2k/3)之和。
这种观点可以在 Pope 2000 书第 88 页找到依据。在 “ The Finite Volume Method in Computational Fluid Dynamics: An Advanced Introduction with OpenFOAM® and Matlab®” 这本书的第 699 页，也提到 k 是被放到压力项里去了，目的在于使动量方程中只含有 νt 这一个跟湍流有关的未知量——来源于xiaopingqiu
functions
{
writeFields // name of the function object
{
type writeObjects;
libs ( "libutilityFunctionObjects.so" );
objects
(
T U rho // list of fields/variables to be written
);
// E.g. write every 1e-5 seconds of simulation time only the specified fields
writeControl runTime;
writeInterval 1e-5; // write every 1e-5 seconds
}
}
