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



  • 论文中的形式如下图:
    微信图片_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数分别对应微信图片_20191012173142.png
      但对于epsilon的Schmidt数来说,论文和代码互为倒数。这里也不一致。

    • 疑问3:CES论文提到对于单相kEpsilon模型,有微信图片_20191012173414.png
      这说明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


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2