L'ajout de particules aux liquides modifie les propriétés optiques et physiques des liquides, p. ex. la couleur, la masse volumique et la viscoélasticité. Les solvants, tels que l'eau ou l'acétone, sont généralement des liquides idéalement visqueux, tandis que les liquides contenant des particules affichent un comportement rhéologique plus complexe. Si les particules sont solubles dans le solvant, le produit final est appelé une solution. Les particules sont dissoutes sous forme de molécules ou d'ions (<1 nm). Dans le cas où les particules sont insolubles, le produit est un mélange biphasé (ou multiphasé) qui est appelé une suspension si les particules sont plus grandes que 1 µm ou appelé un colloïde si les particules mesurent entre 1 nm et 1 µm. De nombreux produits liquides comme les peintures, les encres, les boissons, les médicaments, les boues ou les gels douche contiennent divers types de particules pour obtenir le produit final souhaité ou ajuster les propriétés de transformation. Une faible fraction volumique de particules (appelée "fraction solide" par la suite) induit généralement un comportement de fluage, tandis que des concentrations élevées de particules peuvent entraîner un épaississement au cisaillement. Cependant, la viscoélasticité des matériaux dépend non seulement de la concentration des particules, mais aussi de la forme et de la taille des particules. Dans la plupart des suspensions, les particules diffèrent en taille, ce qui peut être exprimé comme la distribution de taille des particules. De plus, la charge électrique de la surface des particules (potentiel zêta) influence le comportement rhéologique de la suspension. Cet article vise à séparer les effets de la concentration des particules, de la forme des particules, de la distribution de la taille des particules et de la charge de surface des particules, et à nommer les méthodes courantes pour mesurer ces propriétés.
Influence des particules sur la rhéologie en suspension
Particules – la fraction solide dans les suspensions
Les particules dans un liquide agissent comme des obstacles, entravant l'écoulement du liquide et augmentant ainsi la résistance à l'écoulement, c'est-à-dire la viscosité. Un modèle simple pour décrire la viscosité accrue est d'assumer un champ d'écoulement autour d'une seule sphère, représentant des concentrations de particules très faibles, et une viscosité résultante accrue simplement calculée à partir de la viscosité du solvant et de la fraction volumique du solide dispersé (Figure 1, Équation 1, Einstein 1906, 1911). La formule a été modifiée et adaptée pour plusieurs sphères de même taille ("monodisperse") pour des fractions solides maximales de 0,15 à 0,2 (Équation 2) par Batchelor (1977), avec η représentant la viscosité en vrac de la suspension, η0 la viscosité du liquide, et φ la fraction solide dans la suspension.
Équation 1: $\eta = \eta_0(1+\frac{5}{2}\varphi)$ Équation 2: $\eta = \eta_0(1+2.5 \varphi + 6.2\varphi^2)$ Les suspensions avec une fraction solide encore plus élevée peuvent être décrites par le modèle de Krieger & Dougherty (1959, Équation 3), valable pour des taux de cisaillement faibles et élevés. Dans ce modèle, deux paramètres d'entrée supplémentaires sont requis. Φmax est la fraction solide maximale dans la suspension, permettant encore à la suspension de s'écouler. La viscosité intrinsèque [η] est une mesure de la forme des particules (2,5 pour les sphères). Équation 3 : $\eta = \eta_0(1- \frac{\varphi}{\varphi_{max}})^{-[\eta]*\varphi_{max}}$
En supposant des particules sphériques, la densité d'emballage maximale théorique peut être de 0,52 à 0,74, selon le schéma d'empilement des sphères (Figure 2). Il est bien établi que pour l'empilement aléatoire serré de sphères, la valeur est d'environ 0,64 (par exemple, Scott & Kilgour 1969). Pour de faibles fractions solides, la viscosité de la suspension est augmentée. Cependant, aucune interaction significative entre les particules n'a lieu et le comportement d'écoulement de la suspension est toujours contrôlé par le comportement d'écoulement newtonien de la phase solvant. Avec l'augmentation de la fraction solide, il est plus probable que les particules entrent en collision. Ils agissent maintenant comme des obstacles et une force de cisaillement supplémentaire est requise lorsque la friction se produit lors de la collision des particules (Figure 3). Lorsque la fraction solide approche φmax, les interactions entre particules en cours augmentent considérablement la force requise pour ciseler l'échantillon, dominant désormais sur la contrainte de cisaillement requise pour ciseler le solvant.
La réponse rhéologique de la suspension dépend maintenant également du taux de cisaillement, le début du comportement d'écoulement non newtonien, par exemple le caractère de cisaillement-diluant de la suspension, est très variable pour différentes suspensions mais peut commencer à φ = 0,1 à 0,5 (par exemple, Stickel & Powell 2005). De plus, à des taux de cisaillement plus élevés, les particules forment des agrégats et un blocage se produit, entraînant un épaississement par cisaillement (Figure 4). La transition de la fluidité à la viscosité est spécifique à l'échantillon et contrôlée par des facteurs tels que la taille (distribution) des particules et la forme des particules.
Taille des particules et distribution de la taille des particules
Pour une fraction solide donnée, une taille de particule plus petite entraîne un plus grand nombre de particules dans un volume donné (Figure 5). Par conséquent, la surface des particules peut être plusieurs ordres de grandeur plus élevée pour des particules plus petites. De très petites particules affichent un mouvement brownien agissant contre une force de cisaillement appliquée (Perrin 1910). Par conséquent, la charge de surface, l'adsorption et l'hydratation ont un effet potentiellement fort sur la taille effective des particules hydrodynamiques, en particulier pour les petites particules.
Nous avons déjà montré qu'à une certaine valeur de la fraction solide, lorsque l'interaction entre particules devient significative, un comportement de fluage est observé. La différence de viscosité entre les petites et grandes particules est diminuée à des taux de cisaillement plus élevés, car les particules sont réarrangées de manière favorable par rapport à la direction d'écoulement (Figure 6). Cependant, pour la majorité de toutes les suspensions existantes, la taille des particules n'est pas uniforme et est mieux décrite comme une distribution de taille des particules (PSD, Figure 7). La densité d'emballage maximale, et donc la fraction solide maximale, bénéficie de particules de taille variable ("polydispersité") car l'espace peut être rempli plus efficacement.
En ce qui concerne l'équation de Krieger & Dougherty (Équation 3), φmax augmente (il peut être supérieur à 0,74 !) et réduit ainsi la viscosité mesurée pour une concentration de particules donnée. Plus la distribution de la taille des particules (PSD) est large, plus l'effet positif sur la densité d'emballage est important (Figure 7). Nous avons déjà appris que des particules plus petites entraînent une viscosité plus élevée pour un taux de cisaillement et une fraction solide donnés. Si nous comparons maintenant la viscosité de deux PSD monodisperses et d'un PSD polydisperse, un minimum de viscosité est évident pour le mélange polydisperse (Figure 8). Cependant, de telles suspensions polydisperses peuvent encore avoir une contrainte de rendement élevée (Ancey & Jorrot 2001) et pour les suspensions colloïdales avec des distributions de taille de particules bimodales, les plus petites particules peuvent agir comme un lubrifiant (Ancey 2001).

