顶盖驱动流的格子Boltzmann模拟代码.pdfVIP

#includeiostream #includecmath #includecstdlib #includeiomanip #includefstream #includesstream #includestring using namespace std; const int Q=9; const int NX=256; const int NY=256; const double 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,dx,dy,Lx,Ly,dt,rho0,P0,tau_f,niu,error; void init(); double feq(int k,double rho,double u[2]); void evo