init_geo.cpp

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

#include "global.h"
void init_geo()
{
	int x,y;
	/* Driven Cavity 方腔驱动流*/
/*
 *                        BC  (SET_U)(速度边界条件)
 *  y=max_y-1 -|       --------
 *             |       * 	  *
 *             | solid * FLUID* solid  (WALL)
 *             |       *      *
 *        y=0 -|       ********
 *             |         solid
 *             |______________________
 *     		        |         |
 * 				   x=0     x=max_x-1
 */
	

	for(y=0;y<max_y;y++)
	for(x=0;x<max_x;x++)
		wall[y][x]=FLUID;
	for(y=0;y<max_y-1;y++) 
		wall[y][0]=wall[y][max_x-1]=WALL;
	for(x=0;x<max_x;x++) wall[0][x]=WALL;
	for(x=0;x<max_x;x++) wall[max_y-1][x]=SET_U;


}