拉格朗日粒子如何计算一个面上的SMD？

[yfclark]

各位老师们好，请较一下如何获得一个面上的SMD呢，我现在的思路是：
（1）根据粒子与网格的空间位置关系，将粒子直径分布转化为欧拉场，插值到面上求平均
（2）直接提取这个面穿过的网格，把这些网格内的粒子直径求平均
不知道各位是否有更好的思路？

[星星星星晴]

@yfclark cloud function 的 particleCollector？

[yfclark]

谢谢，我看了一下不太一样，还是准备自己写个python小脚本处理

[upc_ngh]

@yfclark 请问关于拉格朗日粒子直径的分布您现在有什么好的建议吗？已经处理出来直径在面上的分布了吗？

[yfclark]

就是用python写的小程序，统计空间位置在区间x和x+deltax里面的液滴，Delta x我用的是一个网格的宽度

[upc_ngh]

@yfclark谢谢回答，我还是有点困惑， 第一步要知道空间粒子的直角坐标系坐标(这个地方用您在另外一个帖下回答的来解决进行转换？)，第二步再自己写个小程序进行转换，是这样的吧。问题是第一步的转换我刚接触，不太明白去怎么操作，有参考资料看一下，或者稍微再讲明白一点吗？万分感谢

[yfclark]

bluecfd的小程序：https://github.com/blueCFD/lagrangianExtraFunctionObjects
我没用这个，但是你可以试一试
我写的脚本可以看看smd处理：https://github.com/graceyangfan/openfoamPostProcess/blob/master/cloudpost.py

[upc_ngh]

@yfclark 谢谢，我认真学习一下

[upc_ngh]

@yfclark 您好，看了您提供的网址决定还是退回到了of4，因为他这个bluecfd提供的小程序我看了一下好像没有在READMD中提到怎么去使用，就安装了bluecfd4可以直接得到直角坐标系下的position，还想再叨扰一下，您这个脚本怎么在paraview下加载，我找了一晚上资料也没太看懂，怎么加载

[yfclark]

你换一下github的版本选项，可以从master换成of5x系列，另外脚本是在Python环境下运行的，不用paraview

[yfclark]

脚本是自己编写的只是提供一个参考，另外你可以学习一门脚本语言，Python和MATLAB是比较好的选择，学了不亏

[upc_ngh]

@yfclark 非常感谢您的回答，MATLAB是会用的，但是脚本这个东西怎么去用以前没接触过，我琢磨一下，谢谢。

[upc_ngh]

@yfclark 真是太不好意思，又来打扰您，今天我又尝试了bluecfd的方法，按照他给的步骤，wmake完成后把库也链接到算例下的controlDict了，按照他的描述应该在每个时间步下/lagrangian/kinematicCloud`内出现一个position.orig但是无论我怎么尝试，也都出现不了，终端显示也是读取了这个命令了

writeCloudOldStyle writeCloudOldStyle1: applying to clouds:
 kinematicCloud

这中间哪里出了问题呢？？？