Re: Imposed velocity profile
Posted: Tue Apr 22, 2014 10:55 pm
Hello,
I'm using Code_Saturne 3.0.3.
I want to impose a non-standard boundary condition : a suction velocity v(x) profile on the upper face of a parallelepiped:
v=-0.12/7*Re^(-1/7)
This topic and the user guide strongly helped me write the following code:
Unfortunately, it doesn't work. Somebody can help me? And tell me what's wrong with the code or if I forgot something? I attach the compile.log file.
Thanks in advance for your help.
Best regards,
Loïc.
I'm using Code_Saturne 3.0.3.
I want to impose a non-standard boundary condition : a suction velocity v(x) profile on the upper face of a parallelepiped:
v=-0.12/7*Re^(-1/7)
This topic and the user guide strongly helped me write the following code:
Code: Select all
call getfbr ('Farfield',nlelt,lstelt)
!==========
do ilelt = 1,nlelt
ifac = lstelt(ilelt)
iel = ifabor(ifac) !cell adjacent to boundary face
do iphas = 1, nphas
itypfb(ifac,iphas) = iindef
icodcl(ifac,iv(iphas)) = 1
rcodcl(ifac,iv(iphas),1) = (0.12d0/7.d0)*((cdgfbo(1,ifac))/((0.0000333d0)**(-1.d0/7.d0)))
rcodcl(ifac,iu(iphas),1) = 1.d0
rcodcl(ifac,iw(iphas),1) = 0.d0
rcodcl(ifac,iu(iphas),2) = rinfin
rcodcl(ifac,iv(iphas),2) = rinfin
rcodcl(ifac,iw(iphas),2) = 0.d0
rcodcl(ifac,ik(iphas),1) = 0.000833d0
rcodcl(ifac,iomg(iphas),1) = 0.1318d0
enddo
enddo
Unfortunately, it doesn't work. Somebody can help me? And tell me what's wrong with the code or if I forgot something? I attach the compile.log file.
Thanks in advance for your help.
Best regards,
Loïc.