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



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



  • 想问作者一个问题,请问是如何在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;
                }
            );
        }
    }
    


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



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



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



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



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



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



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



  • @mechanicsdog 谢谢您,万分感谢



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

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



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


Log in to reply