Euler-Lagrange的一些解析，sprayFoam

星星星星晴

这部分你可以删掉，particlecollector总是有些问题，我也遇到过，后来不必要的东西全删掉了


    forAll(hitFaceIDs_, i)
    {
        label faceI = hitFaceIDs_[i];
        scalar m = p.nParticle()*p.mass();

        if (negateParcelsOppositeNormal_)
        {
            vector Uhat = p.U();
            Uhat /= mag(Uhat) + ROOTVSMALL;
            if ((Uhat & normal_[faceI]) < 0)
            {
                m *= -1.0;
            }
        }

        // add mass contribution
        mass_[faceI] += m;

        if (nSector_ == 1)
        {
            mass_[faceI + 1] += m;
            mass_[faceI + 2] += m;
            mass_[faceI + 3] += m;
        }

        if (removeCollected_)
        {
            keepParticle = false;
        }
    }

vbcwl

@星星星星晴前辈您好，我也是参考了这个教程(GOFUN2018)在pimpleFoam的基础上添加了求解粒子的部分，想请教您一下，这个教程里的方法是不是流体对粒子，粒子对流体，粒子之间，粒子与壁面的相互作用都可以考虑，只需要在kinematicCloudProperties调整相关的设置就可以了？

星星星星晴

@vbcwl 这个也要看你用的是什么of版本，我记得好像在4之前 turbulence dispersion model是有问题的，kinematic parcel好像不做这步，忘了不确定，需要你自己查一下，粒子之间好像kinematic parcel也没有吧，需要你确认，因为有的时候是有这个模型，但是他只是一个空函数，壁面就是简单的 stick,rebound, escape。

然后drag 还有lift的cd cl好像也是有问题的，具体你需要根据你的需求来看，都得查code

vbcwl

@星星星星晴感谢前辈

星星星星晴

@vbcwl
在kinematicparcel.C中

流体对粒子判断受力是通过Fcp.Sp()和Fncp.Su()这两项，显性隐性的区别，calcVelocity方程计算新速度。
calcDispersion方程是dispersionmodel里计算新流场速度，OF自带是基于k-e的

粒子对流体：UTrans这一项传递的

粒子之间：kinematicCloud中的 stochasticCollision().update(td, solution_.trackTime());这一项

粒子与壁面：hitFace， hitpatch这个方程等，具体调用哪个这部分我不是很熟悉
OpenFOAM-8/src/lagrangian/intermediate/submodels/Kinematic/PatchInteractionModel/ 具体模型是在这个部分

具体你用了某个模型的话一定要看好code，有的模型kicematicparcel是没有的

vbcwl

@星星星星晴多谢提醒，一定在了解code以及公式之后在使用模型

chengan.wang

@星星星星晴您好，谢谢回复

[2] file: /home/chengan/OpenFOAM/chengan-2.4.0/run/sprayFoam/watermist/processor2/0/uniform/lagrangian/sprayCloud/sprayCloudOutputProperties.cloudFunctionObject.particleCollector2.massFlowRate at line 0.
[2] 
[2]     From function operator>>(Istream&, Scalar&)
[2]     in file lnInclude/Scalar.C at line 93.
[2] 
FOAM parallel run exiting
[2] 

particleCollector output:
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 2 in communicator MPI_COMM_WORLD
with errorcode 1.

这个问题我大概知道怎么回事了，当resetOnWrite no;如果我采用单核计算，这个错误问题就没有出现；resetOnWrite yes;才能用多核并行，但我不知道有没有啥影响。

第二个问题，我直接在OpenFOAM-2.4.0下面的ParticleCollector.C文件修改，然后在intermediate文件夹下编译,wclean; wmake libso。其它没有动。但是运行结束后在postProcessing文件夹里面没有找到output.txt文件，感觉是没有调用新编译的库文件函数？您的编译方法参考教程能方便发一下么？
Screenshot from 2022-08-27 16-26-52.png

第三个问题，我用polygon做了一个切片，

    particleCollector2
    {
        type            particleCollector;

        mode            polygon;
        polygons
        (
            (
                (-0.01 0.06 -0.01)
                (-0.01 0.06 0.01)
                (0.01 0.06 0.01)
                (0.01 0.06 -0.01)
            )
        );
        normal          (0 0 1);

        negateParcelsOppositeNormal no;
        removeCollected no;
        surfaceFormat   vtk;
        resetOnWrite    yes; 
        log             yes;
    }

Screenshot from 2022-08-27 16-43-54.png

很明显，粒子已经通过我设计的多边形，但是质量流量等于0，这又是什么原因呢？

星星星星晴

@chengan-wang

