twoPhaseEulerFoam pEqn.H bug?

[Jasper 0]

大家好，本人在学习twoPhaseEulerFoam pEqn.H 代码时发现，code中定义的pressure "diffusivity"

 "rAUf",
mag(alphaf1*alpharAUf1 + alphaf2*alpharAUf2)

alpharAUf1在UEqn.H中给出了定义：

fvc::interpolate(max(alpha1, phase1.residualAlpha())*rAU1)

与公式中给出的形式并不一样（下图右边项），

这是代码出现的bug还是我理解的不到位？还请大佬们解惑

[Jasper 0]

@李东岳 现在这个问题还没有搞清楚，东岳大佬能否帮忙解答一下？

[李东岳]

跟下面这个一致么

$$\begin{array}{}\text{(1)}& \mathrm{\nabla }\cdot \left(({\alpha }_{\mathrm{d},\mathrm{P}}\frac{{\alpha }_{\mathrm{d},\mathrm{P}}}{A_{\mathrm{d},\mathrm{P}}}+{\alpha }_{\mathrm{d},\mathrm{P}}\frac{{\alpha }_{\mathrm{d},\mathrm{P}}}{A_{\mathrm{d},\mathrm{P}}})\mathrm{\nabla }{p}_{\mathrm{r}\mathrm{g}\mathrm{h},\mathrm{P}}\right),\end{array}$$

[Jasper 0]

@李东岳 一致，但是我实在看不出来是怎么推出来的这个形式。
我是看的这个网站上的推导http://xiaopingqiu.github.io/2015/05/17/twoPhaseEulerFoam2/
请问李老师您这个公式有相应的文献或者网站吗？万分感谢！困扰好久了

[李东岳]

@jasper-0 https://www.jianguoyun.com/p/DbNqqy8Q9s3ZBhjj5tYCIAA

[Jasper 0]

@李东岳 感谢大佬！我好好研读一下