CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    bunrs湍流耗散力

    OpenFOAM
    1
    2
    293
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • A
      anubis 最后由 编辑

      com版本和org版本的burns湍流耗散力好像不太一样

      OpenFOAM8里的

      Foam::turbulentDispersionModels::Burns::D() const
      {
          const fvMesh& mesh(pair_.phase1().mesh());
          const dragModel& drag =
              mesh.lookupObject<dragModel>
              (
                  IOobject::groupName(dragModel::typeName, pair_.name())
              );
      
          return
              drag.Ki()
             *continuousTurbulence().nut()
             /sigma_
             *pair_.dispersed()
             *sqr(pair_.dispersed() + pair_.continuous())
             /(
                  max(pair_.dispersed(), pair_.dispersed().residualAlpha())
                 *max(pair_.continuous(), pair_.continuous().residualAlpha())
              );
      }
      

      OpenFOAM2012里的

      Foam::turbulentDispersionModels::Burns::D() const
      {
          const fvMesh& mesh(pair_.phase1().mesh());
          const dragModel&
              drag
              (
                  mesh.lookupObject<dragModel>
                  (
                      IOobject::groupName(dragModel::typeName, pair_.name())
                  )
              );
      
          return
              0.75
             *drag.CdRe()
             *pair_.continuous().nu()
             *pair_.continuous().turbulence().nut()
             /(
                  sigma_
                 *sqr(pair_.dispersed().d())
              )
             *pair_.continuous().rho()
             *pair_.dispersed()
             *(
                 1.0/max(pair_.dispersed(), residualAlpha_)
               + 1.0/max(pair_.continuous(), residualAlpha_)
              );
      }
      
      

      公式应该是
      6bd2f00f-aa3e-4c3e-bc02-ba62f01dfd1b-image.png
      按照这个公式,感觉2012版本里的比较接近,但是还不太一样,相对速度Ur,在2012的代码里好像变成了连续相的粘度除以离散相直径,虽然量纲和速度是一样的,但是不太明白为什么会这样?

      A 1 条回复 最后回复 回复 引用
      • A
        anubis @anubis 最后由 编辑

        又仔细看了一下,两个版本应该是差不多的。
        相对速度是因为CdRe不是曳力系数 而是曳力系数×雷诺数 所以可以用粘度、直径还有密度换算成相对速度。

        1 条回复 最后回复 回复 引用
        • First post
          Last post