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

rotating mesh with cyclic boundaries

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 16, 2013, 12:53
Question rotating mesh with cyclic boundaries
  #1
New Member
 
Christian
Join Date: Nov 2013
Posts: 7
Rep Power: 12
ChristianR1988 is on a distinguished road
Hello everybody,

first of all i have to admit, that this is my first bigger project with OpenFoam. I wish to simulate a radial compressor stage using rhoPimpleDyMFoam. In order to get a quick solution i would like to use cyclic patches and rotate the whole mesh with a constant rotational speed between an inlet- and outlet-mesh. So, my question is simple : Is it possible to use cyclic patches instead of meshing the whole stage in this case?
For example: would it be possible to change the tutorial " annularThermalMixer " that way, so that there is just one single blade and a piece of the outer mesh?
I would be very happy, if someone could help me with this question. I already read some stuff about GGI and dynamic Meshes, but until now i didnt found a solution

Thanks for support.
Attached Files
File Type: pdf rotating_system.pdf (93.9 KB, 54 views)
ChristianR1988 is offline   Reply With Quote

Old   November 19, 2013, 06:00
Default fixed Wall
  #2
New Member
 
Christian
Join Date: Nov 2013
Posts: 7
Rep Power: 12
ChristianR1988 is on a distinguished road
OK, I solved the problem by editing an "dynamicMeshDict" -file.Now my mesh rotates, as it should . But now I have a new problem: the upper wall have to stand still. I think that would be no problem, if there is a big gap between the shovel and the channel-wall (like in the annularThermalMixer-case), but in this case the gap is around 0,5 mm thin. So there are just a few number of cells in it. I dont exactly now how the tutorial "annularThermalMixer" work, but in my oppinion it would be the simplest way to do it similar. But I have a little doubt of it, because its so near to the wall and may disturb the solution of the boundary-layer. Does anybody now a more simple way???Is it possible for example to write a boundarycondition in the manner U_w=-Omega*R(z), where the velocity at the wall is negative to the direction of the rotation but not constant in place?

Thanks for help and best regards!
ChristianR1988 is offline   Reply With Quote

Old   November 20, 2013, 03:12
Default
  #3
New Member
 
Thomas Reviol
Join Date: Jul 2011
Location: Germany, Kaiserslautern
Posts: 27
Rep Power: 14
TReviol is on a distinguished road
Hi Chrisitan,

you could choose the BC "rotatingWallVelocity".

Regards

Thomas
TReviol is offline   Reply With Quote

Old   November 26, 2013, 07:51
Default thx
  #4
New Member
 
Christian
Join Date: Nov 2013
Posts: 7
Rep Power: 12
ChristianR1988 is on a distinguished road
Hi Thomas!
Thank you very much for your reply!
Wouldnt be me, if there would be no more question at this point :
Is rotatingWallVelocity-BC at the Rotor on static mesh the same like a moving mesh with uniformValue (0 0 0) -BC on the Rotor?

Best Regards,

Christian!
ChristianR1988 is offline   Reply With Quote

Old   November 26, 2013, 10:00
Default
  #5
New Member
 
Thomas Reviol
Join Date: Jul 2011
Location: Germany, Kaiserslautern
Posts: 27
Rep Power: 14
TReviol is on a distinguished road
Hi Christian,

for simulating moving mesh, I define the moving regions within the rotating mesh with the BC "movingWallVelocity" (value is set to uniform (0 0 0)). For rotating faces within the non-rotating mesh, I set the BC "rotatingWallVelocity".

Regards

Thomas
TReviol is offline   Reply With Quote

Old   November 27, 2013, 14:32
Default
  #6
New Member
 
Christian
Join Date: Nov 2013
Posts: 7
Rep Power: 12
ChristianR1988 is on a distinguished road
Haudi Thomas!
Thanks for your answer!
I have tried for innerwall AND blade (rotor) the following (non-moving mesh) for U:

type rotatingWallVelocity;
origin (0 0 0);
omega 100;
value uniform (0 0 0);

and the outerwall

type fixedValue;
value uniform(0 0 0);

The strange thing is, that, when I plot the solution, the velocity at the innerwall is , as expected not 0, but on the blade it is 0... how is that possible? there must be something wrong...
Im also running the dynamic-mesh case, which seems to do a better job.
I'll post my results as I get them.

Best Regards,

Christian.
ChristianR1988 is offline   Reply With Quote

Reply


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
[snappyHexMesh] Add Mesh Layers doesnt work on the whole surface Kryo OpenFOAM Meshing & Mesh Conversion 13 February 17, 2022 07:34
how to set periodic boundary conditions Ganesh FLUENT 15 November 18, 2020 06:09
3D Hybrid Mesh Errors DarrenC ANSYS Meshing & Geometry 11 August 5, 2013 06:42
[ICEM] Creating 3D mesh in ICEM with defined boundaries martyn88 ANSYS Meshing & Geometry 4 September 9, 2012 15:42
Moving Mesh & Not Rotating Mesh AB Siemens 1 October 25, 2004 03:10


All times are GMT -4. The time now is 21:36.