Mohan
Das
a,
Lucille
Chambon
a,
Zsigmond
Varga
b,
Maria
Vamvakaki
a,
James W.
Swan
b and
George
Petekidis
*a
aIESL – FORTH and Department of Material Science and Technology, University of Crete, GR – 71110, Heraklion, Greece. E-mail: georgp@iesl.forth.gr
bDepartment of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
First published on 13th November 2020
A combination of rheology, optical microscopy and computer simulations was used to investigate the microstructural changes of a semi-dilute suspension of attractive rigid rods in an imposed shear flow. The aim is to understand the relation of the microstructure with the viscoelastic response, and the yielding and flow behaviour in different shear regimes of gels built from rodlike colloids. A semi-dilute suspension of micron sized, rodlike silica particles suspended in 11 M CsCl salt solution was used as a model system for attractive rods’ gel. Upon application of steady shear the gel microstructure rearranges in different states and exhibits flow instabilities depending on shear rate, attraction strength, volume fraction and geometrical confinement. At low rod volume fractions, the suspension forms large, vorticity aligned, particle rich flocs that roll in the flow-vorticity plane, an effect that is due to an interplay between hydrodynamic interactions and geometrical confinement as suggested by computer simulations. Experimental data allow the creation of a state diagram, as a function of volume fraction and shear rates, identifying regimes of stable (or unstable) floc formation and of homogeneous gel or broken clusters. The transition is related to dimensionless Mason number, defined as the ratio of shear forces to interparticle attractive force.
Colloidal gels formed by aggregation of colloidal particles when their attraction strength U, is larger than thermal energy kBT exhibit time dependent phenomena such as ageing25–27 and thixotropy28 and possess solid-like response and elasticity.29,30 At low particle volume fraction, colloidal gels subjected to external stresses can undergo shear driven aggregation,31,32 and other shear induced structural changes in the flow-vorticity plane,33,34 breakage of clusters,35,36 exhibit delayed yielding,37,38 collapse of the particle network and sedimentation.39,40 Due to their thixotropy, shear history can be used as a way to tune viscoelastic response of colloidal gels by controlling the structural heterogeneity by varying pre-shear rate and strain amplitude in steady shear or oscillatory shear respectively.41,42 These properties render model colloidal gels valuable in simulating more complex industrial formulations and processes.
Low particle volume fraction colloidal gels exposed to shear flow in a confined space can exhibit large heterogeneities in the form of anisotropic, particle-rich flocs shaped like “logs” rolling with axis perpendicular to the flow-gradient plane. This so-called log rolling phenomenon has been observed experimentally in a variety of colloidal systems such as carbon black particles, cellulose, carbon nanotubes, silica and alumina particles as well as in attractive emulsions.43–52 Understanding of this phenomenon is critical in applications where attractive colloidal suspensions flow through narrow channels such as a slurry of carbon black or graphene particles used in capacitive energy storage systems.53–55 Initial observations of vorticity aligned flocs were attributed to “hoop stresses” arising from elastic instabilities due to the soft viscoelastic domains within a fluid of lower viscoelasticity. Such stresses were caused by normal forces exerted by a growing particle floc in the flow-gradient plane that is restricted within the rheometer gap. This confinement eventually results in the growth of the particle floc in the flow-vorticity plane.43,44 However, it is not clear how particle flocs remain stable while undergoing shear flow and what is the relation between the particle concentration, attraction strength and gap in the flow geometry, remaining a subject of interest.52 Recently Brownian dynamics simulations showed that this phenomenon arises due to a coupling between hydrodynamic interactions (HI) and the confining boundaries for a system of attractive particles,56 whereas instead in the absence of HI particles arrange themselves into sheets. The dimensionless ratio of hydrodynamic drag force due to shear to the interparticle attraction force (known as Mason number) can be used to define the boundaries of a stability diagram that relates shear rates and gaps necessary to observe stable vorticity aligned flocs. This response is also demonstrated in recent experiments with attractive spherical particles of different sizes and attraction strength.57
Dispersions of attractive rod-shaped colloids with high aspect ratio have a lower percolation threshold than spherical colloids and possess the additional control parameter of orientation which depends on the aspect ratio.58 Hard, attractive rods of length L and diameter D, exhibit Diffusion Limited Cluster Aggregation (DLCA) at significantly lower particle volume fractions compared to spheres and can form stable gels in the semi-dilute regime (Vrod/L3 < φ < Vrod/L2D).58–63 Furthermore, at high volume fractions (φ ≫ Vrod/L2D), rod shaped colloids can form glasses and liquid crystalline phases with short or long range orientational as well as positional order following Onsager's prediction of structural phase transitions.13,64–67 These properties make rod-shaped colloids useful in the field of advanced materials and for various optical applications.68–70 Considering the rich microstructural and rheological information available for colloidal gels made of spherical particles that have considerable scientific and practical implications, it is of interest to explore the dynamics and structure of colloidal gels made up of shape anisotropic particles such as rods.
Microscopic silica rods have been used as model systems to study the dynamics and phase behaviour of rod-like colloidal dispersions.71 Due to their relatively low dispersity in size, high aspect ratio and spherocylindrical shape, they exhibit all the phases associated with hard, rod-like systems.72 Furthermore, due to their microscopic size, they are useful in quantitative 3D real-space studies in concentrated phases at a single particle level73 not easy with other model rod-like particles such as FD-virus74 where scattering techniques are more widely used.75
In this paper we study the effects of external shear (imposing different shear rates) on the structure and the subsequent mechanical response of an attractive rod suspension consisting of micron sized silica rods dispersed in a concentrated salt solution. We mainly focus on the investigation of the conditions and mechanisms related with the formation of vorticity aligned clusters, as a consequence of shear melting (yielding) of rodlike colloidal gels in confined geometries. To this end rods of intermediate aspect ratio (= 10) are utilized while effort has been made to experimentally probe structures at different length-scales in order to most efficiently and comprehensively reveal any shear induced structural formations. For this we use a combination of rheometry and optical techniques such as dark-field microscopy and fluorescent confocal microscopy. The mechanisms of yielding and flow of the attractive rod suspension and the structural relaxation after shear cessation are explored through these methods. Such attractive rod suspensions exhibit flow instabilities at low shear rates arising from a competition between attraction forces and shear under confinement leading to the formation of large heterogeneities and vorticity aligned “log-rolling” clusters. Based on these quantities a state diagram was created describing how rod volume fraction and Mason number affect the evolution of the microstructure in shear flow.
The rods have an average length, L = 4.36 μm with a length polydispersity of 18% and an average aspect ratio of 10 (Fig. S1, ESI†). We assume the rods to be cylinders with spherical end caps and calculate the volume of the rod, Vp = (π/4)LD2 + (π/6)D3. From this we calculate the theoretical overlap volume fraction, φ* (= Vp/L3) which was found to be 0.008. The density of the rods was found to be 1.98 g cm−3 (Fig. S2, ESI†).
(1) |
We should also note that the sedimentation time over the length (≈5 μm) for a single rod parallel to its major axis in dilute solution is approximately 100 s, calculated based on the parallel friction coefficient for a rod,78
We also use a rotational Peclet number, Per (= /Dr) to non-dimensionalize the applied shear rates, with the rotational diffusion coefficient of a rod defined by,84
A second set of rheo-optical measurements (Type A) were performed using a glass plate–plate geometry (diameter: 40 mm) at different gaps.
In a given simulation Nrods rod like particles are modeled by connecting Nbeads spherical domains together in series with stiff and torsionally resistant springs to represent an individual stiff rod. The beads have mean particle radius, a, with an imposed polydispersity of σ = 5% to account for the slight roughness observed in the experimental rod system. Neighboring beads in a rod just touch at the extent of their hydrodynamic radius. The rods are made attractive through a short-ranged inter-particle attraction between beads of different rods and the beads exert a hard sphere-like repulsion when their hydrodynamic radii touch. This attraction is modeled using an Asakura–Oosawa form87 of width δ/a = 0.2 and characteristic depth U. For short-ranged attractions, of this form the potential is essentially linear between the point of inter-bead contact and the finite range of the attractive well. For the dispersions in the simulations we set the thermal energy scale to kBT = 0 as previous studies have observed log-rolling for both thermal and athermal systems. The strength of shear is varied by adjusting the strain rate and the dispersion is strained for up to γ = 300 strain units starting from a randomly dispersed configuration. As the system is athermal, the Mason number, Mn is used to control the shear force exerted on the particles. The spring stiffness holding beads together in a rod is always at least ten times higher than the shear strength to avoid any bending of the rods under flow. In the steady shear experiments, the flow Peclet number is larger than one, while the attractive forces are many time stronger than the thermal energy kBT (≈10kBT) so that the athermal approximation is a good one.
The dispersion is sheared employing Lees–Edwards boundary conditions in the flow–flow gradient direction.88 The aspect ratio of the simulated box in the gradient-vorticity direction is kept constant at 1 and the flow-gradient aspect ratio Axy is varied between 2.5 and 7.5 in order to resolve multiple, periodically spaced flocs. Confining walls along the flow-vorticity plane are represented by a hexagonal close packed lattice of particles connected to each other through rigid bonds. The walls interact hydrodynamically with particles of the colloidal dispersion and a stiff repulsion at contact identical to that between rods. This approach is a rather simplistic representation of confining walls, yet it effectively captures many features of the complex hydrodynamic interactions near walls.89 There are three layers of particles in each wall to ensure sufficient bending stiffness to resist buckling under shear. Two walls in the flow-vorticity plane define the boundaries of the sheared cell and their separation d the height of the confining gap. The specified volume fraction in simulations is given in terms of the total volume of beads relative to the sheared cell.
A similar picture is revealed by the computer simulations (Fig. 2(e)) performed at a rod volume fraction φsim/φ* = 1.25 and 10kBT attraction strength per bead. Simulations are reproducing qualitatively the experimental picture, with the main finding that the rods structure so that they do not exhibit a high degree of alignment. The clusters of particles created by Diffusion Limited Aggregation (DLA) are quite fluffy and easily form a space spanning network.
Below we present the mechanical properties of such structures formed at quiescent conditions as well as the effect of steady shear rate in the formation of new shear induced metastable states.
The measurement time in the DTS measurement was normalized by the rotational relaxation time, τR (= 1/6Dr) and the shear modulus was scaled by the thermal energy density, kBT/(rod volume). All the suspensions exhibit ageing after rejuvenation (Fig. 3(a)). In the case of suspensions at φ ≤ φ*, a clear liquid to solid like transition was observed over time. Since the solvent is essentially Newtonian, the time evolution of the viscoelasticity (ageing) can be associated purely with the formation of a percolated structure aided by diffusion of the attractive rods. As more and more rods form connections over time, the structure evolves after shear cessation. At higher rod volume fractions this behaviour is weaker as rejuvenation does not result in significant change in the suspension microstructure and void volume fraction (or free volume between clusters) and hence the structure reforms faster and reaches a quasi-steady state in a shorter time.41 This picture is revealed during the DTS tests using three dimensional confocal microscope imaging where individual rod aggregates were observed to percolate over time (Fig. 3(b)) (check Video, ESI†). Therefore, the combined linear rheology and confocal microscopy suggests that as shear rejuvenation breaks down the structure only into smaller aggregates at low φ's these re-percolate slower after shear cessation. However, at higher φ's, due to the dense nature of the suspension, the initial re-percolation process is much quicker and is followed by a slower strengthening of the structure as a result of more local rearrangements and contact-driven ageing that was observed in systems with similar interactions.90
Frequency sweep measurements show that both G′ and G′′ are nearly independent of frequency, ω. The frequency of the measurement was rescaled with τR (as ω·τR) indicating we are accessing time scales below the single rod rotational relaxation time (Fig. 4). The measurement suggests that there are no relaxation processes associated with structural changes or rod rearrangements, at different time and length scales. This is consistent with the formation of a strong gel even at very low rod volume fractions.
Fig. 4 Dynamic frequency sweep measurements of the rod suspension at different rod volume fractions. |
Amplitude sweep measurements performed at an angular frequency of 1 rad s−1 show that the rod gel network yields fully beyond a strain amplitude defined conventionally at the crossover of G′ and G′′. This yield point (G′(γy) = G′′(γy)) which is volume fraction dependent, is observed at strain amplitudes lower than 10% for all the rod volume fractions (Fig. 5), with the yield strain, γy decreasing with increasing rod volume fraction. Moreover, at low rod volume fractions the shear modulus, G′, shows a clear dependence on strain amplitude even before yielding, while for the higher rod volume fractions; a smaller yield strain indicates the formation of a more brittle gel. Studies on yielding in colloid polymer gels under oscillatory shear report that at high strain amplitudes the gel structure is broken and liquefied. However, at intermediate strain amplitudes, the gel structural heterogeneity increases, which eventually upon shear cessation results in a more heterogeneous and weaker gel.42 Note that the volume fraction dependence of G′ exhibits a power law slope which is almost half to that found for spherical particle gel as found by Shih et al.,29 while a qualitative agreement is also found in the decrease of the yield strain with φ as the system becomes more concentrated and more brittle (see Fig. S5, ESI†).
The investigation using confocal microscopy for the rod suspension undergoing oscillatory shear flow (Fig. S4, ESI†) shows that for low volume fractions (φ/φ* = 1) at low strain amplitudes (γ < γy), the percolated rod network undergoes reversible bending deformation, where clusters and gel strands bend back and forth, while internally individual rods rearrange reversibly under oscillatory shear (see confocal microscopy Video, ESI†). At large strain amplitudes (γ > γy) however, the percolated structure breaks and new larger and denser aggregates (as appear from microscopy images) are formed as a result of new bond formation between smaller neighbouring rod aggregates. In this way free volume (or void volume) between clusters increases. In the case of high rod volume fraction (φ/φ* = 5), at low strain amplitudes it was not possible to identify such local particle rearrangements. However, at high strain amplitudes individual rods within the aggregates were initially observed to undergo local angular displacements and at extremely large strain amplitudes, the gel structure gets ruptured (see Video, ESI†). Moreover, even at these large strain amplitudes, there is no flow-orientation of the rods.
The bending elasticity in colloidal aggregates formed by screening Coulomb interactions has been demonstrated in detail for spherical particles.90–93 Here however, in addition to contact bending stiffness witnessed in interparticle bonds formed by van der Waals attractions, particle orientation also plays a significant role in determining the bending (or flexure) modulus of the percolated structure and consequently in yielding of the rod gel.
The yielding mechanism can be contrasted with that reported by Mohraz and Solomon91 on gels of both spherical and rod-like colloids. For the higher volume fractions measured here, i.e. 2.5 and 5φ* (φ = 0.04 and 0.08), the structure of the gel network is such that there are no soft points on the backbone allowing angular rotation due to crowding by neighboring rods resulting in increased contact points per rod. The presence of multiple contact points per rod within the gel structure imparts a high bending elasticity of the gel network resulting in a brittle structure with a low yield strain.
At lower volume fractions, near the critical gelation point, the mechanism of yielding involves bending of the flexible gel backbone before rupturing. This flexibility arises mainly due to lower restrictions on individual rod motion in a less crowded environment allowing rotation of individual rods around specific soft points (as shown in the confocal Video in ESI†). Thus similarities with spherical particles DLCA gels maybe identified where relative angular rotation of a singly connected pair is allowed leading to unwinding of the backbone of a fractal cluster before rupturing.
A similar start-up shear response is found at a volume fraction of 2.5φ*. However, at 5φ* this response was different (Fig. S5, ESI†). From macroscopic observations it was established that at this volume fraction there was no formation of log-rolling flocs (Fig. 8). These findings indicate the role of particle volume fraction in creating such heterogeneities. Similar volume fraction dependence of this phenomenon has been reported in other attractive systems of spherical colloids.49,50,52 Furthermore, we should keep in mind that the floc formation is initiated at the center of the cone–plate geometry where the gap is smaller and rods could experience confinement for the specific rod size.
To this end we used parallel plate geometry to quantify the effect of confinement on the formation of vorticity aligned flocs in a low volume fraction suspension (φ/φ* = 1) of the attractive rods by changing the gap h between the plates (Fig. 9). We observed that beyond h = 300 μm (≈66Lrod), there was no formation of the vorticity aligned flocs. Furthermore, the floc diameter varies linearly with h. At larger gaps, we observe a response similar to shear banding (see Video, ESI†). However, with the existing imaging setup this could only be verified qualitatively. Interestingly, in the case of cone–plate geometry, the largest h at the edge is 122Lrod which is still able to induce vorticity aligned flocs. This is due to the initiation of the floc formation within the confinement zone close to the center of the cone–plate geometry. Since h increases from the center towards the edge of the cone–plate geometry, the cluster growth is unhindered at low Mn and the confinement effect is extended to larger values of h. This phenomenon puts into question the relation between the cluster length scale and confinement. It may rather be that once the floc formation is initiated, it grows in size even towards larger gaps until the volume fraction of particles within the floc matches the volume fraction of the gel at quiescent state. This also means that there should be a gradient in particle volume fraction within the flocs that are formed as a result of varying gap as observed in cone–plate geometry.
After discussing the macroscopic effects of steady shear flow on attractive rod suspension under confinement, we explore the microscopic changes using a combination of rheology and confocal microscopy. We capture the suspension microstructure before and after application of a steady shear rate of Per = 1 for the rod volume fraction at which we obtain log rolling flocs. A close inspection of the vorticity aligned flocs reveals the isotropic configuration of the rods within the particle rich floc. It is known that maximum packing for isotropic non-Brownian rods is expected to be around φmax(L/D) ≈ 5–6101 beyond which colloidal rods form ordered liquid crystal structures allowing mobility along their axis.13,102 However, attractions between rods can considerably slow down or arrest the diffusion driven ordering depending on the strength of attraction. Hence, our observations suggest that densely packed attractive rods show little thermal motion and form an isotropic microstructure.
We further need to explore how the heterogeneities arising from steady shear under confinement leads to the formation of log-rolling flocs that appear stable and equally spaced. In addition a quantitative determination of a possible liquid crystalline ordering by extracting the order parameter for the densely packed rods within the flocs would be valuable. As experimental observations are not able to easily provide such detailed structural information, we turn to Brownian dynamics simulations with hydrodynamic interactions in order to complementarily probe these phenomena and deduce answers to the above questions.
We found that for the weakest imposed shear rates (Mn = 0.001), log-like aggregates formed only for rod volume fractions smaller than 0.04 (φsim/φ* = 2). However, in these simulations stronger shear flows were able to break down the rod gel at higher volume fractions and also form vorticity aligned flocs. Fig. 10 shows the simulation results with L/D = 5 at a Mason number of 0.01 across different particle concentrations. Here the different colours indicate the different population size, in a 5% polydisperse size distribution, along the length of the rod which as discussed above simulates some roughness in the experimental rodlike particle. For the four cases that formed log-like aggregates, φsim = 0.02, 0.04, 0.06, 0.08 (or φsim/φ* = 1, 2, 3, 4) orientational order parameter was measured as S = 0.178, 0.085, 0.039 and 0.139, which indicates very little global alignment of the rods within the flocs.
In simulations of L/D = 10, which are closer in dimensions to the experimental suspensions, we also observe the formation of vorticity aligned flocs at volume fractions of up to 0.04 (φ/φ* = 5), when high enough Mason numbers were imposed. At large Mason numbers, Mn > 0.1, the aligned flocs were always destabilized in these simulations. Fig. 11 shows how multiple flocs form simultaneously in simulations and distribute themselves uniformly in the simulation cell with characteristic spacing between the flocs of about 2.7 times the channel height. From prior work we know that the Mason number characterizes the balance of hydrodynamic forces and interparticle attractions, and that this balance sets the preferred size of the aligned flocs. When the preferred floc size is larger than the size of the gap, alignment results. But when the preferred floc size is smaller than the gap, the alignment is disrupted. Fig. 12 depicts a slice through log-like flocs formed in experiment as well as in simulations with L/D = 10 and volume fraction of φ/φ* = 1 (exp.) and φsim/φ* = 1.25 (sim.) at Mn = 0.001 (exp.) and 0.01 (sim.) respectively. We find that there is no global alignment of rods within the floc, though in simulation results neighboring rods show some coordination in their alignment.
Fig. 11 Simulations for a φsim = 0.04 at Mn = 0.01 in which 3 log-like flocs form. The spacing between the flocs matches the usual experimental observation: ∼2.7h. |
There are some differences between the experimental and the simulation conditions. Two key differences are the ratio of the gap to the particle size and the interactions between the particles and the boundaries. In the simulations, gaps bigger than five particle lengths are inaccessible. The largest simulations run here used tens of thousands of rods, or hundreds of thousands of beads. Doubling the gap in the simulation box requires doubling the other dimensions of the simulation cell to maintain its aspect ratio. Doubling the dimensions would thus require eight times as many particles. Future work is needed to achieve small gaps in experiments.
In the experiments the maximum gap at the edge of the cone–plate geometry is 531 μm which is equal to 122Lrod and flocs covered the entire diameter of the geometry. In the center of the cone, the gap is closer to 15Lrod. However in the plate–plate geometry, flocs did not form above a gap of 66Lrod. The simulations probe systematically smaller gaps than in either experimental flow geometries, which means the break-up of aligned flocs should occur at systematically lower Mason numbers in the experiments. From recent work,57 the relationship between the critical Mason number and the gap in a parallel plate geometry is predicted to scale as Mn ∼ h−2.
While the simulations use this parallel plate geometry, experimental geometries differ in significant ways that might alter this scaling prediction. In the cone–plate geometry, the shear rate is the same at all radial distances while the gap grows radially. Thus, the Mason number does not vary throughout the sample, but the degree of confinement does. In the plate–plate geometry, the shear rate grows with radial position, while the gap remains fixed at all positions. Thus, the Mason number varies throughout the sample and the degree of confinement does not. These differences might lead to slightly different predictions for the conditions under which aligned flocs can be observed in each geometry. Further work is needed to understand why this transition is so robust to changes in flow geometry.
Based on the flow instabilities in attractive rod suspensions undergoing shear flow that were confirmed visually, we propose a state diagram (Fig. 13) relating the rod volume fraction and Mn. This is similar to the stability diagram proposed previously for dilute gels made up of spherical particles.56 At Mn > 0.01, the suspension appears to be homogeneous, however, this is due to the reduction in size of the heterogeneities to microscopic scale as confirmed using dark-field microscopy (Type B imaging), rather than complete breaking of the clusters to single particle level, which takes place at even higher rates. The large polydispersity in the floc size and time evolution of their size under steady shear rate flow conditions make it difficult to quantify the confinement (h/Lrod) associated with individual Mn and to propose a stability diagram. One may achieve this by using a parallel plate shear cell geometry where the shear rate and gap between the plates can be maintained constant throughout the duration of the measurement; a task that will be pursued in the future.
Initial yielding of rod gel (γ < 100%) at low rod volume fractions under oscillatory shear may be attributed to the bending and small scale restructuring of the percolated network rather than large scale rupture (Fig. 14(a)). Rupturing of the rod network occurs at larger strains (γ ≥ 200%). At low rod volume fractions, an individual rod in a percolated network has less number of contact points with the neighbouring rods. This allows an individual rod to re-orient itself in the direction of the imposed shear. However, a collective movement under shear of all the individual rods within a network from their initial position leads to a change in conformation of the network that mimics bending or flexing. This bending of the network before rupture pushes the yield strain to higher values.
Fig. 14 Proposed mechanism of yielding of rod gels under externally imposed stresses at low rod volume fraction (a) and at high rod volume fraction (b). |
A similar response has been reported for DLCA clusters of colloidal gels of spherical particles undergoing yielding.91 It is shown that these clusters undergo considerable deformation, such as unwinding, before rupturing. Here the mechanism proposed is the relative angular rotation of a singly connected particle pair. Under imposed stress, the large difference in stretching and bending moments in a singly connected backbone leads to full unwinding of the cluster before rupturing. Other mechanisms reported for a classic DLCA structure, explain the possibility of supporting large stresses if the structure is dominated by single-bond interactions formed via van der Waals forces. Here the bonds get rigid over time through contact-driven ageing and prevent rolling of contacts.90,92 Rod gels at low particle volume fractions are dominated by single-bond interactions and are susceptible to re-orientation at the contact point under shear. In addition they also exhibit contact – driven ageing by van der Waals forces.
At higher rod volume fractions, the percolated structure is rigid due to increased interparticle contact points leading to lower degree of freedom for individual rods to re-orient under externally imposed shear flow. This structure acquires larger G′ (more elastic), but on the other hand, is more brittle leading to fracturing at lower strains (Fig. 14(b)). Our study shows that even in the case of rod gels, mechanism of yielding exhibits a volume fraction dependence that was not reported earlier.91
The rigid rod suspension undergoing flow at steady shear rates, exhibit flow instabilities due to the competition between interparticle attraction and shear induced drag forces. Under confinement in narrow gaps, Fdrag < Fattr. (or Mn < 1), rod aggregates grow in size and form larger flocs. Since the rod aggregates are more porous compared to their spherical counterparts, the effect of confinement is more prominent, and the resultant floc resembles a “tumbleweed”. The continuous growth of particle rich flocs at low Per or Mn < 1 is restricted by the narrow gap in the flow-gradient plane. This leads firstly, to the densification of the porous aggregates and secondly, to their growth in the flow – vorticity plane exhibiting a “log-rolling” response at steady state. The fact that there is no local ordering of individual rods due to the increased rod density within the flocs corroborates the formation of low internal concentration, porous aggregates. Alignment or higher degree of local ordering may then be achieved only if the rod volume fraction is increased substantially.
Simulation results concur with the experimental observations, where at Mn < 1 rod aggregation occurs and at Mn < 0.1, the vorticity aligned flocs are stable. Furthermore, for a gap of 5Lrod, an inter-floc spacing (i.e. the distance between successive “logs”) of 2.7h was obtained at steady state. This reinforces the fact that under similar conditions, the mechanism of vorticity aligned floc formation is independent of particle aspect ratio and is strictly a function of interparticle attraction strength, the gap length scale and the particle volume fraction.
Nonlinear rheology measurements reveal that at low rod volume fractions (φ/φ* ≤ 2.5) and low steady shear rates (Per < 100), attractive rod suspensions exhibit flow instabilities and form large heterogeneities along the vorticity direction. Similar results were obtained using Brownian dynamics simulations. For Mason number, Mn < 1, the interparticle attraction forces dominate and particle aggregates remain and restructure under shear. When the aggregate growth is restricted due to confinement in the flow-gradient plane, cluster grows in the flow-vorticity plane eventually creating log rolling flocs. We also show that large gaps lead to other phenomena such as slip or shear banding as shown in ESI,† for an experiment with h = 500 μm. Simulations show that hydrodynamic forces prevent the flocs from colliding with each other while maintaining a distance between them of about three times (actually 2.7) the channel (or gap) height. Both in experiments and simulations, rods do not form any locally oriented structure within the densely packed flocs. A state diagram identifies the regimes of vorticity aligned floc formation as a function of rod volume fraction and Mason number.
The effects introduced by particle anisotropy on yielding and the log-rolling phenomena may be multiple. Firstly, particle shape anisotropy shifts the volume fractions associated with the transition to lower values compared to spheres. Due to this, shear induced clusters along the vorticity direction of rods are less dense than those of spheres which due to higher bending ability and intracluster rearrangements may densify more than their rod-like counterparts. Note that the latter have the additional degree of orientational freedom where nematic clusters would have a much higher density than isotropic ones. In the present system however, there is direct experimental evidence, backed by computer simulations, that local bending rigidity restricts rod reorientation thereby freezing the rotational degrees of freedom and hindering the internal structures from evolving towards clusters of nematic buddles. As was shown DLCA clusters of colloidal spheres undergo considerable deformation, such as unwinding, before rupturing whereas the same was not observed in the corresponding rods.91 This reflects the inextensibility (or lack of internal degrees of freedom) of a rod in comparison to a linear array of spheres with the same aspect ratio.
Our study shows that flow instabilities such as vorticity aligned floc formation arise when highly heterogeneous attractive particle suspensions flow through confined geometry. The fact that the phenomenon itself is independent of particle aspect ratio proves its universality and needs to be given due consideration when performing fundamental and application oriented studies in formulations including attractive colloids.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d0sm01576h |
This journal is © The Royal Society of Chemistry 2021 |