Structure and rheological behavior of highly charged colloidal particles in a cylindrical pore I. Effect of pore size

In this work we performed nonequilibrium Brownian dynamics (NEBD) computer simulations of highly charged colloidal particles in diluted suspension under a parabolic flow in cylindrical pores. The influence of charged and neutral cylindrical pores on the structure and rheology of suspensions is analyzed. A shear-induced disorder-order-disorder-like transition was monitored for low shear rates and small pore diameters. We calculate the concentration profiles, axial distribution functions, and axial-angular pair correlation functions to determine the structural properties at steady state for a constant shear flow for different pore sizes and flow strengths. Similar behavior has been observed in a planar narrow channel in the case of charged interacting colloidal particles (M.A. Valdez, O. Manero, J. Colloid Interface Sci. 190 (1997) 81). The mobility of the particles in the radial direction decreases rapidly with the flow and becomes practically frozen. The flow exhibits non-Newtonian shear thinning behavior due to interparticle interactions and particle-wall interaction; the apparent viscosity is lower as the pore diameter decreases, giving rise to an apparent slip in the colloidal suspension. The calculated slip velocity was higher than that obtained in a rectangular slit under shear flow.