拉格朗日粒子如何计算一个面上的SMD?
-
各位老师们好,请较一下如何获得一个面上的SMD呢,我现在的思路是:
(1)根据粒子与网格的空间位置关系,将粒子直径分布转化为欧拉场,插值到面上求平均
(2)直接提取这个面穿过的网格,把这些网格内的粒子直径求平均
不知道各位是否有更好的思路? -
@yfclark cloud function 的 particleCollector?
-
谢谢,我看了一下不太一样,还是准备自己写个python小脚本处理
-
@yfclark 请问关于拉格朗日粒子直径的分布您现在有什么好的建议吗?已经处理出来直径在面上的分布了吗?
-
就是用python写的小程序,统计空间位置在区间x和x+deltax里面的液滴,Delta x我用的是一个网格的宽度
-
@yfclark谢谢回答,我还是有点困惑, 第一步要知道空间粒子的直角坐标系坐标(这个地方用您在另外一个帖下回答的来解决进行转换?),第二步再自己写个小程序进行转换,是这样的吧。问题是第一步的转换我刚接触,不太明白去怎么操作,有参考资料看一下,或者稍微再讲明白一点吗?万分感谢
-
bluecfd的小程序:https://github.com/blueCFD/lagrangianExtraFunctionObjects
我没用这个,但是你可以试一试
我写的脚本可以看看smd处理:https://github.com/graceyangfan/openfoamPostProcess/blob/master/cloudpost.py -
@yfclark 谢谢,我认真学习一下
-
@yfclark 您好,看了您提供的网址决定还是退回到了of4,因为他这个bluecfd提供的小程序我看了一下好像没有在READMD中提到怎么去使用,就安装了bluecfd4可以直接得到直角坐标系下的position,还想再叨扰一下,您这个脚本怎么在paraview下加载,我找了一晚上资料也没太看懂,怎么加载
-
你换一下github的版本选项,可以从master换成of5x系列,另外脚本是在Python环境下运行的,不用paraview
-
脚本是自己编写的只是提供一个参考,另外你可以学习一门脚本语言,Python和MATLAB是比较好的选择,学了不亏
-
@yfclark 非常感谢您的回答,MATLAB是会用的,但是脚本这个东西怎么去用以前没接触过,我琢磨一下,谢谢。
-
@yfclark 真是太不好意思,又来打扰您,今天我又尝试了bluecfd的方法,按照他给的步骤,wmake完成后把库也链接到算例下的controlDict了,按照他的描述应该在每个时间步下/lagrangian/kinematicCloud`内出现一个position.orig但是无论我怎么尝试,也都出现不了,终端显示也是读取了这个命令了
writeCloudOldStyle writeCloudOldStyle1: applying to clouds: kinematicCloud
这中间哪里出了问题呢???