At first from reality model to compute model
dt = (2/omega-1)dxdx/6.0/mu;
c = dx/dt;
c_s = c/sqrt(3);
U0 = mu*Re/L0;
the distribution function code is:
void CMainFrame::Feq(double VU,double VV,double rho,double *ff)
{
int i;
double temp,tempp,cs;
cs=1./3.;
VU/=c;
VV/=c;
tempp=VUVU+VVVV;
for(i=0;i<9;i++)
{
temp=ei[i][0]*VU+ei[i][1]*VV;
if(i==4)
ff[i] = 4./9.rho * (1-3./2.tempp);
else
{
if(i%2==1)
ff[i] = 1./9.rho * (1+3.temp + 4.5temptemp-1.5tempp);
else
ff[i] = 1./36.rho * (1+3.temp + 4.5temptemp-1.5tempp);
}
}
}
Can someone tell me why
VU/=c;
VV/=c;
do need?