水、油、水銀によるダム・ブレイク(高分解メッシュ)

OpenFOAM 4.x

ケース

$FOAM_TUTORIALS/multiphase/multiphaseInterFoam/laminar/damBreak4phaseFine

概要

解析モデル 解析モデル

水(water)、油(oil)、水銀(mercury)の3種類の液体によるダム・ブレイク問題を6秒分、解きます。Z 方向は1メッシュとして2次元問題として計算を行います。

各相の物性、相互作用はファイル constant/transportProperties で以下の様に指定しています。ここで nu は動粘性係数(m2/s)、rho は密度(kg/m3)、sigmas は各相の間の表面張力です。

phases
(
     water
     {
         transportModel Newtonian;
         nu [0 2 -1 0 0 0 0] 1e-06;
         rho [1 -3 0 0 0 0 0] 1000;
     }

     oil
     {
         transportModel Newtonian;
         nu [0 2 -1 0 0 0 0] 1e-06;
         rho [1 -3 0 0 0 0 0] 500;
     }

     mercury
     {
         transportModel Newtonian;
         nu [0 2 -1 0 0 0 0] 1.125e-07;
         rho [1 -3 0 0 0 0 0] 13529;
     }

     air
     {
         transportModel Newtonian;
         nu [0 2 -1 0 0 0 0] 1.48e-05;
         rho [1 -3 0 0 0 0 0] 1;
     }
);

sigmas
(
    (air water) 0.07
    (air oil) 0.07
    (air mercury) 0.07
    (water oil) 0.07
    (water mercury) 0.07
    (oil mercury) 0.07
);

メッシュは以下の通りで、メッシュ数は34200です。

メッシュ メッシュ
メッシュ(正面) メッシュ(正面)

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

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

  • 水:0.0
  • 油:1.0
  • 水銀:2.0
  • 空気:3.0
各相(alphas)
初期時刻での各相(alphas) 初期時刻での各相(alphas)
時刻0.75秒での各相(alphas) 時刻0.75秒での各相(alphas)
最終時刻での各相(alphas) 最終時刻での各相(alphas)

最終的に各相が密度ごとに分かれて層を作ることが確認できます。

実行コマンド

cp -r $FOAM_TUTORIALS/multiphase/multiphaseInterFoam/laminar/damBreak4phaseFine damBreak4phaseFine
cd damBreak4phaseFine

cp -r 0.orig 0

blockMesh

setFields
decomposePar
mpirun -np 4 multiphaseInterFoam -parallel
reconstructPar

paraFoam

計算時間

48分52.74秒 ※4並列、Inter(R) Core(TM) i7-2600 CPU @ 3.40GHz 3.40GHz