LBM笔记.docx

Study notes for Lattice Boltzmann MethodD2Q9 lattice and velocities[1]The macroscopic fluid density:where f is the single-particle distribution function.The macroscopic velocity:The macroscopic pressure:Unit conversion[2]Single Relaxation Time BGK[1]Streaming and collision: is the streaming part is the collision termThe D2Q9 equilibrium distributionfunction is defined aswhere the weights are 4/9 for the rest particles (a = 0), 1/9 for a = 1, 2,3, 4, and 1/36 for a = 5, 6, 7, 8, and c is the basic speed on the lattice (1 in the simplest implementation). Note that if the macroscopic velocityu = 0, the equilibrium are simply the weights times the fluid density.The other form is[3]Whereand Coding for BGKmacroscopic fluid density and velocity:// Computing macroscopic density, rho, and velocity, u=(ux,uy).for( j=0; jLY; j++){for( i=0; iLX; i++){u_x[j][i] = 0.0;u_y[j][i] = 0.0;rho[j][i] = 0.0;if( !is_solid_node[j][i]){for( a=0; a9; a++){rho[j][i] += f[j][i][a];u_x[j][i] += ex[a]*f[j][i][a];u_y[j][i] += ey[a]*f[j][i][a];} u_x[j][i] /= rho[j][i];u_y[j][i] /= rho[j][i];}}}Streaming:// Streaming step.for( j=0; jLY; j++){jn = (j0 )?(j-1):(LY-1)jp = (jLY-1)?(j+1):(0 )for( i=0; iLX; i++){if( !is_interior_solid_node[j][i]){in = (i0 )?(i-1):(LX-1)ip = (iLX-1)?(i+1):(0 )ftemp[j ][i ][0] = f[j][i][0];ftemp[j ][ip][1] = f[j][i][1];ftemp[jp][i ][2] = f[j][i][2];ftemp[j ][in][3] = f[j][i][3];ftemp[jn][i ][4] = f[j][i][4];ftemp[jp][ip][5] = f[j][i][5];ftemp[jp][in][6] = f[j][i][6];ftemp[jn][in][7] = f[j][i][7];ftemp[jn][ip][8] = f[j][i][8];}}}Equilibrium Distribution Function:// Compute the equilibrium distribution function, feq.f1=3.;f2=9./2.;f3=3./2.;for( j=0; jLY; j++){for( i=0; iLX; i++){if( !is_solid_node[j][i]){rt0 = (4./9. )*rhoij;rt1 = (1./9. )*rhoij;rt2 = (1./36.)*rhoij;ueqxij = uxij;ueqyij = uyij;uxsq = ueqxij * ueqxij;uysq = ueqyij * ueqyij;uxuy5 = uxeqij + uyeqij;uxuy6 = -uxeqij + uyeqij;uxuy7 = -uxeqij + -uyeqij;uxuy8 = uxeqij + -uyeqij;usq = uxsq + uysq;feq