Asymptotic numerical hypocoercivity of the space-time discontinuous Galerkin method for the Kolmogorov equation

技术摘要：时空不连续伽辽金法在柯尔莫哥洛夫方程中的渐近数值超强制性

问题陈述
本文研究了经典柯尔莫哥洛夫方程（Kolmogorov equation）的数值离散化，该方程是表现出退化扩散特性的动力学类型方程的代表模型。该方程定义在有界空间区域 $\Omega \subset \mathbb{R}^2$ 上，形式为 $u_t + Lu = 0$，其中 $L = -\partial_{xx} + x\partial_y$。该算子的一个核心特征是，扩散仅发生在 $x$ 方向，而输运项 $x\partial_y$ 则在 $y$ 方向起作用。尽管缺乏 $y$ 方向的显式扩散，该系统仍表现出“超强制性”（hypocoercivity）——即输运与退化耗散之间的相互作用将衰减效应传播到整个空间域，从而实现向平衡态的指数收敛。

核心挑战在于证明标准的数值格式，特别是经典的空时不连续伽辽金（dG）方法，能够保留这种离散层面的超强制结构。虽然超强制性在连续设定下已得到深入理解（特别是通过 Villani、Hérau、Nier 等人的工作），但在伽辽金方法中保持该性质的研究仍处于探索阶段，以往的数值结果主要局限于需要专门非标准形式的有限差分、有限体积或非协调有限元方法。

方法论
作者分析了应用于柯尔哥洛夫方程的标准空时 dG 方法。空间离散采用针对输运项的上风 dG 格式以及针对扩散项的对称内罚 dG 格式。时间离散采用标准的 dG 方法，多项式次数为 $q$。

其方法论的核心在于构建一种改进的稳定性分析，以模拟连续超强制框架。在连续设定中，证明超强制性通常通过修改内积（引入混合导数或特定的梯度项）来实现，从而有效地创建一个具有更强控制力的“强”范数，使其具备谱间隙。

为了在不改变 dG 方法的标准双线性型的条件下实现类似结果，作者采用了一种特定的测试函数策略。他们不是仅通过测试离散方程相对于解 $U$ 的情况，而是测试其相对于复合函数 $V = V_1 + V_2 + V_3$ 的情况：

1. $V_1 = U$：标准的能量测试项。
2. $V_2 = \tau(xU_y + U_t)$：模仿输运算子的项，由依赖于网格尺寸 $h$ 和多项式次数 $p$ 的参数 $\tau$ 进行缩放。
3. $V_3 = -\text{div}(A\nabla U)$：涉及修改后的矩阵 $A$（依赖于 $h$ 和 $p$）以及梯度散度项，旨在恢复对缺失的 $y$ 导数的控制。

至关重要的是，作者观察到对于 dG 使用的不连续多项式空间，函数 $V$ 仍然属于与 $U$ 相同的离散空间。这使得作者可以使用 $V$ 作为可容许的测试函数，而无需引入以往使用连续单元时所需的非标准 Petrov-Galerkin 项。

主要贡献与结果

1. 标准 dG 的离散超强制性：本文证明了经典的空时 dG 方法具有渐近数值超强制性。通过在“增强型超强制范数”下证明单侧 inf-sup 条件，作者展示了数值解在长时间内向平衡态呈指数衰减。这是首个针对标准伽辽金方法的此类结果。

2. 增强范数与谱间隙：分析引入了一系列包含数值解全梯度的范数族。作者在这些范数下证明了网格依赖的谱间隙（Poincaré 型不等式）。衰减率 $\kappa$ 取决于离散参数，其标度关系为 $\kappa \sim h_{\min}^4 p^{-8}$。虽然随着网格细化（$h \to 0$）或多项式次数增加（$p \to \infty$）衰减率会发生退化，但该结构确保了在任何固定离散下的稳定性与衰减。

3. Inf-Sup 稳定性：除了超强制性外，本文还证明了空时 dG 方法在增强范数下具有 inf-sup 稳定性。这一结果比现有的针对平流-扩散问题的空时 dG 方法的稳定性结果更强，因为它提供了对全梯度而非仅对平流方向的控制。

4. 先验误差估计：利用已建立的 inf-sup 稳定性，作者在增强的空时范数下推导出了对于足够正则解的最优先验误差界。在满足 $h \approx O(k)$ 的条件下，误差估计相对于空间网格尺寸 $h$ 和时间步长 $k$ 均是优化的。

5. 数值验证：数值实验证实了理论上的收敛速率，并展示了解向平衡态的指数衰减。实验表明，尽管理论衰减率取决于 $h$ 和 $p$，但观察到的向平衡态的收敛过程是稳健的。

意义与主张
作者声称，这项工作代表了对已知经典伽辽金方法实现（渐近）数值超强制性的首次证明。其意义在于，该方法无需像以往的非协调有限元方法那样通过添加额外的 Petrov-Galerkin 项来修改双线性型，即可保留超强制结构。

本文强调，其证明技术具有通用性。通过构造模拟连续超强制理论中修改内积效应的测试函数，该方法有望推广到其他动力学类型方程及退化抛物型问题。这些结果为使用标准高阶空时 dG 方法进行退化扩散问题的长时间模拟提供了严谨的基础，确保了数值格式不会人为地抑制连续模型中固有的物理耗散机制。