mesh.c

这是一个用SIMPLE法模拟的矩形流场的程序

#include "global.h"
void mesh()			//划分网格
{
	x=1;			
	y=1;
	l=12;			
	m=12;
	dx=x/(l-2);	//dx=0.1		
	dy=y/(m-2); //dy=0.1
	
	//设定u节点坐标
	xu[1]=xu[2]=0;
	for(i=3;i<=l;i++)
		xu[i]=xu[i-1]+dx;

	yu[1]=0;
	yu[2]=dy/2;
	for(j=3;j<=m-1;j++)
		yu[j]=yu[j-1]+dy;
	yu[m]=y;
	
	//设定v坐标
	xv[1]=0;
	xv[2]=dx/2;
	for(i=3;i<=l-1;i++)
		xv[i]=xv[i-1]+dx;
	xv[l]=x;

	yv[1]=yv[2]=0;
	for(j=3;j<=m;j++)
		yv[j]=yv[j-1]+dy;

	//设定压力节点坐标
	for(i=1;i<=l;i++)
		xp[i]=xv[i];
	for(j=1;j<=m;j++)
		yp[j]=yu[j];
	

}