OpenFOAM在并行运算结束后 能否并行reconstructPar？或者有没有其他方法能够加快reconstructPar？

[Stan同学]

最近跑了几个网格总数比较大的案例（大概500w~1000w），并行计算结束后，执行了reconstructPar进行重构，但是重构时间特别长。

想请问各位老师，进行重构时能否像运算时那样并行呢？或者有没有其他方法可以加快重构的速度呢？
（听同学说可以给服务器指定一部分线程用于重构某一部分，再指定另外一部分线程去重构另一部分，请问这种办法能够实现吗？）

[李东岳]

mpirun -np 2 redistributePar -reconstruct -parallel

试试这个

[Stan同学]

@东岳 东岳老师，昨天我试了一下 好像是不行

终端显示如下：

FOAM FATAL ERROR:
Invalid option: -reconstruct

FOAM parallel run exiting

MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD with errorcode 1.

[Zephyr]

用paraview的话是不用重构也可以直接后处理的

[Stan同学]

@Zephyr 大佬 我一般只是用paraview来确定一下网格和各种场的设置是否无误 并且paraview后处理的话一般是绘图吧（用paraview用的比较少 不知道对不对 ）

但我这边是想要把计算的某些结果提出来（例如速度和压力） 应该还是要先重组一下的吧？

[Zephyr]

https://www.cfd-online.com/Forums/openfoam-post-processing/64652-bash-script-pseudo-parallel-usage-reconstructpar.html

你一定要并行reconstructPar的话可以试一试这个

[Samuel-Tu]

印象中以前看一个培训材料说分块和重构都只能单核

[Stan同学]

@Zephyr 谢谢！我再用这个尝试一下

[Stan同学]

@Samuel-Tu 前些天邮件咨询了一下Jozsef Nagy教授 他的回复也是说重构只能单核

[李东岳]

https://www.openfoam.com/releases/openfoam-v3.0+/parallel.php

你看看这个并测试一下，感谢反馈

[Zephyr]

reconstructPar的速度真的挺迷的，也不知道到底怎么实现出来的，慢的离谱

[hurricane007]

@Samuel-Tu 我记得我也仔细研究过，重构的确是只能单线程。。。

[Stan同学]

@东岳
东岳老师您好，我用一个tutorial中的小案例测试了一下，如果用mpirun -np 2 redistributePar -parallel，则在进行并行计算时会报错，报错内容是找不到processor中的初始文件p_rgh。

如果并行计算前不用mpirun -np 2 redistributePar -parallel的话，照常decomposePar然后并行计算，最后用mpirun -np 2 redistributePar -reconstruct -parallel还是会报错，报错内容依然是Invalid option: -reconstruct。

[Stan同学]

@Zephyr
刚刚试了一下，确实感觉速度变快了

[strive_DUT]

@李东岳 李老师，您知道Zephyr同学分享的这个脚本如何使用吗？
链接文本

[ZY_Feng]

@strive_DUT 就是执行脚本，第一个参数输入核心数，可以看我写的知乎文章。https://zhuanlan.zhihu.com/p/636421703

[李东岳]

openfoam的decomposePar、reconstructPar是个老大难。这个应该需要大幅度的fund他们才会解决。目前他们对这方面还没有迫切的处理需求。目前openfoam针对reconstructPar这面的处理方式是后处理也是并行来搞，paraview直接看decompose的数据，各种postProcess也都是并行来。针对decomposePar慢的问题，openfoam这面是画网格snappyHexMesh直接并行来。所以我这面尝试处理2个亿网格，从生成到计算，也还可以。

但是如果是第三方生成的2个亿网格，需要decomposePar，那可就慢了。

[Tens]

感觉OpenFOAM的并行计算真的不太友好，decomposePar和reconstructPar速度慢是一方面，并行还会成倍的增加文件数量，经常因为小文件过多被管理员警告

[李东岳]

@Tens 商软那面是怎么处理的。我看商软好像就是一个大文件。我没用过商软就难以理解咋回事。难道是打包压缩了么。。

[Tens]

@李东岳 是一个大文件，不过这样又不好对数据后处理。OF这面如果不需要decomposePar成processor和reconstructPar，并行和单核一样计算就挺好的