Dear Community
In this computeEquilibrium , the first order is used to calculate the equilibrium state distribution function,However, the bgkCollision called by the calculation of uSqr is a second-order equilibrium distribution function.The code below is in ParticleAdvectionDiffusionBGKdynamics. thank you
{ V uSqr = lbm<DESCRIPTOR>::bgkCollision(cell, rho, adjusted_u, omega);
return {rho, uSqr};
};
T computeEquilibrium(int iPop, T rho, const T adjusted_u[DESCRIPTOR::d]) const override
{
return equilibrium<DESCRIPTOR>::template firstOrder(iPop, rho, adjusted_u);
};