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

技術要約：コルモゴロフ方程式に対する時空不連続ガラーキン法の漸近的数値的ハイポコアシビティ（弱双共鳴性）

問題提起
本論文は、退化した拡散を示すキネティック型方程式の代表的なモデルである、古典的なコルモゴロフ方程式の数値離散化について扱っている。この方程式は、有界な空間領域 $\Omega \subset \mathbb{R}^2$ 上で定義され、$u_t + Lu = 0$（ここで $L = -\partial_{xx} + x\partial_y$）という形式をとる。この演算子の決定的な特徴は、拡散が $x$ 方向のみに発生し、輸送項 $x\partial_y$ が $y$ 方向において作用することである。$y$ 方向には明示的な拡散が存在しないにもかかわらず、このシステムは「ハイポコアシビティ（弱双共鳴性）」という性質を示し、輸送と退化した散逸の相互作用によって、空間領域全体に減衰が伝播し、平衡状態への指数関数的な収束をもたらす。

中心となる課題は、標準的な時空不連続ガラーキン（dG）法、特に古典的な空間・時間dG法が、このハイポコアシビティ構造を離散レベルで保持できることを証明することである。ハイポコアシビティは連続系の設定においては十分に理解されているが（主に Villani, Hérau, Nier らの研究による）、ガラーラーキン法におけるその保存に関する研究はほとんど進んでおらず、これまでの数値的な結果は、特殊で非標準的な定式化を必要とする非適合有限要素法や有限差分法、有限体積法に限定されていた。

手法
著者らは、コルモゴロフ方程式に適用された標準的な時空dG法を分析している。空間離散化には、輸送項に対してアップウィンドdGスキームを、拡散項に対して対称内部ペナルティdGスキームを用いている。時間離散化には、多項式次数 $q$ を用いた標準的なdGアプローチを用いている。

手法の核心は、連続的なハイポコアシビティの枠組みを模倣した修正安定性解析を構築することにある。連続系におけるハイポコアシビティは、内積を混合微分や特定の勾配項を含む形に修正することで、実質的に「より強い」ノルムを構築し、スペクトルギャップを得ることで証明されることが多い。

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)$$: 欠落している $y$ 微分の制御を回復するように設計された、修正行列 $A$（$h$ および $p$ に依存）と勾配のダイバージェンスを含む項。

極めて重要な点として、著者らは、dGで使用される不連続多項式空間において、$V$ が $U$ と同じ離散空間内に留まることを指摘している。これにより、従来のペトロフ・ガラーキン項（連続要素を用いた手法で必要とされるもの）を導入することなく、標準的なガラーキン法において $V$ を許容可能なテスト関数として使用することが可能となる。

主要な貢献と結果

1. 標準的dGにおける離散的ハイポコアシビティ: 本論文は、古典的な時空dG法が漸近的数値的ハイポコアシビティ特性を持つことを確立している。 「ハイポコアシビティ強化」されたノルムにおける一方向のinf-sup条件を証明することにより、著者らは、数値解が長時間において平衡状態へ指数関数的に減衰することを実証している。これは、標準的なガラーキン法に関する初の成果である。

2. 強化ノルムとスペクトルギャップ: 全勾配を含む一族のノルムを導入している。著者らは、このノルムにおけるメッシュ依存のスペクトルギャップ（ポアンカレ型不等式）を証明している。減衰率 $\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)$ の条件下での）時間ステップ $k$ に関して最適である。

5. 数値検証: 数値実験により、理論的な収束率を確認し、解が平衡状態へと指数関数的に減衰することを実証した。実験の結果、理論的な減衰率は $h$ と $p$ に依存するものの、平衡への収束の様子は堅牢であることが示された。

意義と主張
著者らは、本研究が、既知の古典的なガラーキン法に対する（漸近的）数値的ハイポコアシビティの最初の実証であると主張している。その意義は、以前の非適合有限要素アプローチで必要とされた、追加のペトロフ・ガラーキン項による双線形形式の修正を必要とせずに、ハイポコアシビティブな構造を保持している点にある。

本論文は、証明手法が汎用的であることを強調している。連続的なハイポコアシビティ理論で使用される修正内積の効果を模倣するテスト関数を構築することで、このアプローチは他のキネティック型方程式や退化した放物型問題にも拡張できる可能性がある。これらの結果は、退化した拡散問題を長期間シミュレーションするための標準的な高次時空dG法を使用するための厳密な基礎を提供するものであり、数値スキームが連続モデルに内在する物理的な散逸メカニズムを人工的に抑制しないことを保証している。