Reply To: Not clear which version of Shan Chen model is implemented
OpenLB – Open Source Lattice Boltzmann Code › Forums › on OpenLB › Bug Reports › Not clear which version of Shan Chen model is implemented › Reply To: Not clear which version of Shan Chen model is implemented
that is another potential resolution. I had assumed that, given that Guo forcing was chosen, implementing something along the lines of https://journals.aps.org/pre/abstract/10.1103/PhysRevE.86.036701 had been the intention. To do that correctly though, in accordance with formulas (15) and (16) in that paper, the force/2 addition needs to occur before the averaging step.
In any case though, if there was no intention to implement that model, it would likely make sense to simply switch to a ShanChen forcing scheme to avoid future confusions, which would bring the model in line with the original Shan Chen paper.