东岳老师CES论文中两相kEpsilon模型实现的疑问
-
论文中的形式如下图:
- 疑问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_)
- 疑问1:
-
@Samuel-Tu 在 东岳老师CES论文中两相kEpsilon模型实现的疑问 中说:
这说明a相的粘度是靠b相的粘度计算来的,但是我没有找到相关代码,希望大家能告诉我一下文件路径,让我学习一下如何实现。
你可以参考下Gosman 1992年那篇文章
对于公式(32)中C3项(绿框部分),是否少了一个epsilon。因为源码C3项中含有epsilon:
应该是少了个epsilon