CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM

Cylindrical coordinate system - tangential velocity

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   December 6, 2010, 09:18
Default Cylindrical coordinate system - tangential velocity
Member's Avatar
Join Date: Apr 2010
Posts: 48
Rep Power: 15 is on a distinguished road

I am running a case in which I have to set up the velocity in cylindrical coordinate system. My case has the axial velocity of 10 m/s and tangential velocity (which is homogeneous in all inlet) of 7 m/s (radial component is zero). With the script for the inlet velocity in the BC, I have the following possibility:

type cylindricalInletVelocity;
axialVelocity 10;
centre (0 0 0); //by default (0 0 0)
axis (0 0 1);
rpm 22.29;
radialVelocity 0;
value uniform (0 0 10);

but this script gives me only an outside tangential velocity (given by the rpm and the radius of the inlet tube) which is not the homogeneous tangential component I need for the simulation (in this case it is proportional to the radius, but this is not my case).

If I set


is not recognized!

does anyone know how to solve this problem?

many thanks is offline   Reply With Quote

Old   December 6, 2010, 20:39
Stefano Wahono
Join Date: Aug 2010
Location: Melbourne, Australia
Posts: 42
Rep Power: 15
swahono is on a distinguished road
I think you better of using groovyBC.

Download instruction from:

The 0/U dict entry should look something like this:

inlet    {        type            groovyBC;        variables       "rpm=8000.0;Un=55.8;c=sum(pos()*mag(Sf()))/sum(mag(Sf()));n=sum(normal())/mag(sum(normal()));p=pos()-c;r=mag(p)+1.0e-10;R=max(r);xt=vector(n.y,-n.x,0);xT=xt/mag(xt);yt=vector(-n.x*n.z,-n.y*n.z,n.x*n.x+n.y*n.y);yT=yt/mag(yt);";        valueExpression "-Un*normal() + (rpm*pi/30)*((p & yT)*xT - (p & xT)*yT)";
c is the centre of the patch, n is the averaged patch normal. You can easily extend it to have a variation on the patch normal component, for instance a parabolic inlet would only require the addition of *(1-pow(r/R,2)) to Un*normal().

swahono is offline   Reply With Quote

Old   December 7, 2010, 05:07
Default Re:
Member's Avatar
Join Date: Apr 2010
Posts: 48
Rep Power: 15 is on a distinguished road
Thank you very much for your answer. I didn't know anything about groovy BC, which is actually always a usefull tool since I am new on OF.

My problem, unfortunately, is not regarding a swirling with rpm.

Imagine that I have a tube in which I want to set up a velocity in cylindrical coordinate with two components: an axial one and a tangential one. If I use rpm, the tangential velocity is proportional to radius, but in my case all the intel patch needs the same tangential component magnitude (see the inlet zone on the attachment). This is due to simulate the swirling motion.

Have a look to the picture of the inlet: setting up the axial and tangential component, inside the tube I will have a swirl.

Many thanks for cooperation

Attached Images
File Type: png sim.png (12.7 KB, 89 views) is offline   Reply With Quote


Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
cylindrical coordinate -Gambit noussa ANSYS Meshing & Geometry 2 June 12, 2010 09:07
velocity of massless particles and velocity of flow field in same coordinate payam_IUST FLUENT 0 October 18, 2009 23:24
Star-CD Cylindrical coordinate hamoudi STAR-CD 4 April 29, 2009 11:11
tangential velocity on a velocity inlet in 2D bretti FLUENT 0 July 3, 2007 08:42
Star-cd:Refine mesh in cylindrical coordinate system Moussa Diarra Main CFD Forum 0 May 11, 2000 08:17

All times are GMT -4. The time now is 07:35.