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

Inlet Velocity Profile BC - Floating Point exception during solution initialization

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 12, 2012, 05:51
Default Inlet Velocity Profile BC - Floating Point exception during solution initialization
  #1
New Member
 
Join Date: Aug 2010
Posts: 11
Rep Power: 7
Janshi is on a distinguished road
Hello,

Ive encountered a very weird problem during the implementation of a field function to describe an inlet velocity problem. I wish to use a velocity profil as inlet bc in this form: v = [k(zref)*(z/zref)^n , 0 , 0]

Ive defined a vector field function in this form:

Function Name: Velocity Profile
Definition: [30 * pow(($$Position[2]/15),0.22), 0 , 0]

Then Ive defined a local coordinate system for my simulation case, called cartesian 1. For the initial conditions and the inlet bc i selected for the velocity nodes for the coordinate System "Laboratory->Cartesian 1" and for the method "Field Function" and so on...
Once I try to initialize the simulation I get following error:

"Error

A floating point exception has occured: floating point exception [invalid operation]. The specific cause cannot be identified. "

Can anyone imagine what the problem could be?

Best regards
Janshi is offline   Reply With Quote

Old   March 14, 2012, 09:44
Default
  #2
Member
 
Ryan Coe
Join Date: Jun 2010
Location: Albuquerque, NM
Posts: 98
Rep Power: 7
ryancoe is on a distinguished road
That's tough...

Have you experimented with simplifying your field function ( to something like [$$Position[2], 0 , 0] ) and see if that resolves the issue. I know that doesn't fix your problem, but it may help you to track down the bug.
__________________
Ryan
ryancoe is offline   Reply With Quote

Old   March 14, 2012, 10:47
Default
  #3
New Member
 
Join Date: Aug 2010
Posts: 11
Rep Power: 7
Janshi is on a distinguished road
Thanks for your answer and your tip ryan.
Okay now things becomes really strange. The simplified vectorfunction [$$Position[2], 0 , 0] operates also a more complicated function like [30*($$Position[2]/15), 0, 0] operates too. It does not work with an exponential syntax like pow(($$Position[2]/15),0.22)...

I am I correct that the expression pow($x,$y) is the syntax Ive to use if i want to raise x to the power of y?


Edit: Okay I think Ive isolated the problem:

It seems that CCM+ does not accept decimal digits for the exponent n in the expression pow($x,n)...
So in my case the exponent n=0.22 is not accepted by CCM+. I checked it with n=1,5 ; 2,5 and so, same errors. When I change the exponent to 1, 2, 3 and so on, no problem occures and he intializes the solution correctly without errors.
Anyone knows the reason for this and how I can fix it?
Janshi is offline   Reply With Quote

Old   March 14, 2012, 11:14
Default
  #4
Member
 
Ryan Coe
Join Date: Jun 2010
Location: Albuquerque, NM
Posts: 98
Rep Power: 7
ryancoe is on a distinguished road
huh! this was so weird I just had to mess with it myself...

From my experimenting, it looks like pow doesn't accept entries to raise a number to a value less than 1. Is this what you see?

Not sure why that would be the case (you should send an email to the CD-adapco support poeple). Maybe as a work around you could create a table and use that to set you boundary values instead...?
__________________
Ryan
ryancoe is offline   Reply With Quote

Old   March 14, 2012, 11:21
Default
  #5
New Member
 
Join Date: Aug 2010
Posts: 11
Rep Power: 7
Janshi is on a distinguished road
Actually it wont work with any decimal digit. For example try to use 1.5 or 2.4 for the exponent, it wont work either.

I will contact the CD-Adapco support. The idea with the table as a workaround is good idea, thanks for your help Ryan.
Janshi 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
Export & Import Velocity profile as Inlet eRzBeNgEl STAR-CCM+ 6 March 26, 2012 05:16
Finished simulation doesn't start: floating point exception [Divide by zero] MaxCFD STAR-CCM+ 3 June 26, 2011 10:31
what is linear solution and a floating point exception? dada1204 CFX 1 November 5, 2010 03:42
Scale discrete inlet velocity profile with groovyBC cboss OpenFOAM 1 June 20, 2010 13:02
NON-UNIFORM VELOCITY INLET BOUNDARY PROFILE MURAT Main CFD Forum 0 December 30, 2003 22:08


All times are GMT -4. The time now is 19:54.