init.c

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

#include "global.h"
void init()
{	
	
	vin1=10;
	vin2=0;
	tin1=30;
	tin2=30;
	vin=vin1+vin2;
	inlen1=2;
	inlen2=0;
	inlen=inlen1+inlen2;
	outlen1=0;
	outlen2=0;
	outlen3=2;
	outlen=outlen1+outlen2+outlen3; 
	flowin=(vin1*inlen1+abs(vin2)*inlen2)*dy;
	for(i=0;i<=l;i++)
		for(j=0;j<=m;j++)
			{
				u0(i,j)=u(i,j)=0;
				v0(i,j)=v(i,j)=0;
				dp(i,j)=0;
				p[i][j]=0;
				dif[1]=0.1;
				dif[4]=0.05;
				den=1;
				tem0(i,j)=tem(i,j)=10;
			}
	for(j=2;j<=inlen1+1;j++)  //入口1速度
		{
			u0(2,j)=u(2,j)=vin1;
			tem0(1,j)=tem(1,j)=tin1;
		}		
	for(j=2;j<=inlen2+1;j++)  //入口2速度
		{
			u0(l,j)=u(l,j)=vin2;
			tem0(l,j)=tem(l,j)=tin2;
		}
	for(i=2;i<=1+outlen1;i++)//出口速度
		v0(i,m)=v(i,m)=flowin/outlen/dx;
	for(i=l-1;i>=l-outlen2;i--)
		v0(i,m)=v(i,m)=flowin/outlen/dx;
	for(j=2;j<=1+outlen3;j++)
		u0(l,j)=u(l,j)=flowin/outlen/dy;
		
	dt=0.1;
	t=0;
}