Figure 8 : Pour une fraction solide constante et un taux de cisaillement constant, les petites sphères montrent des viscosités plus élevées par rapport aux grandes sphères (voir Figure 6). Les suspensions polydisperses ont un φmaxplus élevé. Pour une fraction solide constante (ici : 0,63), φ/φmax diminue et entraîne un minimum de viscosité pour de telles suspensions polydisperses.
Forme des particules et propriétés de surface
Jusqu'à présent, nous avons considéré des particules sphériques idéalisées. Évidemment, les particules affichent souvent une convexité inégale ou une forme non sphérique (allongement). La viscosité intrinsèque [η] est un facteur pour décrire la forme des particules. Par exemple, les sphères sont représentées par une valeur de 2,5 et les ellipsoïdes avec un rapport d'aspect de 5:1, 10:1, 50:1 par des valeurs de 2,94, 3,43, 7,00 (Figure 9a). Pour des rapports d'aspect (R) de 1:1 à 20:1, la viscosité intrinsèque peut être exprimée par la fonction suivante (Brenner 1974, Pabst et al. 2006, Équation 4) : Équation 4 : $[\eta] = 2.5 + 0.123(R-1)^{0.925}$ Les irrégularités de la surface des particules entraînent une viscosité plus élevée pour deux raisons : La déviation des lignes de flux du solvant s'écoulant autour des particules est plus forte que pour des sphères de même taille (Figure 9b), et l'augmentation de la surface spécifique des particules accroît la probabilité des interactions entre particules (Figure 9c). L'effet est plus fort pour des concentrations de particules élevées et des taux de cisaillement élevés.

