东岳老师CES论文中两相kEpsilon模型实现的疑问

Samuel-Tu

论文中的形式如下图：
微信图片_20191012170512.png

疑问1：va,t指湍流粘度，且k和epsilon方程中用的同样的湍流粘度，而在代码中确是用的有效粘度，且k和epsilon方程中的有效粘度不同，分别为DkEff()和DepsilonEff()，其定义如下：

         tmp<volScalarField> DkEff() const
         {
             return tmp<volScalarField>
             (
                 new volScalarField
                 (
                     "DkEff",
                     (this->nut_/sigmak_ + this->nu())
                 )
             );
         }
         tmp<volScalarField> DepsilonEff() const
         {
             return tmp<volScalarField>
             (
                 new volScalarField
                 (
                     "DepsilonEff",
                     (this->nut_/sigmaEps_ + this->nu())
                 )
             );
         }

可以看出是湍流粘度除以一个数(疑似是Schmidt数）+层流粘度。因此对于粘度计算，代码和公式不一致。

疑问2：代码中的sigmak_=1.0, sigmaEps_=1.3，刚好和CES论文中的Schmidt数分别对应
但对于epsilon的Schmidt数来说，论文和代码互为倒数。这里也不一致。
疑问3：CES论文提到对于单相kEpsilon模型，有
这说明a相的粘度是靠b相的粘度计算来的，但是我没有找到相关代码，希望大家能告诉我一下文件路径，让我学习一下如何实现。另外，对于双相kEpsilon模型，是否还是通过b相的粘度计算a相的粘度呢？？
疑问4：对于公式（32）中C3项（绿框部分），是否少了一个epsilon。因为源码C3项中含有epsilon：

- fvm::SuSp(((2.0/3.0)*C1_ - C3_)*alpha()*rho()*divU, epsilon_)

李东岳

@Samuel-Tu 在东岳老师CES论文中两相kEpsilon模型实现的疑问中说：

这说明a相的粘度是靠b相的粘度计算来的，但是我没有找到相关代码，希望大家能告诉我一下文件路径，让我学习一下如何实现。

你可以参考下Gosman 1992年那篇文章

对于公式（32）中C3项（绿框部分），是否少了一个epsilon。因为源码C3项中含有epsilon：

应该是少了个epsilon