CLBM顶盖驱动流.docxVIP

这一段代码源自何雅玲的格子BOLTZMANN方法的理论及应用其后面源代码，经过调试可用，输出文件可以导入到TECPLOT中。使用时，在VISUAL STUDIO中新建一个C++工程，代码赋值进去调试即可。#includeiostream#includecmath#includecstdlib#includeiomanip#includefstream#includesstream#includestringusingnamespace std; constint Q=9; //D2Q9模型constint NX=256; //X方向constint NY=256; //Y方向constdouble U=0.1; //顶盖速度int e[Q][2]={{0,0},{1,0},{0,1},{-1,0},{0,-1},{1,1},{-1,1},{-1,-1},{1,-1}};double w[Q]={4.0/9,1.0/9,1.0/9,1.0/9,1.0/9,1.0/36,1.0/36,1.0/36,1.0/36};double rho[NX+1][NY+1],u[NX+1][NY+1][2],u0[NX+1][NY+1][2],f[NX+1][NY+1][Q],F[NX+1][NY+1][Q];int i,j,k,ip,jp,n;double c,Re,