# Inlet Velocity Profile BC - Floating Point exception during solution initialization

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

 March 12, 2012, 05:51 Inlet Velocity Profile BC - Floating Point exception during solution initialization #1 New Member   Join Date: Aug 2010 Posts: 11 Rep Power: 6 Hello, I´ve 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] I´ve defined a vector field function in this form: Function Name: Velocity Profile Definition: [30 * pow((\$\$Position[2]/15),0.22), 0 , 0] Then I´ve 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

 March 14, 2012, 09:44 #2 Member   Ryan Coe Join Date: Jun 2010 Location: Albuquerque, NM Posts: 98 Rep Power: 7 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

 March 14, 2012, 10:47 #3 New Member   Join Date: Aug 2010 Posts: 11 Rep Power: 6 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 I´ve to use if i want to raise x to the power of y? Edit: Okay I think I´ve 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?

 March 14, 2012, 11:14 #4 Member   Ryan Coe Join Date: Jun 2010 Location: Albuquerque, NM Posts: 98 Rep Power: 7 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

 March 14, 2012, 11:21 #5 New Member   Join Date: Aug 2010 Posts: 11 Rep Power: 6 Actually it won´t work with any decimal digit. For example try to use 1.5 or 2.4 for the exponent, it won´t 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.

 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 eRzBeNgEl STAR-CCM+ 6 March 26, 2012 05:16 MaxCFD STAR-CCM+ 3 June 26, 2011 10:31 dada1204 CFX 1 November 5, 2010 03:42 cboss OpenFOAM 1 June 20, 2010 13:02 MURAT Main CFD Forum 0 December 30, 2003 22:08

All times are GMT -4. The time now is 21:44.