Embedded Trefftz DG method for steady Navier-Stokes flow. Part I: Oseen linearization

技術要約：定常Navier–Stokes流のための埋め込み型Trefftz DG法（第I部：Oseen線形化）

問題設定
本論文は、既知の対流場 $w$ の周りでの定常非圧縮Navier–Stokes方程式の線形化として機能する、定常Oseen問題の数値近似を取り扱う。支配方程式は、有界領域 $\Omega \subset \mathbb{R}^d$ ($d=2,3$) において、境界条件として均一なディリクレ境界条件を満たす速度 $u$ および圧力 $p$ を求めるものである：
$$-\nu\Delta u + (w \cdot \nabla)u + \nabla p = f, \quad \text{div } u = 0$$
主な課題は、対流項 $(w \cdot \nabla)u$ によって導入される変数係数を扱いながら、減少した自由度（DoF）で高い近似精度を維持する、埋め込み型Trefftz不連続ガラーキン（DG）法の開発である。Stokes問題（$w=0$ の場合）とは異なり、Oseen作用素は変数係数を含むため、標準的な定数係数のTrefftz構成では不十分である。

手法
著者らは、局所的なTrefftz基底を明示的に構築することを回避する、埋め込み型Trefftz-DGアプローチを提案している。この手法は、標準的な多項式DG空間 $X$ から出発し、射影制約を通じて局所的にPDEを課すものである。

1. 局所作用素とTrefftz空間：
   各要素 $T$ 上の局所作用素 $K_{T,w}$ を、解空間を双対テスト空間 $Q(T)'$ へ写すものとして定義する。Trefftz空間 $T_w$ は、$X$ の部分空間であり、局所的な残差が $Q(T)$ に対して弱の意味で消滅する関数からなる。具体的には、$T_w$ 内の関数は、低次の多項式（速度については次数 $k-2$、圧力については $k-1$）に対してテストされた場合にのみ、PDEを満たす。

2. 局所・全体分割と補空間：
   重要な理論的ステップは、 $X(T) = T_w(T) \oplus L(T)$ となる局所補空間 $L(T)$ の構築である。

   • 速度成分 ($L_u$): サブボール上で消滅するバブル関数を用い、ファセットのフラックスを制御するための特定の多項式項で強化され、さらに $d$ 個のファセット上で平均フラックスがゼロになるよう制約されたものを用いて構成される。
   • 圧力成分 ($L_p$): 平均ゼロの次数 $k-1$ の多項式として定義される。
   • 可逆性: 本論文では、対流場 $w$ が解像条件（具体的には、粘性に対する離散ノルム $|w|_{h,d}$ に関する微小性の条件）を満たす場合、局所作用素 $K_{T,w}$ が $L(T)$ 上で有界に可逆であることを証明している。この可逆性により、解の局所的な部分を各要素上で独立して計算することが可能となる。

3. 縮退定式化：
   （局所的な作用素によって速度と結び付けられている）高次の圧力成分を局所的に排除することにより、著者らは速度未知数のみに関する**縮退定式化（reduced formulation）**を導出する。この定式化は以下の特性を持つ：

   • 局所レベルでの鞍点構造を排除する。
   • Stokesの設定から縮退Oseen-Trefftz空間へと強凸性を転送する。
   • 離散的な発散制約を課すために、定数圧力成分に対してのみグローバルな結合を保持する。

4. 安定性解析：
   解析では以下を確立している：

   • 双線形形式の連続性。
   • Trefftz部分空間における縮退形式の強凸性。
   • 標準的なDGノルムにおける完全な結合系のinf-sup安定性。
     これらの結果は、抽象的な局所・全体フレームワークと $L(T)$ の具体的な構成に依拠しており、ウェルポーズドネス（適切解の存在）と準最適性を保証している。

主な貢献
本論文は、埋め込み型Trefftz-DG手法をStokes問題から定常Oseen方程式へと拡張し、非線形Navier–Stokes解析（第II部）のための線形的な基礎を提供している。具体的な貢献は以下の通りである：

• 局所補空間の構築: 定数係数のテイラー多項式の議論に依存していたStokesの場合とは異なり、本研究では対流項から生じる変数係数の局所作用素を扱うための $L(T)$ を構築している。
• 非斉次右辺項の取り扱い: 修正項を用いた議論ではなく、非斉次項を局所問題の定式化に直接組み込んでおり、抽象的な局所・全体分割フレームワークに自然に適合させている。
• 縮退定式化と安定性: 縮退した速度のみの定式化を導出し、完全なOseen作用素に対する縮退inf-sup安定性を証明している。これにより、局所的な鞍点構造が排除され、対流場の適切な微小条件の下で強凸性の転送が可能となる。

結果

• 理論的側面: 本論文は、標準的なDGノルムにおける完全な安定性と準最適性を証明している。手法がウェルポーズドであり、誤差が離散空間における最良近似誤差によって（対流場の解像度に依存する定数でスケールされて）抑えられることを確立している。
• 数値的側面:
  • Oseen–Kovasznayテスト: 規定されたOseen–Kovasznay流を用いた数値実験により、Trefftz-DG法が期待通りの漸近収束率（速度に対して $O(h^{k+1})$、圧力に対して $O(h^k)$）を達成することを示している。
  • 効率性: 標準的なDG法は固定された多項式次数において誤差定数がわずかに小さいものの、Trefftz-DG法は、特に高次の多項式を用いる場合において、与えられた誤差レベルに達するために必要な自由度が大幅に少ない。
  • 計算コスト: タイミングの結果は、Trefftz-DG法が一貫して標準的なDG法よりも高速であることを示しており、その優位性は細かいメッシュや高次化において増大する。
  • 限界: 本手法は対流支配領域（低粘性）において感度を示し、標準的なDGよりも早く誤差が悪化する傾向がある。これは、局所的な埋め込み型Trefftz構成の解像度とコンディショニングに起因すると考えられる。

意義
本論文は、第II部で展開される非線形Navier–Stokes理論に必要な線形解析的基盤を提供すると主張している。変数係数を持つOseen問題への埋め込み型Trefftzアプローチの拡張に成功することで、著者らは、明示的な基底構築なしに、変数係数の不均質な作用素に対してTrefftz手法が適用可能であることを示した。縮退定式化は、局所的な圧力をグローバルな系から分離し、反復解法の安定性解析を簡素化するという点で、非線形問題の解析における極めて重要なステップとして強調されている。本研究は、Oseen問題を古典的なDG離散化とTrefftz型のモデル低減の間の架け橋として位置づけ、粘性非圧縮流の高次・低自由度シミュレーションへの道筋を示している。