Setting pressure on an outlet boundary by specifying rcodcl
Posted: Tue May 31, 2011 6:41 pm
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
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