main.c

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

#include "global.h"
main()
{	
	int inter=0;
	relax=0.1;//设置松弛因子
	relaxp=0.5;
	filehandle1=fopen("result/resul.dat","w");
	filehandle2=fopen("result/relax.dat","w");
	mesh();
	t=0;
	init();
	fprintf(filehandle1," TITLE = Velocity Vectors\nVARIABLES =X\nY\nU\nV\nP\nT\n");
	fprintf(filehandle2,"time inter remain\n");
	output();
	for(time=0;time<=300;time++)
	{	inter=0;//迭代的次数
		remain=0;
		do
			{	
				remain=0;
				simple();
				bound();	//计算出口边界速度
				inter=inter+1;
				
			}while(fabs(remain/flowin)>1e-4);//迭代次数
		//温度求解
	//	temperature();
		output();
		fprintf(filehandle2,"%d,%d,%f\n",time,inter,remain);
		for(i=1;i<=l;i++)
			for(j=1;j<=m;j++)
			{
				u0(i,j)=u(i,j);
				v0(i,j)=v(i,j);
				tem0(i,j)=tem(i,j);
				p[i][j]=0;
			}
			printf("\n%d,%f,",inter,remain);

	}
	fclose(filehandle1);
	fclose(filehandle2);
	return 0;
}