The longest-edge bisection algorithm may produce degenerating tetrahedra

技术摘要：“最长边二分算法可能会产生退化四面体”

问题陈述
本文探讨了四面体网格细化理论中的一个关键空白。虽然最长边二分算法是构建嵌套单纯复形（simplicial meshes）的一种标准技术，但其在三维空间中的非退化行为尚未得到充分理解。在二维空间中，通过重复最长边二分生成的三角形的收敛性和形状正则性是已知的。然而，在三维空间中，不受限制的几何指令——即必须二分当前最长的边，但在多个边的长度相等时可以任选其一——缺乏非退化的保证。以往的研究主要集中在标记边方案（确保有限相似类）或特定的四面体族，而对于具有任意平局处理（tie-breaking）规则的一般性、不受限制规则的行为，仍是一个悬而未决的问题。

方法论
作者构建了一个显式的、精确的反例，以证明不受限制的规则会导致退化。其方法论包括：

1. 定义参数族： 定义了一个特定的四面体 $E(a)$，其顶点取决于参数 $0 < a \le 1$。
2. 两步递归过程： 本文分析了一个两步二分过程：
   • 第一步： 二分 $E(a)$ 的唯一最长边 $AD$，并保留子四面体 $O(a)$。
   • 第二步： 在 $O(a)$ 中，边 $BD$ 和 $DC$ 的长度完全相等且均为最长。作者选择 $BD$ 进行二分，保留生成的子四面体 $F(a)$，并对其顶点进行重新标记。
3. 精确全等证明： 文中证明了所得出的四面体 $F(a)$ 与原始参数族成员 $E(a/2)$ 全等。这建立了一个递归关系，即每经过两个二分步骤，形状参数 $a$ 就会减半。
4. 迭代构造： 从 $E(1)$ 开始，该过程生成了一个无限序列：$E(1) \to O(1) \to E(1/2) \to O(1/2) \to \dots$。
5. 正则性分析： 作者计算了当 $a \to 0$（其中 $a_k = 2^{-k}$）时，该序列的归一化体积比 ($|T|/h_T^3$) 和特定的二面角，以测试是否符合标准的正则性准则。

核心贡献与结果
本文提供了严谨的证明，证明了三维空间中不受限制的最长边二分算法可能会产生一系列退化的四面体。具体发现如下：

• 违反形状正则性： 归一化体积比 $|E_k|/h_{E_k}^3$ 随着 $k \to \infty$ 趋于零。具体而言，该比例渐近衰减为 $2^{-k/2}$，证明该族不存在统一的正下界。
• 违反最小角条件： $E(a)$ 在边 $CD$ 处的内二面角 $\delta(a)$ 趋于零。文中显示 $\delta(a) \sim \sqrt{8a/7}$，这意味着该序列包含任意“扁平”的角。
• 违反最大角条件： $O(a)$ 在边 $MC$ 处的内二面角 $\Theta(a)$ 趋于 $\pi$。其余弦值趋于 $-1$，表明四面体在另一种构型下变得任意“片状”（sliver-like）或扁平。
• 平局处理的作用： 退化是由最长边长度的循环平局（$BD = DC$）驱动的。本文证明，一个简单的确定性平局规则（选择与较长边相对的边）足以选中这一“坏”分支。

意义与主张
本文的主要主张是适度且明确的：它证明了在不受限制的最长边规则下，存在一个退化的可容许轨道（admissible orbit）。标题中使用“可能”（may）一词至关重要；该构造并不断言所有平局处理惯例都会导致退化，也不声称所有轨道都是退化的。相反，它确立了如果不将显式的平局解决机制纳入算法并作为正则性定理的一部分进行分析，该算法就无法保证生成非退化网格。

这项工作是对“仅靠几何选择规则本身足以保证三维正则性”这一假设的反例。它强调了在三维空间中，细化的行为关键取决于标记和平局处理规则，并且在相等的最长边之间进行任意选择会导致同时失去形状正则性、最小角界限和最大角界限。