算例topoSet、refineMesh后，无法decomposePar

[hy1112006]

最近在试运行一个算例，整体计算区域是标准六面体，结构化网格（网格单元亦为六面体）：
（1）串行运算：blockMesh --> checkMesh --> XXXFoam(运行求解器)，没有问题；
（2）并行运算：blockMesh --> topoSet --> refinMesh --> checkMesh 至此没有问题，但是开始分解网格时decomposePar，出现错误，过程如下；

--- blockMesh如下（无问题）：

--- topoSet -constant 如下 （无问题）：

--- refineMesh -overwrite 如下（无问题）：

--- decompossePar 如下，报错！！：

请问各位专家，这个有可能是什么原因呐，感谢！

[wsxfyy]

@hy1112006 我之前这种做法分块的时候并没有报错，也能正常运行。你要不换一下分块方式，不要用这种分块方式，用scotch试一下

[hy1112006]

@wsxfyy 我看了下，这个算例就是用scotch 分解的 只有看看以后有没有啥方法了

[xpqiu]

@hy1112006
先 refineMesh，后 topoSet 试试。
topoSet 生成的 cellZone 或者 cellSet 等，保存的网格ID是你refineMesh 之前的， refineMesh 之后，网格数量都变了，所以保存在 cellZone 或者 cellSet 里面的编号跟refine之后的网格都不匹配了。

[hy1112006]

@xpqiu 刚才试了下，如果先refineMesh的话，会报错：

这又是什么原因呐，疑惑~~

[xpqiu]

@hy1112006
哦，你的 refineMeshDict 里面需要一个 cellSet 来定义需要refine 的网格。这个 cellSet 也只是用来起这个作用吧。
所以，你需要在 refineMesh 之前，先 topoSet 把 cellSet 生成出来，然后 refineMesh
但是，在 decomposePar 的时候要排除对 cellSet 进行 decompose，因为我上一条回复说的原因。可以给 decomposePar 加一个选项，-noSets，这样在 decomposePar 的时候就不读取 cellSet 了，也就不会再触发你主楼遇到的错误了。

[hy1112006]

@xpqiu 好的，已经尝试成功，感谢老师指点啊！