回転するローターを含む自由表面解析

界面の様子

粘度が高いケース1の方が抵抗が大きくなり、ローター内に残りやすくなっています。

また、メッシュ解像度を上げることで、より細かい液滴まで再現可能となります。その分、計算負荷が大きくなりますが、OpenFOAMはライセンスフリーのため、並列計算によって計算時間の短縮が見込めます。

ローター側面のせん断応力分布

図2に流体によってローター側面に作用するせん断応力分布を示します。粘度が高いケース1の方が壁面に発生するせん断応力が大きくなっています。

 

図2　ローター側面のせん断応力分布（左：ケース1（高粘度）、右：ケース2（低粘度））

断面における流速分布

図3の断面での流速分布を図4に示します。

流動抵抗は粘度に依存するため、速度分布が変化します。粘度が低い方が抵抗が小さく、大きな流速が発生しやすくなっています。

図3　出力断面

図4　流速分布（左：ケース1（高粘度）、右：ケース2（低粘度））

• 流入、流出がある回転するローターの自由表面解析が可能
• 粘度の違いによる流れの様子の変化を確認できる
• 壁面せん断応力分布を出力可能
• 任意の断面における流速分布、流速ベクトル等を出力可能
• 細かい液滴を再現する高いメッシュ解像度の場合でも、ライセンスフリーで大並列計算が可能
• 着色機能を使うことで、均一に処理できるかを調べることも可能

以上により、粘度や形状の変更によって変化するせん断や流速等を調べることで、効率的な分散処理方法の検討にOpenFOAMを利用できます。