Colloids are ubiquitous in the food, medical, cosmetic, polymer, water purification and pharmaceutical industries. Colloids thermal, mechanical and storage properties are highly dependent on their interface morphology and their rheological behavior.

Numerical methods provide a cheap and reliable virtual laboratory for the study of colloids. However efficiency is a major concern to address when using numerical methods for practical applications.

This work introduces the main building-blocks for an improved lattice Boltzmann-based numerical tool geared towards the study of colloidal rheology and interface morphology.

The efficiency of the proposed model is enhanced by using recently developed and validated migrating multi-block algorithms for the lattice Boltzmann method (LBM).The migrating multi-block was used to simulate single component, multi-component, multiphase and single component multiphase flows. Results were validated by

experimental, numerical and analytical solutions.

The contamination of the fluid-fluid interface influences the colloids morphology. This issue was addressed by the introduction of the hybrid LBM for surfactant-covered droplets. The module was used for the simulation of surfactant-covered droplet deformation under shear and uniaxial extensional flows respectively and under

buoyancy. Validation with experimental and theoretical results was provided.

Colloids are non-Newtonian fluids which exhibit rich rheological behavior. The suppression of coalescence module is the part of the proposed model which facilitates the study of colloids rheology. The model results for the relative viscosity were in agreement with some analytical results.

Biological suspensions such as blood are micro-colloids by nature. The study of the blood flow in the microvasculature was heuristically approached by assuming the red blood cells as surfactant covered droplets. The effects of interfacial tension on the flow velocity and the droplet exclusion from the walls in parabolic flows were in qualitative agreement with some experimental and numerical results. The Fahraeus and the Fahraeus-Lindqvist effects were reproduced.

The proposed modules could be used separately or in combination for the study of a variety of colloids and biological suspensions problems as this was demonstrated throughout this work.