技術要約：高品質制約付きの予測的異方的グリッド・ワーピングを用いたFLIP液体シミュレーション

1. 問題提起

流体粒子法（FLIP）シミュレーションは、質量と運動量の追跡にラグランジュ粒子を用い、圧力投影にオイラー格子を用いる手法を組み合わせたものである。この手法は安定かつ効率的であるが、標準的なFLIPは固定された一様な格子に依存しているため、薄い液膜、急速に移動する自由表面、衝突によって生じる噴流、および強い局所的変形といった重要な流れの特徴に対して、十分な解像度を割り当てることができない。

動的なグリッド・ワーピング（Grid Warping）は、格子ノードを関心領域へと移動させることでこの問題に対処する。しかし、従来のスカラベースのワーピング（Ibayashiら [15] など）は、空間的な重要度（自由表面への距離など）に基づいてノードを集中させるが、方向への意識が欠如している。この不一致は、特定の軸に沿って引き伸ばしや圧縮が発生する急速に変形する流れにおいて、解像度の不足を招く。さらに、過度な異方的ワーピングは、しばしば深刻な局所的セルの歪み（圧縮、せん断、反転）を引き起こし、粒子・格子間の転送および圧力投影の安定性を低下させる。核心となる課題は、数値的不安定性を防ぐために局所的なセルの品質を明示的に制約しながら、大幅な適応的グリッド移動と方向的な解像度再分配を実現することである。

2. 手法

提案手法は、予測的変形推定、行列式正規化された異方的メトリック、および構造化ワープト・グリッドFLIPソルバー内でのヤコビアンに基づく品質制約を統合したものである。パイプラインは以下の通り進行する。

2.1 予測的自由表面の重要度および変形

グリッドの応答遅延を回避するため、本手法は現在の速度と重力を用いて、近未来の粒子位置 ($\hat{x}_p$) を予測する。

• 重要度: 予測された粒子に対してガウス・スプラッティング（Gaussian splatting）を行うことで、予測密度場を構築する。この場の勾配が正規化された表面重要度重み ($q_i$) を定義し、高解像度を必要とする領域を特定する。
• 変形方向: 近傍の流体セルに対する加重最小二乗法により、局所的な速度勾配 ($G_i$) を推定する。変形勾配テンソル ($F_i = I + \tau G_i$) を計算し、その関連するコーシー・グリーン・テンソル ($C_i = F_i^T F_i$) を解析する。最も小さい主伸長（principal stretch, $e_1$）に対応する固有ベクトルを、薄化や圧縮が発生する軸として、好ましい細分化方向として選択する。

2.2 行列式正規化された異方的メトリック

グリッドの総体的な体積を縮小させることなく、方向的に解像度を再分配するために、方向性メトリック $M$ を構築する。

• メトリックは、表面の重要度と変形の方向性を組み合わせる。
• 行列式正規化戦略により、アスペクト比 ($r$) が変形に合わせて変化する場合でも、メトリックの行列式が一定 ($\det(M) = h^{-6}$) であることを保証する。これにより、異方的な引き伸ばしを許容しつつ、グリッドの全体的なスケールが崩壊するのを防ぐ。
• 最終的なメトリックは、表面ベースの成分と変形ベースの成分をブレンドしたものであり、連続性を確保するためにスムージングされる。

2.3 品質制約付きグリッド・ワーピング

格子ノードは、ターゲットとなるメトリックを満たしつつ、厳格な品質制約に従って更新される。

• 最適化: ノードは、メトリック空間における距離と規則的なグリッド間の距離の差 ($E_{metric}$) を最小化するように更新される。
• 正則化: 以下の2つの追加エネルギー項が導入される：
  • $E_{vol}$: セルのヤコビアン行列式の行列式 ($\det(J)$) の1からの偏差にペナルティを課し、体積変化を制御する。
  • $E_{shape}$: 条件数 ($\kappa(J)$) を用いて、ターゲットメトリックからのセルの形状の偏差にペナルティを課す。
• 受理制約: 各増分更新の後、ラインサーチが行われる。候補となるグリッドは、以下を満たす場合にのみ受理される：
  • 体積比: $V_c / V^0_c \geq r_V$ (例: 0.30)。
  • 行列式の境界: $J_{min} \leq \det(J_c) \leq J_{max}$。
  • 条件数: $\kappa(J_c) \leq \kappa_{max}$。
    もし制約に違反した場合、ステップサイズは半分にされ、有効なステップが見つからない場合は反復が終了し、直前の有効な状態が保持される。

