[推荐]CFD基础知识PDF（包括SIMPLE算法编程实例）对初学者很有帮助

bfo_yyw

谢谢~

Hellofriend

十分感谢！

shirazbj

下载的人还真不少，可能都学成老鸟了吧？

我刚刚下载，用g95编译，运行，没有停下来的意思。忍到1000，终于停了。仔细一看代码，就是迭代1000次。

pdf写的挺清楚，似乎明白多一点主网格和次网格的关系和用处了。

c small iterations
             do its=1,30
                do i=2,31
                   pp(i)=(ae(i)*pp(i+1)+aw(i)*pp(i-1)+s(i))/ap(i)
                end do
             end do

这个its的循环不知道是干什么的?its也没在别处引用.空转?

谢谢

smileskyli

多谢楼主！

liwei620715

顶一下！

shensln

多谢分享！

sartin

顶！

shirazbj

程序里有两个u的定义,如下:

c   u --- velocity on main cell center
c   uf --- velocity on main cell face

这比较奇怪. 既然是交错网格,一个网格存U,另一个网格存P. 为什么U要存在两个网格里呢？

还有这段：

c solving momentum equation on 2nd grid
             do i=2,32
                aue(i)=amax1(abs(den*u(i)*AA(i)/2.),viscos*aa(i)/dx)
     *               -den*u(i)*AA(i)/2.
                auw(i)=amax1(abs(den*u(i-1)*AA(i-1)/2.),viscos*aa(i-1)/dx)
     *               +den*u(i-1)*AA(i-1)/2.
                su(i)=a(i)*(p(i-1)-p(i))
                aup(i)=auw(i)+aue(i)
             end do
c          write(6,*)aue(10),auw(10),aup(10),su(10),' aue,auw,aup,su'
c small iterations
             do i=2,32
                uf(i)=(aue(i)*uf(i+1)+auw(i)*uf(i-1)+su(i))/aup(i)
             end do

系数是用的u的，解出来的却是uf

最后这个插值：
c interpolation
             do i=1,32
                u(i)=.5*(uf(i+1)+uf(i))
             end do
愣把uf给插回到u, 这是要负负得正的效果么?

我觉得在第2个网格上解动量方程时,是把u表示成已知的第一个网格上p的函数。解出来的还是u. 在第1个网格上解连续方程时，是把p表示成第2个网格上的u的函数.Uf没有用.

云泉

谢谢

mitmitmit

多谢分享

chaster0118

感謝，對我幫助很大

echun1226

俺下载了，所以要顶一下！非常感谢呢……

sumao

谢谢了

gxm712

谢了

nedu529

挺好的