We study the overall dissipation rate of highly concentrated non-colloidal suspensions of rigid neutrally buoyant particles in a Newtonian fluid. This suspension is confined to a finite size container, subject to shear or extensional boundary conditions at the walls of the container. The corresponding dissipation rates determine the effective shear viscosity μ* and the extensional effective viscosity λ*. We use recently developed discrete network approximation techniques to obtain discrete forms for the overall dissipation rates, and analyse their asymptotics in the limit when the characteristic interparticle distance goes to zero. The focus is on the finite size and particle wall effects in spatially disordered arrays. Use of the network approximation allows us to study the dependence of μ* and λ* on variable distances between neighbouring particles in such arrays.
Our analysis, carried out for a two-dimensional model, can be characterized as global because it goes beyond the local analysis of flow in a single gap between two particles and takes into account hydrodynamic interactions in the entire particle array. The principal conclusion in the paper is that, in general, asymptotic formulae for μ* and λ* obtained by global analysis are different from the formulae obtained from local analysis. In particular, we show that the leading term in the asymptotics of μ* is of lower order than suggested by the local analysis (weak blow-up), while the order of the leading term in the asymptotics of λ* depends on the geometry of the particle array (either weak or strong blow-up). We obtain geometric conditions on a random particle array under which the asymptotic order of λ* coincides with the order of the local dissipation in a gap between two neighbouring particles, and show that these conditions are generic. We also provide an example of a uniformly closely packed particle array for which the leading term in the asymptotics of λ* degenerates (weak blow-up).