技術要約：ミメティック・スペクトル要素法による移流

問題設定

本論文は、有界領域における微分形式の線形移流方程式の数値離散化を取り扱う。古典的なスカラー移流問題 $\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$ は与えられたリプシッツ連続な速度場、$\mathcal{L}_v$ はリー微分である。課題は、物理法則の根底にある幾何学的構造を保存し、局所的な質量保存を保証し、かつ、特定の時刻に評価される量と時間間隔にわたって積分される量を区別しながら、スペクトル精度を実現する離散化手法を開発することにある。

手法

提案手法である**ミメティック・スペクトル要素法（Mimetic Spectral Element Advection）**は、カルタンのホモトピー公式を用いてリー微分を組み込むことにより、先行研究 [4] で確立されたミメティック・フレームワークを拡張したものである。このアプローチは、以下の3つの核心的な柱に基づいている：

1. 微分幾何学的基礎
本手法は微分形式の言語を利用する。リー微分 $\mathcal{L}_v$ は、カルタンのホモトピー公式を用いて次のように表される：
$$\mathcal{L}_v \alpha^{(k)} = d\iota_v \alpha^{(k)} + \iota_v d\alpha^{(k)}$$
ここで、$d$ は外微分、$\iota_v$ は内部積である。定式化における極めて重要な側面は、内積 $\iota_v$ とウェッジ積（外積）の間の随伴関係であり、これは $L^2$ 内積を通じて以下のように定義される：
$$(\iota_v \alpha^{(k)}, \beta^{(k-1)})_{L^2} = (\alpha^{(k)}, v^\flat \wedge \beta^{(k-1)})_{L^2}$$
この双対性により、内部積によって表される物理量は、その双対な微分1形式へと写像することが可能となる。

2. 空間離散化（ミメティック・スペクトル要素）
空間領域は、点、線分、面、および体積からなるセル複体を用いて離散化される。滑らかな $k$-形式の空間は、基底形式 $\epsilon_i^{(k)}$ によって張られる有限次元空間 $\Lambda_h^k$ によって近似される。

• 基底の構成: 基底関数は、四角形要素上の区分的多項式展開（テンソル積）を用いて構築される。1次元において、0形式はガウス・ロバット・ノードで評価されるラグランジュ多項式によって表現され、1形式は特定の積分特性を満たすように導出された「エッジ多項式」によって表現される。
• 可換射影: 外微分と可換である射影演算子 $\pi_h$ （すなわち $\pi_h d = d \pi_h$）が定義される。これにより、離散外微分 $d_h$ が、$\{0, 1, -1\}$ の値のみを含むインシデンス行列として表現されることが保証され、連続作用素のトポロジカルな構造が保存される。

3. 時間離散化（カノニカル・ミメティック・インテグレータ）
時間積分には、カノニカル・ガウス・コロケーション・インテグレータ [5] から導出された任意の次数のシンプレクティック作用素を用いる。

• 交互配置された時間ノード: 本手法は、2組の時間ノードを区別する：
  • ガウス・ロバット・ノード ($t_k$): 解の変数（例：$\rho$）が定義される点。
  • ガウス・ノード ($\tilde{t}_q$): 時間微分（フラックス）が評価される点。
• 離散進化: この交互配置により、ある区間における解の変化が、内部のガウス・ノードで評価されたフラックスに等しいとされる離散インテグレータが得られる。この構造は、リープフロッグ法やインプリシット・ミッドポイント則を反映しており、シンプレクティックな性質を保証する。

4. 離散内部積
内部積 $\iota_v$ は、離散設定において双対ペアリング（式16）を強制することによって離散化される。これにより、基底形式と速度場の内積を含む線形システムを解くことで、離散フラックスが計算されるシステムが構築される。

主な貢献

• 幾何学的整合性: 本スキームは、カルタンの公式を用いてリー微分を明示的に組み込んでおり、微分作用素およびそのヒルベルト随伴の計量に依存しない表現を尊重した離散化を実現している。
• スペクトル精度: 高次の多項式基底関数を利用することで、空間におけるスペクトル収束を可能にしている。
• 局所的質量保存: インシデンス行列の使用と離散外微分の特定の構造により、局所的な質量保存が保証される。
• 時間ステップの区別: 本フレームワークは、時刻に評価される量と時間間隔にわたって積分される量を厳密に分離しており、高次のシンプレクティック積分を実現するために交互配置された時間格子を利用している。

数値実験結果

著者らは、定数およびルッドマン・ボルテックススの速度場における正弦波および正弦ベルの移流をテストするため、周期境界条件を持つ2次元領域を用いた数値実験を提示している。

• 収束性: 本手法は、時間積分の誤差が空間誤差を支配しない限り、$(p+1)$ 次の代数的 $h$-収束およびスペクトル $p$-収束を示す。
• 時間積分の精度: 解の誤差は、時間積分スキームの次数 ($p_t$) に依存することが示されている。時間スキームが十分に正確であれば、初期の離散化誤差が保持されるが、そうでなければ人工的な分散が増大する。
• 質量保存: 総質量誤差は最初の $10^3$ ステップまでマシンゼロ（機械精度）に留まり、$2 \times 10^4$ ステップ後でも $10^{-12}$ を下回っており、優れた保存特性を示している。
• 可逆性: ルッドマン・ボルテックス内での正弦波の移流に続き、流れの方向を反転させるテストにより、本手法が初期解を回復できることが示され、積分法の可逆性が証明された。
• 人工分散: 本論文では、人工分散は時間積分の次数に依存する旨を述べている。低次の時間スキームは、移流される波の周波数に依存する分散誤差を導入する。

重要性と主張

本論文は、導出されたスキームが、物理適合的（ミメティック）なフレームワーク内で、スペクトル精度と局所的質量保存を成功裏に組み合わせていると主張している。カルタンのホモトピー公式を用いてリー微分をミメティック・フレームワークに組み込むことにより、著者らは、移流方程式の幾何学的構造を明確にする離散化手法を提供している。本手法は、物理法則の根底にある幾何学的およびトポロジカルな性質を保存することが不可欠な移流問題を解くための、堅牢なアプローチとして提示されている。著者らは、空間離散化は非常に高精度であるものの、全体の精度と分散特性は使用される時間積分の次数に依存することを控えめに注記している。