Page 2 of 2

Re: Opening like boundary condition

Posted: Wed Jul 15, 2015 9:31 am
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.

Re: Opening like boundary condition

Posted: Wed Jul 15, 2015 4:13 pm
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.