[Abolfazl_cfd]

Hi everybody.

I have already validated my code with compressible flow in a channel with normal boundary condition.
But as I tried to model a periodic boundary condition, I'm getting strange results. The more the code proceeds in time, the more velocity in the center line of the channel decreases. Is it normal??
In my code, there are two columns of ghost cells (two columns at entry on channel and two columns at exit). Here is a the part of code for periodic boundary condition.
I would be so grateful, if you would could help me.

Thank you all.


i=-1
do j=-1,jm+2
!
p (i,j) =p(im-1,j)
t (i,j) =t(im-1,j)
rho(i,j) =rho(im-1,j)
u (i,j) =u(im-1,j)
v (i,j) =v(im-1,j)
!
end do
!
i=0
do j=-1,jm+2
!
p (i,j) =p(im,j)
t (i,j) =t(im,j)
rho(i,j) =rho(im,j)
u (i,j) =u(im,j)
v (i,j) =v(im,j)
!
end do
!
i=im+1
do j=-1,jm+2
!
p (i,j) =p(1,j)
t (i,j) =t(1,j)
rho(i,j) =rho(1,j)
u (i,j) =u(1,j)
v (i,j) =v(1,j)
!
end do
!
i=im+2
do j=-1,jm+2
!
p (i,j) =p(2,j)
t (i,j) =t(2,j)
rho(i,j) =rho(2,j)
u (i,j) =u(2,j)
v (i,j) =v(2,j)
!
end do