Figure 9 : (a) Rapports d'aspect des ellipsoïdes réguliers, (b) champ d'écoulement autour des sphères par rapport au champ d'écoulement autour des particules irrégulières, (c) interactions entre particules des sphères et des particules irrégulières. Les particules irrégulières montrent plus de points de contact et donc une influence plus élevée de la friction inter-particulaire est attendue.
Si nous comparons maintenant la dépendance du taux de cisaillement des sphères et des particules allongées, nous récapitulons que les sphères montrent un cisaillement-diluant à mesure que les interactions entre particules deviennent moins prononcées à des taux de cisaillement élevés. L'augmentation de la surface des particules allongées entraîne une viscosité à faible cisaillement plus élevée par rapport aux sphères. Cependant, les particules non sphériques ont la possibilité d'adapter leur orientation à une direction d'écoulement. Au repos, ils sont orientés de manière aléatoire mais s'alignent dans la direction du flux lorsqu'un taux de cisaillement est appliqué. Cela conduit à une viscosité inférieure par rapport aux sphères à des taux de cisaillement plus élevés car la viscosité intrinsèque est alors plus faible – résultant de la densité d'emballage maximale φmax (Figure 10).

Figure 10 : Viscosité dépendante du taux de cisaillement des sphères par rapport aux ellipsoïdes. Les deux montrent un amincissement par cisaillement. Cependant, les suspensions avec des ellipsoïdes ont des viscosités plus élevées à faibles taux de cisaillement et les suspensions avec des sphères ont des viscosités plus élevées à des taux de cisaillement plus élevés.
L'effet de la charge des particules électriques
Selon l'application de la suspension, une grande variété de substances est utilisée à la fois pour les solvants et les particules, pouvant éventuellement avoir une charge électrique. Le moment électrostatique, qui peut représenter l'attraction ou la répulsion au sein d'une suspension, est quantifié comme le potentiel zêta. Si les particules ont une charge de surface positive ou négative, des molécules ou des ions avec une charge négative ou positive occupent la surface (Figure 11) et forment ensemble une particule avec un potentiel de surface non nul.

Figure 11 : Différents types d'interactions particule-solvant pour une particule avec une charge de surface négative : (a) pas d'adsorption (potentiel zêta positif élevé), (b) adsorption et formation de la couche de Stern (potentiel zêta positif faible) et (c) couche de Stern et adsorption chimique supplémentaire (potentiel zêta négatif faible)
Ce potentiel électrique affecte l'interaction avec le solvant ainsi que le comportement inter-particulaire. Si la charge de surface est élevée, la répulsion entre les particules domine sur l'attraction par les forces de van-der-Waals (Figure 12a) et une barrière énergétique empêche la proximité des particules. En revanche, pour les particules avec une charge de surface presque nulle, la répulsion est plus faible et un minimum d'énergie pour une certaine distance entre particules facilite un agencement ordonné des particules (Figure 12b).

