Topology-Guided Schur Recycling for Parametric Incompressible CFD of Separated Flows

技術要約：剥離流のパラメトリック非圧縮CFDにおけるトポロジー誘導型シュア・リサイクル

問題提起
幾何学的形状、強制力、またはレイノルズ数の変化に伴い、非圧縮流れの巨大な結合速度・圧力系を数百回、あるいは数千回にわたって解く必要があるパラメトリック計算流体力学（CFD）では、計算上の大きなボトルネックが生じる。剥離流においては、これらのサドルポイント系の再構築と求解の繰り返しが大きな負担となる。ソルバーの再利用戦略（低次元モデルやクリロフ部分空間のリサイクルなど）は、これらの計算を加速できる可能性があるが、科学的なリスクを伴う。すなわち、積極的な再利用は、再循環領域の生成、消滅、または分裂といった重要な物理的遷移を平滑化したり、見逃したりする可能性がある。ソルバーは数値的に妥当な残差を生み出す一方で、基礎となる流れのトポロジー（例：剥離構造）が根本的に変化している場合があり、その結果、物理的に誤った結果を導く恐れがある。現在の手法には、再利用された解を受け入れる前に、代数的な整合性と物理的なトポロジーの安定性を同時に検証するメカニズムが欠けていることが多い。

手法：CETERフレームワーク
本論文は、局所的に摂動を与えられた幾何形状におけるパラメトリック非圧縮CFDのために設計された、CETER（Topological Event-Triggered Schur Recycling：トポロジー事象駆動型シュア・リサイクル）フレームワークを提案している。この手法は、再利用されたソルバーのステップを受け入れてよいかを判断するために、4つの異なるコンポーネントを統合している。

1. 局所ランク演算子圧縮： 本フレームワークは、局所的な幾何学的摂動（例：壁面の変形や狭窄）が離散演算子の限定的な部分集合にのみ影響を与えることを前提としている。現在の演算子 $K(\mu)$ と参照演算子 $K(\mu_0)$ の差は、影響を受けるセルパッチに対する特異値分解（SVD）を用いて、低ランク更新（$U M V^T$）へと圧縮される。
2. Woodbury型のシュア・リサイクル： 圧力シュア補完 $S(\mu) = -C - B A^{-1} B^T$ をゼロから再構築する代わりに、CETERはWoodburyの行列恒等式を用いて、逆運動量ブロック $A^{-1}$ の作用を更新する。これにより、高密度なシュア行列を明示的に形成することなく、低次元の補正を通じてシュア補完の作用を修正することが可能になる。
3. 事後検証（A Posteriori Certification）： すべての再利用された求解は、3つの代数チェックを受ける。
   • 残差チェック： 相対的なサドルポイント残差 $\eta_{lin}$。
   • 発散チェック： 離散的な非圧縮性の欠損 $\eta_{div}$。
   • 格子収束性： （オプション）監査実行のための格子収束指数（GCI）。
4. パーシステント・ホモロジーによるトポロジー・ゲート： これが本手法の革新的な核となる部分である。スカラーの分離指標場（例：符号付き軸方向速度、壁面せん断応力、または渦度大きさ）が、パーシステント・ホモロジーを用いて解析される。本手法は、その場の超レベル集合（superlevel sets）に対するパーシステンス図を計算する。「トポロジカルな分離欠損」$\eta_{top}$ は、再利用された解のパーシステンス図と、信頼できる比較対象（完全に再構築された解、またはより細かい格子解）との間の2-ワッサーシュタイン距離を用いて算出される。

決定ロジック
トポロジー安全受容指数（$C_{top}$）が定義され、代数的な残差、発散の欠損、およびトポロジーの欠損を組み合わせる。パラメータのステップは、$C_{top} \leq 1$ かつすべての個別指標が規定の許容範囲内である場合にのみ受け入れられる。もし指数が1を超えた場合、ステップは拒絶され、パラメータの増分が減少され、あるいは近傍での完全な再構築またはメッシュ細分化がトリガーされる。

