Page 1 of 1

Setting pressure on an outlet boundary by specifying rcodcl

Posted: Tue May 31, 2011 6:41 pm
by Gareth Gretton
Hello,

Hopefully this is a very quick question!

I'd like to explicitly set the pressure on all faces of an outlet boundary by specifying the value in the array rcodcl using usclim.f90.

This seems to be pretty simple, as per:

https://code-saturne.info/products/code-saturne/forums/general-usage/906800148
http://code-saturne.org/forum/viewtopic.php?f=9&t=640
http://code-saturne.org/forum/viewtopic.php?f=9&t=627

What I'm a little confused about is exactly which pressure is stored in rcodcl. In the first thread referenced above, the implication is that it is the total pressure, whereas in the final forum thread the implication is that it is the reduced pressure.

Looking at condli.f90 (section 9) suggests that it is the total pressure. Obviously it would be helpful to have this confirmed.

It's worth asking the same about the array rtp as I also wish to initialise the pressure field in usiniv.f90. As I understand, this is the reduced pressure.

Finally, in the post:

http://code-saturne.org/forum/viewtopic ... t=10#p3553

David states as a post script:

"In the formula for the total pressure, I've intentionnally forgotten the reference of the reduced pressure Pred0."

I'm afraid I'm a little confused by this comment - are you able to explain a little?

Thanks,

Gareth