顶盖驱动流.cpp

#include
#include
#include
#include
#include
#include
#include

using namespace std;	///省略名为std的函数书写 std::
const int Q= 9;			//D2Q9 model. 规则网格模型
const int NX= 256;		//x 方向
const int NY= 256;		//y 方向
const double U= 0.1;	//顶盖速度
//c，网格步长，时间步长都为1
int e[Q][2]= {{00}{10}{01}{-10}{0-1}{11}{-11}{-1-1}{1-1}};			//根据模型进行速度配置，二维数组，每个点两个方向，9行2列的数组。
double w[Q]= {4.0/91.0/91.0/91.0/91.0/91.0/361.0/361.0/361.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];		//网格为256×256，节点数为257*257.
int ijkipjpn;
double cRedxdyLxLydtrho0P0tau_fniuerror;

void init（）;
double feq（int kdouble rhodouble u[2]）;
void evolution（）;
void output（int m）;