mpi_d2q9.h

格子Boltzmann方法 格子Boltzmann方法是为了保留格子气自动机方法的优点

// mpi_D2Q9.h
#ifndef MPI_D2Q9_H
#define MPI_D2Q9_H

extern int D, Q, Nx, Ny, n;
extern char e[9][2], e1[9], **boundary, bounceback, filename[80];
extern double ***f, ***f1, **rho, **u, **v, *f0;
extern double mu, Re, omega, dx, dt, t, c, c_s;
extern double Lx, Ly, L0, U0, V0, Rho0;

double CPU_time();
void ShowBdy(int nx, int ny);
int MyUpper(char *str);					// 临时使用的处理字符串函数

void ReadParaFile(char *filename);
void InitData(int numprocs, int myid);
void init_macro(int numprocs, int myid);		// 根据宏观初始流场进行初始化
void init_micro(int numprocs, int myid);		// 初始化微观场
void Finish();						// 释放堆内存，程序结束时调用

void feq(double u, double v, double rho, double *f0);	// 平衡态函数计算

void collision(int numprocs, int myid);			// 碰撞

void flow(int numprocs, int myid);			// 流动
void Jacobi(int numprocs, int myid);
void evolution(int numprocs, int myid);			// 演化（执行碰撞、流动操作）
void CollectData(int numprocs, int myid);
void SaveFile();					// 将速度(u,v)、密度 rho 写入相应的文件

#endif