MOOSE框架中网格转换工具保留边界面错误分析与修复

问题背景

在MOOSE多物理场仿真框架中，网格处理是数值模拟的重要前置步骤。开发者在处理复杂网格转换时发现了一个关于边界面保留的缺陷：当使用ElementsToTetrahedronsConverter工具将六面体网格转换为四面体网格时，原始网格中定义的边界面信息会出现错误——部分边界面被错误保留或丢失。

问题现象

具体表现在使用以下网格处理流程时：

1. 首先生成一个半圆形二维网格
2. 删除右侧边界
3. 将二维网格旋转生成三维网格
4. 将六面体元素转换为四面体元素

转换后，"outer"边界面包含了额外的错误面片，这会导致后续边界条件施加出现偏差，严重影响仿真结果的准确性。

技术分析

网格转换原理

在MOOSE框架中，ElementsToTetrahedronsConverter工具负责将高阶元素（如六面体）分解为低阶元素（如四面体）。这个过程需要正确处理以下信息：

• 节点坐标
• 元素连接关系
• 边界面信息

问题根源

经过代码分析，发现问题的核心在于MooseMeshElementConversionUtils工具类中边界面保留逻辑存在缺陷。当原始六面体元素被分割为多个四面体时，原有的边界面信息没有被正确映射到新生成的四面体面上。

具体表现为：

1. 边界面识别算法没有考虑元素分割后新生成的面
2. 原有边界面与分割后面之间的对应关系建立不完整
3. 边界条件传递过程中缺乏严格的几何验证

解决方案

开发团队通过以下步骤修复了该问题：

1. 完善面映射算法：重新设计六面体到四面体转换过程中边界面映射的逻辑，确保每个原始边界面都能正确对应到新生成的四面体面上。

2. 添加几何验证：在转换过程中增加几何一致性检查，确保新生成的边界面与原始边界面的几何特征（如法向量、面积等）保持一致。

3. 优化数据结构：改进边界面信息的存储方式，使其能够更好地处理元素分割后的复杂情况。

影响评估

该修复确保了以下关键功能：

• 边界条件能够正确施加在转换后的网格上
• 复杂几何形状的网格转换结果更加可靠
• 多物理场耦合仿真中的边界处理更加精确

最佳实践建议

对于MOOSE用户，在处理网格转换时应注意：

1. 在转换复杂网格后，建议使用MeshChecker工具验证边界面完整性
2. 对于关键边界，可在转换前后分别输出网格进行可视化对比
3. 考虑使用更高阶的网格转换算法以获得更好的边界保持效果

该修复已合并到MOOSE主分支，用户可以通过更新代码库获取这一改进。