[論文レビュー] A Hybridizable Discontinuous Galerkin Method for the non--local Camassa--Holm--Kadomtsev--Petviashvili equation

Mukul Dwivedi, Ruben Gutendorf+1 more|arXiv (Cornell University)|Jan 20, 2026

Nonlinear Waves and Solitons被引用数 0

ひとこと要約

論文はカーテシアン格子上で2D CH-KP方程式に対して高次ハイブリダイズド不連続有限要素法（HDG）を用い、非局所作用素を補助変数で局所化し、エネルギー安定性と部分最適 convergence 率を証明する。数値検証とピークオン分解能の結果を含む。

ABSTRACT

This paper develops a hybridizable discontinuous Galerkin method for the two-dimensional Camassa--Holm--Kadomtsev--Petviashvili equation. The method employs Cartesian meshes with tensor-product polynomial spaces, enabling separate treatment of \(x\) and \(y\) derivatives. The non-local operator \(\partial_{x}^{-1}u_{y}\) is localized through an auxiliary variable \(v\) satisfying \(v_x = u_y\), allowing efficient element-by-element computations. We prove energy stability of the semi-discrete scheme and derive \(\mathcal{O}(h^{k+1/2})\) convergence in space. Numerical experiments validate the theoretical results and demonstrate the method's capability to accurately resolve smooth solutions and peaked solitary waves (peakons).

研究の動機と目的

• 非局所項と高次非線形性を含む2次元CH-KPモデルを動機づけ、定式化する。
• テンソル積構造を保持し、要素毎の計算を可能にするHDG離散化を開発する。
• 半離散系のエネルギー安定性を確立し、空間的収束率を導出する。
• 非局所演算子を補助再構成変数で局所化し、計算を容易にする。
• 数値実験を通じて滑らかな解とピークオンの分解能能力を示す。

提案手法

• x-およびy-微分を分離するため、テンソル積のQ_k空間を用いたカーテシアンメッシュを使用する。
• 非局所演算子と高階微分を分解する補助変数を導入し、特に v に対して v_x = u_y などの関連変数 r, z, p, q, s を導入する。
• 要素を結合するグリッド殻上の3つのトレース変数を用いるHDGスキームを構成する。
• 非局所項の適切性を確保するためのプリミティブの再構成にアンカー条件を課す。
• グローバルな安定化パラメータ仮定（Assumption 3.1）下で半離散系のエネルギー安定性を証明する。
• 正則性仮定のもと、L2安定性結果を導出し、k ≥ 2 について O(h^{k+1/2}) の空間収束率を示す。

Figure 1 . Cartesian grid structure showing the element $K_{ij}$ with its boundary faces.

実験結果

リサーチクエスチョン

• RQ1 カーテシアン格子上の高次HDG離散化は2次元の非局所CH-KP方程式を正確に解けるか？
• RQ2 HDGフレームワーク内で非局所演算子 ∂_x^{-1}u_y を安定性を保ちながら局所化できるか？
• RQ3 CH-KP方程式に対する半離散HDGスキームのエネルギー安定性はどうなるか？
• RQ4 この非局所分散系に対する提案HDG法の空間的収束次数はどれくらいになるか？
• RQ5 数値実験で方法はピークオンと滑らかな解を高忠実度で解決できるか？

主な発見

• HDG離散化は指定された安定化仮定の下、半離散設定でエネルギー安定性を達成する。
• k ≥ 2 の場合、合理的な正則性仮定のもとで空間収束率が O(h^{k+1/2}) を達成する。
• 数値実験は理論的収束率を検証し、滑らかな解とピークオンを正しく解決できることを示す。
• 境界アンカーを用いた v = ∂_x^{-1}u_y の再構成法により一意性と要素毎計算の効率性を確保する。
• 本構成はマルチトレースLDG法と比べてグローバル自由度を低減しつつ安定性を維持し、収束性を改善する。

(a) Cross-section $x=0$