rotating domain

mach000 · April 16, 2010, 12:28

Hi to everybody,
I'm simulating a multifase rotating domain in wich the axis of rotation doesn't match with one of the global ones.

So, I defined a local Coord and specified one of them as rotation axis.

The simulation doesn't work because:

This simulation involves a gravity vector which is not aligned with the rotation axis. Because the solver currently interprets the gravity vector as being in the relative frame, the gravity vector must be made a function of time so that it is constant in the absolute frame (ie, counterrotating in the relative frame).

Could someone help me in order to undestand the problem?

I specify that I have experiences with multiphase simulation with rotating domain. I never had problems because the rotation axis was coincident with one of the global ones.

Thank you in advance

ghorrocks · April 18, 2010, 07:42

Looks like you might have to define gravity as a body force yourself and make it a function which makes it point in the direction you want. Then you can do whatever you like to the gravity vector - but you will have to do the maths to make it point in the correct direction, including the rotating frames of reference stuff.

mach000 · April 19, 2010, 10:52

Hi,
the rotating velocity is "omega". Vertical global axis is Z. So, in a stationary domain the gravity vector as theese components:

g_x=0
g_y=0
g_z=-g

Now, if I have to define a time function in order to have the g vector always alligned to Z axis, I think I have to made zero g_x and g_y defining expression such as

g_x(t)=g_y(t)=g*sin(omega*t).

Is it correct?

bharath · October 24, 2016, 02:04

Thanks for thread.

As you told in the topic that expressions are working.

ghorrocks · October 24, 2016, 02:12

I suspect Ciro's equation has an error - I would not expect the x and y components to be equal in general, and for x and y I would expect one term to be a sine and one to be a cosine (ie: https://en.wikipedia.org/wiki/Rotating_reference_frame).

Also note that I suspect there are high order terms which arise from using a moving acceleration vector like this. If your rotational speed is slow they are likely to be negligible, but as the rotational speed gets faster they will be come important. So be aware this simple approach is probably not universally applicable.

April 16, 2010, 12:28	rotating domain	#1
mach000 Member Ciro Cannavacciuolo Join Date: Mar 2009 Posts: 35 Rep Power: 17	Hi to everybody, I'm simulating a multifase rotating domain in wich the axis of rotation doesn't match with one of the global ones. So, I defined a local Coord and specified one of them as rotation axis. The simulation doesn't work because: This simulation involves a gravity vector which is not aligned with the rotation axis. Because the solver currently interprets the gravity vector as being in the relative frame, the gravity vector must be made a function of time so that it is constant in the absolute frame (ie, counterrotating in the relative frame). Could someone help me in order to undestand the problem? I specify that I have experiences with multiphase simulation with rotating domain. I never had problems because the rotation axis was coincident with one of the global ones. Thank you in advance

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
transient simulation of a rotating rectangle	icesniffer	CFX	1	August 8, 2009 07:25
Single Rotating Domain	fluentnoob	FLUENT	4	June 23, 2009 10:48
Difference b/t rotating wall and rotating domain?	Juno Kim	CFX	2	February 20, 2007 20:08
rotating domain in rotating domain, different axis	Robert Stringer	CFX	3	December 4, 2006 07:04
Size of rotating domain	Adam	CFX	4	November 28, 2006 10:28

April 18, 2010, 07:42		#2
ghorrocks Super Moderator Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 17,701 Rep Power: 143	Looks like you might have to define gravity as a body force yourself and make it a function which makes it point in the direction you want. Then you can do whatever you like to the gravity vector - but you will have to do the maths to make it point in the correct direction, including the rotating frames of reference stuff.

April 19, 2010, 10:52		#3
mach000 Member Ciro Cannavacciuolo Join Date: Mar 2009 Posts: 35 Rep Power: 17	Hi, the rotating velocity is "omega". Vertical global axis is Z. So, in a stationary domain the gravity vector as theese components: g_x=0 g_y=0 g_z=-g Now, if I have to define a time function in order to have the g vector always alligned to Z axis, I think I have to made zero g_x and g_y defining expression such as g_x(t)=g_y(t)=gsin(omegat). Is it correct?

October 24, 2016, 02:04		#4
bharath Senior Member Bharath kumar Join Date: Apr 2009 Posts: 169 Rep Power: 17	Thanks for thread. As you told in the topic that expressions are working.

October 24, 2016, 02:12		#5
ghorrocks Super Moderator Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 17,701 Rep Power: 143	I suspect Ciro's equation has an error - I would not expect the x and y components to be equal in general, and for x and y I would expect one term to be a sine and one to be a cosine (ie: https://en.wikipedia.org/wiki/Rotating_reference_frame). Also note that I suspect there are high order terms which arise from using a moving acceleration vector like this. If your rotational speed is slow they are likely to be negligible, but as the rotational speed gets faster they will be come important. So be aware this simple approach is probably not universally applicable.