Liquid agitation in a rotating and oscillating prism

OpenFOAM 4.x

Case directory

$FOAM_TUTORIALS/multiphase/interDyMFoam/ras/testTubeMixer

Summary

Model geometry (Prism and liquid surface) Model geometry (Prism and liquid surface)

We calculate the water flow in a shaking prism for one second using the mesh movement function. The motion of the prism is defined by superimposing two types of motion.

The properties of the water and air are specified in the file "constant/transportProperties".

Mesh movement is specified in the file "constant/dynamicMeshDict" as follows.

dynamicFvMesh   solidBodyMotionFvMesh;
solidBodyMotionFvMeshCoeffs
{
    solidBodyMotionFunction multiMotion;
    multiMotionCoeffs
    {
        // Table rotating in z axis
        rotatingTable
        {
            solidBodyMotionFunction rotatingMotion;
            rotatingMotionCoeffs
            {
                origin          (0 0.1 0);
                axis            (0 0 1);
                omega           6.2832; // rad/s
            }
        }
        // Tube rocking on rotating table
        rotatingBox
        {
            solidBodyMotionFunction oscillatingRotatingMotion;
            oscillatingRotatingMotionCoeffs
            {
                origin          (0 0 0);
                omega           40;         // rad/s
                amplitude       (45 0 0);   // 45 degrees max tilt
            }
        }
    }
}

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

Mesh Meshes

The calculation result is as follows.

Surface boundary of water
Surface boundary of water at initial time Surface boundary of water at initial time
Surface boundary of water at final time Surface boundary of water at final time

Commands

cp -r $FOAM_TUTORIALS/multiphase/interDyMFoam/ras/testTubeMixer testTubeMixer
cd testTubeMixer

blockMesh
cp 0/alpha.water.orig 0/alpha.water
setFields
interDyMFoam

# Display in 3D
paraFoam

Calculation time

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

Reference