Oscillating channel

Update: June 1, 2017
OpenFOAM 4.x

Case directory

$FOAM_TUTORIALS/incompressible/pimpleDyMFoam/oscillatingInletACMI2D

Summary

We woll calculate the flow in a channel with inlet and outlet, where an inlet part of the channel oscillates in the vertical direction of the main flow.

Model Model

The fluid flows in from the region "inlet" at (1, 0, 0) m/s and out from the region "outlet". The inlet channel oscillates in the Y-axis direction with an amplitude of 0.5 m and a frequency of 3.14 rad/s (0.5 rps). The Oscillating region, direction, amplitude and frequency are specified in the file constant/dynamicMeshDict as follows.

dynamicFvMesh   solidBodyMotionFvMesh;
motionSolverLibs ( "libfvMotionSolvers.so" );
solidBodyMotionFvMeshCoeffs
{
    cellZone        inletChannel;

    solidBodyMotionFunction  oscillatingLinearMotion;

    oscillatingLinearMotionCoeffs
    {
        amplitude       (0 0.5 0);
        omega           3.14;
    }
}

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

Mesh at time 0 Mesh at time 0

The calculation result is as follows.

Mesh movement
Flow velocity (U)
Flow velocity at final time (U) Flow velocity at final time (U)
Pressure at final time (p) Pressure at final time (p)

Commands

cp -r $FOAM_TUTORIALS/incompressible/pimpleDyMFoam/oscillatingInletACMI2D oscillatingInletACMI2D
cd oscillatingInletACMI2D

blockMesh
topoSet -constant
createBaffles -overwrite
cp -rf 0.orig 0

decomposePar
mpirun -np 4 pimpleDyMFoam -parallel
reconstructPar

paraFoam

Calculation time

  • 1 minute 54.28 seconds *4 parallel, Inter(R) Core(TM) i7-2600 CPU @ 3.40GHz 3.40GHz