A pressure-robust HHO method for the solution of the incompressible Navier-Stokes equations on general meshes

技術要約：一般的なメッシュ上における非圧縮ナビエ・ストークス方程式のための、圧力頑健なHHO法

問題設定
本論文は、非圧縮ナビエ・ストークス方程式の数値近似、特に圧力頑健性（pressure robustness）の達成という課題に取り組んでいる。標準的な離散化では、速度場の誤差評価が圧力に依存することがあり、これは物体力が大きな回転なし成分（すなわち、スカラーポテンシャルの勾配）を含む場合に性能が悪化するという問題を引き起こす。標準的な有限要素法や一致する単体メッシュに対しては圧力頑健な手法が開発されているが、これらを一般的な多面体メッシュ（任意の多面体要素や非一致界面をサポートするもの）へと拡張することは、依然として大きな課題である。著者らは、彼らの先行研究 [11] における未解決の問い、すなわち「一般的な多面体メッシュに対して圧力頑健なハイブリッド高次（HHO）法を構築できるか」に対し、回答を与えることを目的としている。

手法
提案手法は、2次元および3次元で有効な、完全な圧力頑健性を備えたHHOスキームである。核心となる革新は、各メッシュ要素 $T$ のサブトリギュレーション（部分細分化）上で定義された、局所的な発散保存型速度再構成 ($R^k_T$) の構築にある。

1. メッシュとサブメッシュ: 本手法は、多面体メッシュ $\mathcal{M}_h$ 上で作動する。極めて重要な点として、各要素 $T$ 内に一致する単体サブメッシュ $\mathfrak{T}_T$ が存在することを仮定している。特定の幾何学的仮定として、$\mathfrak{T}_T$ 内のすべての単体は共通の頂点 ($x_T$) を共有しており、これにより錐体状または非錐体状のサブメッシュが可能となる。
2. 発散保存型再構成: 与えられた離散速度未知数 $\mathbf{v}_T$ に対し、再構成 $R^k_T \mathbf{v}_T$ は、サブメッシュ $\mathfrak{T}_T$ 上の局所的な混合問題の解として得られる。この問題は、以下の条件を満たす Raviart–Thomas–Nédélec 空間 $RT_k(\mathfrak{T}_T)$ 内の関数を求めるものである：
   • 要素境界における離散速度の法成分と一致すること。
   • 離散的な発散を保存すること（離散発散作用素 $D^k_T$ と可換であること）。
   • Koszul 空間 $G_{c, k-1}(T)$ および要素の内部に関して直交条件を満たすこと。
     この構成は、制限的な面の幾何形状を必要としたり、勾配の有界性に関する理論的な未解決問題を抱えていたりした従来のアプローチ（例：Wachspress 座標への依存）の限界を回避している。
3. 離散定式化:
   • 粘性項: サブメッシュ上の標準的なHHO勾配再構成を使用する。
   • 物体力: 力 $\mathbf{f}$ は、生の離散未知数ではなく、再構成された速度 $R^k_h \mathbf{v}_h$ に対してテストすることで離散化される。これは、回転なしの力（勾配）が速度解に影響を与えないことを保証するために、発散保存特性を利用している。
   • 対流項: 再構成された速度 $R^k_h$ と拡張された勾配再構成 $G^{2(k+1)}_{\mathfrak{T}_T}$ を用いた、新しい三線形形式が導入されている。この定式化は、圧力頑健性を維持しながら、非散逸性と一貫性を確保する。

主な貢献

• 一般的なメッシュへの対応: 本論文は、圧力頑健なHHO法を、単体や特定の多面体幾何形状への制限を克服し、一般的な多面体要素からなるメッシュへと拡張できるという肯定的な回答を提供している。
• 新規の再構成: サブトリギュレーション上の局所的な混合問題に基づく、発散保存型速度再構成の導入。このアプローチは、サブメッシュが共通頂点仮定を満たす限り、任意の多項式次数 $k \ge 0$ および一般的な多面体に対してウェルポーズ（適正）であることが示されている。
• 理論的解析: 著者らは厳密な誤差評価を確立している。速度誤差が圧力および物体力の回転なし成分（ホッジ分解におけるパラメータ $\lambda$）に依存しないことを証明している。収束率は、標準的な正則性仮定の下で、エネルギーノルムにおいて $O(h^{k+1})$、速度の $L^2$ ノルムにおいて $O(h^{k+2})$ である。
• 頑健性の証明: 理論的結果は、本手法が大きな回転なしの物体力に対して頑健であることを示しており、これは標準的なHHOスキームには備わっていない特性である。

数値結果
理論的な知見は、2次元における包括的な数値テストを通じて検証されている。

1. Kovasznay 流: Cartesian（直交格子）、hexagonal（六角格子）、および Kershaw メッシュにおいて収束率が検証されている。結果は、提案手法と標準的なHHOの両方について、最適な収束次数（エネルギーノルムで $h^{k+1}$、$L^2$ 速度誤差で $h^{k+2}$）を確認している。
2. 頑健性テスト: 大きな回転なしの物体力（$\lambda = 10^6$）を持つ問題が解かれている。提案手法は、$\lambda$ の大きさに関わらず、厳密な解を再現するか（$k=1$ の場合）、急速に収束する。対照的に、[8] の標準的なHHO法は、メッシュ細分化を行っても減少しない大きな誤差を示しており、圧力頑健性の欠如を裏付けている。
3. Lid-Driven Cavity（蓋駆動空洞流）: レイノルズ数 $Re=1000$ において、古典的な蓋駆動空洞問題に適用されている。結果は文献の参照解と一致している。さらに、回転なしの物体力を加えた場合でも、速度プロファイルは変化しないのに対し、標準的な手法は異なる誤った解へと収束する。

意義と主張
本論文は、完全に圧力頑健なHHO法が、多面体要素からなる一般的なメッシュ上で動作することを主張している。その意義は、多面体手法の柔軟性（複雑な幾何形状や非一致界面の処理）と、圧力頑健なスキームの物理的忠実度（回転なしの力を正しく扱う能力）の間の溝を埋めることにある。著者らは、本構成がサブトリギュレーション上の局所的な混合問題に依拠していることを強調しており、これにより、以前の多面体アプローチを悩ませていた幾何学的制限（三角形の面のみなど）や理論的な未解決問題（Wachspress 座標の有界性など）なしに、必要な連続性と近似の有界性を導出できるとしている。本研究は、任意の次数の精度を維持しながら、圧力頑健な離散化の適用範囲をより広範なクラスのメッシュへと拡張するものである。