setup0.c

数值传热学NHT C语言源代码及解释 陶文铨院士 的经典例题中3个题目的解答

#include "global_var.h" 
void default_setting()
{ 
  int i;
  nfmax=10;              
  np=11;
  nrho=12;
  ngam=13;
    /****************
	     nf[i][j][0]~到f[i][j][nfmax]是这样分配的:
		 nf[i][j][1]: x方向的速度u(sorry,现在我的认识水平只停留在直角坐标系上） 
         nf[i][j][2]: y方向的速度v
         nf[i][j][3]: 压力修正值p'
         nf[i][j][4]: 温度t
		 nf[i][j][5]~nf[i][j][10]: 其它φ变量
         nf[i][j][11]: 压力p
		 nf[i][j][12]: 密度ρ
         nf[i][j][13]: 扩散系数Γ
   *********************/
  lstop=0; //bool .false.
  for (i=1;i<=nfmax;i++) 
  {
	  lsolve[i]=0; //不求解nf[i][j][0]~到f[i][j][nfmax]            
	  lblk[i]=1;   //进行块修正                                    
	  ntimes[i]=1; //在solve子程序中求解变量f[i][j][nf]的扫描一次  
  }
  for (i=1;i<=13;i++)
  {
	  lprint[i]=0;  //不打印中间结果
	  relax[i]=1;   //relax=1其实就是不松弛了                      
  }
  mode=1;           //直角坐标系
  last=5;           //迭代5次                                      
  time=0;           //开始计算的时间点		 
  dt=1.0e10;        //时间步长                                     	 
  ipref=1;
  jpref=1;          //压力参考点的(i,j)为(1,1)
  rhocon=1;         //不可压缩流的密度为1
}