Shan-Chen multicomponent fluid

Hello, I am using the Shan-Chen model implemented in Palabos to model the flow of two fluids. I compared the numerical solution with the analytical one for the Poiseuille flow. I am not able to create the interface when the viscosities are different (tau1 = 1.0, tau2 = 3.0, for example) or when the densities are different (not greatly different). Therefore, the numerical solution is quite different from the analytical one. In the paper of Chin et al., 2002, the authors found that the Shan-Chen model still works well for viscosity ratio < 10. I am not able to obtain this result. Did any one find the same issue? Please let me know if you know how to solve this issue.

The model of Porter et al., 2012 is capable of simulate flows with higher density ratio and viscosity ratio. Did any one use this model and implement it in Palabos? If yes, please tell me how to implement this model.

I am looking forward to your answer.