关于compressibleInterFoam中压力方程的疑问

[ZX]

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))
    );

根据压力修正方程
$$\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{\partial p_{rgh}}{\partial t} = \frac{\alpha_1 }{\rho_1} \frac{\partial \rho_1}{\partial t}$$
这样的话岂不是重复了？

[shrine]

@zx 请教一下，thermophysicalProperties三个文件该如何设置，看不懂各参数什么意思

[ZX]

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

[shrine]

@zx 我想计算管道内的水锤，因为管道有调压井，所以又有自由液面

看到tutorials里面只有层流的例子，这个不能计算紊流吗？

那三个文件中
phases (water air);

pMin 10000;

sigma 0.07;
这两个变量不知道什么意思
还有水的文件中
equationOfState
{
C (0.001278 -2.1055e-06 3.9689e-09 4.3772e-13 -2.0225e-16);
}

这个也不知道是啥，我这个算例是不是不用变