Hi friends

in documement "choice of unit in LBM’ sec 2.4 mentioned that dt=2*10**(-4) and dx=1/100. in every text there is this formula: c=dx/dt. and we use c in equilbrium distribution: feq(i)=w(i)**rho*(1+3/(c**2)(e(i)*u)**2+… in this example : c=dx/dt=50!!! and we must set c=50 in program for feq(i)???

Thanks a lot.

Hello,

when implementing the LBM algorithm one usually takes dt=dx=1. With these “lb” units c=1/sqrt(3) (for the D2Q9 lattice).

People that are using dx/dt=c are just taking this 1/sqrt(3) are the reference velocity of the physical setup. In the case of the document you are mentioning the reference velocity is an arbitrary value U=0.02.

Thanks Malaspin

in other words we can say :notice the difference between dx,dt(physical) and dx,dt(lb) :c which is used in algorithm is c=dx(lb)/dt(lb) not in physical unit. and in lb unit always we use: dt=dx=1.

Is it right?

Hi guys,

So does this mean that in the algorithm and where c is used in calculation of Equilibrium functions, we should always use c=1 ??? I have seen a lot of literature that have adopted various values of c=dx/dt, some even as high as 500. I am still not sure to find the appropriate c=dx/dt that should be used in the feq formulas. Thanks in advance for your help.