# Applying Rotational Boundary Condition on Surface with Radial Wall

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

July 17, 2012, 13:15
Applying Rotational Boundary Condition on Surface with Radial Wall
#1
New Member

Join Date: Jul 2012
Posts: 11
Rep Power: 6
Hello,

I am new to the world of CFD and I have to say that this forum has been a great asset in over coming many of the issues I had faced, many thanks are due.

Up to now I have been stuck in a problem and cannot find the solution to here, maybe someone has had previous experience with this. I am attempting to model the flow in a certain domain, where the inner wall should be rotating and the outer wall should be stationary as seen in Figure 1.

Thus I applied a rotating boundary condition at the rpm I want to the parts being rotated, but I get an error of the tangential velocity being normal to a surface, which is true as can be seen in the cross section of the part in Figure 2 and a close up in Figure 3.

To try to fix this, I made the whole domain rotate and applied a counter rotating wall condition on the outer wall, thus making it stationary in the absolute frame. In CFX Post this seems to work, as when I plot an isosurface of zero velocity with respect to the stationary frame I get what I expect that the wall has zero velocity (Figure 4 and Figure 5 show part before and after isosurface plot, respectively).

Will this give me an accurate model?
Is there another way of going around doing this?

Thank you
Attached Files
 Figure 1.pdf (39.3 KB, 48 views) Figure 2.pdf (48.2 KB, 31 views) Figure 3.pdf (33.3 KB, 22 views) Figure 4.pdf (41.6 KB, 20 views) Figure 5.pdf (42.8 KB, 19 views)

July 17, 2012, 13:30
#2
New Member

Join Date: Jul 2012
Posts: 11
Rep Power: 6
To validate what I have done I am taking a simple case of two concentric cylinders and see what happens where when I rotate the domain (Figure A) and apply a counter rotating condition on the outer wall. I get the same results with the isosurface (Figure B)

I plot the velocity along the line shown in Figure C, and I get a solution that makes sense with the expected tangential velocity at the inner surface and almost zero at the outer surface as seen in Figure D.

However, when I increase the diameter of the cylinders keeping the gaps constant and plot the same line as in Figure C I get the following plot as seen in Figure E.
Attached Files
 Figure A.pdf (52.8 KB, 8 views) Figure C.pdf (57.3 KB, 7 views) Figure B.pdf (62.7 KB, 5 views) Figure D.pdf (39.0 KB, 8 views) Figure E.pdf (36.2 KB, 7 views)

Last edited by tango711; July 17, 2012 at 15:32.

 July 17, 2012, 14:24 #3 Member   Luis Filipe Fabiani Join Date: Apr 2009 Posts: 43 Rep Power: 9 Tango 771, Isn't this a case of a rotating DOMAIN? What boundary conditions are you using? Best regards

July 17, 2012, 14:34
#4
New Member

Join Date: Jul 2012
Posts: 11
Rep Power: 6

No in reality I this domain is not rotating, only the inner wall sufraces. This is a sub domain of a larger domain which makes up a pump seal.

My inlet and outlet are actually on another subdomin, this part is connected to the rest of the domains through interfaces (frozen rotor interface when needed) above and below it. Figure 6 shows the whole computational domain with this specific part highlighted in green. All surfaces except the outer wall should be rotating.

Please let me know if this helps.

Thanks
Attached Files
 Figure 6.pdf (89.5 KB, 21 views) Figure 6 part 2.pdf (69.3 KB, 12 views)

 July 17, 2012, 14:34 Help me for a model #5 New Member   Josept David Revuelta Acosta Join Date: Jul 2012 Posts: 1 Rep Power: 0 Hello Im David, I am studing my master degrees and I am relatively new in CFD tool, my case is the following: I need to simulate a struture with a fluid impacting it, I am new in AUTODYN and my english isnt very clear, but i will try to explain me: In a Natural Channel Condition, I have a concrete matress on the Channel, I need doing to flow a water into the channel and I want to observe the structural displacement of the matress and its displacements. Is that possible with AUTODYN? and if that is possible i'd appreciate your help. Thanks and regards

 July 17, 2012, 17:18 #6 Member   Luis Filipe Fabiani Join Date: Apr 2009 Posts: 43 Rep Power: 9 tango711, Just to be clear, what is the error? Does CFX crashes or you are saying that the simulation is not correct? Can you plot a cross-section plane with the tangential velocity and vectors? Best regards

July 18, 2012, 10:08
#7
New Member

Join Date: Jul 2012
Posts: 11
Rep Power: 6
Quote:
 Originally Posted by lffabiani tango711, Just to be clear, what is the error? Does CFX crashes or you are saying that the simulation is not correct? Can you plot a cross-section plane with the tangential velocity and vectors? Best regards
