Palabos: 关于rho和rhoBar的探讨

提问者Arpon表示

在AdvectionDiffusionDynamics代码里，他感到困惑，因为rhoBar是f’s之和，rho为rhoBar+1，他在文献里看BGK equilibrium项应该是这样的：
f_i^eq = t_i * rhoBar * (1 + (u*v_i)/c_s²)

然而代码里却是这样的：
f_i^eq = t_i * (rhoBar + (rho * u*v_i)/c_s²)

jonas1回复：

在Palabos的定义中f_i^bar = f_i-t_i，目的是为了精确度，同样的，Palabos定义的f_i^eq,bar = f_i^eq - t_i

所以通常的equilibrium如下：
f_i^eq = t_i * rho * (1 + (u*v_i)/c_s²)

代换rhoBar=rho-1则有
f_i^eq,bar = t_i * rho * (1 + (u * v_i)/c_s²) - t_i = t_i * (rhoBar + (rho * u *v_i)/c_s²)

我推导了一下的确是这样的。