Re: Opening like boundary condition
Posted: Wed Jul 15, 2015 9:31 am
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.
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.