固液二相流のシミュレーション（再投稿）

[yuki yuki]

お世話になっております．

大学院の研究において，今年度からopenfoamを使っている学生です．

Multiphase/icoReactingMultiphaseInterFoam/solidMelting2D

のチュートリアルを改変して，「液相流れに球（今回は二次元なので円柱）を固定して，球の体積現象，融解」を目指し，シミュレーションを実装しています．

問題点として

①setfieldsdictを用いて，一部分を固体体積率１にしていますが，固体部分が流れによって移動してしまう．

snappyhex,decomposepardic等を用いたメッシュ作成も考えましたが，円柱内部の体積変化が観察できないと思いました．

②計算が途中で破談してしまう．

添付ファイルとして，以下に添付します．

①可視化結果（中央部に移る円環はcylinder境界条件のものです．無視してください．）②計算破断スクショ③setfieldsdict④controldict

長くなってしまいましたが，一番に知りたいのはsetfiledsの内部場への固定です．

よろしくお願いいたします．