Physically Inspired Particle Framework for Recalescence-Induced Self-Jumping Ice Dynamics

技術要約：再加温に起因するセルフジャンピング氷力学のための、物理学に触発された粒子フレームワーク

問題提起
コンピュータグラフィックスは、氷（つらら、洞窟生成物、表面の不透明度など）の幾何学的成長や静的な外観を広範にモデル化してきたが、既存の手法は、凍結の初期段階で発生する非平衡な過渡現象を捉えることに大きく失敗している。具体的には、再加温に起因するセルフジャンピング（recalescence-induced self-jumping）——すなわち、急速な潜熱放出と非対称な蒸発により、過冷却された液滴が自発的に離脱し、跳ね上がる現象——は、軽量なシミュレーションフレームワークでは未だ対処されていない。この挙動を完全に解明するには、熱伝達、相変化、蒸気流、および剛体力学を結合させた、高コストな多相熱流体ソルバーが必要となるが、これらはリアルタイムのグラフィックスアプリケーションにおいては計算量的に極めて困難である。本研究では、完全な流体方程式を解くことなく、本質的なインパルス生成メカニズムを捉えるための、低次かつ制御可能なモデルが欠如している。

手法
著者は、非対称な蒸発フラックス（再加温フロント付近）による支配的な運動量生成メカニズムを近似することで、セルフジャンピングをシミュレートする低次・粒子ベースの物理モデルを提案する。本フレームワークは、以下の4段階のパイプラインで動作する。

1. 自動核生成（Automatic Nucleation）： システムは、Schrageに基づく蒸発冷却モデル（上向きの面でより強力な方向バイアスを持つ）を用いて表面冷却をシミュレートする。表面粒子の温度が臨界過冷却閾値（$T_{nuc} = T_m - \Delta T_{nuc}$）を下回ったときに、核生成が自動的にトリガーされる。この最初の閾値通過地点が、核生成のシード（種）として定義される。
2. 再加限フロントの伝播（Recalescence-Front Propagation）： 核生成後、凍結フロントはシードから表面に沿って双方向に伝播する。モデルは、温度を融点（$T_m$）付近で安定させるために「マッシュバンド（Mushy band）」熱クランプを採用し、液体分率（$\alpha$）を徐々に減少させることで、液相から固相への相転移をシミュレートする。
3. 弧限定的な方向性インパルスの蓄積（Arc-Limited Directional Impulse Accumulation）： ジャンプを捉えるために、モデルは氷／マッシュ表面（$J_i$）と液体表面（$J_w$）の間の蒸発フラックス差を計算する。極めて重要な点として、全表面にわたって積分すると水平成分が相殺されてしまうため、本フレームワークは、インパルス蓄積をフロントの先端付近の狭い弧限定領域（または3Dにおける測地線バンド）に制限する。この選択的な積分により、非対称な運動量の不均衡が保持され、正味の線形インパルス（$J_{acc}$）および角インパルス（$\tau_{acc}$）が生成される。
4. 剛体離脱（Rigid-Body Detachment）： 蓄積されたインパルスが閾値を超えると、凍結した液滴は剛体運動へと移行する。蓄積されたインパルスはジャンプ力（$F_{jump}$）とトルク（$\tau_{jump}$）に変換され、元の形状を維持したまま、核生成シードとは逆方向に物体を推進させる。

任意の幾何形状への一般化
本フレームワークは、2Dの円形ベースモデルから、任意の2D画像マスクおよび3D三角形メッシュへと拡張される。

• 2Dマスク： 円の角度座標は、表面グラフから導出された境界弧長パラメータ化に置き換えられる。
• 3Dメッシュ： フロントの伝播は、メッシュ表面上の測地線距離を用いてパラメータ化される。インパルス蓄積領域は、フロント先端の周囲の測地線バンドとなる。これにより、モデルは複雑なトポロジー（ドラゴン、バニー、ティーポットなど）を扱いながら、同一のフロントベースのロジックを維持できる。

主な貢献

1. フロント認識型相転移モデル： 自動核生成、双方向フロント伝播、およびマッシュバンド熱安定化を結合した統一的な構造を提供し、完全な熱流体結合なしに再加温段階のシミュレーションを可能にする。
2. 弧限定インパルス積分： フロント先端付近の狭い領域内でのみ方向性インパルスを蓄積する新しいスキーム。これにより、水平方向の相殺を効果的に抑制し、実験的に観察されるジャンプの方向（核生成サイトの反対側）を再現する。
3. 形状に依存しない拡張性： 2D粒子定式化を、測地線パラメータ化を通じて任意の3D三角形メッシュへと正常に一般化し、セルフジャンピングのダイナミクスが表面のトポロジーや曲率に依存しないことを示した。
4. 計算効率： PyCUDAベースの軽量な実装により、完全なマルチフィジックスソルバーを回避しながら、リアルタイムに近い性能（2Dで約1.16 ms/フレーム、3Dで約3.08 ms/フレーム）を達成した。

実験結果
著者は以下の手法を通じてフレームワークの検証を行った：

• 2Dシミュレーション： 上部表面がより速く冷却され、最初に核生成するという非対称な凍結パターンをモデルが再現することを示した。結果として、ジャンプの方向は一貫してシードの反対側であり、その大きさは形状の対称性（例：非対称なスパイラルは対称なハートよりも強いジャンプを生む）に基づいて変化することが示された。
• サンプリングの堅牢性： 本手法は、グリッドサンプリングとポアソン円盤サンプリングの両方において安定しており、局所的な密度不規則によるインパルス量の変動は軽微であった。
• 3Dメッシュへの適用： 複雑なモデル（Dino、Bunny、Dragon、Teapotなど）に適用した結果、システムは局所的な曲率およびフロント伝播経路に沿った幾何学依存のジャンプベクトルを生成した。
• ユーザー制御： 本フレームワークは、自動核生成とユーザー定義のシード位置の両方をサポートしており、アーティストが凍結方向とジャンプ方向をインタラクティブに制御することを可能にする。

意義と主張
本論文は、再加温に関する流体力学の研究とコンピュータグラフィックス・シミュレーションの間の溝を埋めるものであることを主張している。その主要な意義は、計算可能な、低次表現としての再加温駆動型インパルス生成を提供することにある。著者は、提案するフレームワークが以下の特性を持つと断言している：

• 完全な熱流体ソルバーのコストをかけることなく、セルフジャンピングの本質的な物理（熱的不対称性 $\to$ 蒸発反動 $\to$ 離脱）を捉える。
• 視覚効果、ゲーム、および科学的可視化のための、形状に依存しないかつ制御可能なツールを提供し、妥当な相変化ダイナミクスの生成を可能にする。
• 凍結アニメーション・パイプラインを置き換えるものではなく、そこに統合可能な「再加温段階インパルスモジュール」として機能する。

著者は、本モデルが現象論的な近似であることを認め、限界についても明記している。本モデルは、温度や速度場の正確な予測を行うものではなく、経験的に調整されたパラメータ（インパルスゲイン、弧の幅など）に依存し、3Dの体積相変化、ガス相互作用、または不純物輸送を解決するものではない。しかし、視覚的に妥当なリアルタイム・グラフィックス・シミュレーションという目的において、提案されたフレームワークは、再加温によるセルフジャンピングの支配的なメカニズムを正常に再現していると主張している。