decomposePar 的手动分块/手动指定计算区域分块

[hurricane007]

但是我试了下，运行出了个错误，感觉搞笑得一塌糊涂，找不到constant文件？？不是应该找的是constant/cellDist么？

--> FOAM FATAL IO ERROR:
cannot find file

file: /home/he/catalyticFOAM/case/cloneCaseStudy/constant at line 0.

    From function regIOobject::readStream()
    in file db/regIOobject/regIOobjectRead.C at line 72.

FOAM exiting 

[李东岳]

忒复杂了 有点爱莫能助

[浪迹天大]

邱老师的这篇博客也讲了这个方法，我按此操作过，没遇到问题。

[李东岳]

http://www.cfd-china.com/topic/231/openfoam并行计算如何手动划分区域/2
只不过文件名不同？没细看

[bestucan]

包围同样的面积，圆形周长最短。但是圆形不能无缝契合。
所以蜜蜂做蜂巢用的是接近圆形的正六边形。
接触面积最小分区方式应该在正六边形的基础上变。考虑计算区域、区域和边界条件的关系、网格大小不均匀，变成啥样就看经验了。

找不到文件这个错误，可能是当前路径有问题
python脚本不是把celldist写道constant了么，为什么应该找不到它？
ls看一下constant是不是被复写了。python脚本里，写文件那部分，看着挺别扭的

[bestucan]

一步一步来嘛，先手动操作成功了，再写脚本。这样可以定义问题是在操作思路上还是在自动化脚本上。

[hurricane007]

@浪迹天大 对的，邱老师这篇博客讲的是一样的。
@东岳 对，和你这个也是一样的。
@bestucan 包围面积的话，是不是scotch考虑的是通信最少，所以采用的是边界网格数最少的方式？我去重新手动不要用python脚本试试，找不到constant文件我真是惊呆了

[hurricane007]

@bestucan update一下，的确手动复制就可以了，用python代码就不对。我得找找那个代码哪儿的问题

[bestucan]

@hurricane007 我说蜂窝只是不加其他考虑的回答你的那个问题，怎么样让接触边界小。
我觉得通信量的大小倒不是很影响并行速度，主要是区域间的互相影响。

---

比如如果有个长管喷嘴，照着喷嘴中心竖切分成两个区域。这样喷嘴左边和喷嘴右边都可以各自影响不是很大的一起计算。
如果是横切两个区域，一个包含喷嘴，一个不包含喷嘴。那没有喷嘴的区域在初始时间没有什么计算量。因为还没喷过来。

---

这是个很糙的例子，最恶劣的状况，网格分区，犬牙交错。那么计算过程更像是核之间的接力赛，你方唱罢我登台。那并行的优势是一点点都没有了。

所以要避免这样的状况，我觉得要结合边界条件对流域有个定性的预估。至于scotch的工作原理我是不懂的

[hurricane007]

@bestucan 对的对的，我也觉得是这样，那个“很糙”的例子其实很好，当计算区域内计算负载不均衡的时候，均匀的分区其实并不能均匀分配负载在各个线程上，我想做这个手动分区其实就是想搞这个，把计算负载最大的catalytic 边界分在更多的核上算。
BTW，既然blockMeshDist里面能指定blcoks了，意味着在blockMesh里面已经是分块划分网格了，有没有一种方式直接在decomposeParDist里面指定每一个blcoks的分区方式呢？

[hurricane007]

再update一下python的代码，不知道哪个为啥错，但是重新写了下执行对了。但是响了下应该是可以写个script执行的，如果电脑上没python。。。

#read file
f1 = open("0/cellDist", "r")
lines = f1.readlines()
f1.close()

#replace and delete unwanted lines
lines[11] = lines[11].replace("volScalarField", "labelList")
lines[12] = lines[12].replace("0", "constant")
del lines[17:21]
bfline = 0
for i, line in enumerate(lines):
    if (line.find("boundaryField") != -1):
        bfline = i
        break
del lines[i:]

#write file
f2 = open("constant/cellDist", "w")
for line in lines:
    f2.write(line)
f2.close();

[李东岳]

感谢诸位大神分享！

[bestucan]

@hurricane007 比block还简单粗暴，横着分两份，竖着分两份。第三个维度没有。

[hurricane007]

@bestucan 如果我在blockMeshDict里面定义了两个blocks，然后我想把第一个block分成两快，第二个block分成6块，有简单的实现方法么。。

[bestucan]

@hurricane007 我的实际操作经验还是挺少的。没见有可以把block和分区结合一起的方法，毕竟这两个分的依据都不太一样。有个方式是manual，看起来很适合你的样子。直接指定哪个网格属于哪个分区。但没见过相关例子。

[hurricane007]

@bestucan
这个就是我这个帖子说的。。。。

[bestucan]

@hurricane007 那这种方法挺好的
如果手动加工去关联block和分区。可以用写脚本，bash、python都行。读constant/polymesh里的数据，计算生成celllist。
这个工作量并不是很大，就是逻辑理着麻烦。
linux一切皆文本的好处，什么信息都可读，操作方法通用，容易嫁接。

[hurricane007]

@bestucan 其实这个方法是用 setFields 先赋值一个cellLabel 场，每个点的值为想分的区，再把这个场给转化成manualDist，然后再用decomposePar 根据这个manualDist 分区了。
但是我觉得可能自己分区不一定有算法分的好，那么最好就是手动分成几块，然后算法再把每块分成几小块

[bestucan]

@hurricane007 这个思路挺不不错。我没见过手动分区的例子，不懂了 加油。