contact angle treatment in Lee and Lin (2005) model!


I have been trying to implement Briant and Yeomans (2004) in Lee and Lin (2005) multicomponent multiphase (MCMP) model to simulate proper contact angles during flow of immiscible fluids in porous media.

I wonder if anyone has gone this far using Lee and Lin (2005) model. I would appreciate any help in this context.

T. Lee and C.-L. Lin, A stable discretization of the lattice Boltzmann equation for simulation of incompressible two-phase
flows at high density ratio, J. Comput. Phys. 206, 16 (2005)

Briant, A.J., Yeomans, J.M., 2004. Lattice Boltzmann simulations of contact line motion: II. Binary fluids. Phys. Rev. E 69, 031603.