Rheology of viscoelastic suspensions of spheres under small and large amplitude oscillatory shear by numerical simulations

Abstract

The dynamic response of a viscoelastic suspension of spheres under small and large amplitude oscillatory shear is investigated by three-dimensional direct numerical simulations. A sliding triperiodic domain is implemented whereby the computational domain is regarded as the bulk of an infinite suspension. A fictitious domain method is used to manage the particle motion. After the stress field is computed, the bulk properties are recovered by an averaging procedure. The numerical method is validated by comparing the computed linear viscoelastic response of Newtonian and non-Newtonian suspensions with previous theories and simulations. The numerical predictions are in very good quantitative agreement with experimental data for the Newtonian case, whereas deviations are found with respect to some sets of experiments for semidilute and concentrated viscoelastic suspensions. To investigate on such discrepancies, the effect of aggregates in the bulk of the suspension is examined. The simulations show that the presence of structures significantly alters the loss modulus. Such an effect is more pronounced as the volume fraction increases. In this light, the above mentioned disagreement between simulations and data (and among experimental data themselves) can be rationalized, as its origin can be attributed to inhomogeneous particle configurations. For increasing strain amplitudes, both loss and storage moduli depart from the linear viscoelastic values. Although the deviations are qualitatively similar to the large amplitude response of the unfilled suspending matrix, our results for dilute and semidilute suspensions show that the decrease of the moduli is more and more pronounced as the volume fraction is higher. Furthermore, a higher concentration of solid particles reduces the value of strain amplitude such that the nonlinear behavior is observed. Simulations at higher frequencies also correctly capture the overshoot in the loss modulus for intermediate strain amplitudes. Finally, the effect of fluid elasticity on the particle motion is analyzed. The particles are found to move away from their starting positions and the average distance, computed at the beginning of each cycle with respect to the initial configuration, linearly increases with the number of cycles. The change in the microstructure is attributed to the long-range hydrodynamic interactions mediated by fluid viscoelasticit

Rheology of a dilute viscoelastic suspension of spheroids in unconfined shear flow

“…Non-Newtonian liquids, for instance, are responsible for complex particle dynamics (D'Avino and Maffettone 2015) that, in turn, alters the suspension rheology as compared to the Newtonian case (Mewis and Wagner 2012). The bulk rheology of suspensions of spheres in non-Newtonian media has been thoroughly investigated at various concentrations (Chan and Powell 1984;Aral and Kalyon 1997;Schaink et al 2000;Haleem and Nott 2009;D'Avino et al 2013). In the limit of dilute suspensions of spheres, the bulk rheology (viscosity, first and second normal stress differences) has been recently examined through theoretical analysis by assuming a weakly viscoelastic suspending medium (Greco et al 2005;Koch and Subramanian 2006;Greco et al 2007;Housiadas and Tanner 2009).…”

Analysis of dynamic oscillatory rheological properties of PP/EVA/organo-modified LDH ternary hybrids based on generalized Newtonian fluid and generalized linear viscoelastic approaches

“…There are various kinds of standard flows [1][2][3] to study the rheological properties of non-Newtonian fluids which all are similar in basics. In these methods, a material is subjected to defined stress or strain, and then, its response to the imposed stress or strain is measured.…”

Simulations of the start-up of shear flow of 2D particle suspensions in viscoelastic fluids: Structure formation and rheology

“…The situation gets even more difficult if the suspending fluid is viscoelastic, where effects such as particle alignment [10,11] and particle migration [12,13] might play a role in the rheology. Suspensions in viscoelastic fluids have been studied both experimentally [14,15] and computationally [16,17] (an extensive overview was given recently in [18]), yet a full understanding is lacking. For example in a highly elastic Boger fluid, typical particle dynamics observed in viscoelastic suspensions are absent [19].…”