code_saturne User's Forum

Dear Support and Users,
 
I try to set velocity vector direction normal to the inlet using surfbo(...) function. Unfortunately I get an error. Please take a look at the attachments (uscpl.f90 and listing file containing error information). The speed vector directions set by using surfbo(...) are currently commented. When setting speed vector direction without surfbo(...) the calculation runs ok.   Am I am doing something wrong?

Look forward for your advice.

With kind regards

Norbert

Hello,

Why are you some times using only the first 2 components of surfbo, and sometimes with a different sign for those components ?

Also, it is surprising that you only use surfbo, not velocity*surfbo/norm(surfbo).

Note that surfbo is not a unit normal, but the unit normal * face surface.

Best regards,

Yvan