CFD中文网

@李东岳 明白了，感谢李老师！

代码中对网格进行了遍历，如何指导某个网格是否是壁面处的网格呢？或者如何遍历壁面附近的网格呢？

forAll(pCells, cellI)
{
    ...
}

@shrine thermophysicalProperties文件里设置两相的名称还有表面张力系数，表面张力系数如果使用liquidProperties类型就是从指定的液体组分读取例如水的话就是H2O.C那个文件。
另外两个分别定义两相的热物性，具体是哪部分有疑问呢？

p_rghEqnComp1 =
    pos(alpha1)
    *(
        (
            fvc::ddt(alpha1, rho1) + fvc::div(alphaPhi1*rho1f)
          - (fvOptions(alpha1, mixture.thermo1().rho())&rho1)
        )/rho1
      - fvc::ddt(alpha1) - fvc::div(alphaPhi1)
      + (alpha1*psi1/rho1)*correction(fvm::ddt(p_rgh))
    );

根据压力修正方程
$$\mathrm{\nabla }\cdot U+\frac{{\alpha }_1}{{\rho }_1}\frac{D{\rho }_1}{Dt}+\frac{{\alpha }_2}{{\rho }_2}\frac{D{\rho }_2}{Dt}=0$$
这里p_rghEqnComp1应该是
$$\frac{{\alpha }_1}{{\rho }_1}\frac{D{\rho }_1}{Dt}$$
代码的前面部分，即

    pos(alpha1)
    *(
        (
            fvc::ddt(alpha1, rho1) + fvc::div(alphaPhi1*rho1f)
          - (fvOptions(alpha1, mixture.thermo1().rho())&rho1)
        )/rho1
      - fvc::ddt(alpha1) - fvc::div(alphaPhi1)

经过简化就已经是$\frac{{\alpha }_1}{{\rho }_1}\frac{D{\rho }_1}{Dt}$了，所以后面的(alpha1*psi1/rho1)*correction(fvm::ddt(p_rgh))是哪里来的呢，特别是Foam::correction这个函数不太理解，看了fvMatrix.H源码也不甚清楚，我的理解后面这段代码就是
$$\frac{{\alpha }_1{\psi }_1}{{\rho }_1}\frac{\mathrm{\partial }{p}_{rgh}}{\mathrm{\partial }t}=\frac{{\alpha }_1}{{\rho }_1}\frac{\mathrm{\partial }{\rho }_1}{\mathrm{\partial }t}$$
这样的话岂不是重复了？

@gengchen 之前跑les也是，只有中心差分对得上实验，但是震荡。。

@gengchen geng大神

phi一般是流场解出来之后，用速度在网格面心插值得到的surfaceScalarField。也就是说，phi是求解出来的量，不能随意改动吧？如果要改温度方程对流项，应该再叠加一个phir，类比interFoam求解器相分数方程的interface compression term。