请教：LBM 对poiseuille流的模拟

1058343691

各位前辈，我现在在学习编写poiseuille流的模拟，遇到了一些问题：
     我是想模拟poiseuille流从开始到充分发展这个过程，入口处，上下边界都是采用非平衡态外推格式，在出口边界也想用非平衡态外推格式，       
                                //边界处理
                                for(j=1;j<NY;j++)//左右边界
                                        for(k=0;k<Q;k++)
                                        {
                                                rho[NX][j]=rho[NX-1][j];
                        u[NX][j][1]=u[NX-1][j][1];
                       u[NX][j][0]=u[NX-1][j][0];
                               f[NX][j][k]=feq(k,rho[NX][j],u[NX][j])+f[NX-1][j][k]-feq(k,rho[NX-1][j],u[NX-1][j]);
                                                rho[0][j]=rho[1][j];
                           u[0][j][1]=u[1][j][1];
                          u[0][j][0]=u[1][j][0];
                                                f[0][j][k]=feq(k,rho[0][j],u[0][j])+f[1][j][k]-feq(k,rho[1][j],u[1][j]);
                                        }
                                        for(i=0;i<=NX;i++)//上下边界
                                                for(k=0;k<Q;k++)
                                                {
                                                        rho[0]=rho[1];
                          u[0][1]=0;
                          u[0][0]=0;
                        f[0][k]=feq(k,rho[0],u[0])+f[1][k]-feq(k,rho[1],u[1]);
                                rho[NY]=rho[NY-1];
                                u[NY][1]=0;
                           u[NY][0]=0;
        f[NY][k]=feq(k,rho[NY],u[NY])+f[NY-1][k]-feq(k,rho[NY-1],u[NY-1]);
                                                }
     发现程序发散。请问各位高手，这个问题该如何解决

echome

上下边界处有问题，成一维的了，绕流出口一般设置为充分发展流。

1058343691

呵呵，粘的时候不晓得为啥子上下边界所有的设置都变一维的了，谢谢你的指点。
  改了一下，还是发散的，不晓得问题出在哪里了。

echome

那有可能是绕流边界处或者是松弛因子有问题，如果可以的话，把绕流处边界发出来，一起看看是不是有问题。我是将边界分为八段用非平衡外推做的。

[ 本帖最后由 echome 于 2013-5-19 13:51 编辑 ]

1058343691

O(∩_∩)O谢谢你哦，我已经给你发送了短消息了。