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

Problem rotating frame of reference.

Register Blogs Community New Posts Updated Threads Search

Like Tree1Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 13, 2017, 07:41
Smile Problem rotating frame of reference.
  #1
New Member
 
Rocky
Join Date: Jan 2017
Posts: 21
Rep Power: 9
Rockytime is on a distinguished road
Hi All,

CFX solve pop out this error before it starts to run.

The angle between the specified velocity and the element surface is|
| 43.727 degrees at this face. This is considered an error because |
| it implies that the mesh is moving. The following are possible |
| reasons for the error message.

In the CFX setup, I set my fluid wall to counter rotating wall. It works in cylindrical enclosure but it fail to work in a rectangular enclosure. Anyone can tell me how should I setup? in order to make the rectangular enclosure stay stationary?

Fluid Wall


Rockytime is offline   Reply With Quote

Old   February 13, 2017, 18:31
Default
  #2
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,872
Rep Power: 144
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
The counter rotating wall condition applies a tangential velocity to the selected wall such that it opposes the domain rotation and results in a stationary wall in the absolute frame.

But, as the error message says this requires the velocity to be tangential and that means axisymmetric about the rotation axis. Your rectangular box is not axisymmetric, so you cannot model the walls of your rectangular box with counter rotating walls. This approach will work with cylindrical walls as they are axisymmetric.
Red Ember likes this.
ghorrocks is offline   Reply With Quote

Old   February 13, 2017, 20:17
Default
  #3
New Member
 
Rocky
Join Date: Jan 2017
Posts: 21
Rep Power: 9
Rockytime is on a distinguished road
Hi Glen,

Thanks

So did you have anyway to let my rectangular duct stay in stationary frame while I setting the domain as rotating motion?
Rockytime is offline   Reply With Quote

Old   February 13, 2017, 20:43
Default
  #4
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,872
Rep Power: 144
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
The rotor goes in a rotating frame of reference and the stationary rectangular duct goes in a stationary frame of reference. So there will be 2 domains, and they are connected by a GGI with some form of frame change model.
ghorrocks is offline   Reply With Quote

Old   February 13, 2017, 21:18
Default
  #5
New Member
 
Rocky
Join Date: Jan 2017
Posts: 21
Rep Power: 9
Rockytime is on a distinguished road
So for the frame change of model, can I just select general connection with frozen rotor?
Rockytime is offline   Reply With Quote

Old   February 13, 2017, 22:18
Default
  #6
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,872
Rep Power: 144
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
If frozen rotor is an applicable frame change model for you then yes. From what I can see of your model I suspect it is appropriate.
ghorrocks is offline   Reply With Quote

Old   February 13, 2017, 22:29
Default
  #7
New Member
 
Rocky
Join Date: Jan 2017
Posts: 21
Rep Power: 9
Rockytime is on a distinguished road
Another question, I have 1 enclosure and 1 import turbine geometry. Do I need to use boolean to subtract them into 1 body? No right?
Rockytime is offline   Reply With Quote

Old   February 14, 2017, 01:05
Default
  #8
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,872
Rep Power: 144
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
Probably - but I do not know what your enclosure is. I can guess what your turbine is.

This question is best asked on the geometry and meshing forum. But before you do, make sure you do the tutorial examples on turbomachinery for the meshing software and CFX so you understand the basics. You can download these from the ANSYS customer website.
ghorrocks is offline   Reply With Quote

Old   February 14, 2017, 01:18
Default
  #9
New Member
 
Rocky
Join Date: Jan 2017
Posts: 21
Rep Power: 9
Rockytime is on a distinguished road
Ok. Thanks for your advice.

I Have set my turbine in rotary domain with a rotational velocity and setup the wall boundary as well. I'm just wondering do I need to turn on mesh deformation and specific mesh motion at my wall turbine boundary? Is it ok if i specify my mesh motion with surface of revolution? Sorry for asking too much question. I try to understand mesh motion through the cfx tutorial but I still not quite sure how it function.