2.4 時間的トラッキング

計算コストを削減するため、メトリック変化ガイド付きのトラッキング戦略を採用する。フレーム間のターゲットメトリックの変化が閾値を下回る場合、フルな再最適化ではなく、簡略化された軽量な更新、または単一の更新を用いてワーピング最適化を行う。これにより、緩やかに変化する流れの状態における時間的コヒーレンスを利用する。

2.5 FLIPへの統合

ワープト・グリッドは、規則的なデカルト格子（Cartesian grid）の接続性を維持する。粒子・格子間および格子・粒子間の転送には、変形したヘキサヘドラル（六面体）セル上の三線形補間（trilinear mappings）を利用する。圧力投影は、変形した面における有限体積定式化を用いて行われ、発散が依然として高い場合には欠損補正（defect-correction）ソルブが適用される。

3. 主な貢献

1. 予測的異方的メトリック: 近未来の自由表面の重要度と局所的な変形方向を組み合わせ、行列式正規化を用いてグローバルなスケールを維持しながら空間解像度を再分配する手法。
2. 品質制約付きワーピング・スキーム: ヤコビアンによる正則化を、行列式、条件数、およびセル体積に関する明示的な受理テストと結合させたワーピングアルゴリズム。これにより、大規模なグリッド移動を可能にしつつ、局所的なセル歪みを厳格に制限できる。
3. メトリック変化ガイド付きトラッキング: ターゲットメトリックが緩やかに変化する場合に、繰り返しのグリッド最適化を減らす時間的戦略。これにより、グリッドの品質を損なうことなくパフォーマンスを向上させる。

4. 実験結果

本手法は、様々な2Dおよび3Dシーン（ダム崩壊、液滴の飛散、渦パッチ、環状崩壊、スロッシングタンク）において、Uniform FLIP、従来のスカラ動的ワーピング、およびリアクティブ（非予測的）な異方的バリアントと比較して評価された。

• グリッド適応強度: 提案手法は、スカラ・ワーピングと同等またはわずかに高い正規化グリッド変位（平均1.52 対 1.45）を達成し、品質のために適応性を犠牲にしていないことを確認した。
• 局所グリッド品質:
  • ヤコビアン歪み: 平均ヤコビアン行列式の偏差 ($E_J$) は約45.1%減少した (0.204 から 0.112)。
  • 条件数: 平均条件数は約25.1%減少した (1.99 から 1.49)。
  • 高歪みセル: 歪み閾値を超えるセルの割合は約70.9%減少した (9.36% から 2.72%)。
  • 体積保存: 平均セル体積誤差は約47.4%減少した。
• シミュレーションの安定性: 投影後の発散は約30.2%減少し、圧力計算における数値的安定性の向上を示した。
• パフォーマンス: フル手法（トラッキングあり）は、スカラ・ワーピングと比較して約9.3%の時間オーバーヘッドが発生したが (1フレームあたり13.63 ms 対 12.47 ms)、大幅に優れた品質を提供した。トラッキング機構により、非トラッキングの予測版と比較してワーピング段階のコストが約31.9%削減された。
• 3Dの堅牢性: 本手法は、双方向の衝突や局所的な飛散を含む複雑な3Dシナリオにおいても、安定したヤコビアン指標と低い歪み比を維持した。

5. 意義と主張

本論文は、予測的な方向性適応が、構造化ワープト・グリッドFLIPSにおいて、明示的な局所幾何品質制御と正常に組み合わせることができると主張している。主な意義は、「流れの特徴を捉えるために十分に動かすこと」と「セルの品質を維持すること」の間のトレードオフによって、グリッド・ワーピングが制限される必要はないことを示した点にある。適応の強さと局所的な歪みの制約を分離することで、本手法は以下を実現している：

• 方向的効率性: 解像度が単なる空間的近接性ではなく、流れの引き伸ばしや圧縮に沿って配置される。
• 安定性: 明示的なヤコアン制約により、攻撃的な異方的ワーピングが通常陥るセルの崩壊やせん断を防ぐ。
• 時間的コヒーレンス: トラッキング戦略により、ワープト・グリッドの品質を低下させることなく、効率的な計算が可能となる。

著者は、限界として、急激な衝突時における予測誤差の可能性や、極端なシナリオにおいてチューニングを要する固定された品質閾値への依存を挙げている。しかし、結果は、このアプローチが解像度と数値的安定性の両方が極めて重要となる高忠実度液体シミュレーションにおいて、堅牢なフレームワークを提供することを支持している。