I'm using Code_Saturne version 3.0.3.
Currently, I'm trying to simulate the flow on a fat plane. I chose a rectangular domain (2D). The definition of the geometry consists of two joint blocks (using face joinning).
I put a wall condition to the bottom wall, symmetries conditions to the side walls, an inlet and an outlet condition to the extremities. I also wrote the following code for the upper face of the domain:
Code: Select all
call getfbr ('Farfield',nlelt,lstelt)
!==========
do ilelt = 1,nlelt
ifac = lstelt(ilelt)
iel = ifabor(ifac) !cell adjacent to boundary face
itypfb(ifac) = iparoi
icodcl(ifac,iv) = 1
icodcl(ifac,iu) = 3
icodcl(ifac,iw) = 1
icodcl(ifac,ik) = 1
icodcl(ifac,iomg) = 1
rcodcl(ifac,iv,1) = 0.d0
rcodcl(ifac,iu,3) = 0.d0
rcodcl(ifac,iw,1) = 0.d0
rcodcl(ifac,ik,1) = 0.000833d0
rcodcl(ifac,iomg,1) = 0.1318d0
enddo
My simulation run but the results are not good and I don't understand why.
I join a pdf with the graphic of the evoluation of the Cp (Coefficient of pressure) on the upper wall and the graphic with the inlet profil and the oulet profil. We see that the flow has not accelerated between the inlet and the outlet of the domain. The conservation law is not respected. Normally due to the development of boundary layer on the wall, the flow will have to accelerate but it's not the case here.
Anyone has an idea about what's wrong with the simulation? I put my .xlm file and my mesh file if someone wants to test it.
Thanks in advance for your help.
Best regards,
Loïc.