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

问题描述

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

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

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

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

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

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

解决方法

这在 the documentation for Gmsh2D 中讨论。

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

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

然后您的 FiPy 脚本将使用

访问它们

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