d2q9.h

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

// D2Q9.h
#ifndef D2Q9_H
#define 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, rho_max;

void Init(char *filename);	// 初始化(申请堆内存，其中调用其它初始化函数)
void init_macro();			// 根据宏观初始流场进行初始化
void init_micro();			// 初始化微观场
void SetParameters();		// 仅改变参数 re, mu, omega 时，从头开始计算
void Finish();				// 释放堆内存，程序结束时调用

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

void collision();			// 碰撞
void flow();				// 流动

void evolution();			// 演化（执行碰撞、流动操作）

void SaveFile();			// 将速度(u,v)、密度 rho 写入相应的文件

#endif