CFD中文网

看很多文献中都是将计算域的长、宽、高和二分之一的高度$\delta$做比值等于多少倍的$\pi$，这其中的原理是什么？

@fubianhanshu
在particleIO.C下添加一行代码，注释原先的输出可以输出颗粒全局坐标。

os  << position() << token::SPACE << celli_;

先采用最笨的方法，把保存时间步设置小，遍历读取坐标，其他的再慢慢实现，看着lagrangian粒子程序太复杂了

@星星星星晴 谢谢哈，这个昨天我也看到了，我再看看他的这个实现思路。另外应该还有两种思路，一是：src/functioObjects/lagrangian/cloudInfo可以在计算过程总输出颗粒质量随时间的文本信息（这正是我想实现的效果），应该也可以将particle的坐标输出的程序代码到这里新建类，但是对于底层代码理解的不透，尝试了还没编译通。二是：至少v1706版本以后的functionObject/field/中都有particleDistribution库函数，也可以将它改写到of版本里实现，这个工作应该更大些，涉及太多的依赖头文件和关键字写法的不一样。

DPMFoam中将是不是可以实现颗粒位置信息按照时间顺序输出到postProcessing中的一个文本中？正常情况颗粒位置是保存在每个输出的时间步文件夹下有positions文本中的，我计算颗粒比较少，想要输出颗粒连续时间的位置，若是是每个时间步保存文件大小太大，且后期通过shell脚本遍历每个结果路径下的positions内容写入单独文本中比较麻烦。是不可以改写src/lagrangian/basic/IOPosition中的远程序实现，不知道具体怎么改，有前辈可以指点下吗？

@李东岳 这个实现了就能很好的解决关于颗粒的后处理问题了。一种是要想知道某个颗粒归属哪个网格就一直追踪这个颗粒的坐标，然后再确定网格；另一种是要知道某个网格内有哪几个粒子就像你说的遍历一遍网格确定粒子，会很慢。但可能还会遇到下面的问题，我在做CFD-DEM的时候在paraview中选取的粒子发现有两个编号，一个是：ID，一个是id，而且在多核计算的时候发现id会在跨越block时进行重新编号（后处理显示时），DPMFoam中应该不存在。

@sungda 好的，谢谢哈。我使用的5.X版本，跟liggghts耦合用，想要给定随距离变化的速度梯度作为初始场，我安装试下，有问题了向你多请教

@sungda 你好，请问你是在openfoam哪个版本实现的呀

@李梦瑶 粒子可以循环，你试下就知道了。

目前是用的MPPICFOAM

@星星星星晴 我都是合成后处理的，选定某个粒子，整个时间序列ID号应该是不变的，多核分割后怎么处理界面的也不太清楚，我注意到分割后并行计算时，每个块内的粒子是动态变化的

@upc_ngh 选完后显示全域，装上Python后，用find找也很方便，就像楼下道友说的

@bestucan 你是行家呀，厉害厉害，终于知道原因了

要是虚拟机的话，装一个版本临时看看结果，处理下，后期处理我感觉转成VTK格式在Windows版本上流畅很多，速度也快

@fubianhanshu 上面的图没输出好，贴个清晰点的

@fubianhanshu 自己先回复下吧，在paraview中可以处理，我先贴出来一个图，某个固定ID粒子的Z坐标随时间的变化曲线图，具体的还在摸索，不是很清楚，希望得到大家更好的建议

@zhangxc0223 我不清楚你怎么做的，不过我刚发现可以在paraview中可以显示粒子ID，随着时间步推移，好像每个粒子的ID号是不变的，这样就可以通过后处理输出某个固定ID的各个时间步的信息了，可以交流下

楼主你这个问题解决了吗？

欧拉拉格朗日方法的MPPICFoam求解器计算的颗粒输运，该怎么输出或处理粒子的位置信息，比如:给每个粒子编号，追踪某个固定编号（id=100）的粒子在每个时间步的空间位置、速度？纠结纠结······ 跪求懂得前辈指点

@Jacob-Z 在 CFDEM安装 中说：

@fubianhanshu CFDEM®coupling 3.8.0=LIGGGHTS® 3.8.0+OF 5.x

好的，谢谢！我也尝试下，以后一起交流