init.cpp

Lattice Boltzmann用于模拟方腔流的程序。 生成的plt文件需用tecplot软件打开。

#include "global.h"
void init()
{
	int x,y;

	/*u=v=0 时候的粒子分布函数*/

	for(y=0;y<max_y;y++)
	{
		for(x=0;x<max_x;x++)
		{
			f[0][y][x]=rho_0*4.0/9.0;
			f[1][y][x]=f[2][y][x]=f[3][y][x]=f[4][y][x]=rho_0/9.0;
			f[5][y][x]=f[6][y][x]=f[7][y][x]=f[8][y][x]=rho_0/36.0;
		}
	}
}