并行效率疑问

[random_ran]

@xpqiu

非常感谢你的建议， renumberMesh 在一个800万网格 96核的模拟中，时间缩短了40%。

以后多核计算前，我都会renumberMesh 的，也建议大家用上这个好用的功能，加速你的算例。:D

[CFDngu]

@random_ran 你好，我想问一下，我运用scotch方法剖分了网格，并行运算时提示网格出现下述问题是咋回事：

大概提示是找不到法向量，我对应到程序中是这个地方：

   103     if (mag(pn) < VSMALL)
   104     {
   105         FatalErrorInFunction
   106             << "Cannot determine normal vector from patches."
   107             << abort(FatalError);
   108     }
   109     else
   110     {
   111         pn /= mag(pn);
   112     }
   113

大概意思就是分块的面积太小了，找不到法向量？

用simple方法提示的错误不一样，大意几个地方网格的交界面有问题。

可是不管哪个方法，我都检查了网格，checkMesh没问题，单核也可以计算。

我怀疑是不是加密的关系，一般感觉加密的level到3就容易出问题，我这个就到了3，请问你遇到过类似的情况吗？

[random_ran]

@CFDngu

你好,CFDngu:

你这个问题我没有遇到过，decompose的时候没有遇到过报错的问题。可能是我的几何体比较简单的原因？

建议你把：

decomposeParDict 内容
$ checkMesh 结果
O.F. 版本

贴出来，或许有碰到过类似问题的朋友能帮助到你。

[CFDngu]

@random_ran 多谢建议，不过确实都没有任何问题。。。

等过一段我贴上来：）

[CFDngu]

@xpqiu renumber的功能是什么呢？我试着用了一下，发现结果完全不对啊。。。
不renumber的结果就对了

[qjh888]

好贴子，马住。另外推荐一个连接：
http://nscc-gz.cn/newsdetail.html?6151

[xpqiu]

@CFDngu renumberMesh 的作用是减小矩阵的 bandwidth，让非零元素更集中到主对角附近。
https://en.wikipedia.org/wiki/Cuthill–McKee_algorithm

[CFDngu]

@xpqiu 可是我的算例为什么renumber之后就结果就不对了呢？只有初始场是正确的，第二步开始就完全错了。和算例是2D的有关系吗？

[xpqiu]

@CFDngu 2D 也可以renumberMesh，具体是怎么不对，log贴上来看看。不太可能是renumberMesh 导致的，可能是其他原因。

[CFDngu]

好像真是有问题，最下面的。。。。

facesets出问题了吗？

[xpqiu]

@CFDngu renumberMesh运行完以后应该重新提取一下faceSet

[CFDngu]

@xpqiu :expressionless: 不好意思我没用过这个，请问您那里有这方面的资料吗

[bingningmeng45]

你是在哪个超算中心跑的算例？用的是什么规格的机器

[random_ran]

@bingningmeng45

我是在 Université de Sherbrooke 的长毛象2号超算上算的。这篇帖子的一楼最后我详细描述了硬件。

我没有在专门搞CFD的组里搞CFD，机时的分配只有100-CPU年。

[bingningmeng45]

@random_ran 在超算上安装哪个版本的openfoam？是自己安装的么？我现在跑算例要在超算上计算，现在还处在调研阶段，对超算的情况不是很了解。希望你能多给我一些建议

[bingningmeng45]

@random_ran renumberMesh这个命令你用过么？（我看贴上有人用过）用过的话你跑的是什么算例模型，计算效率提高了多少，计算结果准确么？我的算例也是圆柱绕流湍流模型验证 是否可用这个命令？

[赵一铭]

@random_ran
顶长毛象

@bingningmeng45
renumberMesh看你需求吧，锦上添花的东西，不用也可以，用了更好，莫非你并行遇到了瓶颈。

[random_ran]

@赵一铭

哈哈，这下可好了，每次登录长毛象的时候，都会自动脑补这张图了。 :joking:

[bingningmeng45]

@random_ran renumberMesh 这个如何操作运行？直接在终端输入吗

[程迪]

可以试试metis，可能会比scotch要好一些。