Update-Magnitude State Redistribution (UM-SRD): A Shut-off Extension of Weighted SRD for Cut-Cell Methods

技術的概要：更新大きさ状態再分配（UM-SRD）

問題定義
カットセル法は、体積適合メッシュなしで複雑な幾何学形状を有限体積法で扱えるようにするが、任意に小さなカットセルに起因する深刻な安定性の課題を導入する。状態再分配（SRD）法、特に Berger & Giuliani（2024）が開発した重み付き SRD フレームワークは、小さなセルを近傍にマージし、マージされたメッシュ上で平均を進化させ、結果を元のメッシュに再分配することで、これらの安定性問題に対処する。

しかし、既存の重み付き SRD 定式化の決定的な限界は、「シャットオフ」メカニズムの欠如にある。フラックス再分配の古典的手法では、フラックスの不均衡がゼロの場合、ポストプロセッシングは自然に消滅するのに対し、重み付き SRD は有限体積更新の大きさに関わらず、すべての時間ステップで非自明な再分配演算子を適用する。この継続的な活性化は、2 つの実用的な欠陥をもたらす：

1. 定常状態の保存失敗：フラックスのバランスが厳密にゼロであっても、SRD は点ごとの値を体積重み付き平均に置き換えるため、機械精度への収束を妨げ、静的領域に $O(1)$ の誤差を導入する。
2. 不要な散逸：小さく非ゼロの更新を伴う滑らかな領域において、再分配の継続的な適用は、実際の物理的ダイナミクスとは無関係に、近傍の勾配に比例する数値拡散を導入する。

手法
本論文は、有限体積更新の局所的な大きさに基づいて SRD 演算子を単位演算子と混合させるように設計された、重み付き SRD のシャットオフ拡張である**更新大きさ状態再分配（UM-SRD）**を導入する。

本手法は標準的な 3 段階の時間ステップ枠組み内で動作する：

1. 有限体積更新：ベーススキームを用いて、再分配前の状態 $U^*$ を計算する。
2. 近傍集約：マージされた近傍 $M_j$ 上で、重み付き近傍平均 $\hat{Q}^*_j$ を計算する。
3. 混合再分配：SRD 演算子 $S_j$ を直接適用する代わりに、UM-SRD は凸結合を適用する：
   $$R_j(s_j) = (1 - s_j)I + s_j S_j$$
   ここで、$I$ は単位演算子であり、$s_j \in [0, 1]$ は混合パラメータである。

混合パラメータ $s_j$ は、滑らかで局所的に定義されたシャットオフ関数 $h(\eta_j)$ によって決定される：
$$s_j = \frac{\eta_j^p}{\eta_j^p + \tau^p}$$
ここで、$\eta_j$ は近傍 $M_j$ 内の最大絶対更新大きさの正規化指標である（具体的には $\eta_j = \frac{\max \|\Delta U_i\|}{\epsilon + \max \|\Delta U_k\|}$）。

• ゼロ更新領域：近傍内のすべての更新が機械精度レベル（$\Delta U_i \approx 0$）の場合、$\eta_j = 0$、$s_j = 0$ となり、手法は厳密に単位演算子に帰着する（$U^{n+1} = U^*$）。
• 活性領域：重要な更新が発生する場合（例えば、不連続点の近くや、小さなセルの剛性が活性化されている場合）、$\eta_j \approx 1$、$s_j \approx 1$ となり、手法は完全な重み付き SRD として振る舞う。

主要な解析的貢献
本論文は、単一の小さなカットセルを有する 1 次元線形移流モデルに対する UM-SRD の厳密な解析を提供する：

• 保存性：混合演算子は、任意の $s_j \in [0, 1]$ に対して、構成上保存的であることが証明される。
• 全変動減少（TVD）：ベースの 1 次風上スキームと同じ CFL 条件の下で、UM-SRD は元のカットセルメッシュ上で TVD であることが証明される。これは、単位演算子と TVD である SRD 演算子の凸結合に依存する。
• 定常状態の保存：UM-SRD がベーススキームの定常状態を厳密に保存することが証明される。有限体積更新がゼロであれば、$s_j=0$ となり、再分配は発生しない。これにより、標準的な SRD に内在する偽のドリフトを回避する。
• 精度：滑らかな領域において、UM-SRD はベーススキームの摂動であることが示される。混合パラメータがメッシュサイズ（または実務で使用される正規化指標）に応じて適切にスケーリングする場合、局所打ち切り誤差への追加寄与は、ベースの 1 次スキームよりも高次となる。したがって、UM-SRD は、基礎となる風上法の 1 次収束率を維持する。

数値結果
1 次元および 2 次元のモデル問題に対する数値実験は、理論的予測を確認する：

• 収束性：UM-SRD は、滑らかな移流に対して $L_1$ ノルムおよび $L_\infty$ ノルムで 1 次収束を示し、誤差はベーススキームおよび標準的な重み付き SRD と区別がつかない（差異 < 0.01%）。
• シャットオフ挙動：コンパクトなパルスを用いたテストにおいて、混合パラメータ $s_j$ はパルスが存在しない場合に厳密にゼロに低下し、ベーススキームを回復する。パルスが小さなセルと相互作用する際には、$s_j$ は約 0.99 まで上昇する。
• 安定性：UM-SRD は、フルメッシュ時間ステップ（小さなセルにおける局所 CFL が安定性限界の 10 倍である場合）においてベーススキームを安定化させ、ベース風上スキームで観測される指数関数的発散を防止する。
• 非振動性：不連続データを用いたテストは、UM-SRD が小さなセルの近くで標準的な SRD と同様の挙動を示し、非振動的であることを示す。
• 定常状態の保存：実験により、標準的な SRD が静的プロファイルに $O(1)$ の誤差を導入するのに対し、UM-SRD は初期プロファイルを機械精度まで維持することが確認された。
• 2 次元への拡張：本アルゴリズムは、埋め込み境界を有する 2 次元直交格子に成功裏に適用され、混合ロジックの変更なしに 1 次収束を示した。

意義
本論文は、UM-SRD が既存の重み付き SRD インフラに対する実用的で最小限の変更を提供すると主張する。確立された近傍および重み構造を再利用し、局所的な混合メカニズムのみを追加することで、SRD にゼロ更新特性を付与する。これにより、任意に小さなカットセルに対する状態再分配の安定性利点と、フラックス再分配手法の精度保持特性が組み合わされ、特に厳密な定常状態の保存を可能にし、滑らかな領域における不要な散逸を低減する。著者らは、現在の解析は 1 次元スカラー線形移流に限定されているが、アルゴリズム構造は次元に依存せず、生産用カットセルコードへの統合に準備ができていると指摘している。