有理関数を使ったCIP法. CIP法の欠点は関数を移流させるとき,若干のオーバーシュートがでることである. このオーバーシュートは移流が進んでも大きくはならないため問題がないように思えるが, 水と空気の二相流など,計算したい物理量の差が非常に大きい場合に, 小さなオーバーシュートが大きな問題となる. RCIP法 *1, *2 はこの問題を解決する. RCIP法はCIP法の3次関数の代わりに有理関数

eqa_rcip1.gif

を用いる. ここで,eqa_rcip2.gifはパラメータで0か1の値をとる. eqa_rcip3.gifのときはCIP法と同じであり,eqa_rcip4.gifのとき, この関数は単調性と凹凸性を保存する有理関数となる. eqa_rcip5.gifであり, eqa_rcip6.gifである. 上式の各係数は以下である.

eqa_rcip7.gif

ここで,

eqa_rcip8.gif

RCIP法以外にも,Tangent変換を用いることで オーバーシュートを抑える方法も提案されている *3. Tangent変換を用いた手法では物体界面がシャープに保たれるという特徴を持つ. これは数値流体計算上は非常に有益であるが, グラフィックスで用いた場合,グリッド解像度によっては 液体表面があまり滑らかにならず,保存性も保証されない.

RCIP法による密度関数の移流では,界面を表す0から1への変化部分が拡散し, 関数の傾きがなだらかになる. この滑らかな界面はCGでの利用においては美しいレンダリング結果を生むが, シミュレーションでは数値的な不安定性を引き起こす. 密度関数の体積移動によるSTAA手法 *4 を使うことで, 界面拡散を抑制することができる. また,拡散をある程度の幅で制御する方法 *5 もある.


*1 F. Xiao, T. Yabe, and T. Ito, "Constructing oscillation preventing scheme for advection equation by rational function", Comput. Phys. Commun. 93, pp.1-12, 1996.
*2 F. Xiao, T. Yabe, G. Nizam and T. Ito, "Constructing a multi-dimensional oscillation preventing scheme for the advection equation by a rational function", Comput. Phys. Commun. 94, pp.103-118, 1996.
*3 T. Yabe and F. Xiao "Description fo Complex and Sharp Interface during Shock Wave Interaction with Liquid Drop", Journal of the Physical Society of Japan, 62, pp.2537-2540, 1993.
*4 池端昭夫, 肖鋒, "保存型自由表面捕獲スキームと固液気3相流への適用", 日本機械学会2002年度年次大会, 2002
*5 M. Fujisawa and K. T. Miura, "Animation of ice melting phenomenon based on thermodynamics with thermal radiation", Proc. GRAPHITE2007, pp.249-256, 2007.

添付ファイル: fileeqa_rcip2.gif 818件 [詳細] fileeqa_rcip3.gif 836件 [詳細] fileeqa_rcip4.gif 743件 [詳細] fileeqa_rcip5.gif 818件 [詳細] fileeqa_rcip6.gif 766件 [詳細] fileeqa_rcip7.gif 803件 [詳細] fileeqa_rcip8.gif 850件 [詳細] fileeqa_rcip1.gif 809件 [詳細]

トップ   編集 凍結 差分 履歴 添付 複製 名前変更 リロード   新規 一覧 検索 最終更新   ヘルプ   最終更新のRSS
Last-modified: 2024-03-08 (金) 18:06:06