Fully discrete analysis of pressure-robust enriched Galerkin methods for the time-dependent Stokes equations

技術要約：時間依存型ストークス方程式に対する、完全離散化された圧力ロバストな高次補完ガラーキン法の解析

問題提起
本論文は、粘性を持つ非圧縮性流体の流れを支配する、時間依存型ストークス方程式の数値解法について扱っている。古典的な混合有限要素法（テイラー・フドー法やクロイツィー・ラヴィアート法など）における主要な課題は、「圧力ロバスト性（pressure-robustness）」の欠如である。標準的な定式化では、離散化された速度誤差が、連続的な圧力および強制項の非回転成分に直接依存してしまう。その結果、たとえ厳密な速度が滑らかであっても、圧力が大きい場合や急峻な勾配を持つ場合、速度近似の精度が著しく低下する。この問題は、低粘性領域において特に深刻であり、誤差境界に粘性パラメータ（$\mu^{-1}$）の逆数項が含まれることで、非物理的な結果を招くことがある。定常問題や他の様々な定式化（$H(\text{div})$適合空間や仮想要素法など）には圧力ロバストな手法が存在するが、時間微分項に速度再構成を組み込んだ高次補完ガラーキン（Enriched Galerkin; EG）法を用いた、時間依存型ストークス系に対する完全離散化されたパラメータ明示的な解析は、これまで確立されていなかった。

手法
著者らは、連続ガラーキン（CG）法の計算効率と、不連続ガラーキン（DG）法の局所的な質量保存特性を組み合わせた、完全離散化されたEGフレームワークを提案している。

1. 空間離散化: 速度空間 $V_h$ は、連続な1次ラグランジュ空間（$C_h$）を、要素ごとに定義された不連続な1次の平均ゼロベクトル空間（$D_h$）によって高次補完することで構築される。圧力空間 $Q_h$ は、区分定数関数からなる。厳密な局所質量保存を実現し、圧力ロバスト性を達成するために、速度再構成演算子 $R: V_h \to BDM_1(\mathcal{T}_h)$ が用いられる。この演算子は、EG基底の連続成分を保持しつつ、不連続な補完成分を $H(\text{div})$ 内の最低次ラヴィアート・トーマス（$RT_0$）部分空間へと写像する。
2. 時間離散化: 著者らは、後退オイラー法（1次精度）とクランク・ニコルソン法（2次精度）の両方を用いた完全離散化スキームを開発している。重要な手法上の革新は、再構成演算子 $R$ を、強制項だけでなく、離散的な時間微分項にも組み込んだ点にある。具体的には、標準的な差分商 $(u_h^{m+1} - u_h^m)/\tau$ が $(R u_h^{m+1} - R u_h^m)/\tau$ に置き換えられる。
3. 解析: 著者らは、無条件安定性を確立し、先験的な誤差評価式を導出している。解析では、メッシュサイズ（$h$）、時間ステップサイズ（$\tau$）、および粘性（$\mu$）への依存関係を明示的に追跡している。

主な貢献

• 完全な圧力ロバストEG定式化: 不連続な補完成分を $H(\text{div})$ 適合空間へと写像することで、厳密な局所質量保存を実現する、圧力ロバストなEG空間離散化を導入した。
• 時間微分における再構成: 時間依存型問題において圧力ロバスト性を維持するためには、速度再構成を強制項だけでなく、離散的な時間微分項にも適用する必要があることを示した。これにより、定常問題で確立された手法を過渡現象へ拡張する際の制限が解消される。
• 最適次数の誤差評価: 無条件安定性を証明し、最適次数のパラメータ明示的な先験的誤差評価式を導出した。極めて重要な点として、速度の誤差境界は連続的な圧力および強制項の非回転成分に依存せず、粘性の逆数因子（$\mu^{-1}$）を含まない。
• 質量保存: 再構成された速度は、各要素上で点ごとに発散ゼロ（divergence-free）であることが示されており、各タイムステップにおいて厳密な局所質量保存を保証する。

結果
2次元および3次元における数値実験により、理論的知見が検証された。

• 収束率: 本スキームは予測された収束率を達成している。$\tau \sim h$ とした後退オイラー法では、速度に対して1次の収束が観察される。クランク・ニコルソン法において $\tau \sim h^{1/2}$ とした場合、速度誤差は時間に対して2次に近い収束（結合関係に応じて実質的に $O(h)$ またはそれ以上）を示しながら、ロバスト性を維持する。
• 粘性へのロバスト性: 低粘性領域（$\mu = 10^{-6}$ から $10^{-10}$）において、提案手法は精度を維持する。一方で、再構成を強制項にのみ適用した手法（SR-BE）は、最適な収束率を回復できず、誤差の悪化を示す。
• 質量保存: 3次元テストにおいて、再構成された速度は機械精度（$\approx 10^{-14}$）まで発散ゼロであり、一方で生の離散速度は、蓄積することのない非ゼロの発散欠損を保持している。
• 実用的な応用: 本手法は、矩形障害物の周囲の3次元流動に適用され、混合境界条件および低粘性下での安定性と、複雑な流れの特徴（剥離、再結合、圧力降下）を捉える能力が示された。

意義および主張
本論文の主要な貢献は、高次補完ガラーキン法を用いた、時間依存型ストークス方程式に対する完全離散化された圧力ロバストな解析を確立したことにある。著者らは、再構成演算子を離散的な時間微分項に含めることが不可欠であることを強調している。これを行わない限り、定常問題で確立された圧力ロバスト性は過渡現象には拡張されない。導出された誤差評価式は、粘性の逆数因子を含まないという点で重要であり、流体が低粘性になっても数値精度が低下しないことを保証している。本研究は、圧力ロバスト性と厳密な質量保存が極めて重要となる非圧縮性流体シミュレーションにおいて、EG法を使用するための厳密な理論的基礎と実用的な検証を提供するものである。