東京工科大学 メディア学部 ブログ

FLIP は 2005 年に Zhu と Bridson らによって発表された手法 [1] で，流体の運動を記述するナビエ・ストークス方程式に基づき，移流項と圧力項に分離して順々に計算します．

移流項の計算には粒子を使ってラグランジュ的に，圧力項の計算には格子を使ってオイラー的に解きます．したがって，移流項は流速値を持った粒子を流れに沿って移動させるだけで計算できる上，粒子の分布をそのまま液体表面の追跡として利用できるので，レベルセット法などの界面追跡法を別途導入する必要がありません．
また，圧力項は格子法で離散化するので，完全な粒子法ですべてを計算する手法と比べて未知数の行列のバンド幅が小さく，高速に圧力場を求めることができます．

いわば，「粒子（ラグランジュ）」と「格子（オイラー）」のいいとこ取りな手法ですね（笑）．


FLIP シミュレーションの流れを下図に示します．
FLIP では，最初に圧力項を計算します．圧力項は，① 粒子の流速場を格子に転写し，② いったん格子上で圧力を計算します．粒子での圧力勾配は ③ 格子で計算された圧力の勾配を粒子位置で補間して求めます．移流項は圧力項のあとに ④ 粒子の位置で格子の流速を補間し，その流速に従って粒子位置を移動して計算します [2]．

図．FLIP シミュレーションの概要

FLIP シミュレーションを実装した結果を，以下に示します（菊池研究室での実装結果です）．

映像：FLIP シミュレーションの例（菊池研での実装）

FLIP そのものは，1980 年代の CFD（Computational Fluid Dynamics 数値流体力学）分野において，おもに密度が変化するようなダイナミックな流体の動きを算出する手法として活用されていました．それを 2000 年代中盤に CG 分野でも応用するようになり，現在の流体シミュレーションの代表的な手法のひとつとなっています（ただし，FLIP にも欠点はあることを断っておきます）．

このように，異分野で活用されている手法を新しい分野でも改良・拡張しながら利用していくということは，意外と珍しいことではありません．

これからの研究者には，広く様々な領域に目を向け，日々勉強を続けながら精進していくことが求められます．