Fipy 中的特殊网格划分（空心圆柱体的 3d 网格划分）答案

【问题标题】：Special meshing in Fipy(3d meshing of a hollow cylinder)Fipy 中的特殊网格划分（空心圆柱体的 3d 网格划分）
【发布时间】：2020-11-09 21:31:20
【问题描述】：

我想在图片中对圆柱形裂缝进行网格划分，如下所示： enter image description here

网格中的单元格具有 6 个面和 8 个节点的六面体形式。该网格将用于 fipy 求解扩散方程。（径向流）任何人都可以帮助回答我关于使用 gmsh 代码的网格的问题吗？因为它是与 fipy 兼容的唯一方法。我感谢您的帮助。非常感谢。

【问题讨论】：

【解决方案1】：

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

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

如果您的问题不是圆柱对称的，那么正如@RanjeethMahankali 评论的那样，请向我们展示您尝试过的操作。

【讨论】：

@jeguyer 先生，我这里的问题是3D 圆柱形 对称，所以我需要顶点的三个分量坐标（x，y，z）；使用圆柱形网格2D，我只有 x,y
一开始我在 Grid3D 中模拟流动（将 fipy 与名为 PFC3D 的离散元素软件耦合，以在求解流动的扩散方程时考虑粒子：self.mesh = fipy.Grid3D(nx=10, ny=20, nz=10, dx=0.01, dy=0.01, dz=0.01) 现在我想要移动到 3D 径向流...
是否可以在 z 方向上平移圆柱网格 2D 以产生 3D 问题？真的很感谢你的帮助！
我不明白你的问题。 CylindricalGrid2D 定义在 r 和 z 中，关于 r=0 具有圆柱对称性。
Guyer 先生，我认为这种类型的网格无法通过 fipy 网格生成器执行。问题是我需要三个节点（x、y、z）组件来构建我的网格我的代码（Itasca PFC3D）。CFD 模拟只能在 3D 问题上执行。CylindricalGrid2D 只会给我 2 个组件。我在这里发布一个新问题 [stackoverflow.com/questions/64818780/problem-in-loading-msh-file- in-fipy-via-gmsh3d] 因为我从 Gmsh 导入网格时遇到了一些问题。感谢您的帮助。