Last edited by Rockytime; February 14, 2017 at 03:36.
Rockytime is offline   Reply With Quote

Old   February 14, 2017, 17:42
Default
  #10
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,872
Rep Power: 144
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
If the device only rotates then you do not need mesh motion. It can be done as rotating frame of reference.

You only need mesh motion if it does some other motion (eg translates and rotates) or it deforms.
ghorrocks is offline   Reply With Quote

Old   February 15, 2017, 00:52
Thumbs up
  #11
New Member
 
Rocky
Join Date: Jan 2017
Posts: 21
Rep Power: 9
Rockytime is on a distinguished road
Ok. I run my simulation in transient analysis, but it seem the result is not in my expectation. Will it be different if i run my result in steady state analysis?

The reason I run my result in transient is I can vary my inlet velocity with CEL function. So I can see how much my torque increase when the flow increase in one graph.
Rockytime is offline   Reply With Quote

Old   February 15, 2017, 06:02
Default
  #12
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,872
Rep Power: 144
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
Quote:
So I can see how much my torque increase when the flow increase in one graph.
A common misconception. Do not do this - how do you know the transient effects have faded out? Also this approach is much slower than doing a series of steady state simulations.

As for your first sentence: I have no idea what results you got, what results you expect and how you have set your model up. I cannot answer that question with more information.
ghorrocks is offline   Reply With Quote

Old   February 15, 2017, 10:53
Default
  #13
New Member
 
Rocky
Join Date: Jan 2017
Posts: 21
Rep Power: 9
Rockytime is on a distinguished road
Ok. My simulation consists 2 domain. 1 is the Rectangular Duct ( Fluid Body) and another 1 is rotor.

I setup both domain as fluid domain with air @ 25 degree. The Rectangular Duct domain motion is set as stationary domain and the rotor is set as rotating domain and specific with a rotating speed (-4500 rpm). The fluid model is set as Shear Stress Transport (SST) with the advanced turbulence model curvature correction and production limiter.

Rectangular Duct Domain, contain 3 boundary, Inlet, Outlet and Wall Boundary. The face that located at +y axis is inlet, -y axis is outlet, the remaining faces is set as Wall boundary with stationary domain. No slip wall is applied on the wall boundary.

Rotor Domain, I named all the selection ( Shaft, blade, hub) and create them together as Wall boundary with (frame type= Rotating). No slip is applied on the wall boundary.

The result I expect is the turbine should be stall at some point when the flow increase jsut like the what picture below shown.

Now the problem im facing is I keep getting the same torque value even with I set different value of flow for inlet flow velocity.





Bellow is the Cel Code for your better understanding.
Attached Files
File Type: docx CEL Code.docx (15.5 KB, 10 views)
Rockytime is offline   Reply With Quote

Old   February 15, 2017, 11:43
Default
  #14
Senior Member
 
Join Date: Jun 2009
Posts: 174
Rep Power: 17
turbo is on a distinguished road
Your geometry setup is wrong. CFD is only for a negative volume from CAD where fluid is filling. If sticking to your original intent, you should have a single domain. But such an approach is wrong, too, in the square duct enclosure case.
turbo is offline   Reply With Quote

Old   February 15, 2017, 11:55
Default
  #15
New Member
 
Rocky
Join Date: Jan 2017
Posts: 21
Rep Power: 9
Rockytime is on a distinguished road
Quote:
Originally Posted by turbo View Post
Your geometry setup is wrong. CFD is only for a negative volume from CAD where fluid is filling. If sticking to your original intent, you should have a single domain. But such an approach is wrong, too, in the square duct enclosure case.
Hi turbo, I cant get what you mean. Can you further explain what wrong with my geomtery?
Rockytime is offline   Reply With Quote

Old   February 15, 2017, 15:18
Default
  #16
Senior Member
 
