CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    DPMFoam算流化床压降波动过于剧烈

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

      用DPMFoam算二维喷动床流化现象,求解出颗粒相的运动好像一切正常,如下图
      2.JPG 1.JPG
      但是反应器的压降过于剧烈是什么情况?压力求解结果不对吗?
      3.JPG
      压力虽然波动剧烈,但是平均值都没问题,不知道是什么问题,换了几个入口速度都有类似的情况。按照文献压力波动不会这么大,±2000Pa以内比较正常。

      1 条回复 最后回复 回复 引用
      • 上级
        上级 最后由 编辑

        想问作者一个问题,请问是如何在paraview里面把颗粒上色,显示不同的速度的大小的呢。
        此外,是否压力本来就该取平均值,在controlDict里面添加平均库函数来获得平均场,在这里贴上我的用过的例子:

        functions
        {
            fieldAverage
            {
                type          fieldAverage;
                libs          ("libfieldFunctionObjects.so");
                writeControl  writeTime;
                timeStart     1;
                fields
                (
                    UField
                    {
                        mean          on;
                        prime2Mean    on;
                        base          time;
                    }
                             
                    alpha
                    {
                        mean          on;
                        prime2Mean    on;
                        base          time;
                    }
                    
                    alpha.air
                    {
                        mean          on;
                        prime2Mean    on;
                        base          time;
                    }
                    
                    U.air
                    {
                        mean          on;
                        prime2Mean    on;
                        base          time;
                    }
                );
            }
        }
        
        M 1 条回复 最后回复 回复 引用
        • M
          mechanicsdog @上级 最后由 编辑

          @上级 第一个问题的方法是
          1、用foamToVTK命令把他转成VTK格式
          2、用paraview打开这个VTK文件夹中的拉格朗日文件,并用高斯点表示
          3、设置半径就可以了
          第二个问题:确实就是用平均压力来表示该流化速度下的压降,但是这个流化速度下的压降随时间变化也不该这么剧烈,下面是一篇文献中的压降随时间变化的曲线
          捕获.JPG
          所以兄台你的压力波动变化大吗:135:

          上级 1 条回复 最后回复 回复 引用
          • 上级
            上级 最后由 编辑

            我也遇到过相同的问题,平均压力没问题,但是瞬时波动很大,暂时没有找出原因所在,难道是算法的问题么

            1 条回复 最后回复 回复 引用
            • 上级
              上级 @mechanicsdog 最后由 编辑

              @mechanicsdog 您好,我想再请教您一个问题,我按照您的方法实现颗粒上显示速度了,现在我在算一个带颗粒级配的流化床问题,颗粒直径是高斯分布,请问后处理的时候您知道如何显示颗粒的大小么,现在按照高斯点好像只能显示一样大小的颗粒,如果想这么显示您知道如何操作么?谢谢您
              捕获.JPG

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

                压力那个问题 可以跟基金会提bug

                线上CFD课程开始报名:http://www.dyfluid.com/class.html

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

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

                  https://www.sciencedirect.com/science/article/pii/S0010465519303728 我在这篇文章里面,是基于DPMFoam,把拉格朗日的东西换成了GPBE的东西。代码改动非常大,压力方程最后推出来有一些差别。我们这个代码压力没有问题。DPMFoam压力波动那个我也注意到了。一直没仔细debug,有新进展可以交流。

                  线上CFD课程开始报名:http://www.dyfluid.com/class.html

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

                  上级 1 条回复 最后回复 回复 引用
                  • M
                    mechanicsdog @上级 最后由 编辑

                    @上级 您好,这个可以在高斯点中开启高斯阵列,如下图
                    1.JPG

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

                      @李东岳 读了东岳老师的这篇文章,不知道公式(22)的量纲是否有点问题,按照公式(21),弛豫时间的$\tau$量纲应该是$[M^{-1}L^{3}T^{1}]$,那么公式(22)右边第二项的两处量纲都是对不上的,大括号内$\mathbf{g}$那一串不是速度量纲,右边的$\tau$应该是少一个-1次方。另外,请问李老师,您文章onewaycoupledtestcase里面计算颗粒分离时,用的弛豫时间$10d^{2/3}$,请问一下这部分该怎么理解,2/3次方表示的是什么呢?另外,弛豫时间一般理解,量纲不是时间量纲么?谢谢李老师。
                      贴出您文章的公式(21)和(22)
                      捕获.JPG

                      1 条回复 最后回复 回复 引用
                      • 上级
                        上级 @mechanicsdog 最后由 编辑

                        @mechanicsdog 谢谢您,万分感谢

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

                          https://www.jianguoyun.com/p/DTXEFVsQ9s3ZBhj9x8oD 量纲那个你看下这个公式 有详细的推导步骤 现在我一时半会想不起来了 方程11-18

                          $\tau=d^{2/3}$,当时只是数值这么给的,没什么特别的含义。确实应该说明一下。

                          线上CFD课程开始报名:http://www.dyfluid.com/class.html

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

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

                            @李东岳 谢谢李老师,我再详细看一下推导过程

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