Page 1 of 1

Can't impose velocity profile in Code Saturne

Posted: Tue Feb 14, 2017 5:31 pm
by Jundi He
Hi:

In Code Saturne, I use constant velocity in the inlet surface, and the case can run, but when I use the user law and write my own velocity profile, it failed, said that:

SIGFPE signal (floating point exception) intercepted!

This is my setting of the inlet velocity profile, I think there is nothing wrong with the equation. Anyone knows what's wrong? Thank you!

Best Regards!
Jundi He

Re: Can't impose velocity profile in Code Saturne

Posted: Tue Feb 14, 2017 6:25 pm
by Yvan Fournier
Hello,

I tested your expression under Python, and it seems you start having negative numbers when the radius goes from 0.01 to 0.02.

I'll let you double check.

Regards,

Yvan

Re: Can't impose velocity profile in Code Saturne

Posted: Wed Feb 15, 2017 2:50 pm
by Jundi He
Hi:

Many thanks, I need to use the "abs" to make sure that it is positive value.

Best Regards!
Jundi He