两相流通过湍流模型求解粘度项和直接求解粘度项不一致？

[WangQi 0]

最近在尝试改不同的两相流求解器，遇到了一个问题。
当前OF的两相粘度项是通过湍流模型计算的。
turbulence->divDevRhoReff(rho, U)
当我的case不需要湍流使用laminar时，结果和文献结果不匹配：

当我改成直接用拉普拉斯求解：

• fvm::laplacian(nu*rho, U)
  结果：
  

此时与文献结果一致，请问李老师这是什么原因啊？

[李东岳]

即使结果能对上，你植入laplacian(nu,U)是不正确的，nu在VOF里面不是常数，你这么植入缺少项。
http://dyfluid.com/interFoam.html
方程31

[WangQi 0]

@李东岳 李老师，在写fvm::laplacian项之前已经通过体积分数alpha计算过nu了。因为我觉得直接修改turbulence::类里面的nu很麻烦。所以我就自己定义了提及标量场nu然后正常使用alpha计算之后带入laplacian，所以获得了与文献一致的结果。