LBM example code in python has factor of 2 missing

just thought I’d point out that the LBM sample code in python: has a factor of two missing in the viscosity equation.
In the code:
nulb = uLB * r / Re
should be:
nulb = uLB * 2 * r / Re

Because the characteristic length for a cylinder is it’s diameter not it’s radius.

where is this code? (i can’t remember :-))

It’s at:
about half way down the page.

you are completely right. thanks. we’ll correct that.

No problem, I’ve had a lot of run with this code. I added tracer particles to it and made the following video: Maybe someday I’ll upgrade my computer and get Palabos up and running.

1 Like

Very nice! can you share how you make the particle tracing algorithm please

For the 2D simulation, the tracer particles are stored in a NumPy array of size p[N,3]

where N is the number of particles.

p[k,0] is the x position in lattice units for particle k.

p[k,1] is the y position in lattice units for particle k.

p[k,2] is the color map index for particle k.

The modified version of with tracer particles is here:

1 Like