Skip to content
  • 最新
  • 版块
  • 东岳流体
  • 随机看[请狂点我]
皮肤
  • Light
  • Cerulean
  • Cosmo
  • Flatly
  • Journal
  • Litera
  • Lumen
  • Lux
  • Materia
  • Minty
  • Morph
  • Pulse
  • Sandstone
  • Simplex
  • Sketchy
  • Spacelab
  • United
  • Yeti
  • Zephyr
  • Dark
  • Cyborg
  • Darkly
  • Quartz
  • Slate
  • Solar
  • Superhero
  • Vapor

  • 默认(不使用皮肤)
  • 不使用皮肤
折叠
CFD中文网

CFD中文网

  1. CFD中文网
  2. OpenFOAM
  3. 东岳老师CES论文中两相kEpsilon模型实现的疑问

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

已定时 已固定 已锁定 已移动 OpenFOAM
2 帖子 2 发布者 1.4k 浏览
  • 从旧到新
  • 从新到旧
  • 最多赞同
回复
  • 在新帖中回复
登录后回复
此主题已被删除。只有拥有主题管理权限的用户可以查看。
  • S 离线
    S 离线
    Samuel-Tu
    写于 最后由 编辑
    #1

    论文中的形式如下图:
    微信图片_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_)
    
    1 条回复 最后回复
  • 李东岳李 在线
    李东岳李 在线
    李东岳 管理员
    写于 最后由 编辑
    #2

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

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

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

    对于公式(32)中C3项(绿框部分),是否少了一个epsilon。因为源码C3项中含有epsilon:

    应该是少了个epsilon

    http://dyfluid.com/index.html
    需要帮助debug算例的看这个 https://cfd-china.com/topic/8018

    1 条回复 最后回复

  • 登录

  • 登录或注册以进行搜索。
  • 第一个帖子
    最后一个帖子
0
  • 最新
  • 版块
  • 东岳流体
  • 随机看[请狂点我]