Blowing air into water column

Update: July 10, 2017
OpenFOAM 4.x

Case directory

$FOAM_TUTORIALS/multiphase/multiphaseEulerFoam/bubbleColumn

Summary

Model geometry Model geometry

We calculate the mixing of air into the water column for 100 seconds. Air flows in from the region "inlet" at the bottom of the water column at (0, 0.1, 0) m/s and out from the region "outlet". We assume that the ratio of water to air is 1:1 in the region "inlet", and the calculation is performed as a two-dimensional problem with one mesh in the z-direction.

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

phases
(
    air
    {
        rho             1;
        nu              1.6e-05;
        kappa           0;
        Cp              0;

        diameterModel   constant;
        constantCoeffs
        {
            d               3e-3;
        }
    }

    water
    {
        rho             1000;
        nu              1e-06;
        kappa           0;
        Cp              0;

        diameterModel constant;
        constantCoeffs
        {
            d               1e-4;
        }
    }
);

The interaction between water and air is also specified in the file "constant/transportProperties" as follows.

drag
(
    (air water)
    {
        type blended;

        air
        {
            type SchillerNaumann;
            residualPhaseFraction 0;
            residualSlip 0;
        }

        water
        {
            type SchillerNaumann;
            residualPhaseFraction 0;
            residualSlip 0;
        }

        residualPhaseFraction 1e-3;
        residualSlip 1e-3;
    }
);

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

Mesh Meshes

The calculation result is as follows.

Volume ratio of water (alpha.water)
Volume ratio of water at final time (alpha.water) Volume ratio of water at final time (alpha.water)

Commands

cp -r $FOAM_TUTORIALS/multiphase/multiphaseEulerFoam/bubbleColumn bubbleColumn
cd bubbleColumn

blockMesh

setFields
multiphaseEulerFoam

paraFoam

Calculation time

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