Opening like boundary condition

Questions and remarks about code_saturne usage
Forum rules
Please read the forum usage recommendations before posting.
Brian Angel

Re: Opening like boundary condition

Post by Brian Angel »

Hello,

Below is the coding that we have used for the free inlet-outlet boundary condition with a given pressure. You will need to change the name of the group from XXXX to the name of the group that you are using and also comment out the lines for specifying a pressure if you do not want to fix a pressure at your boundary. Can you please try this coding and let me know what happens?

call getfbr('XXXX', nlelt, lstelt)
!==========

do ilelt = 1, nlelt

ifac = lstelt(ilelt)

itypfb(ifac) = ifrent

icodcl(ifac,ipr) = 1
rcodcl(ifac,ipr,1) = P0+11200.d0

icodcl(ifac,iu) = 3
icodcl(ifac,iv) = 3
icodcl(ifac,iw) = 3

rcodcl(ifac,iu,3) = 0
rcodcl(ifac,iv,3) = 0
rcodcl(ifac,iw,3) = 0

enddo


Regards,

Brian Angel.
vgm567
Posts: 17
Joined: Tue Oct 14, 2014 3:51 pm

Re: Opening like boundary condition

Post by vgm567 »

Thanks Brian. This code snippet works.


Brian Angel wrote:Hello,

Below is the coding that we have used for the free inlet-outlet boundary condition with a given pressure. You will need to change the name of the group from XXXX to the name of the group that you are using and also comment out the lines for specifying a pressure if you do not want to fix a pressure at your boundary. Can you please try this coding and let me know what happens?

call getfbr('XXXX', nlelt, lstelt)
!==========

do ilelt = 1, nlelt

ifac = lstelt(ilelt)

itypfb(ifac) = ifrent

icodcl(ifac,ipr) = 1
rcodcl(ifac,ipr,1) = P0+11200.d0

icodcl(ifac,iu) = 3
icodcl(ifac,iv) = 3
icodcl(ifac,iw) = 3

rcodcl(ifac,iu,3) = 0
rcodcl(ifac,iv,3) = 0
rcodcl(ifac,iw,3) = 0

enddo


Regards,

Brian Angel.
Post Reply