Cavity flow (mapping data to different shape meshes)

Update: June 1, 2017
OpenFOAM 4.x

Case directory

$FOAM_TUTORIALS/incompressible/icoFoam/cavity/cavityClipped

Summary

To run this case, you need to run the case "cavity (cavity flow)" at first.

The upper wall "movingWall" is moving in parallel at 1 m/s and the other walls "fixedWalls" are stationary, and we solve the flow without turbulence model. There is no inflow or outflow of fluid, and the calculation is performed as a 2-dimensional problem with a single mesh in the Z direction.

Model geometry Model geometry

The meshes are as follows.

Mesh Meshes

The cavity calculations result will be map to the fields by mapField utility, and used as initial conditions.

Mapped flow velocity at 0.5 seconds (U) Mapped flow velocity at 0.5 seconds (U)

The calculation result is as follows.

Flow velocity at 0.6 seconds (U) Flow velocity at 0.6 seconds (U)
Pressure at 0.6 seconds (p) Pressure at 0.6 seconds (p)

We can see that a new vortex that has generated near the cut off region.

Commands

cp -r $FOAM_TUTORIALS/incompressible/icoFoam/cavity/cavityClipped cavityClipped
cd cavityClipped

blockMesh
cp -r 0 0.5
mapFields ../cavity -sourceTime latestTime
icoFoam

paraFoam

The pre-calculated cavity-flow data directory is assumed to be in the same directory which cavityClipped exists in.

After the meshing, a directory for the data at 0.5 second is created and the cavity-flow data will be mapped to the directory with mapFields. The configuration file for mapFields is system/mapFieldsDict.

Calculation time

--- *Single, Core(TM) i7-2600 CPU @ 3.40GHz 3.40GHz