Equivalence of mixed and nonconforming methods on general polytopal partitions. Part I: Multiscale and projection methods

这是一份关于 Simon Lemaire 所著论文《一般多面体划分下混合与非协调方法的等价性：第一部分：多尺度与投影方法》的详细技术总结。

1. 研究背景与问题定义

核心问题：
本文研究在变扩散系数（variable diffusion）问题中，定义在一般多面体（polytopal）网格上的混合有限元方法（Mixed Methods）与非协调原形方法（Nonconforming Primal Methods）之间的等价性。

数学模型：
考虑定义在有界连通 Lipschitz 区域 $\Omega \subset \mathbb{R}^d$ ($d \in \{2, 3\}$) 上的扩散问题：
$$\begin{aligned} \mathbf{b}^{-1}\boldsymbol{\sigma} - \nabla u &= 0 \quad \text{in } \Omega, \\ -\nabla \cdot \boldsymbol{\sigma} &= f \quad \text{in } \Omega, \\ u &= 0 \quad \text{on } \partial\Omega, \end{aligned}$$
其中 $\boldsymbol{\sigma}$ 为通量，$u$ 为势（标量），$\mathbf{b}$ 为可变的各向异性迁移率张量场，$f$ 为源项。

研究动机：

• 混合方法（Mixed）：直接求解通量 $\boldsymbol{\sigma}$ 和势 $u$，通常具有 $\mathbf{H}(\text{div})$ 连续性，但导致鞍点问题（Saddle-point system）。
• 原形方法（Primal）：仅求解势 $u$，通常是非协调的（Nonconforming），导致对称正定（SPD）系统，但通量重构可能不连续或不准确。
• 混合化（Hybridization）：通过引入面自由度（face unknowns），可以将混合方法转化为非协调原形格式，反之亦然。
• 挑战：在一般多面体网格（Polytopal meshes）和变系数 $\mathbf{b}$ 的情况下，建立这两种方法在离散层面的严格等价性，并确定何时可以构造出实用的“单层”（one-level）算法。

---

2. 方法论与理论框架

作者采用了一种基于虚拟空间（Virtual Spaces）的高层视角，通过“杂交化”（Hybridization）技术连接原形与混合变分公式。

2.1 连续层面的等价性

文章首先回顾了连续层面的四种变分形式：

1. 原形问题 (Problem 1)：标准的弱形式。
2. 混合问题 (Problem 2)：标准的混合弱形式。
3. 杂交原形问题 (Problem 3)：放松 $H^1_0$ 连续性，引入通量迹作为拉格朗日乘子。
4. 杂交混合问题 (Problem 4)：放松 $\mathbf{H}(\text{div})$ 连续性，引入势迹作为拉格朗日乘子。

通过引理 1 和引理 2，证明了：

• 杂交原形解 $\implies$ 混合解（通量由 $\mathbf{b}\nabla u$ 精确重构）。
• 杂交混合解 $\implies$ 原形解。
  这构成了图 1 所示的等价性闭环。

2.2 离散层面的虚拟构造 (Section 3)

为了在离散层面保持上述等价性，作者定义了一组有限维的虚拟空间：

• 通量空间 $\mathbf{\Sigma}_{\ell, k}(K)$：定义为满足特定散度和法向迹多项式约束的 $\mathbf{H}(\text{div})$ 空间。
• 势空间 $U_{\ell, k}(K)$：定义为满足特定扩散方程（$\nabla \cdot (\mathbf{b} \nabla v) \in P_{\ell-1}$）和法向通迹约束的空间。
• 关键性质：
  • 构造满足 $\mathbf{\Sigma}_{\ell, k}(K) = \mathbf{b} \nabla U_{\ell, k}(K)$，即精确保持本构关系。
  • 通量空间保持全局 $\mathbf{H}(\text{div})$ 连续性。
  • 势空间保持弱 $H^1_0$ 连续性（面矩意义下）。
• 等价性定理：证明了基于这些虚拟空间的离散混合问题（Problem 6）与离散原形问题（Problem 5）是等价的。

2.3 投影方法 (Projection Methods, Section 4)

由于虚拟空间 $U_{\ell, k}(K)$ 通常包含非多项式函数（除非 $\mathbf{b}$ 为常数且网格为特定形状），直接计算不可行。作者引入了投影方法的概念：

