twoPhaseEulerFoam pEqn.H bug?
-
大家好，本人在学习twoPhaseEulerFoam pEqn.H 代码时发现，code中定义的pressure "diffusivity"
"rAUf", mag(alphaf1*alpharAUf1 + alphaf2*alpharAUf2)alpharAUf1在UEqn.H中给出了定义：
fvc::interpolate(max(alpha1, phase1.residualAlpha())*rAU1)与公式中给出的形式并不一样（下图右边项），
这是代码出现的bug还是我理解的不到位？还请大佬们解惑
-
@李东岳 现在这个问题还没有搞清楚，东岳大佬能否帮忙解答一下？
-
跟下面这个一致么
\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} -
@李东岳 一致，但是我实在看不出来是怎么推出来的这个形式。
我是看的这个网站上的推导http://xiaopingqiu.github.io/2015/05/17/twoPhaseEulerFoam2/
请问李老师您这个公式有相应的文献或者网站吗？万分感谢！困扰好久了 -
-
@李东岳 感谢大佬！我好好研读一下
