Dam-break problem with water, oil, and mercury

OpenFOAM 4.x

Case directory

$FOAM_TUTORIALS/multiphase/compressibleMultiphaseInterFoam/laminar/damBreak4phase

Summary

Model geometry Model geometry

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

The properties of each fluid are specified in the file "constant/thermophysicalProperties.*", and the surface tensions between the fluids are specified in the file "constant/thermophysicalProperties".

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

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/compressibleMultiphaseInterFoam/laminar/damBreak4phase damBreak4phase
cd damBreak4phase

rm -rf 0
cp -r 0.orig 0

blockMesh
setFields
compressibleMultiphaseInterFoam

paraFoam

Calculation time

2 minutes 33.4 seconds *Single, Inter(R) Core(TM) i7-2600 CPU @ 3.40GHz 3.40GHz

Reference