d2q9.h

是格子boltzmann中的一段代码,希望对做多相流的朋友有所帮助.

#ifndef D2Q9_H_
#define D2Q9_H_

const int Q=9;       //方向数
const int M=260;     //列数
const int N=260;     //行数

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 C, rho[N+1][M+1], f[N+1][M+1][Q],  F[N+1][M+1][Q], u[N+1][M+1][2];
int m, m0, i, j, k, id, jd;
double Re, omega;
double  error_rho, error_u, error_v, duration;
double temp_rho[1000], temp_u[1000], temp_v[1000];
time_t start, finish;

void SetPrameter();
double Computefeq(int k,double rho,double u[2]);
void init_macro();
void init_micro();
void evolution();
void output_screen();
void datadeal();

#endif