Euler-Lagrange的一些解析，sprayFoam
-
@星星星星晴 您好，谢谢回复
[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文件，感觉是没有调用新编译的库文件函数？您的编译方法参考教程能方便发一下么？
第三个问题，我用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; }
很明显，粒子已经通过我设计的多边形，但是质量流量等于0，这又是什么原因呢？
-
-
我没有试过resetOnWrite 这部分，所以我也不知道你这么弄有没有问题。我仅仅是使用output的文件作后处理，套用了particlecollector这个cloudfunction而已，对于其他的东西我并不在意，而且之前总出错，与其用一个无法确定是否正确的东西，不如原始一点自己处理。
-
编译教程你可以在google上找很多，最主要的是你编译了一个东西以后，要对其重命名，然后相关引用到该lib的都要引用你新编译的模型，仔细看了一下，你在solver部分是有问题的，但是你编译了spray这个文件了么？编译的新spray有没有用到新的intermediate？solver没有引用新的spray。
3.如果你像我说的删除掉那部分code的话，自然没有质量的计算啊，所以后处理就是自己写个python script或者你用matlab或者什么后处理文件，因为你文件里面包含了所有hitpatch的数据。
游荡
-
-
- 额，这我就不是很了解了，我也不是code大神，仅仅是使用推荐的方法进行编译，这肯定不会错，正常情况下，建议无论更改了 什么，最好是重新编译一下相关的东西，以免发生错误。有些东西没办法解释。。。就是这么做，也不会思考为什么。
你修改原始包也没问题，但是不用重新命名。
但是如果有些东西你自己编译错了，很可能会引发连锁反应，这也是为什么人家推荐用的是$USER_LIB_SRC$LIB_SRC不是在人家原始的包里面做修改的原因，这样用以区别你自己的东西和原始的东西。
另外并不是所有人都有root的权限的。
我个人觉得还是人家推荐的方式比较好。- 这我就不是很清楚了。
-
@星星星星晴 您好，第三个错误我找到了，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
dsmcFoam吧
我搜了一下，希望有所帮助，如果有条件请尽量使用某G搜材料，百度太废物下面是有人写的解析
https://zhulianhua.github.io/2016/10/07/of-dsmcFoam/这还有一片代尔夫特大学的硕士毕业论文
https://repository.tudelft.nl/islandora/object/uuid:b0922fd5-9472-4037-9b9c-80af46736ad9/datastream/OBJ/download -
-
@星星星星晴 我尝试着调用dsmc的cloud，但是在UEqn.H中加入粒子对流体的作用力时无法成功添加，图片上传一直失败，我手打一下报错信息。如果参考GOFUN那个，直接 +invrhoInfdsmcCloud.SU(U) 报错为：expect primary-expression before'.' token +invrhoInfdsmcCloud.SU(U) ，如果吧dsmcCloud改为dsmc 在UEqn后面加+invrhoInf*dsmc.SU(U)，报错为has no member named 'SU',请问您知道该怎么把粒子力加到UEqn.H里么
-
@vbcwl 添加该项主要是为了2-way coupling，另外可能调用dsmccloud并不是这样的，请参考之前给你的那个解析，每个cloud调用的方法可能不一样，你可以参考一下其他的cloud是怎么被调用的。
看看其他求解器Info<< nl << "Constructing dsmcCloud " << endl; dsmcCloud dsmc("dsmc", mesh); Info<< "\nStarting time loop\n" << endl; while (runTime.loop()) { Info<< "Time = " << runTime.timeName() << nl << endl; dsmc.evolve(); //<<------ 这里是演化步 dsmc.info(); //<<------ dump到屏幕一些基本信息 runTime.write(); //<<------- 这里很关键，到底写的是什么 Info<< nl << "ExecutionTime = " << runTime.elapsedCpuTime() << " s" << " ClockTime = " << runTime.elapsedClockTime() << " s" << nl << endl; } Info<< "End\n" << endl;游荡
-
@星星星星晴 您好，您推荐的2. 输出文件方法我基本弄明白了，特别适合输出一个平面上不同时刻的统计数据，然后用python数据处理。
我遇到的问题主要是输出三维规则网格空间点上的水滴质量，把这个数据导入到别的程序中。我尝试过用paraview输出粒子所在空间点数据，因为很不规则，所以之后用python三维插值，但是效果很差。
所以我觉得您提出的1. 将拉格朗日场转到EULER场比较适合我的情况。但是关于这个代码使用还想请教您几个问题：
1.template<class CloudType> Foam::volScalarField& Foam::IVT8_StochasticCollisionModel_2022<CloudType>::PPC() { if (!PPCPtr_.valid() ) { const fvMesh& mesh = this->owner().mesh(); PPCPtr_.reset ( new volScalarField ( IOobject ( this->owner().name() + ":PPC", mesh.time().timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE ), mesh, dimensionedScalar ("zero", dimless, 0.0), zeroGradientFvPatchScalarField::typeName ) ); } return PPCPtr_(); }这段代码应该是定义了一个函数，应该放在求解器的哪个文件中呢？
或者是自己重新定义一个头文件？-
这个函数中的PPC应该是一个通用模板吧？比如我想输出 p.mass()质量的数据到EULER场，这个PPC直接用mass替换？
-
调用的时候您说用下面的代码，我想知道这段代码又是放在对应的哪个文件中呢？
scalar& PPC = PPCPtr_->primitiveFieldRef() [cellI]; PPC += 1.0; -
-
@chengan-wang
同样可以放到cloudfunction中，对应的改名什么的
第一步写这个function
第二步在constructor里面调用，生成这个field
第三步在你需要的地方调用这个方法比较复杂，需要多个地方定义，我这仅仅是提供一个思路，可以在任何地方生成这个field，只要在计算的时候能够赋值即可，你也可以在你的solver里面生成这个field，只要在拉格朗日部分能调用到这个field就行。
为什么要插值？我没get到你的点。插值也有很多方法的，另外当你的数据量不够的时候自然会存在很大的误差，这是统计学的问题了。
方法一天然的会损失位置信息，因为网格中储存的数据都是在网格中心，如果做复杂的运算的话，这样会损失掉很多信息，anyway这是你的计算，你可以试试。
我个人觉得方法二足矣应对所有的统计问题，无非就是后处理计算复杂一点，无法paraview中直接连续性的可视化
游荡
恨不得送你3个本子