用LBM模拟3d方腔流代码

// 三维方腔.cpp : Defines the entry point for the console application.
//

#include “stdafx.h“


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

using namespace std;                        //变量名定义


const int Q=15;

const int NX=80;//尺度
const int NY=40;
const int NZ=40;
const int Re=400;
const double RHO=1.0;
const double U=0.1;                                  
double niu=U*NX/Re;
double tau=3.0*niu+0.5;


int e[Q][3]={  {000}{100}{-100}{010}{0-10}{001}{00-1}{111}{-1-1-1}{11-1}{-1-11}{1-11}{-11-1}{1-1-1}{-111}  };


double w[Q]={2.0/91.0/91.0/91.0/91.0/91.0/91.0/91.0/721.0/721.0/721.0/721.0/721.0/721.0/721.0/72};      //速度场权函数