我没有试过resetOnWrite 这部分，所以我也不知道你这么弄有没有问题。我仅仅是使用output的文件作后处理，套用了particlecollector这个cloudfunction而已，对于其他的东西我并不在意，而且之前总出错，与其用一个无法确定是否正确的东西，不如原始一点自己处理。
编译教程你可以在google上找很多，最主要的是你编译了一个东西以后，要对其重命名，然后相关引用到该lib的都要引用你新编译的模型，仔细看了一下，你在solver部分是有问题的，但是你编译了spray这个文件了么？编译的新spray有没有用到新的intermediate？solver没有引用新的spray。

3.如果你像我说的删除掉那部分code的话，自然没有质量的计算啊，所以后处理就是自己写个python script或者你用matlab或者什么后处理文件，因为你文件里面包含了所有hitpatch的数据。

chengan.wang

@星星星星晴谢谢

我直接改原始文件内容，然后在原始的安装包编译，也不用改名，我的理解是既然动态链接，不改文件名，就自然调用了，所以spray这个文件没有重新编译。原始的安装包2.4就是一系列文件编译，我改了其中一个文件内容，相当于安装包版本升级了一下编译，这种理解不对么？
我没有删除掉code，没有质量的计算，所以觉得很奇怪。

星星星星晴

@chengan-wang

额，这我就不是很了解了，我也不是code大神，仅仅是使用推荐的方法进行编译，这肯定不会错，正常情况下，建议无论更改了什么，最好是重新编译一下相关的东西，以免发生错误。有些东西没办法解释。。。就是这么做，也不会思考为什么。

你修改原始包也没问题，但是不用重新命名。
但是如果有些东西你自己编译错了，很可能会引发连锁反应，这也是为什么人家推荐用的是

$USER_LIB_SRC

$LIB_SRC

不是在人家原始的包里面做修改的原因，这样用以区别你自己的东西和原始的东西。
另外并不是所有人都有root的权限的。
我个人觉得还是人家推荐的方式比较好。

这我就不是很清楚了。

chengan.wang

@星星星星晴谢谢哈，麻烦您这么久

chengan.wang

@星星星星晴您好，第三个错误我找到了，normal之前设置不对改成(0 1 0);就可以了。搞不懂为啥坐标都给了，自己还要定义面的法线方向。

通过和您交流和自己实践理解，感觉对这个后处理问题最好还是用“2. 输出文件”方法，“1. 将拉格朗日场转到EULER场”这种方法我没有尝试，但是感觉应该属于空间插值，准确行应该是不如直接统计自己指定特定位置上的数据。

还有两个小疑问，如果我统计某一条线上的数据，实际上粒子很难碰到这条线，所以我应该设置一个很细长的面或者是圆柱吧？

另外就是有什么比较快捷的方式在polygons设置很多个平面吗？比如20个切面，或者50个。我可以提前编辑好，累加粘贴，但是感觉太麻烦了

        polygons
        (
            (
                (-0.01 0.06 -0.01)
                (-0.01 0.06 0.01)
                (0.01 0.06 0.01)
                (0.01 0.06 -0.01)
            )
        );

星星星星晴

@chengan-wang 在 Euler-Lagrange的一些解析，sprayFoam 中说：

第三个错误我找到了，normal之前设置不对改成(0 1 0);就可以了。搞不懂为啥坐标都给了，自己还要定义面的法线方向。
通过和您交流和自己实践理解，感觉对这个后处理问题最好还是用“2. 输出文件”方法，“1. 将拉格朗日场转到EULER场”这种方法我没有尝试，但是感觉应该属于空间插值，准确行应该是不如直接统计自己指定特定位置上的数据。

个人认为，定义法向是因为平面可以正面接受parcel 或者背面接受parcel
不属于空间差值，就相当于一个cell不停的累加数据，比如你想要这个cell内平均的拉格朗日颗粒的直径，通过方法1可以这样，cell内不停的积累数据，然后做一个时间平均，也可以统计cell内通过多少parcel，统计数量等，最大的问题在于不能精确到每个parcel但是可以相对简洁的做统计，比后处理要快，而且是可以可视化的。最大的问题在于如何平均，时间平均还是数目平均，有的时候平均会消除一些统计学数据，见仁见智吧，每个方法都有其相对的好处。
我不用polygons 所以没对这个部分看，我直接是circle做一个截面。

vbcwl

@星星星星晴前辈您好，想向您请教一下，openFoam如何调用硬球模型呢？

星星星星晴

@vbcwl 硬球模型？碰撞？

默认是这个弹簧模型
OpenFOAM-8/src/lagrangian/intermediate/submodels/Kinematic/CollisionModel/PairCollision/PairModel/PairSpringSliderDashpot/

默认的话是这个在kinematiccloud里面，不同的cloud调用的一些模型是不一样的，你要具体看你用的是什么cloud，具体看吧，

vbcwl

@星星星星晴前辈您好，我是参考的GOFUN在pimpleFoam上加的粒子。然后这里只有那个PairSpringSliderDashPot软求模型，请问您知道怎么才能调用src/lagrangian/DSMC/submodels/bindaryCollisionModel/LarsenBorgnakkeVariableHardSphere这个可变硬球模型么？