DPMFoam追踪算法失效导致颗粒丢失

[李东岳]

@浪迹天大 个人喜好问题很久没装过paraFoam了， @Alvin 或许你可以装一下paraFoam。另外，重心坐标这个我没研究过，不过很感兴趣。本质上，脱离CFD，也是很有意思的一个东西。如果你有空，可以详细研究研究，分享给我们。我可以给你在CFD界推一把。

[Alvin]

@浪迹天大 谢谢您提供的这个思路。

[Alvin]

@东岳 重心坐标只是第一步，理解of50那一整套新的颗粒追踪算法满有难度的，看能否找到相关文献吧

[dzw05]

@浪迹天大 我用了paraFoam为什么还是不能显示粒子呢？是直接输入paraFoam命令吧？

[浪迹天大]

@dzw05 对呀，直接使用 paraFoam 打开就能看。

• 使用 paraFoam 打开之后，需要勾选需要的信息的，如下图：
  

• 如果使用 touch a.foam && paraview a.foam，则会报错：

ERROR: In /home/ubuntu/OpenFOAM/ThirdParty-dev/ParaView-5.4.0/VTK/IO/Geometry/vtkOpenFOAMReader.cxx, line 8285
vtkOpenFOAMReaderPrivate (0x417d3d0): Error reading line 20 of /root/OpenFOAM/OpenFOAM-dev/tutorials/lagrangian/sprayFoam/aachenBomb/5e-06/lagrangian/sprayCloud/positions: Expected punctuation token ')', found 0

[Alvin]

@dzw05 你用的哪个版本的of？of50下的paraFoam能够读取质心坐标的

[dzw05]

@浪迹天大 可能因为我用的blueCFD，在windows系统中运行，所以即使用paraFoam也会出现报错：(。

[李东岳]

理解of50那一整套新的颗粒追踪算法满有难度的

这是个有意思的研究方向，算法 + 测试 + 应用，可以试试Journal of Computational Physics，Computer Physics Communications，或者可以投个powder technology，Computers Fluids。但这种东西往往需要一定时间去搞。受雇人员（包括硕、博、以及科研高校的课题组成员）往往不能独立的去研究这东西。没时间，雇主也不好交代。

只有某些特殊的人会有时间、有精力去搞，比如我 哈哈。在国内，只要受雇于人，就缺乏独立性，太高调会引火烧身，当时比较年轻，但现在总不能把我的网站名改成“无名流体”吧？哈哈

[dzw05]

@alvin of5.x，是blueCFD版本的。

[Alvin]

@东岳 说的太对了，确实需要花费专门的时间去搞，在企业有限制。做研究就要沉得住气。openfoam-v1712 ....../etc/bashrc文件中有颗粒位置记录方式的选择，也可以调用程序将粒子的质心坐标转换为空间直角坐标：

openfoam1712记录颗粒位置的文件是“coordinates"，通过将”writeLagrangianPositions“ 值置为非零，结果文件夹下就出现了熟悉的”positions“ ：

[Alvin]

@dzw05 blueCFD还没有用过

[Alvin]

@东岳 您可以考虑做培训之外的企业技术咨询服务

[李东岳]

@alvin 在 DPMFoam追踪算法失效导致颗粒丢失 中说：

您可以考虑做培训之外的企业技术咨询服务

一旦公开，基本自断后路 嘿嘿 https://www.zhihu.com/question/33445584

东岳流体logo中间的D表示CFD中的有限控制体，也表示物质导数，左右的弧线表示CFD中的流线。本网站所有内容由东岳流体®出品，不借助（依附）任何组织和单位。各平台恕不承接任何广告及项目。CFD学术合作（非咨询）可联系li.dy@dyfluid.com，问题咨询请前往CFD中文网。

[upc_ngh]

@Alvin 您好，请问有什么好的办法将of5.0下的粒子坐标转化成空间坐标吗？

[yfclark]

@Alvin blue cfd有一个小程序在GitHub开源的，可以用于openfoam5.0,可以去搜索一下，不过不太好用，只能得到各个核的数据，我喜欢reconstruct后再用，它的函数没法重构

[dxl]

可是为什么我在openfoam10中还遇到了颗粒丢失的问题。。。

[sxz0823]

@dxl 您好，在你的 https://www.cfd-china.com/topic/394/如果用dpmfoam求解稀相流会怎么样-误差大么/47 这篇文章中我也遇到了相同的情况可以讨论一下吗