Google OpenFOMAのユーザー様
いつもお世話になります。PIMPLE法のアルゴリズムで分かりかねる部分があったので投稿させていただきました。
SIMPLE法とPISO法に関するアルゴリズムはOpenFOAMwikiで詳しく説明されていますが、PIMPLE法のアルゴリズムではそういった情報が無く、SIMPLEとPISO法の解説とPIMPLEのプログラムコードから以下のようなアルゴリズムと解釈しています。特に不明な部分は、SIMPLEループの処理が1.境界条件の
設定前からなのか
設定後からのループなのかがよく分かりませんでした。
ご存じの方がいらっしゃいましたら以下のアルゴリズムのご確認とご指摘の程よろしくお願い致します。
1.境界条件の設定
2.離散化方程式を解き、中間速度場を計算
3.セル境界面でのフラックスを計算
4.圧力ポアソン方程式を解き、不足緩和による安定化処理?(OpenFOAMでは違うような)
5.セル境界面でのフラックスを修正
6.更新された圧力場から速度を修正
7.境界値を更新
8.3から指定された回数反復計(PISOループ)
9.2から?指定された回数反復計算(SIMPLEループ)
10.時間を増加し、1から繰り返し反復