关于二维模型的计算，这两种snappyHexMesh方法有什么不同，应该采用哪一个？

[DY大世界]

@xpqiu 好的，谢谢老师的指点回复！

[DY大世界]

@史浩 好的，很感谢大佬的回复，我对比看看结果！

[DY大世界]

@xpqiu 你好大神，我想再问您一下，这两天我看了一些文献，有很多文章都是在使用snappyHexMesh后，没有再用extrudeMesh命令，而是直接把前后面的type设置成empty进行计算。但是仿真结果与试验的对比也挺好。

我想问一下这种使用snappyHexMesh后通过直接设置front和back的type为empty的方法，能否表示二维模型，这种网格是算二维的还是三维的。因为他们的计算结果精度挺好，是否表明这种方法对于二维模型是可行的？

[李东岳]

把文章贴上来看看，这应该是错误的

[DY大世界]

@东岳 感谢李老师回复！

例如这篇文章没有提到extrudeMesh，说的是将三维转换成二维只需要将边界指定为empty。

现在我的案例是关于波浪-固体耦合的，如果不用extrudeMesh还好，一使用计算就发散，不知道使用extrudeMesh后还有什么要求，比如setFields什么的。

[xpqiu]

@DY大世界 某个方向设置为 empty，那么这个方向的流场就不会再计算。以前的版本如果empty方向的网格不是固定层数，是无法计算的，checkMesh 也会报一个提示：

 ***Total number of faces on empty patches is not divisible by the number of cells in the mesh. Hence this mesh is not 1D or 2D.

经测试，现在新的版本可以计算，但是结果我没有验证过。从计算的角度来讲，即使这样算出来结果是合理的，也增加了很多没必要的计算，因为这样处理，网格数量会比单层的2维网格多很多。
另外，假设举个例子，一个方形计算区域，中心部分用snappyHexMesh进行了加密，然后x 方向设置为 empty ，那么你计算完的结果，在x 方向上不同位置的截面上流场可能是不一样的，因为网格数量可能不一样，这样显然跟二维计算的概念是不符合的。

[DY大世界]

@xpqiu 非常感谢大佬的回答，那不用extrudeMesh命令所得到的确实是不正确的二维网格。而且我对比了使用extrudeMesh与不使用extrudeMesh的计算案例，得到的结果也是不完全吻合。

[Samuel-Tu]

@xpqiu @DY大世界
我也遇到这个问题了。看网上说的是先createPatch，就是说把empty的一个面先提取出来，然后用extrudeMesh把这个patch沿着垂直方向推一个网格，形成三维的网格。但是我发现，沿着那个方向竟然不是一层网格。原因在于悬挂节点，因为shm的时候用了refine的设置，加密的交界面是悬挂节点，这导致在extrude的时候没法是一层网格，见下图
首先是我按上述说法做出的沿厚度方向只有一层网格的图：

然后提取加密交接面的网格，外表看起来是一层网格

然而细看，竟然是两层网格

我就是通过patch extrude出来的一层网格，但是有的地方（主要是加密变化造成的悬挂节点的地方）还是出现两层网格，很疑惑。。

[cresendo]

@Samuel-Tu 你好，最近自己也刚好尝试过一个二维算例，网格划分参照的是OpenFOAM/OpenFOAM-8/tutorials/incompressible/pimpleFoam/RAS/wingMotion这个自带算例。里面包括如何先用snappy划分，再利用extrudeMesh和createPatch。会不是顺序有哪里出错了呢

[Samuel-Tu]

@cresendo 顺序没有错，那个案例我看过，下午我再去检查一下那个wingMotion案例的网格是否是严格二维的。。因为按理说二维网格的厚度方向只能有一层网格，但是只要shm里用了refine加密，加密的截面上就有悬挂节点，加密截面的那里extrudeMesh生成网格在厚度方向就不止一层网格了

[李东岳]

为何不extrude之后再加密？

[veen]

cfmesh好像能定义不同方向不同的加密

[Samuel-Tu]

@李东岳 在shm里面加密的原因是第一步切割和第二步snap的时候能更好的拟合stl文件。。extrude之后加密对shm是没有影响，但是对拟合stl文件里的形状就没有作用了。。。我现在是现在gmsh里面画好渐变 网格，在shm里面就不refine了

[Samuel-Tu]

@veen cfmesh试过，算的很快，但是遇到复杂的模型发现网格效果没有shm好。。。

[veen]

@Samuel-Tu 好滴，gmsh和salome哪个好一些呢？

[xpqiu]

@Samuel-Tu

有可能是 paraview 显示的问题，你把 decompose polyhedra mesh 关掉再看看网格。
snappyHexMesh 画二维网格的流程可以参考 tutorials 里面 incompressible/pimpleFoam/RAS/wingMotion 这个算例，前面也有人提到了。

[Samuel-Tu]

@veen 我对gmsh熟一点，但是我感觉salome上限更高。。一般我只用salome画stl几何文件，不怎么用它画网格，没太搞明白如何精准设置网格类型和参数，如果有教程可以推给我看看。。

[Samuel-Tu]

@xpqiu 我算了这个算例，结果如下：
首先看全局的网格似乎确实只有一层：

提取一下加密界面的网格发现：

还是有这种两层的网格。

@xpqiu 在 关于二维模型的计算，这两种snappyHexMesh方法有什么不同，应该采用哪一个？ 中说：

有可能是 paraview 显示的问题，你把 decompose polyhedra mesh 关掉再看看网格

我没有找到这个按钮，但是后来我打开了VTK Polyhedra按钮，发现确实是显示的问题，双层网格不见了，加密界面都是单层的多面体了。

[xpqiu]

@Samuel-Tu 在 关于二维模型的计算，这两种snappyHexMesh方法有什么不同，应该采用哪一个？ 中说：

嗯，看来是显示的问题

Decomopose polyhedra 我指的是如下按钮，用 paraview 打开 .foam 后缀读取 openfoam 结果的时候有这个选项。
或者还有一个方法是用 foamToEnsight 转换成 ensight 格式来看

[veen]

@Samuel-Tu 最近几天看了看salome，感觉照着官方的tut操作一遍，使用上应该没啥问题了，除了操作有点繁琐，其他感觉用起来还挺顺手的，我先学学，学会了可能在b站录个系列教程。
salome的第一个tut就可以创建只有一层的网格，楼主有需要可以自己去看看。