check_convergence.cpp

Lattice Boltzmann用于模拟方腔流的程序。 生成的plt文件需用tecplot软件打开。

#include "global.h"

int check_convergence(int t)		//判断是否收敛
{
  int x,y;
  double p, diff=(double)0;
	
  for(y=0;y<max_y;y++)
    {
      for(x=0;x<max_x;x++)
        {		
          p=1.0/3.0*(f[0][y][x]+f[1][y][x]+f[2][y][x]+f[3][y][x]+f[4][y][x]+f[5][y][x]+f[6][y][x]+f[7][y][x]+f[8][y][x]);  //1./3.为当地音速的平方，计算压强pressure=rho*cs*cs							
          diff=max(diff,abs(p-pressure[y][x]));	   //计算压强的最大波动，也可以考察密度
          pressure[y][x]=p;       //开始时pressure[][]=0,随后将p赋给pressure,当压强趋于稳定时，可中止
        }
    }
  printf("最大误差：t=%6d\t%g\n",t,diff);
  
  if((diff)<=eps)
    {
      return(TRUE);
    }
  else
    {
      return(FALSE);
    }							
}