Join Date: Jun 2009
Posts: 174
Rep Power: 17
turbo is on a distinguished road
1) Make a negative volume solid.
2) Make two solids : one is a cylindrical solid immediately surrounding the rotor, the other is the rest volume with the square duct boundary.
3) One is for rotating frame, the other for stationary.
4) Stage interface.
turbo is offline   Reply With Quote

Old   February 16, 2017, 03:17
Default
  #17
New Member
 
Rocky
Join Date: Jan 2017
Posts: 21
Rep Power: 9
Rockytime is on a distinguished road
Quote:
Originally Posted by turbo View Post
1) Make a negative volume solid.
2) Make two solids : one is a cylindrical solid immediately surrounding the rotor, the other is the rest volume with the square duct boundary.
3) One is for rotating frame, the other for stationary.
4) Stage interface.
Hi turbo, do you have any example what is negative volume solid?

Is it like the step below?

1)Create---> Rectangular Enclosure
2) Create --> Boolean ---> Operation: subtract, target body: select recetangular enclosure, tool body: select Rotor. Preserve tool bodies? Yes
3) Create ---> Cylindrical Enclosure around rotor.

Am i right?

Last edited by Rockytime; February 16, 2017 at 04:30.
Rockytime is offline   Reply With Quote

Old   February 16, 2017, 07:33
Default
  #18
Senior Member
 
Join Date: Jun 2009
Posts: 174
Rep Power: 17
turbo is on a distinguished road
Yes, the negative volume is what CAD people say using Boolean.
turbo is offline   Reply With Quote

Old   February 17, 2017, 02:27
Default
  #19
New Member
 
Rocky
Join Date: Jan 2017
Posts: 21
Rep Power: 9
Rockytime is on a distinguished road
HI turbo and Glen,

I had done the negative volume (booelan) in the geometry.

)Create---> Rectangular Enclosure--> target bodies rotor
3) Create ---> Cylindrical Enclosure ---> target bodies rotor.
2) Create --> Boolean ---> Operation: subtract, target body: select cylindrical enclosure , tool body: select Rotor. Preserve tool bodies? No.

And after that I have create 2 domain.

1st Domain - Rectangular Enclosure

Stationary Domain

- Boundary

1)Inlet (Constant Velocity) ,

2)Outlet (Average Static Pressure, Relative Pressure = 0 Pa,

3)Fluid Wall - Location 4 surface around the Rectangular Enclosure, No Slip Wall Condition.

Rotating Domain

- Boundary

1) Wall Boundary - " Whole Rotor " - Frame type Rotating - No Slip Wall

2) Wall Boundary - Rotating Wall - Frame Type Rotating - Wall Velocity -Counter Rotating Wall - No Slip Wall

Interface

Frozon Rotor

I run the simulation as steady state and the torque value I got from the blade does not increase even though I keep increasing the flow. I expect it
will increase like the graph i show above.

Did I do anything wrong on my setup part? I attach my CCL inside the word file if you need any further information.



StreamLine


Rotor rotating


Counter Rotating Wall -
Attached Files
File Type: docx CEL Code.docx (15.0 KB, 4 views)
Rockytime is offline   Reply With Quote

Old   February 17, 2017, 02:28
Default
  #20
New Member
 
Rocky
Join Date: Jan 2017
Posts: 21
Rep Power: 9
Rockytime is on a distinguished road
Hi Guys, I know it is quite long. But I really hope u guys can help me. thank you. Apologize in advance if I did any basic mistake.
Rockytime 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
Second Derivative Zero - Boundary Condition fu-ki-pa OpenFOAM 11 March 27, 2021 05:28
Simple piston movement in cylinder- fluid models arun1994 CFX 4 July 8, 2016 03:54
problem with solving pipe flow in rotating reference frame rpienika Main CFD Forum 0 March 11, 2016 16:25
Different solutions for Rotating Reference Frame fafou09 FLUENT 1 July 24, 2012 11:06
OpenFOAM on MinGW crosscompiler hosted on Linux allenzhao OpenFOAM Installation 127 January 30, 2009 20:08


All times are GMT -4. The time now is 05:47.