不规则多孔介质的网格划分

[hongjiewang]

多孔介质如上图所示，我想给灰蓝色区域也进行网格划分，并将两个网格结果耦合在一个网格内，我需要什么样的操作，有提到用createPatch，是把两个交界面耦合在一起？有没有详细一点的介绍资料或者类似算例可以提供参考。感谢各位~

[李东岳]

你看看传热那面的snappyMuyltiRegionHeater适合你的要求么，这是一个多计算域求解器

[veen]

@hongjiewang 空的区域是多孔介质吗？之前在fluent处理是分成不同的区域，交界面设置interface，或许可以考虑用商软划分网格再导入openfoam。
有个小问题想请教：在openfoam里边多孔介质是直接加fvoption计算还是有自己的求解器啊？

[hongjiewang]

@veen 在Icem里面对不同区域划分，之后导入OpenFOAM会针对不同区域产生多个文件夹，比如会有solid和fluid两个，这种是不是只适用于类似cht的求解器，本身也对应两个区域的求解，而我的求解器是把固体作为一个单独的相，只有一个控制方程进行求解的。
我的求解器是基于针对多孔介质的控制方程进行编写的，相变方面有参考interCondensatingEvaporatingFoam这个求解器。

[hongjiewang]

@李东岳 感谢东岳老师的回复~我的求解器只有一套控制方程，不用单独分出来固体区和流体区，是直接一起求解的，类似porousSimpleFoam求解器对多孔区的处理。分区主要是一些参数设置的不同，对于一些简单几何模型直接采用topoSet是可行的，对于不规则的还是需要借助stl进行区域的设定。我在网上也找了一些方法，感觉下面这个对我的求解器比较可行，我现在困惑的是如何把两个stl文件划分的网格耦合在一起，不清楚有没有这样的算例。

[zousiyu]

你可以看看这个意大利小组的文献，说不定有帮助。他们自己写了个求解器，专门用于求解固定床的非均相催化现象。系统内同时存在自由流动区域和多孔介质区域。
[1] T. Maffei, G. Gentile, S. Rebughini, M. Bracconi, F. Manelli, S. Lipp, A. Cuoci, M. Maestri, A multiregion operator-splitting CFD approach for coupling microkinetic modeling with internal porous transport in heterogeneous catalytic reactors, Chemical Engineering Journal 283 (2016) 1392-1404.

[wwzhao]

@hongjiewang 在 不规则多孔介质的网格划分 中说：

对于一些简单几何模型直接采用topoSet是可行的，对于不规则的还是需要借助stl进行区域的设定。

你在另一个帖子里问我了，我在这里回复你吧。

在 snappyHexMeshDict 里面可以通过 STL 文件（必须封闭）指定 cellZone。例如：

geometry
{
    sphere
    {
        type triSurfaceMesh;
        file "sphere.stl";
    }
};

refinementSurfaces
{
    sphere
    {
        level (0 0);           
        faceZone sphere1;
        cellZone sphere1;
        cellZoneInside inside;
    }
}

[hongjiewang]

@wwzhao 嗯嗯谢谢老师。我需要在不同的case下对骨架区域和流体区域进行相同的操作，之后利用mergeMesh耦合网格。是这样理解的吗

[hongjiewang]

@zousiyu 嗯嗯好的 谢谢老师~我会认真阅读的 有问题再请教

[wwzhao]

@hongjiewang 不需要merge，只需要一个 case 文件夹就行。

[hongjiewang]

@wwzhao 感谢回复~ 这种两个区域都需要进行网格划分也不需要mergeMesh吗，是说把同一个stl添加两次并命名不同，导入进来之后，在snappyHexMesh里面一个取inside，一个取outside？

[wwzhao]

@hongjiewang 不需要对两个区域分别划分，利用 sHM 可以直接同时生成一整套网格。

[hongjiewang]

@wwzhao 嗯嗯好的 明白了 谢谢老师~

[hongjiewang]

@wwzhao 老师 我今天尝试了一下，已经成功生成。但是不确定我的操作是不是对的~麻烦老师帮我看一下。我的stl文件是一个中间的圆，我在refinementSurfaces中设置cellZoneInside选择inside；然后在refinementRegions下面的locationInMesh选择的是圆外一点。我是这样设置的，最后生成两个区域共存的网格，中间圆的网格比较密了。

感谢老师~

[wwzhao]

@hongjiewang 看起来应该没问题，中间的网格密是因为用refinement region加密过了。

[hongjiewang]

@wwzhao 嗯嗯是的 中间是加密之后的~ 非常感谢老师 要不然我要走弯路了。