Yes, it is the second part. I am wondering if the simulation would be correct if instead of applying a rotational wall boundary condition on the inner wall, I have the whole domain rotate and apply to a counter rotating wall boundary condition on the outer wall.

As for my second post, I was checking to see if what I did make sense so I applied the same conditions to two concentric cylinders, but I got strange results. I attached cross-sections with the tangential velocities in Figure F for a small scale which works, and Figure G which is a large scale, and as you can see half way between the two cylinders the vectors change direction, I am not sure why.

Attached Files
 Image F.pdf (38.6 KB, 12 views) Image G.pdf (43.2 KB, 7 views)

 July 18, 2012, 10:23 #8 Member   Luis Filipe Fabiani Join Date: Apr 2009 Posts: 43 Rep Power: 9 My advice would be to: 1- Turbulence model. Is it the right one? 2- Verify your mesh. Check the Y+ values on walls with the recommended for the turbulence model, size sensitivity check. 3- Verify the convergence (http://www.cfd-online.com/Wiki/Ansys_FAQ) Hope it helps!

July 18, 2012, 10:59
#9
New Member

Join Date: Jul 2012
Posts: 11
Rep Power: 6
Quote:
 Originally Posted by lffabiani My advice would be to: 1- Turbulence model. Is it the right one? 2- Verify your mesh. Check the Y+ values on walls with the recommended for the turbulence model, size sensitivity check. 3- Verify the convergence (http://www.cfd-online.com/Wiki/Ansys_FAQ) Hope it helps!

I am checking it in the laminar case. However, for when I try the large scale (the one with funny results) the solver finishes normally within 5 iterations. While the smaller scale one took around 60 iterations.

Is the quick convergence of the large scale case an issue? I cannot find fast convergence problems in the link you provided or in other parts of the forum.

Thank you,

 July 18, 2012, 11:14 #10 Senior Member     Alex Join Date: Jun 2012 Location: Germany Posts: 1,266 Rep Power: 22 5 Iterations until "convergence" is very unusual. How did you check if your simulation is converged? Run the case for a few hundred iterations and see if the results are more plausible.

 July 18, 2012, 11:59 #11 Member   Luis Filipe Fabiani Join Date: Apr 2009 Posts: 43 Rep Power: 9 Alexander is right. Try setting the minimum iterations to 80 and check if the results change

July 18, 2012, 12:27
#12
New Member

Join Date: Jul 2012
Posts: 11
Rep Power: 6
Quote:
 Originally Posted by flotus1 5 Iterations until "convergence" is very unusual. How did you check if your simulation is converged? Run the case for a few hundred iterations and see if the results are more plausible.
Hello Alexander,

I set my minimum iterations to 100 and I ended up with stranger results than previously; now the vectors change direction twice as in Figure H.

Would having large diameter concentric cylinders be affecting this somehow through the shearing of fluid?

What do you guys think?

Many thanks,
Attached Files
 Figure H.pdf (95.0 KB, 11 views)

 July 18, 2012, 13:47 #13 Senior Member     Alex Join Date: Jun 2012 Location: Germany Posts: 1,266 Rep Power: 22 Would you mind showing the residuals? The result in your last picture does not look like a converged solution. If I get you right, the inner and outer wall of your domain are counter-rotating. For this kind of shear driven flow, slow convergence is a typical issue. It is possible that more than thousand Iterations are necessary. BTW: if you upload your pictures in a picture format, then it is much easier to open them.

July 18, 2012, 14:02
#14
New Member

Join Date: Jul 2012
Posts: 11
Rep Power: 6
Quote:
 Originally Posted by flotus1 Would you mind showing the residuals? The result in your last picture does not look like a converged solution. If I get you right, the inner and outer wall of your domain are counter-rotating. For this kind of shear driven flow, slow convergence is a typical issue. It is possible that more than thousand Iterations are necessary. BTW: if you upload your pictures in a picture format, then it is much easier to open them.
Yes, the whole domain is rotating including the inner wall and I applied a counter rotating wall in the outer wall so that I can model it as if it were stationary.

I am doing this as a check, because in my actual domain, I want to have the inner wall rotating, but cannot due to a "tangential velocity normal to a wall" error, which is because I have a radial wall. Thus I am trying to see if using this method, will achieve similar results.

Hope this clears it a bit more.

I have also attached the summary of the run. Sorry for the pdf files, I will attach in JPEGs from now on.

Thanks again,
Attached Files
 Summary.txt (11.3 KB, 17 views)

 Tags counter, isosurface, radial, rotate, rotation

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post mgdenno OpenFOAM Running, Solving & CFD 13 April 7, 2014 10:58 Yr0gErG FLUENT 3 June 12, 2013 02:12 mart.hein OpenFOAM Programming & Development 1 April 14, 2011 12:44 CFD XUE FLUENT 0 July 9, 2010 02:53 CFD XUE FLUENT 0 July 8, 2010 06:49

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