求助：两相VOF求解器中当一相是可压的，另一相不可压时，压力泊松方程如何构建

[lrl3512]

模拟气液两相蒸发燃烧，不考虑高速流动情况，参考了两个求解器，icoReactingMultiphaseFoam和compressibleInterFoam，前者为不可压缩求解器，压力泊松方程为

            fvc::div(phiHbyA)
          - fvm::laplacian(rAUf, p_rgh)

而compressibleInterFoam为

               (max(alpha1, scalar(0))/rho1)*p_rghEqnComp1()
              + (max(alpha2, scalar(0))/rho2)*p_rghEqnComp2()
            )
          + p_rghEqnIncomp

两相均为可压缩流体

那气液两相中，液相为不可压缩，气相为可压的情况下，压力泊松方程如何构建

[李东岳]

在压力方程中，不可压缩项pEqnIncom、可压缩项pEqnComp是分开处理的，求解的时候是两个贡献之和。所以你看一下compressibleInterFoam怎么处理的就可以了。

[学流体的小明]

@lrl3512 compressibleInterFoam可以参考这个
https://zhuanlan.zhihu.com/p/66444729
我也是看这个学的