Can't impose velocity profile in Code Saturne

[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

[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

[Jundi He]

Hi:

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

Best Regards!
Jundi He