主な貢献
本論文は、以下の5つの具体的な貢献を概説している。

• (C1) 離散化された非圧縮ナビエ・ストークス／オーセン（Oseen）系における、幾何学的摂動に起因する局所ランクモデルの定式化。
• (C2) 高密度なシュア補完の明示的な形成を回避する、Woodbury更新を用いたシュア補完のリサイクル戦略の導出。
• (C3) 再循環に関連するスカラー場のパーシステンス図を利用して、構造的変化を検出するトポロジー・ゲートの導入。
• (C4) パラメータステップの受け入れ、削減、再構築、または細分化を制御するための、トポロジー安全受容指数の定義。
• (C5) 製造された演算子チェック、蓋駆動キャビティの監査、チャネル継続のトポロジー・ゲート適用、およびスケール正規化された比較による背面段差型監査を含む、検証経路。

結果と検証
本論文は、予測的な高精度エンジニアリング・シミュレーションではなく、検証指向の数値結果を報告している。主な知見は以下の通りである。

• 代数的整合性： 行列レベルのテストにより、低ランク更新に対してWoodbury更新が完全な直接解を機械精度に近い精度で再現することが確認された。
• 演算子の検証： 製造された微分演算子のチェックにより、有限差分実装において2次収束性が示された。
• トポロジー・ゲーティングの有効性： 局所的な抵抗パッチを持つパラメトリック・チャネル・テストにおいて、残差のみに基づく戦略ではすべてのパラメータステップを受け入れてしまう。しかし、CETERのトポロジー・ゲートは、再循環成分の数とサイズが変化した地点である $\mu = 0.3$ からステップを正しく拒絶した。パーシステンス距離による監査は、$\mu = 0.2$ と $\mu = 0.3$ の間の遷移を、最も強力なトポロジカルな事象として独立して特定した。
• 背面段差監査： レイノルズ数スイープ（$Re = 100$から$400$）において、CETERは再循環構造が安定している場合（$Re \leq 250$）にはステップを受け入れ、構造が変化する場合（$Re \geq 300$）には再構築をトリガーした。再付着長に関する参照トレンドは、Armalyらによるデジタル化データと比較して最大相対誤差が4%未満であり、NASA/Driver–Seegmillerの乱流ベンチマークのオーダーとも一致した。
• パフォーマンス： テストされた密な行列プロトタイプにおいて、トポロジー・ゲート付きのリサイクルは、純粋な代数的リサイクルと比較してオーバーヘッドが発生した。しかし、物理的な変化に関わらずすべてのステップを受け入れてしまう「見逃されたトポロジー事象」を、代数的リサイクルのみを使用した場合と比較して、効果的に回避することに成功した。

意義と主張
本論文の主要な貢献は、代数的なリサイクルツール（Woodbury恒等式、低ランク更新）と、トポロジカル・データ解析（パーシステント・ホモロジー）を、単一のソルバー再利用決定規則へと結合させたことにある。著者らは、トポロジーは単なる後処理の可視化ツールではなく、計算的な事象検出器であることを強調している。

その意義は、剥離流のパラメトリック研究において、代数的な精度だけでは不十分であるという事実を確立した点にある。ソルバーは代数的に一貫している（低い残差を持つ）一方で、物理的には無効（誤った剥離トポロジー）である可能性がある。CETERは、代数的リサイクルを流れのトポロジー情報によってガードするための数学的およびアルゴリズム的な基礎を提供し、物理的な剥離の性質が変化する場合にのみ、ソルバーの再利用を中断させることを保証する。本論文は、これを検証指向のステップとして位置づけており、今後の課題には幾何学的に一致した研究や乱流モデルが含まれるべきであるとしつつ、現在の研究がフレームワークの内部一貫性と決定ロジックを確立したものであると述べている。