tem.c

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

#include "global.h"
void temperature()
{	n=4;
	//ae和aw设定
	for(j=2;j<=m-1;j++)
	{
		for(i=1;i<=l-1;i++)
		{
			d=dif[n]*dy/dx;
			flow=den*dy*u(i+1,j)/2;
			acf();
			aw[i+1][j][4]=d*adf+max(flow,0);
			ae[i][j][4]=aw[i+1][j][4]-flow;
		}
	}

	//an和as设定
	for(i=2;i<=l-1;i++)
	{
		for(j=1;j<=m-1;j++)
		{
			flow=den*dx*v(i,j+1)/2;
			acf();
			as[i][j+1][4]=d*adf+max(flow,0);
			an[i][j][4]=as[i][j+1][4]-flow;
		}
	}
		//边界条件
	for(j=2+inlen1;j<=m-1;j++)
		aw[2][j][4]=0;
	for(j=2+inlen2;j<=m-1;j++)
		ae[l-1][j][4]=0;
	for(i=2;i<=l-1;i++)
		as[i][2][4]=an[i][m-1][4]=0;

	//ap和b设定
	for(i=2;i<=l-1;i++)
		for(j=2;j<=m-1;j++)
		{
			ap0[i][j][4]=den*dx*dy/dt;
			ap[i][j][4]=ae[i][j][4]+aw[i][j][4]+an[i][j][4]+as[i][j][4]+ap0[i][j][4];
			b[i][j][4]=ap0[i][j][4]*tem0(i,j);
		
		}

	//求解温度场
	ist=2;
	jst=2;
	solve();
	n=0;
}