空気、水、油、水銀の遠心分離

OpenFOAM 4.x

ケース

$FOAM_TUTORIALS/multiphase/multiphaseEulerFoam/mixerVessel2D

概要

空気、水、油、水銀が均等に混ざった流体に対して MRF 機能を用いて遠心分離を行います。領域 rotor が Z 軸を回転軸として 10.472 rad/s で回転します。Z 方向は1メッシュとして2次元問題として計算を行います。

モデル形状 モデル形状

回転軸、回転数はファイル constant/MRFProperties 内で以下の様に定義されています。

MRF1
{
    cellZone    rotor;
    active      yes;

    nonRotatingPatches ();

    origin    (0 0 0);
    axis      (0 0 1);
    omega     10.472;
}		

回転領域はセル・ゾーン「rotor」として下図のように定義されています(白色部)。

セル・ゾーン「rotor」 セル・ゾーン rotor

全体のメッシュは以下の通りで、メッシュ数は3072です。

メッシュ メッシュ

計算結果は以下の通りです。

各相と値の対応関係は以下の通りです。

  • 水:0.0
  • 油:1.0
  • 水銀:2.0
  • 空気:3.0
各相(alphas)
最終時刻での空気の体積率(alpha.air) 最終時刻での空気の体積率(alpha.air)
最終時刻での油の体積率(alpha.oil) 最終時刻での油の体積率(alpha.oil)
最終時刻での水の体積率(alpha.water) 最終時刻での水の体積率(alpha.water)
最終時刻での水銀の体積率(alpha.mercury) 最終時刻での水銀の体積率(alpha.mercury)

最終的に密度が小さい順に中心側から層を作ることが確認できます。

実行コマンド

cp -r $FOAM_TUTORIALS/multiphase/multiphaseEulerFoam/mixerVessel2D mixerVessel2D
cd mixerVessel2D

m4 < system/blockMeshDict.m4 > system/blockMeshDict
blockMesh
topoSet
setsToZones -noFlipMap

multiphaseEulerFoam

paraFoam

計算時間

22分55.96秒 ※シングル、Inter(R) Core(TM) i7-2600 CPU @ 3.40GHz 3.40GHz