Skip to content

Wrong radial velocity near the rotating wall (rotatingWallVelocity)

I am using the type rotatingWallVelocity to simulate taylor-couette flow. Strangely, the radial velocity of the rotating wall is equal to 0, while the cells close to the rotating wall, especially the first layer cells, show abnormal radial velocities. In order to observe this phenomenon more clearly, I simulated the case where the inner and outer walls have the same angular velocity. The inlet and outlet are set to be cyclic. The front patch and back patch are set to be empty.

The magnitude of velocity

image

The radial velocity(ur)

image

The radial velocity(ur) is obtained by codes following:

dimensionedVector zNormal("zNormal", dimless, vector(0.0, 0.0, 1.0));volVectorField C = mesh.C(); volVectorField xy = (C - zNormal * (C & zNormal)); ur = U & xy / mag(xy);

The run file is attached below rotatingWallVelocity_test.zip

OpenFOAM version : v2012 Operating system : centos

Edited by Jinghong Su