Rotor wit SRF function

Update: June 1, 2017
OpenFOAM 4.x

Case directory

$FOAM_TUTORIALS/incompressible/SRFPimpleFoam/rotor2D

Summary

A rotor rotating at 60 rpm around the Z-axis is calculated for 2 seconds. The rotation of the computational domain is modeled by adding centrifugal and Coriolis forces (SRF function).

The rotor surface (white part) is set to the no-slip condition, and the problem is solved as 2 dimensional problem with a single mesh in the Z direction.

Model geometry Model geometry

The outer region "freestream" (red part) is set to the velocity (1, 0, 0) m/s in the inertia system. This condition is specified in the file 0/Urel.

freestream
{
	type            SRFFreestreamVelocity;
	UInf            (1 0 0);
	value           uniform (0 0 0);
}

The rotation axis and the rotation speed are specified in the file constant/SRFProperties as shown below.

SRFModel        rpm;
origin          (0 0 0);
axis            (0 0 1);
rpmCoeffs
{
    rpm             60;
}

The standard k-ε model is used as the turbulence model.

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

Mesh Meshes

If you want to visualize the flow velocity in the inertia system, check "U" in the "Properties" tab on ParaView. Note that the data at time 0 does not have "U" values.

Check U Check "U"

The calculation result is as follows.

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/SRFPimpleFoam/rotor2D rotor2D
cd rotor2D

m4 < system/blockMeshDict.m4 > system/blockMeshDict
blockMesh

SRFPimpleFoam

paraFoam

Calculation time

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

Reference

OpenFOAM -回転領域を含む流体計算 (Rotating Geometry)-, Fumiya Nozaki