在 gmsh 中为物理域赋予属性

Question

基于我之前的问题physical entities for gmsh，我正在努力了解是否可以为*.geo 文件中定义的物理实体赋予域特征。

我正在开发 FEM 代码，在这种情况下，我试图定义网格文件中的所有内容，我的意思是：域、边界、域的属性等，以及涉及问题设置的所有内容。请参阅以下*.geo 文件：

Point(1) = {0, 0, 0, 1.0};
Point(2) = {1, 0, 0, 1.0};
Point(3) = {1, 1, 0, 1.0};
Point(4) = {0, 1, 0, 1.0};
Point(5) = {0, 0.5, 0, 1.0};
Point(6) = {1, 0.5, 0, 1.0};
Line(1) = {1, 2};
Line(2) = {2, 6};
Line(3) = {6, 3};
Line(4) = {3, 4};
Line(5) = {4, 5};
Line(6) = {5, 1};
Line(7) = {5, 6};
Line Loop(8) = {4, 5, 7, 3};
Plane Surface(9) = {8};
Line Loop(10) = {7, -2, -1, -6};
Plane Surface(11) = {10};
Physical Surface('top') = {9};
Physical Surface('bottom') = {11};

我打算在*.geo 文件中写入上述定义。这可以吗？

Answer 1

简而言之，you should not（请参阅 GMSH 常见问题解答部分中的问题 12）。根据 GMSH 哲学，并且

这项工作最好留给求解器，GMSH 提供了一种简单的机制来标记元素组

——你已经在使用了。

长版：

我猜，其目的是使用*.msh 文件（由 GMSH 处理相应的*.geo 文件获得）作为 FEM 求解器的输入。因此，我们可以使用GMSH *.msh 的规范。 （由于GMSH最近刚更新到v4，我们也可以看看之前的legacy MSH v2）。

在 v4（v2 提供的更少）中，可以使用以下部分（一些，可选）：

1. MeshFormat - 标题信息
2. PhysicalNames - 物理组编号和名称之间的联系
3. Entities、Partitioned Entities、Nodes、Elements、GhostElements – 几何定义
4. Periodic – 周期信息
5. NodeData、ElementData、ElementNodeData - 后处理-意图数据集。

因此，没有一个部分可以自然地满足您的需求。 你有两个糟糕的选择（我能想到的）：

• 由于标题未知的部分会被忽略，因此您可以手动（在生成 *.msh 之后）在此处添加数据，并根据您的需要为求解器提供所需的输入功能
• 尝试使用 postprocessing-intent ElementData，但这会使物理组标签的使用变得无用（因为您必须根据单个元素 ID 进行操作）

因此，基于此，我将停止正确分配物理组标签和编号，并在 FEM 求解器一侧进行材料属性分配。