通过 FiPy 解决圆柱体外网格表面的最佳方法，

Question

我可以通过任何 FiPy 或 ... 模块处理外部网格面中心吗？对于半径为 'R' 的圆柱体，会在其上创建类似 O'grid 扫描的相关网格。看起来，'R' 大于最外面的网格面中心；所以它们之间是有区别的（图），不能轻易用“R”来表示。

我需要提及围绕圆柱壁（不是顶壁和底壁）的所有外网格面中心。网格系统由“FiPy Gmsh”导入（FiPy version = 3.1 and Python 2.7）。 FiPy 'solver.mesh.getFaceCenters()' 得到一个简单的网格面中心数组。我认为外六边形网格层的圆柱壁相关外表面可以通过以下代码解决：

reduce(numpy.logical_and, (FORMULAs))               # where FORMULAs are specifying mesh face center coordinates limits

我尝试使用以下代码，但如果 FORMULA 适用，则需要对段 '> R' 进行一些更改。

x,y,z = solver.mesh.getFaceCenters()
np.sqrt(x ** 2 + y  ** 2) > R       #  as FORMULAs

如果有人能帮助我解决这个问题，我将不胜感激。我认为最好的方法一定是使用 FiPy 和 numpy 模块。

Answer 1

这在 the documentation for Gmsh2D 中讨论。

虽然可以参数化地定义这样的面，但由于离散化网格的有限分辨率，有必要允许一些位置误差。在抽象几何定义中定义你想要的面并直接引用它们要好得多，例如，你在 Gmsh 中的 GEO 文件会有类似

的内容

Physical Volume("FORMULA") = {5, 6, 7};

然后您的 FiPy 脚本将使用

访问它们

formula_faces = solver.mesh.physicalFaces["FORUMULA"]

通过 FiPy 解决圆柱体外网格表面的最佳方法，

Best way to address outer mesh surfaces of a cylinder by FiPy ,

python-2.7

numpy

mesh

fipy

numpy-ufunc