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

Floating point exception error [High rotating speed tube]

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 5, 2014, 23:49
Default Floating point exception error [High rotating speed tube]
  #1
New Member
 
wen kang
Join Date: Feb 2014
Posts: 9
Rep Power: 3
wkwong is on a distinguished road
Hi,

I have been encountering this Floating point exception error [overflow]/[divide by zero] whenever I run my simulation over 1200 rpm. Below that, its running well. Once over 1200 rpm it will gives me this problem. I understand that this is because the solution diverge, but I cannot identify the reason why, because all of my boundary conditions (mine is really simple project) seems correct. I need to run until 6000 rpm.

My project:
Cooling of rotating tube, rotating at 6000 rpm around an axis outside of itself and convection is allowed everywhere except the base of the tube (cone shape horizontally). Its inital condition is 369k and surrounding is 300k.




Any help PLEASE??

Updates: I managed to get it running at 6000 rpm without getting simulation error immediately by reducing the time step to 1E-14s (unbelievable low), but the residual usually suddenly rise to super huge value after around 300 iterations and eventually I will get the same error again...

my email is masterlancer@hotmail.com

Last edited by wkwong; February 12, 2014 at 11:43.
wkwong is offline   Reply With Quote

Old   February 6, 2014, 01:45
Default
  #2
New Member
 
Join Date: Mar 2009
Location: Belgium
Posts: 13
Rep Power: 8
rabat is on a distinguished road
Hi,

I would try to run the calculation only isothermal.
If it 's converge, than try to modified you mesh, or change the under relaxation factors.

Regards
rabat is offline   Reply With Quote

Old   February 6, 2014, 07:27
Default
  #3
New Member
 
wen kang
Join Date: Feb 2014
Posts: 9
Rep Power: 3
wkwong is on a distinguished road
hi, @rabat, tried and it did not even converge only isothermally, at 6000 rpm.
wkwong is offline   Reply With Quote

Old   February 6, 2014, 20:33
Default
  #4
Senior Member
 
Join Date: Nov 2010
Posts: 418
Rep Power: 8
me3840 is on a distinguished road
Is this transient or steady? A tube rotating outside of its axis sounds like an unsteady problem. Are you changing your timestep to account for the new higher rotation rate?

Does it diverge immediately or after some amount of time?

What kind of mesh are you using? Is it still adequate for the higher rotation rate?
me3840 is offline   Reply With Quote

Old   February 6, 2014, 20:43
Default
  #5
New Member
 
wen kang
Join Date: Feb 2014
Posts: 9
Rep Power: 3
wkwong is on a distinguished road
@me3840,

It is transient. You mean make the time step smaller?

At 6000 rpm, it diverge almost immediately, within 3 iterations.

My model is
Implicit unsteady
Coupled Flow
3D
Laminar
IF97 (water)

My mesh is:
Polyhedral, Prism layer, surface remesher, surface wrapper

I have around 690k cells.. should be fine thou.
wkwong is offline   Reply With Quote

Old   February 8, 2014, 05:46
Default
  #6
New Member
 
wen kang
Join Date: Feb 2014
Posts: 9
Rep Power: 3
wkwong is on a distinguished road
Quote:
Originally Posted by me3840 View Post
Is this transient or steady? A tube rotating outside of its axis sounds like an unsteady problem. Are you changing your timestep to account for the new higher rotation rate?

Does it diverge immediately or after some amount of time?

What kind of mesh are you using? Is it still adequate for the higher rotation rate?
I tried to decrease the time step REALLY small to around 1E-14 s and thanks to you I am able to get the simulation running without the error as mentioned previous. All the parameters's residual fall to somewhere 0.01 however, my energy and continuity residual is at 10^5 region. The simulation can continue, but why is this so?
wkwong is offline   Reply With Quote

Old   February 10, 2014, 07:12
Default
  #7
New Member
 
wen kang
Join Date: Feb 2014
Posts: 9
Rep Power: 3
wkwong is on a distinguished road
Any would be greatly appreciate please. Thank you.
wkwong is offline   Reply With Quote

Old   February 12, 2014, 02:42
Default
  #8
Member
 
Gajendra Gulgulia
Join Date: Apr 2013
Location: India
Posts: 86
Rep Power: 4
ggulgulia is on a distinguished road
Hey Wkong

I think 1e-14 is too small a time step. Your simulation will take a lot of time to complete even 1e-4 second. I suggest check your courant number. Sometimes it's value is 50. If it's so, then change it to 5 and try running the simulation. 1e-5 second with 10-20 inner iteration should be good enough. You should approach the problem in as simple way as possible and then go on to add the complexities as desired when you see your simulation actually works. Therefore I recommend to change the fluid to water constant density or ideal gas liquid H2O

Secondly you have given Laminar flow for the turbulence model. Flow of water at 6000 RPM cannot be laminar. I suggest you change the model to k-e with wall y+ model.

Thirdly your mesh size for this simulation seems very small. Can you post a picture of the domain ? I need to have a look at the grid spacing near the wall boundaries and see if it could be resolved into 2-d axisymmetric problem?
ggulgulia is offline   Reply With Quote

Old   February 12, 2014, 11:32
Default
  #9
New Member
 
wen kang
Join Date: Feb 2014
Posts: 9
Rep Power: 3
wkwong is on a distinguished road
Hi ggulgulia, thanks I have posted my mesh size and my model. It is a fluid region rotating about the Z axis at 6000 rpm. Base of the cone is insulated and the rest is allowed convection at 500 W/m2K
wkwong is offline   Reply With Quote

Old   February 14, 2014, 14:09
Default
  #10
Member
 
Gajendra Gulgulia
Join Date: Apr 2013
Location: India
Posts: 86
Rep Power: 4
ggulgulia is on a distinguished road
Hey Wkwong

I still am wondering if you are using RBM or Overset method for making the tube rotate. Regardless of that you can check the points I have mentioned below...

1. I would suggest that you check with a different turbulence model
2. Try to increase the time step to 1e-5 second
3. Check your initialization value. Sometimes the problem occurs due to improper initialization
4. Try reducing the URF for pressure to 0.5 and then ramp it up to 0.9 after 100-150 iterations.
5. Check the courant number. Sometimes the default value is 50. If it's so then bring it down to 5 or somewhere less than 5 but greater than 1.
6. Go for a first order differencing scheme for 100 iterations and then change the differencing scheme to second order implicit.
ggulgulia is offline   Reply With Quote

Reply

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
Floating point exception with pimpleDyMFoam ebah6 OpenFOAM Running, Solving & CFD 8 May 13, 2013 09:01
Inlet Velocity Profile BC - Floating Point exception during solution initialization Janshi STAR-CCM+ 4 March 14, 2012 11:21
simpleFoam Floating point exception error -help sudhasran OpenFOAM Running, Solving & CFD 3 March 12, 2012 17:23
Pipe flow in settlingFoam floating point exception jochemvandenbosch OpenFOAM Running, Solving & CFD 4 February 16, 2012 04:24
block-structured mesh for t-junction Robert@cfd ANSYS Meshing & Geometry 20 November 11, 2011 05:59


All times are GMT -4. The time now is 01:39.