Figure 12 : Interaction particule-particule en fonction de la distance inter-particulaire. (a) Pour un potentiel zêta élevé, la répulsion augmente avec la diminution de la distance inter-particulaire jusqu'à un maximum, représenté par une barrière énergétique. (b) Pour un potentiel zêta faible, l'agglomération est favorisée par un minimum énergétique à des distances inter-particulaires relativement petites.
En conséquence, pour un échantillon avec un faible potentiel zêta, l'accumulation ou la floculation des particules est favorisée, tandis qu'à un potentiel zêta élevé, les particules peuvent se repousser et éviter efficacement l'agglomération. Cependant, pour comprendre le comportement rhéologique des suspensions, nous devons également tenir compte de la taille des particules et de la fraction solide. Les dispersions avec de petites particules (<1 µm) sont particulièrement sensibles aux changements de potentiel zêta, car la répulsion des particules induite et l'augmentation du volume hydrodynamique (Maranzano & Wagner 2001, Figure 13) augmenteront la viscosité, en particulier à faibles taux de cisaillement. Ceci n'est possible que parce que les forces gravitationnelles agissent de manière subordonnée à de si petites et légères particules, tandis que le mouvement brownien dominera. En supposant des fractions solides relativement faibles, ce type de dispersion maintiendra un plateau de viscosité à cisaillement nul.
En revanche, pour les particules plus grandes, l'effet de sédimentation est généralement plus fort et, avec le potentiel zêta, il contrôlera le comportement d'agglomération. À un faible potentiel zêta, les particules s'agglomèrent et en augmentant progressivement la taille de l'agglomérat, la gravitation finira par provoquer la sédimentation. Une telle agglomération peut induire un point de rendement si la fraction solide est suffisamment élevée et que les particules sont capables de former un réseau tridimensionnel. Notez le comportement contrasté à des taux de cisaillement élevés, où un épaississement par cisaillement peut se produire, car la répulsion inter-particulaire entraîne des viscosités plus élevées pour les échantillons avec un potentiel zêta élevé (Figure 14).
Méthodes pour caractériser les suspensions
Les propriétés physiques en vrac de la suspension et du solvant, comme la viscosité et le point de rendement, peuvent être mesurées avec des viscosimètres ou des rhéomètres lors de tests rotationnels ou oscillatoires. La caractérisation de la fraction solide nécessite la connaissance de la taille des particules, de la distribution de la taille des particules (DSP) et de la forme des particules. Ces propriétés nécessitent des méthodes supplémentaires telles que la diffusion dynamique de la lumière ou la diffraction laser (pour la taille des particules et la PSD) et l'imagerie haute résolution (également pour la forme des particules) à mesurer. De plus, la diffusion de la lumière électrophorétique est utilisée pour déterminer le potentiel zêta. La diffraction laser est capable d'analyser des particules allant de la sous-µm à la taille mm. Le principe de mesure est basé sur l'orientation d'un laser vers l'échantillon et plusieurs détecteurs pour enregistrer la diffusion à différents angles. Le comportement de diffusion de la lumière diffère en fonction de la taille des particules : les particules plus petites diffusent plus fortement que les particules plus grandes. Cependant, la mesure est basée sur l'hypothèse de particules sphériques et la connaissance des propriétés optiques du solvant ainsi que des particules est requise. La diffusion dynamique de la lumière permet de mesurer la taille des particules jusqu'à la plage sub-nm. Le principe de mesure utilise des fluctuations d'intensité sur des échelles de temps très courtes dues au mouvement brownien. Ainsi, la vitesse des particules peut être corrélée à la taille des particules. Pour la diffusion de la lumière par électrophorèse, deux électrodes sont connectées à l'échantillon et un champ électrique est appliqué. Les particules chargées se déplaceront maintenant vers l'électrode positive ou négative, la vitesse est mesurée avec la technique Doppler laser comme un décalage de fréquence ou comme un décalage de phase. Les valeurs obtenues peuvent être converties en distribution de potentiel zêta ou en un potentiel zêta moyen, respectivement. L'imagerie automatisée permet de déterminer la forme et la taille des particules si les particules ont un diamètre d'au moins 1 µm. Des milliers d'images sont prises pour garantir des valeurs statistiquement solides. Cette méthode est particulièrement utile si les particules montrent une forte déviation par rapport à la sphéricité. À l'exception de la diffusion dynamique de la lumière (uniquement applicable aux échantillons humides), ces méthodes peuvent être utilisées pour des échantillons à la fois secs et humides.
Conclusion
De nombreux produits liquides contiennent une certaine fraction de particules et la compréhension du comportement viscoélastique de ces suspensions est d'une importance vitale pour un large éventail de produits et de processus industriels, naturels et biologiques. La fraction de particules, la forme, la taille (distribution) et la charge de surface influencent le comportement rhéologique d'une suspension, qui est typiquement thixotrope mais peut également être newtonien pour de faibles fractions de particules, ou est cisaillement-épaississement pour de fortes fractions de particules à des taux de cisaillement élevés. Une combinaison de méthodes pour déterminer les propriétés rhéologiques et les propriétés des particules est nécessaire pour une caractérisation complète des suspensions.
Références
Ancey C (2001) Rôle des contacts lubrifiés dans des suspensions polydisperses concentrées. Journal de Rhéologie 45 : 1421–1439. doi : 10.1122/1.1413504 Ancey C, Jorrot H (2001) Tension de rupture pour les suspensions de particules dans une dispersion d'argile. Journal de rhéologie 45 : 297–319. doi : 10.1122/1.1343879 Batchelor KG (1977) L'effet du mouvement brownien sur la contrainte de masse dans une suspension de particules sphériques. Journal de Mécanique des Fluides 83 : 97–117. doi : 10.1017/S002211207700106 Brenner H (1974) Rhéologie d'une suspension diluée de particules brunes axisymétriques. Revue internationale d'écoulement multiphasique 1 : 195–341. doi : 10.1016/0301-9322(74)90018-4 Einstein A (1906) Une nouvelle détermination des dimensions moléculaires. Annalen der Physik 324: 230–247. doi: 10.1002/andp.19063240204 Einstein A (1911) Correction de mon travail : Une nouvelle détermination des dimensions moléculaires. Annalen der Physik 339: 591–592. doi: 10.1002/andp.19113390313 Krieger IM, Dougherty TJ (1959) Un mécanisme pour l'écoulement non newtonien dans des suspensions de sphères rigides. Transactions de la Société de Rhéologie 3 : 137–152. doi : 10.1122/1.548848 Maranzano BJ, Wagner NJ (2001) Les effets des interactions interparticulaires et de la taille des particules sur l'épaississement réversible par cisaillement : dispersions colloïdales à sphères dures. Journal de Rhéologie 45 : 1205–1222. doi : 10.1122/1.1392295 Pabst W, Gregorova E, Berthold C (2006) Forme des particules et rhéologie des suspensions des systèmes à fibres courtes. Journal de la Société Européenne de Céramique 26 : 149–160. doi : 10.1016/j.jeurceramsoc.2004.10.016 Perrin J (1910) Mouvement brownien et réalité moléculaire. Taylor & Francis, 93 p. ISBN : 9780486174723 Scott GD, Kilgour DM (1969) La densité de l'emballage aléatoire serré de sphères. Journal de la physique D : physique appliquée 2 : 863. doi : 10.1088/0022-3727/2/6/311 Stickel JJ, Powell RL (2005) Mécanique des fluides et rhéologie des suspensions denses. Revue annuelle de mécanique des fluides 37 : 129–149. doi : 10.1146/annurev.fluid.36.050802.122132