code_saturne User's Forum

Hi,

I've been trying to define in cs_user_boundary_conditions.f90 a free inlet/outlet condition with external head loss (just like the GUI). I've found so far that the itypfb for the free inlet/outlet is ifrent.

However, I'm not able to find how the external head losses must be defined. Is there any way to define them? Or there is any example available?

Thanks!

Hello,

I'm not sure. Did you check examples in src/user_examples or the Doxygen documentation ? For Doxygen, use the documentation with the code distribution; the one on the web site might not be up to date.

Regards,

Yvan

Hi,

Yes, I did check the examples and the documentation. I also searched for posts and I couldn't find any clear reference. I found a head loss example but (correct me if I am wrong) I understand this external head loss is something different.

Maybe if you can point me where the free inlet/outlet is being computed I can find a reference for the external head losses.

Regards,

Hi,

I have found what I was asking. I will comment here just in case it is useful to other people:

The answer is in the file "typecl.f90" on line 941 where it says:

b_head_loss(ifac) = rcodcl(ifac,ivar,2)

Therefore, the external head loss constant "K" (found in the GUI) can be set through "cs_user_boundary_conditions.f90" by setting the value of rcodcl(ifac,ivar,2).

Thanks for all the help and support!