24小时热门版块排行榜    

查看: 878  |  回复: 1

花开时节0931

铁虫 (初入文坛)

[求助] 格子boltzmann中求出速度分布函数,如何求压力? 已有1人参与

格子boltzmann中求出速度分布函数,如何求压力?
下面是打算计算宏观量的 Psum是想要计算的压力,rho是密度,其实只要算出密度就好了 。但是我算出来的 最后高达10的10级次方,和迭代步数有关
        for(i=1;i<NX;i++)//计算宏观量
                                for(j=1;j<NY;j++)
                                {
                                        u0[j][0]=u[j][0];
                                        u0[j][1]=u[j][1];
                                        rho[j]=0;
                                        u[j][0]=0;
                                        u[j][1]=0;
                                               
                                        for(k=0;k<Q;k++)
                                        {
                                                f[j][k]=F[j][k];
                                                rho[j]+=f[j][k];
                                                u[j][0]+=e[k][0]*f[j][k];
                                                u[j][1]+=e[k][1]*f[j][k];
                                          }
                                        Psum[j] = rho[j] * c * c / 3;
                                        u[j][0]/=rho[j];
                                        u[j][1]/=rho[j];
                                }
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

孤独的漫步

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
花开时节0931: 金币+10, ★★★很有帮助, 谢谢 2016-02-26 12:59:50
宏观量这一段没问题,应该是碰撞迁移部分或者边界部分出错了
你必须非常努力,才能看起来毫不费力
2楼2016-01-12 12:35:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 花开时节0931 的主题更新
信息提示
请填处理意见