MFEM中处理带圆形孔洞的方形网格生成与细化技术

概述

在MFEM有限元计算框架中，处理复杂几何形状的网格生成是一个常见需求。本文将详细介绍如何在MFEM中创建包含圆形孔洞的方形计算域网格，并实现高质量的网格细化。

几何建模基础

在MFEM中，通常需要借助外部网格生成工具如Gmsh来创建复杂几何形状的网格。对于包含圆形孔洞的方形域，我们需要：

1. 定义方形主域
2. 定义圆形子域
3. 进行布尔运算得到最终几何

Gmsh的几何描述语言(.geo文件)是实现这一过程的关键。正确的几何定义应该如下：

Rectangle(1) = {x0, y0, 0, width, height, 0};
Circle(2) = {xc, yc, 0, radius, 0, 2*Pi};

其中参数含义为：

• (x0,y0): 矩形左下角坐标
• width,height: 矩形宽度和高度
• (xc,yc): 圆心坐标
• radius: 圆半径

布尔运算与物理分组

通过布尔差运算可以得到带孔洞的方形域：

BooleanDifference{ Surface{1}; Delete; }{ Surface{2}; }

为了在MFEM中区分不同区域和边界，我们需要为几何元素分配物理属性：

Physical Surface("MainDomain", 1) = {1};
Physical Surface("Hole", 2) = {2};
Physical Curve("OuterBoundary", 1) = { /* 矩形边界曲线 */ };
Physical Curve("InnerBoundary", 2) = { /* 圆形边界曲线 */ };

高阶网格生成

MFEM支持高阶有限元计算，因此我们需要在Gmsh中生成高阶网格：

Mesh.ElementOrder = 3;  // 生成三次单元
RecombineMesh;         // 重组三角形为四边形

网格细化策略

MFEM中的网格细化会自动考虑几何形状，但需要注意：

1. 初始网格质量直接影响细化结果
2. 圆形边界需要足够多的初始节点以保证几何精度
3. 可通过Gmsh控制局部网格尺寸：

Mesh.MeshSizeMin = 0.01;  // 最小网格尺寸
Mesh.MeshSizeMax = 0.1;   // 最大网格尺寸

MFEM中的边界条件处理

在MFEM代码中，可以通过属性标记来区分不同边界：

// 获取圆形边界标记
auto circle_bdr = mesh.bdr_attribute_sets.GetAttributeSetMarker("InnerBoundary");

// 在圆形边界上施加边界条件
x.ProjectBdrCoefficient(bc_func, circle_bdr);

常见问题解决

1. 几何位置异常：确保矩形和圆的坐标参数正确理解，矩形后三个参数是尺寸而非坐标

2. 网格质量差：适当调整MeshSizeMin/Max参数，增加RecombineSurface操作

3. 边界识别错误：检查物理组的定义是否正确

4. 高阶单元变形：确保初始网格足够精细以支持高阶插值

结论

在MFEM中处理带圆形孔洞的方形网格需要结合Gmsh的几何建模能力和MFEM的网格处理功能。通过合理设置几何参数、物理属性和网格参数，可以生成高质量的计算网格，为后续有限元分析奠定基础。高阶网格和适当的细化策略能够有效提高计算精度，特别是在曲边附近的解精度。