[论文解读] 3D Volume Calculation For the Marching Cubes Algorithm in Cartesian Coordinates
本文提出了一种基于散度定理的新方法，用于在三维笛卡尔网格中计算部分体积，替代传统马奇立方算法中的四面体剖分。通过选择一个散度恒定的矢量场（例如 f = (x,y,z)），体积通过等值面的曲面积分计算，相比四面体分解，该方法更加鲁棒、高效且便于编码。
From a scalar field defined at the corner of a cube, an isosurface can be extracted using the Marching Cube algorithm. The isosurface separates the cell into two or more partial cells. A similar situation arises when an material interface in the Front Tracking method cuts through the computational cells. A popular method to calculate the volumes of the partial cells is to first partition the cells into tetrahedra and then sum together the volumes of the tetrahedra for the corresponding partial cells. In this paper, the divergence theorem is used to calculate the volumes of the partial cells generated by the Marching Cubes algorithm. This method is both more robust and efficient compared with the tetrahedralization approach.
研究动机与目标
- 解决马奇立方算法中计算部分单元体积缺乏高效且鲁棒方法的问题。
- 克服在复杂界面构型下基于四面体剖分的体积计算所面临的复杂性和编码困难。
- 提供一种基于散度定理的数值稳定且计算高效的体积积分替代方法。
- 使嵌入边界法和前缘追踪法在多相流和界面问题中实现精确的体积计算成为可能。
- 通过使用统一的、带标签的模板系统，简化所有23种唯一立方体构型的实现。
提出的方法
- 应用散度定理将体积计算为曲面积分：体积 = (1/3) ∫∂Ω n·f dS，其中 f = (x,y,z) 确保 ∇·f = 3。
- 该方法使用矢量场 f = (x,y,z)，使得 ∇·f = 3，从而可通过在封闭界面进行曲面积分实现精确的体积计算。
- 界面被三角剖分为平面面片，曲面积分通过二阶求积公式进行计算，可在三角形网格上获得精确结果。
- 采用双字母编码方案（例如 v0, e0, W, I）系统性地标识顶点、边和面，便于代码实现。
- 利用对称性将256种可能的立方体构型简化为23种唯一情况，每种情况均预计算了三角剖分和曲面积分贡献。
- 对每种情况，通过曲面积分计算一个分量的体积，另一分量则通过从立方体总体积中减去得到。
实验结果
研究问题
- RQ1散度定理能否有效应用于马奇立方算法中三维笛卡尔网格的部分体积计算？
- RQ2与基于四面体剖分的体积计算相比，基于散度定理的方法在鲁棒性和效率方面表现如何？
- RQ3能否为所有23种唯一立方体构型开发出一致且便于编码的标签与三角剖分系统？
- RQ4该方法是否适用于需要精确界面体积计算的嵌入边界法和前缘追踪法？
- RQ5在该情境下，选择何种矢量场可确保通过曲面积分实现精确的体积计算？
主要发现
- 基于散度定理的方法在三角形界面上使用二阶求积公式，可精确计算体积，避免了四面体剖分带来的数值误差。
- 该方法比四面体剖分更具鲁棒性和效率，尤其在处理具有多个边和面交叉的复杂界面拓扑结构时表现更优。
- 该方法通过使用统一的标签系统（如 v0, e0, I, W）标识顶点、边和面/界面三角形，显著简化了编码实现。
- 23种唯一立方体构型已完整枚举，所有情况均预计算了三角剖分和曲面积分贡献。
- 该方法可应用于椭圆/抛物型界面问题、两相不可压缩流体流动以及磁流体动力学流动等场景，实现精确的体积计算。
- 通过使用立方体的总体积计算互补分量的体积，减少了计算冗余。
更好的研究，从现在开始
从阅读论文到最终审阅，大幅缩短您的研究时间。
无需绑定信用卡
本解读由 AI 生成，并经人工编辑审核。