Head Loss Region Syntax

[Vance-Turner]

I am trying to define a cylindrical head loss region (to approximate an actuator disk). The disk has some thickness in x; but is cylindrical in y and z (the axis of the cylinder is in the x-direction). The dependence on y would be something like: y <= sqrt(1-z^2); but how do I specify this in code-saturne? Is this even possible? I've searched for such syntax; but found none.

[Brian Angel]

Hello Vance,

You should be able to use the source terms subroutine "cs_user_source_terms.f90" in the SRC/REFERENCE directory to define the head losses for the zone that represents the actuator disc.

Best regards,

Brian Angel.

[Jacques Fontaine]

Vance,

You may use the subroutine uskpdc (uskpdc.f90) to impose head losses. You also can use ustvnv (cs_user_source_terms.f90) as said Brian to impose momentum source terms.

Regards,

[Vance-Turner]

Thank-you very much. I will look into it.