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.
Opening like boundary condition
Forum rules
Please read the forum usage recommendations before posting.
Please read the forum usage recommendations before posting.
Re: Opening like boundary condition
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.