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

mechanicsdog

用DPMFoam算二维喷动床流化现象，求解出颗粒相的运动好像一切正常，如下图

但是反应器的压降过于剧烈是什么情况？压力求解结果不对吗？

压力虽然波动剧烈，但是平均值都没问题，不知道是什么问题，换了几个入口速度都有类似的情况。按照文献压力波动不会这么大，±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;
            }
        );
    }
}

mechanicsdog

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

上级

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

上级

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

李东岳

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

李东岳

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

mechanicsdog

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

上级

@李东岳读了东岳老师的这篇文章，不知道公式（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}$，当时只是数值这么给的，没什么特别的含义。确实应该说明一下。

上级

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