CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    颗粒曳力的孔隙率alphac是怎么读取的(插值计算还是网格中心值)?

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

      ErgunWenYuDragForce.C中,计算颗粒曳力的代码如下:

           scalar alphac(alphac_[p.cell()]);
        
           if (alphac < 0.8)
           {
               return forceSuSp
               (
                   Zero,
                   (mass/p.rho())
                  *(150.0*(1.0 - alphac)/alphac + 1.75*Re)*muc/(alphac*sqr(p.d()))
               );
           }
           else
           {
               return forceSuSp
               (
                   Zero,
                   (mass/p.rho())
                  *0.75*CdRe(alphac*Re)*muc*pow(alphac, -2.65)/(alphac*sqr(p.d()))
               );
           }
      

      这里第一行的alphac(alphac_[p.cell()]),是插值计算网格的alphac值(cellPoint方式),还是只读取网格中心值(cell方式)?因为没有查到关于这个[p.cell()]的解释,所以想请教一下论坛里的前辈

      李东岳 1 条回复 最后回复 回复 引用
      • 李东岳
        李东岳 管理员 @hotairballon1997 最后由 编辑

        @hotairballon1997 p.cell()是拉格朗日粒子所在单元的标号,那就是你说的只读取网格中心值(cell方式)

        CFD高性能服务器 http://dyfluid.com/servers.html

        H 1 条回复 最后回复 回复 引用
        • H
          hotairballon1997 @李东岳 最后由 编辑

          @李东岳 谢谢李老师解答。

          那是不是可以理解为,在比如下图中(打开颗粒对流体的影响),颗粒的中心点在网格1,但是颗粒的实际大小有如图显示那么大。现在就只有网格1受到颗粒的影响,而2,3,4号网格都完全感受不到颗粒的存在?

          不知道这对总体结果会造成多大的影响?

          mesh_and_particle.jpg

          1 条回复 最后回复 回复 引用
          • 李东岳
            李东岳 管理员 最后由 编辑

            应该影响不大。这种颗粒应该不多,不会所有的颗粒都恰好那样。所以对alpha的影响也很小。除非颗粒直径非常大。不过直径非常大的话,那就会引起收敛性的问题。

            CFD高性能服务器 http://dyfluid.com/servers.html

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