A DEM-based analysis of the influence of aggregate structure on suspension shear yield stress

This study theoretically examined the effect of aggregate structure on the suspension shear yield stress. The aggregation process of colloidal particles was simulated using the discrete element model (DEM) combined with the well-known DLVO theory. The predicted aggregate structural characteristics, namely the coordination number and inter-particle forces were then used in a modified version of the Flatt and Bowen mechanistic model [6] to calculate the corresponding suspension yield stress. The effect of key parameters such as solid volume fraction, suspension pH and ionic strength on the aggregate structure and hence the yield stress of the suspension was investigated. The results showed that the yield stress increased significantly under conditions that were favourable for formation of complex net-like aggregate structures, such as high solid volume fractions, pH values near the iso-electric point, and high ionic strengths. In such cases, the mean coordination number reached a maximum value which was considered to be dependent on the particle size and size distribution. The suspension yield stress exhibited a power law dependency on the solid volume fraction. The interconnected network structure developed at high solid volume fractions was found to be the major contributing factor to the observed high suspension yield stress. As the particle–particle repulsion became significant, a decrease in both the number of bonds and the mechanical bonding strength of the aggregate structure was observed. That was considered to be responsible for the reduction in the suspension yield stress. The suspension yield stress became independent of the suspension ionic strength when the ionic strength exceeded the critical coagulation concentration. Satisfactory agreements were obtained between simulation results and the published experimental data