[論文レビュー] Simulating the Dynamic Behavior of Shear Thickening Fluids

Oktar Özgen, Marcelo Kallmann+1 more|arXiv (Cornell University)|Oct 28, 2015

Rheology and Fluid Dynamics Studies参考文献 12被引用数 3

ひとこと要約

本稿では、分数階微積分に基づく履歴依存剛性モデルを用いて、視覚的に正確なせん断強化流体（STF）のシミュレーションを初めて実現した。分数導関数を用いて粘弾性の記憶効果を捉えることで、衝撃抵抗、固体的反発、遅延沈下といった不連続せん断強化現象を、コンピュータグラフィックス分野でこれまで達成されていなかった形で再現した。

ABSTRACT

While significant research has been dedicated to the simulation of fluids, not much attention has been given to exploring new interesting behavior that can be generated with the different types of non-Newtonian fluids with non-constant viscosity. Going in this direction, this paper introduces a computational model for simulating the interesting phenomena observed in non-Newtonian shear thickening fluids, which are fluids where the viscosity increases with increased stress. These fluids have unique and unconventional behavior, and they often appear in real world scenarios such as when sinking in quicksand or when experimenting with popular cornstarch and water mixtures. While interesting behavior of shear thickening fluids can be easily observed in the real world, the most interesting phenomena of these fluids have not been simulated before in computer graphics. The fluid exhibits unique phase changes between solid and liquid states, great impact resistance in its solid state and strong hysteresis effects. Our proposed approach builds on existing non-Newtonian fluid models in computer graphics and introduces an efficient history-based stiffness term that is essential to produce the most interesting shear thickening phenomena. The history-based stiffness is formulated through the use of fractional derivatives, leveraging the fractional calculus ability to depict both the viscoelastic behavior and the history effects of history-dependent systems. Simulations produced by our method are compared against real experiments and the results demonstrate that the proposed model successfully captures key phenomena observed in shear thickening fluids.

研究の動機と目的

• コンピュータグラフィックス分野において、せん断強化流体（STF）の完全な動的挙動をシミュレートする計算モデルの不足に対処すること。
• 衝撃抵抗、応力下での固体的応答、流体状態への遅延回復といった、STFの非直感的物理現象を再現すること。
• 不連続せん断強化に不可欠な履歴依存挙動を捉える、計算的に効率的な手法を導入すること。
• 分数微積分が、安定なSPHフレームワーク内での非ニュートン流体の粘弾性および記憶効果を効果的にモデル化できることを示すこと。

提案手法

• 方法は、1/2階のカプート分数導関数に基づく履歴依存剛性項を滑らかさ粒子法（SPH）に拡張したものである。
• 分数導関数の定式化により、過去の速度状態を記録し、STFに内在する記憶効果を再現できる。
• 剛性項は、事前に計算された重みを用いて過去の粒子速度の重み付き和として計算され、効率性を確保する。
• この履歴依存力は、標準的なSPH力式に統合され、計算複雑性を維持しつつ物理的リアリズムを向上させた。
• 方法は標準的なSPHフレームワークに実装されており、履歴長に比例する線形メモリオーバーヘッドのみを要する。
• 実世界の実験（ボウリングボールの衝撃や振動駆動による指状形成）との比較により、妥当性が検証された。

実験結果

リサーチクエスチョン

• RQ1 計算的に効率的なSPHベースのモデルが、実際のSTFで観察される不連続せん断強化挙動をシミュレートできるか？
• RQ2 ヒステリシスや遅延的固体化といった履歴依存粘弾性効果を、物理的に妥当かつ視覚的に正確にモデル化する方法は何か？
• RQ3 従来の粘弾性モデルに比べ、分数微積分がSTFの記憶効果をどれほど効果的に表現できるか？
• RQ4 本手法が、コンピュータグラフィックス分野で先行シミュレーションなしに、衝撃反発や遅延沈下といった複雑な現象を再現できるか？
• RQ5 SPHシミュレーションに履歴依存剛性項を組み込む際の性能コストはどの程度か？

主な発見

• 提案手法は、ボウリングボールの衝撃時の反発やエネルギー散逸後の遅延沈下を含む、STFの主な視覚的現象を成功裏に再現した。
• 20,000粒子のシミュレーションで、Intel i7-2600Kで1フレームあたり約1.6秒の実行速度が得られ、履歴ステップ数500の場合、わずか40msのオーバーヘッドにとどまった。
• 分数微積分を用いて計算された履歴依存剛性項により、外部力が終了してもヒステリシスや固体的挙動のシミュレーションが可能になった。
• 視覚的に妥当な結果が得られ、特に固体的状態から流動的状態への遷移を実世界の実験とよく一致させた。
• 分数導関数の定式化により、真の1/2階微分の滑らかで安定した近似が可能となり、最小限の計算コストで現実的な記憶効果を実現した。
• 強力な視覚的忠実度にもかかわらず、振動実験における指状形成の高さが低く見積もられ、微細スケールの不安定性を捉える限界が示された。