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

Problems with Small Values in the Rigid Body Solver

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

Like Tree1Likes
  • 1 Post By ghorrocks

Reply
 
LinkBack Thread Tools Display Modes
Old   September 24, 2012, 22:41
Post Problems with Small Values in the Rigid Body Solver
  #1
New Member
 
Josh Yen
Join Date: Jul 2011
Posts: 14
Rep Power: 6
dae4god is on a distinguished road
My simulation involves rotating a rigid body at very small angular accelerations of the order of \alpha = 4.18\times 10^{-3} rad s^{-2}. However, when I look at the Angular Acceleration Z in one of the CFD-Solver monitors (see screenshot), it shows \alpha = 0 rad s^{-2}.

Does anyone know why this happens, or how I can rotate my rigid body at small angular accelerations?
Attached Images
File Type: jpg Angular Acceleration.JPG (42.4 KB, 13 views)
dae4god is offline   Reply With Quote

Old   September 25, 2012, 08:11
Default
  #2
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,951
Rep Power: 85
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
Have you zoomed in on the graph axes? You cannot see the small number with such a wide range.
ghorrocks is offline   Reply With Quote

Old   September 25, 2012, 08:58
Default
  #3
New Member
 
Josh Yen
Join Date: Jul 2011
Posts: 14
Rep Power: 6
dae4god is on a distinguished road
Great suggestion!

I just looked at the CFX-Solver monitor again and I'm not sure how to zoom in. Instead, I exported the plot data and looked at the individual numerical values. Unfortunately, the Angular Acceleration still turns out to be \alpha=0.00E0 rad s^{-2}.

It still seems as if the Rigid Body Solver won't accept small values of Angular Acceleration.
dae4god is offline   Reply With Quote

Old   September 25, 2012, 19:00
Default
  #4
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,951
Rep Power: 85
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
Right click on the chart and set the properties to zoom in.

I would be very surprised if there is a minimum angular acceleration which it accepts. I suspect your issue is more that you are not extracting the acceleration from the solver correctly.

Can you post your CCL?
ghorrocks is offline   Reply With Quote

Old   September 25, 2012, 21:07
Default
  #5
New Member
 
Josh Yen
Join Date: Jul 2011
Posts: 14
Rep Power: 6
dae4god is on a distinguished road
Ah, that's a clever zoom feature!

I've zoomed in and the Angular Acceleration still looks like it is sitting on 0 rad/s^2. I've attached another screenshot.

I wasn't sure if you wanted the actual CFX-Pre .CCL file, or just wanted to see the settings I used. I've done the latter and attached the CFX-Solver Output in the .txt file which shows the settings I used. If you would like the actual CFX-Pre .CCL file, let me know.

Thanks for all your suggestions and help!
Attached Images
File Type: jpg Angular Acceleration Zoom.JPG (39.6 KB, 11 views)
Attached Files
File Type: txt CFX-Solver Output.txt (34.5 KB, 10 views)
dae4god is offline   Reply With Quote

Old   September 25, 2012, 22:02
Default
  #6
Member
 
Hamed Abdul Majeed
Join Date: May 2012
Posts: 72
Rep Power: 5
hamed.majeed is on a distinguished road
Hi,
give the complete output file, also attach a pic showing rigid body convergence with force, torque, and mesh motion convergence.
hamed.majeed is offline   Reply With Quote

Old   September 25, 2012, 22:23
Default
  #7
New Member
 
Josh Yen
Join Date: Jul 2011
Posts: 14
Rep Power: 6
dae4god is on a distinguished road
Hi Hamed,

I've attached the complete CFX-Solver Output and a screenshot of the Torque and Motion convergence (there are no forces applied on the Rigid Body so there is no Force convergence plot).

Also, I have tried running the same simulation using different Rigid Body Control settings including:
  1. Every Timestep
  2. Every Coefficient Loop
  3. General Coupling Control
  4. With and without Simo Wong Angular Momentum control

As you can see from the Output File, the screenshots in this thread are for the General Coupling Control case. Only a couple of timesteps were solved, but I can restart the simulation to test anyone's theories.
Attached Images
File Type: jpg Rigid Body Convergence.JPG (47.5 KB, 6 views)
Attached Files
File Type: zip CFX-Solver Output Complete.zip (57.0 KB, 5 views)
dae4god is offline   Reply With Quote

Old   September 25, 2012, 22:43
Default
  #8
Member
 
Hamed Abdul Majeed
Join Date: May 2012
Posts: 72
Rep Power: 5
hamed.majeed is on a distinguished road
Hi,
I had my own issues regarding rigid body convergence. I have also tried all options but rigid body does not show any convergence.
rigid body convergence issue
hamed.majeed is offline   Reply With Quote

Old   September 25, 2012, 22:50
Default
  #9
New Member
 
Josh Yen
Join Date: Jul 2011
Posts: 14
Rep Power: 6
dae4god is on a distinguished road
Hi Hamed,

Thanks for your reply and I had a look at your post. I believe you can set the convergence target for the motion, force and torque; similar to setting target residuals for mass and momentum.

Under the Rigid Body Control:
  1. Select either Every Coefficient Loop or General Coupling Control
  2. Check Internal Coupling Data Transfer Control
  3. Check Mesh Motion/Force/Torque Data Transfer Control
  4. Check Convergence Target
  5. Set Convergence Target
