CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

rotatingWallVelocity bc

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 9, 2010, 07:42
Default rotatingWallVelocity bc
  #1
New Member
 
Join Date: Jul 2009
Posts: 11
Rep Power: 16
Lodda is on a distinguished road
Hi Foamers,

i try to calculate a flow between two cylinders. The inner cylinder is rotating so i used the rotatingWallVelocity bc for the inner wall.

My 0/U file:

dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
movingWall
{
type rotatingWallVelocity;
origin (0 0 0);
axis (0 0 1);
omega 25;
}

fixedWall
{
type fixedValue;
value uniform (0 0 0);
}

frontandback
{
type empty;
}
}

Paraview shows that the rotating wall velocity is exact twice bigger then expected.
The diameter of the inner cylinder is 0.035 m OpenFoam wall velocity 0.875 m/s.
What is the unit of omega in the U ?

Thanks an best regards

Lodda
Lodda is offline   Reply With Quote

Old   September 9, 2010, 08:30
Default
  #2
Senior Member
 
Anonymous
Join Date: Mar 2009
Posts: 110
Rep Power: 17
madad2005 is on a distinguished road
velocity (m/s) = omega (1/s) * radius (m)

Are you using the cylinder diameter or the radius to compute the rotational velocity? That could be why your velocity is twice as big as you were expecting.
madad2005 is offline   Reply With Quote

Old   September 9, 2010, 08:54
Default
  #3
New Member
 
Join Date: Jul 2009
Posts: 11
Rep Power: 16
Lodda is on a distinguished road
The OpenFoam result for the rotating wall velocity is 0.875 m/s. I expected according to v = omega * r = 25 * 0.0175 = 0.4375 m/s. I dont know what is wrong.

boundaryField
{
movingWall
{
type rotatingWallVelocity;
origin (0 0 0);
axis (0 0 1);
omega 25;
value nonuniform List<vector>
200
(
(0.0137422 -0.874784 0)
(0.0412131 -0.873921 0)
(0.0686432 -0.872195 1.86596e-27)
(0.0960057 -0.869609 -0)
(0.123273 -0.866164 0)
(0.150419 -0.861864 0)

This are the components of the calculated velocity. The magnitude at the rotating wall is always 0.875 m/s.



Best regards

Lodda
Lodda is offline   Reply With Quote

Old   September 9, 2010, 09:23
Default
  #4
Senior Member
 
Anonymous
Join Date: Mar 2009
Posts: 110
Rep Power: 17
madad2005 is on a distinguished road
Ok, one thing to check (which I am sure you have done) is that the diameter you are using is correct. Second is to make sure your origin is in the right place. What am I looking at in the image? Do you have a cylinder within a cylinder? Apart from this, I'm a bit stumped tbh. I have it working without an issue on a more complex problem that this one.
madad2005 is offline   Reply With Quote

Old   September 9, 2010, 09:33
Default
  #5
New Member
 
Join Date: Jul 2009
Posts: 11
Rep Power: 16
Lodda is on a distinguished road
This are two cylinders and the inner cylinder is rotating with omega = 25 as explained above and thats it.

Sorry, i checked my blockMeshdict the diameter is wrong!

Thanks for answering
Lodda is offline   Reply With Quote

Old   September 9, 2010, 09:47
Default
  #6
Senior Member
 
Anonymous
Join Date: Mar 2009
Posts: 110
Rep Power: 17
madad2005 is on a distinguished road
Heh, no worries. It's a surprisingly easy mistake to make!
madad2005 is offline   Reply With Quote

Old   September 12, 2010, 03:41
Default
  #7
Member
 
David
Join Date: Aug 2010
Location: North Ossetia, Vladikavkaz
Posts: 30
Rep Power: 15
davidmd is on a distinguished road
Send a message via ICQ to davidmd
Lodda, I'm trying to simulate cylinder with rotating bottom, filled with a viscous laminar incompressible fluid. Your case is very similar to mine. Maybe you could send me your files, because I'm new to OpenFOAM, and can not figure out how I do it.
Sorry for my English
davidmd is offline   Reply With Quote

Old   February 14, 2012, 08:53
Default train moving in tunnel
  #8
Member
 
Join Date: Jan 2012
Posts: 58
Rep Power: 14
sheikh nasir is on a distinguished road
hello ,
can any body help me as i am getting floating error: invalid number. my email sheikhnasir39@gmail.com
with regards
sheikh nasir is offline   Reply With Quote

Old   June 18, 2018, 12:31
Default
  #9
Member
 
Viraj Belekar
Join Date: Jun 2016
Posts: 68
Rep Power: 9
viraj20feb is on a distinguished road
Quote:
Originally Posted by madad2005 View Post
Ok, one thing to check (which I am sure you have done) is that the diameter you are using is correct. Second is to make sure your origin is in the right place. What am I looking at in the image? Do you have a cylinder within a cylinder? Apart from this, I'm a bit stumped tbh. I have it working without an issue on a more complex problem that this one.
Greetings!

I am wanted ask to what do you mean by 'make sure your origin is in the right place'? I also have an annulus-type geometry in which the inner cylinder is rotating. Also, if we look at the source code:
const vectorField Up
(
(-om)*((patch().Cf() - origin_) ^ (axis_/mag(axis_)))
);

Shouldn't the distance be the perpendicular distance from the axis rather than the distance between the face center and the origin?

Thanks
Viraj
viraj20feb is offline   Reply With Quote

Reply

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
How to use rotatingWallVelocity boundary condition terrybarnaby OpenFOAM 26 February 7, 2018 19:53
timeVaryingUniformFixedValue andrea.pasquali OpenFOAM 18 October 22, 2012 08:30
rotatingWallVelocity BC milleniumrider OpenFOAM Running, Solving & CFD 1 September 16, 2010 12:00
Rotating cylinder ole OpenFOAM 6 February 10, 2010 05:45


All times are GMT -4. The time now is 04:32.