• 核心思想：当通量空间 $\mathbf{\Sigma}_{\ell, k}(K)$ 是多项式空间（如 Raviart-Thomas 或 BDM 空间），但本构关系 $\mathbf{b}\nabla u \in \mathbf{\Sigma}$ 无法精确满足时，通过投影算子 $\boldsymbol{\pi}^{\mathbf{b}^{-1}}_{\mathbf{\Sigma}}$ 将 $\mathbf{b}\nabla u$ 投影到通量空间。
• 投影公式：
  $$(\mathbf{b}^{-1} \boldsymbol{\pi}^{\mathbf{b}^{-1}}_{\mathbf{\Sigma}}(\mathbf{b} \nabla u), \boldsymbol{\tau}) = -(\nabla \cdot \boldsymbol{\tau}, \pi u) + \langle \boldsymbol{\tau} \cdot \mathbf{n}, \pi u|_{\partial K} \rangle$$
• 实用判据：提出了一个简单且可验证的满射性条件（Surjectivity of the flux map），即通量自由度映射到内部梯度和面迹的映射必须是满射且有界右逆。这是保证混合/原形问题适定性和等价性的充分条件。

---

3. 主要贡献与结果

3.1 多尺度方法的等价性 (Multiscale Methods)

• 结果：建立了四种不同的“无过采样”（oversampling-free）多尺度方法在离散层面的第一层等价性。
• 具体对应：
  • 证明了 MsHHO (Multiscale Hybrid-High Order) 方法与 MHM (Multiscale Hybrid-Mixed) 方法的等价性（推广了文献 [10] 的结果）。
  • 在混合设置下，证明了新的混合方法等价于基于虚拟势空间的原形方法。
• 意义：统一了不同多尺度方法的理论框架，表明它们本质上是同一数学结构的不同表现形式。

3.2 投影方法与自稳定混合方法 (Projection & Self-stabilized Methods)

• 结果：将投影有限元方法的概念扩展到虚拟势空间。
• 重新解释：证明了某些自稳定混合高阶方法（Self-stabilized Hybrid High-Order, HHO）和弱伽辽金方法（Weak Galerkin, WG）实际上是投影方法。
• 关键发现：
  • 对于一般的混合方法（如高阶 Raviart-Thomas 或 BDM 空间），如果本构关系不能被多项式空间精确保持，则其等价的原形形式必须包含一个投影算子。
  • 这为这些自稳定方法提供了混合有限元的等价形式，解释了其数值稳定性来源。

3.3 特殊情况的“银弹” (The Silver Bullet)

• 简单情形：当网格为单纯形（simplices）或超矩形（hyperrectangles），且迁移率 $\mathbf{b}$ 为单元常数时：
  • 最低阶 Raviart-Thomas (RT1) 混合方法等价于气泡增强（bubble-enriched）的 Crouzeix-Raviart (CR) 原形方法（在单纯形上）或 Rannacher-Turek 方法（在超矩形上）。
  • 新发现：证明了该等价性在任意形状的单纯形上均成立，无需假设网格正则性（此前文献通常假设正则单纯形）。
  • 静态凝聚：展示了如何通过静态凝聚将混合系统的鞍点问题转化为缩减规模的 SPD 系统，并给出了显式的重构公式。

3.4 一般多面体网格上的适定性判据

• 提出了一个基于通量映射满射性的判据，用于验证一般多面体网格上混合/原形方法的适定性和等价性。该判据是文献 [3] 中定理 4 的简化重述，但更易于在实际算法中检查。

---

4. 结论与意义

理论意义：

1. 统一框架：在一般多面体网格和变系数背景下，建立了混合方法与非协调原形方法之间严格的等价性理论框架。
2. 桥梁作用：连接了多尺度方法（MsHHO, MHM）与投影方法，揭示了自稳定混合方法（HHO, WG）背后的混合有限元本质。
3. 虚拟空间视角：通过引入虚拟空间，阐明了何时可以保持本构关系的精确性，以及何时必须引入投影算子。

实际意义：

1. 算法设计：为设计新的多面体网格数值方法提供了指导。如果通量空间满足满射条件，则可以直接构造等价的混合或原形格式。
2. 计算策略：对于满足特定条件（如 RT1 在单纯形上）的情况，提供了高效的静态凝聚策略，将混合问题转化为更易求解的 SPD 系统。
3. 鲁棒性：证明了在任意形状单纯形上的等价性，消除了对网格几何正则性的依赖，增强了方法在复杂几何网格上的适用性。

后续工作：
本文（Part I）主要关注多尺度方法和投影方法。第二部分（Part II）将专门讨论一般多面体单元方法（Polytopal Element Methods），特别是当通量空间不满足 $\mathbf{H}(\text{div})$ 全局连续性时的情况。

总结：
这篇文章通过引入虚拟空间和投影算子，深刻揭示了变扩散问题中混合与非协调方法之间的内在联系。它不仅推广了现有的多尺度方法理论，还为理解现代自稳定混合高阶方法提供了新的数学视角，是计算数学领域在多面体网格方法研究中的重要理论进展。