Hello everyone,

I’m developing a code using the Shan & Chen model to perform multiphase simulations. I tried to evaluate the forces using the Momentum Exchange Algorithm (MAE) and I realized the populations (f0, f1, …f8) in my simulations are not the ones the LBM theory mentions.

I thought the populations (distribution functions) stood for the probability of molecules moving towards another node in the grid and the sum of them should be close to 1 (since this is the density). However, my results show very low values for the population (f_i values are between 0 and 0.1, but the sum never reaches 1). Besides, I also get negative values and I don’t know why.

Could someone help me out and clarify this to me. Thanks!

By the way, I get the populations using:

std::auto_ptr<MultiScalarField2D > population_0_light=computePopulation(lightFluid,Box2D(0, nx, 0, ny ),0);

More info about the function:

computePopulation(lattice, domain, iPop)

Extract the population iPop on each cell, and deposit the result in a scalar-field. The value iPop ranges from 0 to q-1.

(http://www.palabos.org/documentation/userguide/appendix-functions.html)