Was this what you were asking in your other thread? Have you tried using these settings?
dae4god is offline   Reply With Quote

Old   September 25, 2012, 23:09
Default
  #10
Member
 
Hamed Abdul Majeed
Join Date: May 2012
Posts: 72
Rep Power: 5
hamed.majeed is on a distinguished road
I applied residuals for torque, force, and mesh motion. Still divergent.

hamed.majeed is offline   Reply With Quote

Old   September 25, 2012, 23:17
Default
  #11
Member
 
Hamed Abdul Majeed
Join Date: May 2012
Posts: 72
Rep Power: 5
hamed.majeed is on a distinguished road
Also I tried.

hamed.majeed is offline   Reply With Quote

Old   September 25, 2012, 23:22
Default
  #12
Member
 
Hamed Abdul Majeed
Join Date: May 2012
Posts: 72
Rep Power: 5
hamed.majeed is on a distinguished road
I also tried.

hamed.majeed is offline   Reply With Quote

Old   September 26, 2012, 01:37
Default
  #13
New Member
 
Josh Yen
Join Date: Jul 2011
Posts: 14
Rep Power: 6
dae4god is on a distinguished road
For the torque, force and motion convergence, try bumping up the maximum number of iterations to say 100. Does this change the results?
dae4god is offline   Reply With Quote

Old   September 26, 2012, 14:38
Default
  #14
Member
 
Hamed Abdul Majeed
Join Date: May 2012
Posts: 72
Rep Power: 5
hamed.majeed is on a distinguished road
Hi,

I have tried that also, results have improved yet same behavior.
I haven't explained my case. A rigid body falls under gravity in a column of fluid.

As I have shown that the rigid body convergence diverges. However, in ANSYS help I found this.

Force Convergence

It is defined as ratio of difference of forces for previous and present time step divided by the greater of the two forces.
Now, in my case the rigid body falls under gravity.
At each time step the force acting on the body differs as the rigid body rotates and translates. Is that the explanation of the curve behavior???

Also torque convergence is defined similar to force convergence but for torque.
hamed.majeed is offline   Reply With Quote

Old   October 1, 2012, 19:57
Default
  #15
New Member
 
Josh Yen
Join Date: Jul 2011
Posts: 14
Rep Power: 6
dae4god is on a distinguished road
Hi Hamed,

This is a screenshot of my Rigid Body Convergence for the simulation, run over a much longer time. You can see that the motion convergence is quite smooth, but the torque convergence is a bit more scattered.

Does this help?
Attached Images
File Type: jpg Rigid Body Convergence 2.jpg (90.4 KB, 11 views)
dae4god is offline   Reply With Quote

Old   October 4, 2012, 00:35
Default
  #16
Member
 
Hamed Abdul Majeed
Join Date: May 2012
Posts: 72
Rep Power: 5
hamed.majeed is on a distinguished road
Hi,
Thanks for the reply. Is this your final curve. I mean is your solution convergent for the curves shown?
How can we say the torque of rigid body converges in this scenario!
What were the rigid body settings you used.
Thank you

Regards
Hamed
hamed.majeed is offline   Reply With Quote

Old   October 4, 2012, 01:41
Default
  #17
New Member
 
Josh Yen
Join Date: Jul 2011
Posts: 14
Rep Power: 6
dae4god is on a distinguished road
Hi Hamed,

The image in my last post shows the convergence levels for Motion and Torque, once all the other (mass and momentum) residuals had converged to 1E-06.

Not exactly sure if we can say that the torque and motion have converged.

You can see my Rigid Body settings from the CFX-Solver Output I posted earlier in the thread.
dae4god is offline   Reply With Quote

Old   October 4, 2012, 02:34
Default
  #18
Member
 
Hamed Abdul Majeed
Join Date: May 2012
Posts: 72
Rep Power: 5
hamed.majeed is on a distinguished road
Hi,
That means the solution is not acceptable. As rigid body has not converged... shouldn't it converge for solution to be acceptable?
Thanks
hamed.majeed is offline   Reply With Quote

Old   October 4, 2012, 07:44
Default
  #19
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,951
Rep Power: 85
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
But how tight convergence is required? It may already be sufficiently converged.

A sensitivity study would help here, but as you do not appear to be able to converge tighter that might be difficult.

Of course I suggested ages ago you do a test with an initial condition of cross flow to see if this is a low Reynolds number affect but you appear to have ignored that suggestion.
hamed.majeed likes this.
ghorrocks is offline   Reply With Quote

Old   October 5, 2012, 07:40
Default
  #20
Member
 
Hamed Abdul Majeed
Join Date: May 2012
Posts: 72
Rep Power: 5
hamed.majeed is on a distinguished road
I will try that, thnx ghorrocks
hamed.majeed is offline   Reply With Quote

Reply

Tags
acceleration, cfx, rigid bodies

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Difficulty in calculating angular velocity of Savonius turbine simulation alfaruk CFX 8 December 3, 2013 17:51
Rigid body - pivoting about a fixed point Catenary CFX 0 April 26, 2012 02:07
Numerical errors in nested domain with pre-calculated boundary values Arnoldinho OpenFOAM Running, Solving & CFD 3 April 4, 2012 10:31
Problems with SUPG body force term FEM question Main CFD Forum 0 January 21, 2006 18:51
CFX 5.5 Roued CFX 1 October 2, 2001 16:49


All times are GMT -4. The time now is 08:53.