Dam-break problem with water, oil, and mercury (High resolution mesh)

Update: July 10, 2017
OpenFOAM 4.x

Case directory

$FOAM_TUTORIALS/multiphase/multiphaseInterFoam/laminar/damBreak4phaseFine

Summary

Model geometry Model geometry

We calculate a dam-break problem with three liquids: water, oil, and mercury, for 6 seconds. The calculation is performed as a two-dimensional problem with a single mesh in the Z direction.

The properties and interactions of each fluid are specified in the file "constant/transportProperties" as follows where nu is the kinematic viscosity (m2/s), rho is the density (kg/m3), and sigmas is the surface tension between the fluids.

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
);

The meshes are as follows, and the number of mesh is 34200.

Mesh Meshes
Meshes (front view) Meshes (front view)

The calculation result is as follows.

The ID values of each fluid are as follows.

  • Water: 0.0
  • Oil: 1.0
  • Mercury: 2.0
  • Air: 3.0
Fluids (alphas)
Fluids at initial time (alphas) Fluids at initial time (alphas)
Fluids at 0.75 sec (alphas) Fluids at 0.75 sec (alphas)
Fluids at final time (alphas) Fluids at final time (alphas)

At final time, we can see that each fluid is divided into layers by density.

Commands

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

Calculation time

48 minutes 52.74 seconds *4 parallel, Inter(R) Core(TM) i7-2600 CPU @ 3.40GHz 3.40GHz