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{equation}\label{incompressiblep}
\nabla\cdot\left(\left(\alpha_{\rd,\rP}\frac{\alpha_{\rd,\rP}}{A_{\rd,\rP}}+\alpha_{\rd,\rP}\frac{\alpha_{\rd,\rP}}{A_{\rd,\rP}}
\right)\nabla p_{\mathrm{rgh},\rP}\right),
\end{equation}

[Jasper 0]

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

[李东岳]

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

[Jasper 0]

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