Cylindrical coordinate system - tangential velocity
hi,
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: inlet { 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 "tangentialVelocity" is not recognized! does anyone know how to solve this problem? many thanks |
I think you better of using groovyBC.
Download instruction from: http://openfoamwiki.net/index.php/Contrib_groovyBC The 0/U dict entry should look something like this: Code:
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)"; Cheers, Stefano |
Re:
1 Attachment(s)
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 regards |
All times are GMT -4. The time now is 11:27. |