我已经生成了一个 gmsh 网格,并通过 Gmsh2D 导入到 FiPy 中。我想解决表面问题以设置边界条件,但我不知道如何操作。
在我在文档中发现的 FiPy 示例中,建议命名某些行以便稍后处理它们。将网格导入 fipy 后如何执行此操作?
// note: if you do not use any labels, all Cells will be included.
Physical Surface("Outer") = {1};
Physical Surface("Middle") = {2};
Physical Surface("Inner") = {3};
// label the "north-west" part of the exterior boundary
// note: you only need to label the Face elements
// (Physical Line in 2D and Physical Surface in 3D) that correspond
// to boundaries you are interested in. FiPy does not need them to
// construct the Mesh.
Physical Line("NW") = {5};
-----------------
编辑:
对于简单的曲面,这将起作用:
我忽略了mesh.exteriorFaces
。
对于一个简单的圆,这导致了简单的解决方案:
xfc, yfc = mesh.faceCenters() # returns x faceCenters coordinate, ...
xcc, ycc = mesh.cellCenters()
plot(xfc[where(mesh.exteriorFaces == False)],yfc[where(mesh.exteriorFaces == False)],'ro', label='inside')
plot(xfc[where(mesh.exteriorFaces == True)],yfc[where(mesh.exteriorFaces == True)],'ko', label='surface')
legend()
尽管如此,我仍在寻找有关如何从外部访问 gmsh 代码的答案,但这也许对其他人有帮助:)
最佳答案
“外”不是一张物理面孔,而是一组物理细胞。要访问“Outer”的边界面,您可以添加
Physical Line("Outer Boundary") = {1, 2, 3, 4, 5, 6, 7, 8};
到您的 Gmsh 脚本,然后使用
应用约束var.constrain(value, where=squaredCircle.physicalFaces["Outer Boundary"])
正如您所发现的,您始终可以访问mesh.exteriorFaces
来获取定义整个网格边界的面(即,只有一个单元格的面)一边)。使用Physical Surface
定义的Gmsh域不一定受mesh.exteriorFaces
限制。
关于python - 将网格导入 Fipy 后如何从 python 访问 gmsh 代码?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/59612531/