Fipy中的特殊网格划分(空心圆柱体的三维网格划分)

Question

我想在下面的图片中对柱状骨折进行网格处理：enter image description here

网格中的单元具有6个面的六面体形式，8个nodes.This网格将在fipy中用于求解扩散方程。(径向流动)有人能帮助回答我关于使用gmsh代码的网格的问题吗？因为这是与fipy兼容的唯一方法。我很感谢你的帮助。非常感谢。

Answer 1

如果需要圆柱对称，可以使用CylindricalGrid1D或CylindricalGrid2D生成网格。偏移网格以在中间创建孔，例如：

mesh = CylindricalGrid2D(dr, dz, Lr, Lz) + [[r0], [0]]

如果你的问题不是圆柱对称的，那么就像@RanjeethMahankali评论的那样，请告诉我们你试图做什么。