技术摘要：模仿谱元平流法

问题陈述

本文研究了有界区域内微分形式线性平流方程的数值离散化问题。经典的标量平流方程 $\partial_t \rho + \nabla \cdot (v\rho) = 0$ 在微分几何框架下被推广为：
$$\frac{\partial \alpha^{(k)}}{\partial t} + \mathcal{L}_v \alpha^{(k)} = 0$$
其中 $\alpha^{(k)}$ 是一个 $k$-微分形式（根据 $k$ 的不同，代表标量、向量或更高维度的量），$v$ 是给定的 Lipschitz 连续速度场，$\mathcal{L}_v$ 是李导数（Lie derivative）。其挑战在于开发一种既能保持物理定律底层几何结构，又能确保局部质量守恒，同时能区分在特定时刻评估的量与在时间间隔内积分的量的离散化方法，并实现谱精度。

方法论

所提出的方法——模仿谱元平流法（Mimetic Spectral Element Advection），通过引入 Cartan 同伦公式将李导数纳入先前的模仿框架 [4] 中。该方法建立在三个核心支柱之上：

1. 微分几何基础
该方法使用微分形式语言。李导数 $\mathcal{L}_v$ 通过 Cartan 同伦公式表示为：
$$\mathcal{L}_v \alpha^{(k)} = d\iota_v \alpha^{(k)} + \iota_v d\alpha^{(k)}$$
其中 $d$ 是外微分，$\iota_v$ 是内乘（interior product）。公式的一个关键方面是内乘与外积之间的伴随关系，该关系通过 $L^2$ 内积定义：
$$(\iota_v \alpha^{(k)}, \beta^{(k-1)})_{L^2} = (\alpha^{(k)}, v^\flat \wedge \beta^{(k-1)})_{L^2}$$
这种对偶性允许将由内乘表示的物理量映射到其对偶微分 1-形式。

2. 空间离散化（模仿谱元）
空间域通过由点、线段、曲面和体积组成的胞腔复形（cell complex）进行离散化。光滑 $k$-形式的空间由由基底形式 $\epsilon_i^{(k)}$ 张成的有限维空间 $\Lambda_h^k$ 来近似。

• 基函数构建： 基函数使用在四边形单元上的分片多项式展开（张量积）进行构建。在 1D 中，0-形式由在 Gauss-Lobatto 节点上评估的 Lagrange 多项式表示，而 1-形式则由旨在满足特定积分性质的“边多项式”表示。
• 交换投影： 定义了一个投影算子 $\pi_h$，使得该算子与外微分交换（$\pi_h d = d \pi_h$）。这确保了离散外微分 $d_h$ 由仅包含 $\{0, 1, -1\}$ 值的关联矩阵表示，从而保留了连续算子的拓扑结构。

3. 时间离散化（规范模仿积分器）
时间积分采用源自规范 Gauss 配点积分器 [5] 的任意阶辛算子（symplectic operator）。

• 交错时间节点： 该方法区分了两组时间节点：
  • Gauss-Lobatto 节点 ($t_k$)： 解变量（如 $\rho$）在此处定义。
  • Gauss 节点 ($\tilde{t}_q$)： 时间导数（通量）在此处评估。
• 离散演化： 这种交错结构产生了一个离散积分器，其中在一个区间内的解的变化量等于在内部 Gauss 节点处评估的通量。这种结构模仿了 leap-frog 方法和隐式中点规则，确保了辛特性。

4. 离散内乘
通过在离散设置中强制执行对偶配对（公式 16）来对内乘 $\iota_v$ 进行离散化。这导致了一个通过求解涉及基底形式与速度场内积的线性系统来计算离散通量的系统。

主要贡献

• 几何一致性： 该方案通过 Cartan 公式显式地结合了李导数，确保了离散化尊重微分算子及其 Hilbert 伴随的无度量表示。
• 谱精度： 该方法利用高阶多项式基函数，实现了空间上的谱收敛。
• 局部质量守恒： 使用关联矩阵和离散外微分的特定结构保证了局部质量守恒。
• 时间步长区分： 该框架严格区分了在时间瞬时评估的量与在时间间隔内积分的量，利用交错时间网格来实现高阶辛积分。

数值结果

作者在具有周期性边界条件的 2D 域上进行了数值实验，测试了在常数速度场和 Rudman 涡旋速度场下正弦波和正弦钟形波（sine bells）的平流情况。

• 收敛性： 该方法展示了 $(p+1)$ 阶的代数 $h$-收敛和谱 $p$-收敛，前提是时间积分误差不会占据主导地位。
• 时间积分精度： 解的误差表明其取决于时间积分方案的阶数 ($p_t$)。如果时间方案足够精确，则初始离散误差将被保留；否则，人工色散会随时间增加。
• 质量守恒： 总质量误差在最初的 $10^3$ 个时间步内保持在机器零附近，并在 $2 \times 10^4$ 步后仍低于 $10^{-12}$，展示了极佳的守恒特性。
• 可逆性： 一个涉及在 Rudman 涡旋中平流正弦波，随后反转流向的测试表明，该方法可以恢复初始解，证明了积分方法的可逆性。
• 人工色散： 文中指出，人工色散是时间积分阶数的函数；较低阶的时间方案会引入取决于被平流波频率的色散误差。

重要性与主张

本文声称，所推导的方案成功地在物理兼容（模仿）框架内结合了谱精度与局部质量守恒。通过利用 Cartan 同伦公式将李导数扩展到模仿框架中，作者提供了一种能够阐明平流方程几何结构的离散化方法。该方法被呈现为一种稳健的方法，用于解决那些需要保留物理定律底层几何和拓扑属性的平流问题。作者也谦虚地指出，虽然空间离散化非常精确，但整体精度和色散特性取决于所使用的时间积分方案的阶数。