We combine the Shan-Chen multicomponent lattice Boltzmann model and the link-based bounce-back particle suspension model to simulate particle motion in binary immiscible fluids. The impact of the slightly mixing nature of the Shan-Chen model and the fluid density variations near the solid surface caused by the fluid-solid interaction, on the particle motion in binary fluids is comprehensively studied. Our simulations show that existing models suffer significant fluid mass drift as the particle moves across nodes, and the obtained particle trajectories deviate away from the correct ones. A modified wetting model is then proposed to reduce the non-physical effects, and its effectiveness is validated by comparison with existing wetting models. Furthermore, the first-order refill method for the newly created lattice node combined with the new wetting model significantly improves mass conservation and accuracy.

A robust and efficient solver coupling computational fluid dynamics (CFD) with discrete element method (DEM) is developed to simulate particle-laden flows in various physical settings. An interpolation algorithm suitable for unstructured meshes is proposed to translate between mesh-based Eulerian fields and particle-based La-grangian quantities. The interpolation scheme reduces the mesh-dependence of the averaging and interpolation procedures. In addition, the fluid-particle interaction terms are treated semi-implicitly in this algorithm to improve stability and to maintain accuracy. Finally, it is demonstrated that sub-stepping is desirable for fluid-particle systems with small Stokes numbers. A momentum-conserving sub-stepping technique is introduced into the fluid-particle coupling procedure, so that problems with a wide range of time scales can be solved without resorting to excessively small time steps in the CFD solver. Several numerical examples are presented to demonstrate the capabilities of the solver and the merits of the algorithm.