如何解决通过 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"]
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。