MPPICFoam中如何输出某个颗粒在各个时间步的位置信息？

[fubianhanshu]

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

[fubianhanshu]

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

[fubianhanshu]

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

[upc_ngh]

您好，请问上面那个固定ID粒子的纵坐标图怎么做出来的？可以说的在详细一点吗？

[upc_ngh]

@fubianhanshu 前辈您好请问这个粒子高度变化图可以给指点一下吗？跪求 QQ:1514969704希望得到前辈的指点

[星星星星晴]

hi 关于粒子编号什么的 在of中会有origid 和origProcID两个东西来显示该粒子的编号，但是问题是，如果用并行计算，这个东西是怎么走的，目前我是不太了解。 如果你有特定的parcel需要追踪，在kinematicparcel 的move中 单独写一个输出，选择好特定的parcel，然后输出就好了，但是并行这个真的很头大。。希望给你一点思路

 /*   std::ofstream outfile;
    outfile.setf (ios_base::fixed, ios_base::floatfield);
    outfile.precision (7);
    outfile.open ("Injectior_Magnitude_kinematic_parcel", ios_base::app);
            << p.position().component (0)           << tab
            << p.position().component (1)           << tab
            << p.position().component (2)           << tab
            << mag(p.U())                  << tab
            << p.U().component (0)                          << tab
            << p.U().component (1)                               << tab
            << p.U().component (2)                              << tab
            << p.d()*1e6                   << tab            

            << nl;*/

[upc_ngh]

好的得到数据之后再自己作图就好了，差不多明白了，谢谢前辈！

[bestucan]

paraview 工具栏 edit 里有个 find data，可以用”与或“逻辑添加条件过滤目标信息。

[fubianhanshu]

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

[fubianhanshu]

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

[星星星星晴]

@fubianhanshu so, 这就是个问题，因为不太懂具体是怎么追踪的，我的想法是在injection model中，直接给parcel 添加一个maker， 然后在paraview中 单独输出这一个粒子就能得到你想要的信息了，就是费劲点，或者就是在move中添加一个输出这样做。

[upc_ngh]

@fubianhanshu 好的我自己试一下，谢谢各位前辈 @星星星星晴 @fubianhanshu