Published online by Cambridge University Press: 18 October 2019
We study the rheology of a non-colloidal suspension of rigid spherical particles interacting with a viscoelastic matrix. Three-dimensional numerical simulations under shear flow are performed using the smoothed particle hydrodynamics method and compared with experimental data available in the literature using different constant-viscosity elastic Boger fluids. The rheological properties of the Boger matrices are matched in simulation under viscometric flow conditions. Suspension rheology under dilute to semi-concentrated conditions (i.e. up to solid volume fraction $\unicode[STIX]{x1D719}=0.3$) is explored. It is found that at small Deborah numbers
$De$ (based on the macroscopic imposed shear rate), relative suspension viscosities
$\unicode[STIX]{x1D702}_{r}$ exhibit a plateau at every concentration investigated. By increasing
$De$, shear thickening is observed, which is related to the extensional thickening of the underlying viscoelastic matrix. Under dilute conditions (
$\unicode[STIX]{x1D719}=0.05$), numerical results for
$\unicode[STIX]{x1D702}_{r}$ agree quantitatively with experimental data in both the
$De$ and
$\unicode[STIX]{x1D719}$ dependences. Even under dilute conditions, simulations of full many-particle systems with no a priori specification of their spatial distribution need to be considered to recover precisely experimental values. By increasing the solid volume fraction towards
$\unicode[STIX]{x1D719}=0.3$, despite the fact that the trend is well captured, the agreement remains qualitative with discrepancies arising in the absolute values of
$\unicode[STIX]{x1D702}_{r}$ obtained from simulations and experiments but also with large deviations existing among different experiments. With regard to the specific mechanism of elastic thickening, the microstructural analysis shows that elastic thickening correlates well with the average viscoelastic dissipation function
$\unicode[STIX]{x1D703}^{elast}$, requiring a scaling as
$\langle \unicode[STIX]{x1D703}^{elast}\rangle \sim De^{\unicode[STIX]{x1D6FC}}$ with
$\unicode[STIX]{x1D6FC}\geqslant 2$ to take place. Locally, despite the fact that regions of large polymer stretching (and viscoelastic dissipation) can occur everywhere in the domain, flow regions uniquely responsible for the elastic thickening are well correlated to areas with significant extensional component.
To send this article to your Kindle, first ensure no-reply@cambridge.org is added to your Approved Personal Document E-mail List under your Personal Document Settings on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part of your Kindle email address below. Find out more about sending to your Kindle. Find out more about saving to your Kindle.
Note you can select to save to either the @free.kindle.com or @kindle.com variations. ‘@free.kindle.com’ emails are free but can only be saved to your device when it is connected to wi-fi. ‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.
Find out more about the Kindle Personal Document Service.
To save this article to your Dropbox account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you used this feature, you will be asked to authorise Cambridge Core to connect with your Dropbox account. Find out more about saving content to Dropbox.
To save this article to your Google Drive account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you used this feature, you will be asked to authorise Cambridge Core to connect with your Google Drive account. Find out more about saving content to Google Drive.