Mixing air, water, oil and mercury in a cylindrical mixer

Update: July 10, 2017
OpenFOAM 4.x

Case directory

$FOAM_TUTORIALS/multiphase/multiphaseInterFoam/laminar/mixerVessel2D

Summary

We calculate the mixing of air, water, oil, and mercury in a mixer with MRF function. The region "rotor" rotates at 6.2831853 rad/s with the Z-axis as the rotation axis, and the calculation is performed as a two-dimensional problem with one mesh in the Z direction.

Model geometry Model geometry

The rotation axis and the rotation speed are defined in the file "constant/MRFProperties" as follows.

MRF1
{
    cellZone    rotor;
    active      yes;

    nonRotatingPatches ();

    origin    (0 0 0);
    axis      (0 0 1);
    omega     constant 6.2831853;
}

The rotation area is defined as the cell zone "rotor" as shown in the figure below (white part).

Cell zone rotor Cell zone "rotor"

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

Mesh Meshes

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 final time (alphas) Fluids at final time (alphas)

Commands

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

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

multiphaseInterFoam

paraFoam

Calculation time

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