A dichotomy of finite element spaces and its application to an energy-conservative scheme for the regularized long wave equation

技術要約：有限要素空間の二分性と、正則化長波方程式に対するエネルギー保存スキームへの適用

問題提起
本論文は、エネルギー保存型ガラーキン離散化において観測される特定の現象、すなわち、有限要素空間の多項式次数に基づく収束挙動における特異な二分性について論じている。具体的には、連続ラグランジュ有限要素空間を用いる場合、奇数次の多項式次数では最適な収束率が得られる一方で、偶数次の多項式では収束率が低下するという現象である。著者らは、この挙動が特定の数値スキームによるアーティファクトではなく、有限要素空間自体の構造に内在するものであることを証明することを目的としている。さらに、本論文では、質量、インパルス（運動量）、およびエネルギーの保存則を持つ、長波表面波の基本モデルである正則化長波（RLW）方程式（またはベンジャミン・ボナ・バミー（BBM）方程式）へのこの特性の影響を調査している。

手法
解析は、一様な周期格子上で定義された一次元連続ラグランジュ有限要素空間（$V^k_h$）の枠組み内で行われる。核心となる手法は、有限要素空間 $V^k_h = S^1_h \oplus S^2_h$ の直和分解である。ここで、$S^1_h$ は頂点ベースの関数からなり、$S^2_h$ はバブル関数からなる。

1. $L^2$ 射影の解析： 著者らは、導関数の $L^2$ 射影 $P$ を厳密に解析している。彼らは、十分に滑らかな周期関数に対して、誤差項 $\|P[(Pu - u)_x]\|$ が奇数次多項式次数（$k = 2\kappa + 1$）においてのみ「超近似（super-approximation）」特性を示し、$O(h^{k+1})$ のオーダーを達成することを証明している。対照的に、偶数次（$k = 2\kappa$）の場合、誤差は標準的な最適次数 $O(h^k)$ に留まる。この区別は、直交多項式および特定の基底関数（$\psi$ や $\phi_j$ など）の性質を利用して、頂点部分空間とバブル部分空間への射影に関する詳細な推定を通じて導出されている。
2. 半離散定式化： RLW方程式のための修正ガラーキン有限要素法が提案されている。このスキームは、質量と三次エネルギー汎関数を共に保存するように設計されている。半離散定式化は、非線形項および分散項の $L^2$ 射影を含む常微分方程式系として表現される。
3. 誤差解析： 半離散定式化に対する理論的なア・プリオリ誤差評価が確立されている。この証明は、第2節で確立された導関数の $L^2$ 射影の二分性に基づいている。解析では、奇数次と偶数次の $k$ を区別することで、解 $u$ およびその導関数 $u_x$ の具体的な収束率を導出している。
4. インパルス近似： 提案されたスキームは質量とエネルギーを保存するが、インパルス汎関数を厳密には保存しない。著者らはインパルスの近似精度を分析し、それが多項式次数に依存した高い精度で保持されることを示している。
5. 完全離散スキーム： エネルギー汎関数を保存しながら時間方向にこの半離散系を積分するために、リラクゼーション・ルンゲ＝クッタ法を採用している。これらの手法は、離散エネルギーの正確な保存を保証するために、リラクゼーション・パラメータ $\gamma_n$ を介して時間ステップを調整する。

主な貢献と結果

• 理論的二分性： 主要な貢献は、導関数の $L^2$ 射影の超近似特性が、奇数次のラグランジュ有限要素空間に特有のものであるという証明である。奇数次 $k$ に対しては $\|P[(Pu - u)_x]\| = O(h^{k+1})$ であり、偶数次 $k$ に対しては $O(h^k)$ である。
• RLW の収束率： この特性をRLW方程式に対するエネルギー保存スキームに適用することで、著者らは以下を確立した：
  • 奇数 $k$ の場合：数値解は最適次数 $O(h^{k+1})$ で収束し、導関数近似 $u_x$ もまた $O(h^{k+1})$ のオーダーで収束する（超収束）。
  • 偶数 $k$ の場合：収束は劣最適であり、解は $O(h^k)$ で、その導関数は $O(h^{k-1})$ で収束する。
• インパルス保存： 本論文は、修正スキームがインパルス汎関数を厳密には保存しないものの、インパルス近似の誤差が、奇数 $k$ に対しては $O(h^{2k})$、偶数 $k$ に対しては $O(h^{2k-2})$ のオーダーであることを示している。
• 数値検証： 広範な数値実験により、理論的予測が裏付けられた。誤差 $\|P[(Pu - u)_x]\|$ およびRLW解の収束率は、導出された境界と一致している。さらに、シミュレーションにより、質量とエネルギーはマシン精度まで保存される一方で、インパルス誤差は小さく、予測された収束パターンに従うことが示された。

意義
本論文は、エネルギー保存型スキームにおける収束挙動の二分性が、数値的な異常ではなく、使用される有限要素空間の根本的な特性であることを主張している。導関数の $L^2$ 射影を介して、多項式次数のパリティ（奇偶）に直接収束率を結びつけることで、本研究は、先行文献（例：[28]）における過去の経験的な観察に対する理論的な説明を提供している。これらの結果は、三次エネルギー汎関数と導体射影を伴う問題において、奇数次の多項式を選択することは単なるヒューリスティック（経験則）ではなく、最適な精度を達成するための必然であることを示唆している。提案されたRLW方程式に対するエネルギー保存スキームは、適切な有限要素空間が選択されている限り、高次の精度と基本的な物理量（質量およびエネルギー）の保存を同時に達成できることを示す、具体的な応用例となっている。