Effects of salinity on flows of dense colloidal suspensions

,


I. INTRODUCTION
Rheology of colloidal suspensions has been an active field for years, and it is well known that interactions (van der Waals, electrostatics, polymer brushes, etc.) between particles play an important role in the suspension behavior [1][2][3].In the last decade, numerous experimental and numerical works have led to a complete theoretical framework to describe the rheology of dense non-Brownian suspensions [4].The main difference with colloidal suspensions is that bigger particles are not sensitive to thermal motion, and that surface forces are usually negligible.However, surface interactions between grains have recently been proposed as the key ingredient to explain several macroscopic rheologic behavior, such as shearthickening [5,6] or shear-thinning [7] in non-Brownian suspensions.It is therefore interesting to study how changes in particles interactions can affect the flowing properties of the suspensions.In colloidal suspensions, salt concentration has been historically used as a way to tune the interactions between particles [8,9].In non-Brownian suspensions, salt has also been used recently as a way to effectively enhance the friction between particles [10,11], but the exact microscopic effect of the salt remains unclear.
The cross-over between thermal "colloids" and athermal "granular" suspensions is controlled by the gravita-tional Péclet number: with d the diameter of the particles, m = ∆ρ 6 πd 3 the mass the particles corrected by the buoyancy (∆ρ = ρ silica − ρ fluid is the difference of density between the particle and its surrounding fluid), g the gravitational acceleration, k B the Boltzmann constant, and T the temperature.In this work, we place ourselves in the intermediate regime of "dense colloidal suspensions" (P e ≳ 1), where the particles are fully sedimented, inducing a high concentration in particles, but the thermal agitation and the surface interactions cannot be neglected [12][13][14].In particular, we use silica micro-particles that are negatively charged and show a repulsive interaction in water, that we can tune by adding ions in solutions [15].In a previous work [14] we have shown that such dense colloidal suspensions show peculiar flow properties when inclined in rotating drum experiments.Above a threshold angle θ c , those suspensions exhibit a "fast avalanche" regime, that is similar to the one observed in non-Brownian ones.Below θ c they show a "slow creep" regime, that is thermally activated and depends heavily on the Péclet number.
In this study, we explore the flowing behavior of such dense colloidal suspensions when the repulsive interactions between the particles is progressively screened by salt added in the suspension.We focus in particular on the transition region between totally repulsive "frictionless" particles, and totally adhesive "non-flowing" particles.We interpret our results with the theoretical framework developed for non-Brownian suspensions, and we perform additional microscopic measurements of surface roughness, repulsion forces, and piles compacity, in order to connect the rheology of the suspensions to the microscopic interactions between the particles.

A. Microfluidic drums
Suspensions are made of silica particles from microParticles GmbH, with diameter d = 2.12±0.06µm, dispersed in solutions of NaCl (Sigma-Aldrich) in deionized water (ELGA Purelab®Flex, 18.2 MΩ • cm) with concentration C ranging from 0 to 2.5 × 10 −2 mol • L −1 .The gravitational Péclet number of those particles in suspension is P e ≈ 21.The suspensions are held in polydimethylsiloxane (PDMS) microfluidic drums, made with standard soft-lithography techniques.The drums have a diameter of 100 µm and a depth of 45 µm.
The microfluidic drums are filled using the following protocol: a PDMS sample is made with an array of thousands circular holes with the desired diameter and depth (once sealed, these holes will become the drums which contain the colloidal suspension).The PDMS sample is carefully washed and rinsed, first with isopropyl alcohol, then with deionized water.It is then cleaned for 15 min in deionized water in an ultrasonic bath.Next, it is immersed in a beaker containing the saline solution with the desired NaCl concentration C, and is let to degas for 15 min in the ultrasonic bath.The PDMS sample is removed from the ultrasonic bath and placed on a sample holder, the drums facing up.At this stage, the drums are only filled with the saline solution, and a drop (200 µL) of this solution is added on top of the sample, to avoid bubble formation due to evaporation.Then a droplet (30 µL) of a concentrated microparticles suspension is injected with a micropipette on top of the microdrums.The particles are let to sediment for 1 min in the drums.Finally, the drums are closed by placing a clean glass slide 1 on top of the PDMS sample, and by pressing it against the PDMS.The glass slide is maintained in position by six screws in the sample holder, which guarantees that the drums remains correctly sealed during the whole experiment.The particles typically fill ∼25 % of the drums volume.
The observation is made with the custom-made experimental set-up shown in Figure 1 a).It is a horizontal video-microscopy apparatus, made of a CCD camera (Basler acA2440-75um) linked to a microscope turret with long working distance microscope objectives (Olympus MPLFLN x10, and LUCPLFLN x20 ) through a lens tube (InfiniTube™Standard ), in front of a motorized rotation stage (Newport URB100CC ), with a manual 2D translation (Owis KT 90-D56-EP ) for the sample holder.To guarantee correct visualization of the sample, the rotation axis of the rotation stage is aligned with the optical axis of the video-microscopy system with a very high precision (up to a few microns).This axis is horizontal, so 1 The cleaning procedure for the glass slide is the same as the one for the PDMS sample.that the field of view contains the vertical gravity vector.To avoid external vibration, the whole set-up in installed on an optical table with passive isolation mounts (Thorlabs PWA075 ).
Before each measurement, the sample is shaken so that the particles are suspended, then let to sediment for 8 min, ensuring that the initial horizontal state of the pile is the same for each experiment.Then, the drums are rotated by an initial angle θ S = 30 • and images are taken with a logarithmic framerate for 24 h while the pile relaxes toward horizontal (at the beginning of the experiment the frame-rate is 20 images per second, at the end it is 10 −3 images per second).An example of experimental image is shown in Figure 1 b).Thanks to the use of low magnification microscope objectives, we are able to record simultaneously the flows in 20 different drums for each experiment.Images are then analyzed using contrast difference (contour finding algorithm from scikit-image) to automatically detect the top surface of the pile, and extract its angle with respect to the horizontal θ as a function of the time t.The typical dispersion between the 20 different angles θ i (t) obtained for the 20 different drums in the same experiment is of a few degrees (see fine brown curves in Figure 3).This leads to an accuracy of the average angle better than 0.5 • .AFM force spectroscopy studies are performed on a MFP-3D (from Asylum Research, Oxford Instrument) with a homemade colloidal probe [16].The probe is a silica bead of 10 µm, from the same manufacturer (mi-croParticles GmbH), glued to the end of a silicon nitride cantilever (DNP cantilever -Bruker ) with a bicomponent adhesive (Loctite EA3430 ).We used the thermal noise method [17,18] to know the nominal spring constant of the tooled cantilever for quantitative measurements.
The force curves are recorded in deionized water and in three different NaCl solutions: forces between the silica probe and a flat silica substrate are measured during the movement of the probe at constant velocity of 1 µm • s −1 towards the surface until contact is made (the maximum applied force is 6 nN) and then during the return.
Surface imaging of 2 µm silica particles is also performed in tapping mode (PPP NCHR AFM tip from NanoAndMore).Typical images have a spatial resolution of 2 nm • pixels −1 , and the total field of view is about 1 × 1 µm2 (512 × 512 pixels).To determine their RMS roughness value, the curvature of the spherical cap of the bead was subtracted from the image by a flatten of order 2 along the X and Y axis.

C. Confocal Microscopy
Stacks of images of dense colloidal suspensions at rest are obtained with a confocal microscope (Leica SP5, excitation wavelength 488 nm).The measurement is made with two different salt concentration ("pure" deionized water and 1 × 10 −2 mol • L −1 NaCl).Each suspension is introduced in a container, and a small amount of Rhodamine B is added (concentration 6 × 10 −6 mol • L −1 ).Before measurement, the suspension is let to sediment for 10 min.Scanns are performed at 400 Hz, with an oilimmersion microscope objective (HCX PL FLUOTAR 63.0x 1.25 OIL).The images stacks have a spatial resolution of 0.2 µm • pixels −1 in all three directions (XYZ), and the total field of view is about 200 × 200 × 10 µm3 .Before tracking the particles' coordinates, the contrast of each image is corrected.An example of the image obtained is shown in Figure 2. Due to the relatively high monodisperisty of the particles, and the high Péclet number, the sediment shows a poly-crystalline structure, as it has been observed in other experimental [19,20] or numerical [21] works.
We use the TrackPy software [22] to obtain the 3D coordinates of the particles present in the stack.After analysis, about 40000 particles distributed on up to 5 successive layers of sediment are found.

µm
FIG. 2. Image of the bottom layer of the suspension of 2.12 µm particles in deionized water, observed with the confocal microscope.The contrast has been inverted so that the particles appear light on a dark background.Poly-crystalline structure is visible, due to the relatively high Péclet number and high monodisperisty of the particles.

A. Microfluidic drums
The effect of the salt concentration on the average flow curves of dense colloidal suspensions in micro-fluidic drums is presented in Figure 3.For the lowest ionic strength 2 (bottom curve), the typical flow behavior is retrieved, with a "fast avalanche" regime at high angle, and a "slow creep" regime, which is logarithmic in time, below a threshold angle θ c [14].However, when the ionic strength is increased, the flow behavior progressively changes, as the repulsive force between the particles is progressively screened.The most noticeable changes are: the increase of a waiting time plateau at the beginning of the experiment before the fast avalanche regime starts, the increase of the time needed for the fast avalanche regime to reach the threshold angle θ c , the increase of the threshold angle θ c itself, and the increase of the final angle at the end of the experiment.There might also be a small change in the slow creep regime, however, as discussed later, this effect is less clear.For the highest ionic strength of 5 × 10 −2 mol • L −1 , almost no flow is observed when the drums are rotated 3 , and particles agglomerates are visible in the pile.This is consistent with the fact that ionic strengths above 5 × 10 −2 mol • L −1 can be used to generate floculated suspensions [23] or colloidal gels [24] from silica particles dispersed in water.Those agglomerates of particles also lead to irreproducible avalanche curves (as shown in Figure 3), because they can sustain non-flat pile interfaces where the angle θ is not properly defined anymore.To better quantify the effect of the salt concentration, we define a few experimental quantities, schematically presented in Figure 4. We call τ s the "starting time" of the fast avalanche, that is the time required for the pile to reach 95 % of its initial angle θ S .We fit both the end of the fast avalanche regime, and the slow creep regime by a linear function in the semilogarithmic plot (i.e.θ = A log(t) + B with A and B two constants).We define the threshold angle θ c as the crossing point between the two fitted regimes.We call "avalanche speed" ∆θ/∆t, the average flowing rate of the fast avalanche regime.We call S the slope of the slow creep regime in the semilogarithmic time-scale.Note that those four quantities are defined in different temporal regions of the flow curve, and are mathematically independent one from another.
The measured values are presented in Figure 5.Both τ S and θ c increase slightly faster than linearly with the ionic strength, with a steep increase when the ionic FIG. 4. Schematical description of the quantities of interest measured for each flow curve.θS is the initial starting angle ; τS is the starting time needed to reach 0.95 θS ; ∆θ/∆t is the average speed of the avalanche ; θc is the threshold angle between the two regimes ; and S is the slope of the creep regime.
strength is close to 5 × 10 −2 mol • L −1 .The avalanche speed ∆θ/∆t seems to decrease almost linearly with the ionic strength.Finally, the slope of the creep regime S seems to first increase up to a maximum when the ionic strength is close to 1 × 10 −2 mol • L −1 , and then slightly decreases.The colloids surface roughness is measured by AFM imaging.The RMS roughness of a 2 µm silica particle is less than 1 nm over 1 µm × 1 µm.A typical image of the surface imaging is shown in Figure 6 a FIG. 6. a) AFM imaging of a 2 µm silica particle's surface roughness (the lateral resolution is 2 nm • pixels −1 , and total field of view is 1 × 1 µm 2 ).The topographic profiles extracted from the image along a vertical and horizontal lines are respectively shown to the left and on top of the image.b) Repulsive forces measured between a 10 µm silica particle and a flat silica surface, for different concentrations of the saline solution.The black dashed-lines correspond to numerical fits with the theoretical formula for the electrostatic double-layer interaction force between a sphere and a flat surface (see eq.

and table I for best fitting parameters values).
Typical repulsive forces F are shown in Figure 6 b) as a function of the distance D between the particle's surface and the flat silica surface, for different salt concentrations.At large enough separation distance, they show a good agreement with the theoretical electrostatic double-layer forces between surfaces in liquids at constant potential [15]: where, d is the particle diameter, λ D is the Debye length, and Z is a constant equal to 9.22 × 10 −11 tanh 2 (ψ 0 /103) J • m −1 at 25 • C, with ψ 0 the surface potential in mV.Dashed lines in Figure 6 b) are the numerical fits of the forces curves F (D) with eq. 2, using two free parameters (ψ 0 and λ D ).The average best fitting values obtained for the surface potential ψ 0 and the Debye length λ D are presented in table I.For each salt concentration, at least 4 independent force curves have been measured.The errorbars are estimated from the fits accuracy and data dispersion.Values of the surface potential ψ 0 are expected to be found between −70 mV and −20 mV in pure water [25][26][27][28], and to decrease with the salt concentration [25,26].Note that the exact value depends on the surface state (cleanness, roughness) of the silica [27], and that theoretical values can be hard to determine [29].On the contrary, the expected values of the Debye lengths can easily be computed from the ionic strength of the solutions [15]: at 25 • C, λ D = 0.304/ √ C nm, with C the monovalent salt concentration in mol•L −1 .The expected values are presented in table I and show a good agreement with the ones measured experimentally.Note that the ionic strength of the suspension made with "pure" deionized water is unknown, but is expected to be about a few 10 −5 mol • L −1 due to water contamination from dissolved carbon dioxyde [26] and from the colloids themselves.The value λ D ≈ 65 nm that we find corresponds to a ionic strength of 2 × 10 −5 mol • L −1 .To estimate the typical volume that is occupied by one particle on the bulk of each suspension, we compute the 3D Voronoï tessellation (scipy.spatiallibrary based on the Qhull library) of the particles 3D coordinates obtained by confocal microscopy.We first remove the particles on the side/edges of the sample, then compute the distribution of volumes of the Voronoï cells of the remaining particles.The probability density functions (PDF) of the volume occupied by a single particle in the bulk are shown in figure 7, for two different salt concentration.
Finally, the mode of the distribution is taken as the typical volume around each particle.This value is used to compute the initial packing fraction Φ 0 of the suspension, given that the actual volume of one particle is know (πd 3 /6).The measured packing fractions are presented in table II for the two tested salt concentrations.

IV. DISCUSSION
In this section, we provide physical explanations for the changes in flowing behaviors of the dense colloidal suspensions when the salt concentration is increased.

A. Critical angle
Salt added to screen repulsive interactions between silica particles has been used as a way to modify the friction between the grains in non-Brownian suspensions [10,11].The main idea is that the double-layer repulsive force, which size is typically given by the Debye length λ D , prevents the particles from having frictional contacts as long as λ D is bigger than the surface roughness of the grains When salt is added, the Debye length becomes smaller, up to a point where the particles can touch each others.This induces a frictionless to frictional transition for the suspension around λ D ≈ r, which leads to shear-thickening behavior [10] or hysteretic flows [11].
This idea can be used to simply explain the observed increase of the critical angle θ c when the ionic strength increases (Figure 5 b)).Indeed, θ c corresponds roughly to the "angle of repose" of the granular suspension: it's the angle below which no flow should be observed if the pile was non-Brownian.Both numerical [30][31][32] and experimental [33] studies have shown that the angle of repose of a dry granular pile increases when the microscopic friction coefficient between the grains µ p is increased.Therefore, one can expect that the increase in salt concentration, increases the effective friction between the particles, which then increases the angle of repose of the pile.Notably, we see that the measured critical angle θ c in deionized water is about 4.6 • , which is close4 to the angle of repose 5.76 • that is observed in numerical simulation for frictionless particles [34].Moreover, θ c in the solution with the highest salt concentration (4 × 10 −2 mol • L −1 ) is about 17.2 • , which is not too far from the repose angle of 25 • for macroscopic glass beads [35].
Our AFM measurements support this hypothesis.As shown in Fig. 6, the typical peak surface roughness r of our particles is about 2 nm, and the repulsive force F (D) between one particle and a flat silica surface is well described by the theoretical double-layer electrostatic theory (Equation 2).Therefore one can estimate that the critical salt concentration where the Debye length becomes equal to the surface roughness (λ D ≈ r) is about C c = 2.3 × 10 −2 mol • L −1 .This is consistent with the critical ionic strength at which we see a transition from colloidal piles completely flowing back to horizontal C ≤ 1 × 10 −2 mol • L −1 , to completely arrested colloidal piles C ≥ 5 × 10 −2 mol • L −1 (see Figure 3).Note that we cannot directly measure the microscopic friction coefficient µ p between particles with our experimental set-up.However, values gathered in the literature can be found in Lee et al. [36]: silica microparticles have a typical friction coefficient 0.03 ≤ µ p ≤ 0.1 in milli-Q water, µ p ≈ 0.3 in NaCl solution with concentration C = 1 × 10 −3 mol • L −1 , and µ p ≈ 0.9 in alkaline solution with ionic strength 16 × 10 −2 mol • L −1 .

B. Starting time
The increase of the delay before the start of the "fast avalanche" regime (Figure 5 a)) is reminiscent of the dilatancy effects that are observed in macroscopic granular suspensions [4,37,38].When a pile of grains is fully immersed in a Newtonian fluid, its flowing properties strongly depend on its initial packing fraction.After a sudden inclination, loosely packed piles tend to flow almost immediately, while densely packed piles show a time delay before the initiation of the flow.This phenomenon is explained by a pore pressure feedback scenario: a densely packed pile has to dilate before being able to flow, and during this dilation, the surrounding fluid is sucked into the granular layer, which tends to stabilize the pile and delay the start of the flow [38].
Since the ionic strength reduces the double-layer repulsive force between the particles, one can expect that it reduces the mean distance between particles, hence increasing the initial packing fraction Φ 0 of the pile.Therefore, we can expect that the ionic stress increases the time delay τ S before the start of the flow, due to dilatancy effects.Our set-up does not allow us to measure the packing fraction Φ of the pile during the flow, to directly observe dilatancy effects.However, our confocal microscopy measurements support the fact that the initial packing fraction Φ 0 of the sedimented pile increases with the ionic strength of the suspension.As shown in table II, the packing fraction is about 51 % in deionized water and increases to about 61 % in a solution with NaCl concentration 1 × 10 −2 mol • L −1 .Notably, the critical packing fraction Φ 0C above which dilatancy effect are observed in macroscopic granular suspensions [38] is about 58 • .This is consistent with the fact that we observe almost immediate flow in deionized water (τ S = 0.87 s), while we observe significant start delay with high ionic strength suspensions (τ S = 40.9s for C = 4 × 10 −2 mol • L −1 ).
Nevertheless, two points must be noted.First, the increase of the initial packing fraction Φ 0 that is observed with the increase of the salt concentration might seem surprising.Indeed, for macroscopic granular materials, it is known that the packing fraction obtained after sedimentation of the granular medium (random loose packing) decreases when the friction between the particles increases [39][40][41].Since we have already shown that the effective friction between the particles µ p increases with the salt concentration, one could expect that the packing fraction would rather be lower when the ionic strength of the suspension is higher.The solution to this apparent contradiction comes from the fact that our suspensions are Brownian: we believe that the thermal agitation helps the suspension to always reach the highest accessible packing fraction (random close packing).Second, the fact that we observe dilantancy effects is itself a proof that the friction between the grains increases with the salt concentration.Indeed, numerical simulations have shown that frictionless grains do not show dilatancy effects [34].

C. Avalanche speed
The fact that the salt concentration increases the effective friction between the particles can also explain the observed decrease of the avalanche speed ∆θ/∆t (Figure 5 c)).Indeed, both numerical simulations [42][43][44][45][46] and experimental works [47,48] have shown that the rheology of dense non-Brownian suspensions depends on the microscopic friction coefficient between the grains.A recent review can be found in Lemaire et al. [7], and we only recall here a few key results.For example, in volumeimposed simulations, the viscosity of the suspension η S increases with the microscopic friction coefficient µ p .In pressure-imposed simulations, the stress ratio µ (which can be seen as the macroscopic friction coefficient) increases with µ p , while the volume fraction Φ decreases with µ p , at fixed viscous number J 5 .
In general, it is expected that the flow rate of the suspension Q decreases when the viscosity η S increases, when the stress ratio µ increases, and when the volume fraction Φ decreases.Thus, an increase of the microscopic friction coefficient is expected to lead to a decrease of the avalanche speed ∆θ/∆t.Direct comparison are difficult to achieve, since it is non-trivial to compute the theoretical flow rate Q in the rotating drum geometry 6 .But the orders of magnitude are reasonable.In our experiment we observe that ∆θ/∆t decreases by a factor of ∼ 10 between pure water (∆θ/∆t = 0.55 • • s −1 ) and high ionic strength suspensions (∆θ/∆t = 0.03 • • s −1 for C = 4 × 10 −2 mol • L −1 ).In simulations [45], the viscosity of the suspension increases by a factor of 10 when the microscopic friction coefficient µ p increases from 0 to 1 at volume fraction Φ = 55 %.

D. Slope of the creep regime
Following previous work [14,50], the time evolution of the pile angle θ during the creep regime can be described with a simple model where particles in the top layers are considered blocked by their neighbors, and the creep occurs when they jump above those neighbors thanks to thermal agitation.This model gives the following mathematical expression: 2 αP e arcoth exp t τ αP e e −αPeθc coth (αP e θ S /2) (3) where: α is a dimensionless geometric parameter, P e is the gravitationnal Péclet number, τ is a characteristic time depending on the fluid's properties and drum geometry, θ c is the critical angle, θ S is the initial inclination angle (if θ S ≤ θ c ).
If P e ≫ 1, and θ c ≪ θ ≪ 0, the equation 3 can be approximated by: Equation 4 directly gives the slope of the creep regime: S = 1/(α P e ).Knowing that P e ≈ 21.4 for the 2.12 µm particles and that α ≈ 2.6 was found in previous experiments [14], we can expect S ≈ 0.0177 rad = 1.01 • .This is consistent with the values we measured (0.6 • ≤ S ≤ 1.6 • , see Figure 5 d)).However, the model does not predict a significant variation of S with the salt concentration C. Indeed, when salt is added to the suspension, P e only slightly varies because the density ρ fluid of the salted water varies.Even with the most concentrated solution (C = 5 × 10 −2 mol • L −1 ) the density only increases by ∼ 3 %, which leads to a small decrease P e ≈ 20.7.As for α, it corresponds to the "height" of the barrier that one particle has to cross to jump over its neighbors.One can assume that this value slightly decreases when the Debye length λ D decreases because the particles has to jump above a particle of effective diameter d + 2λ D .For example, if λ D goes from 50 nm (deionized water) to 1 nm (high salt concentration), this would predict that α decreases by 0.1 2.12 ≈ 5 %.In the end, following the model, the slope of the creep regime S should monotonically increase when the salt concentration increases, and should not vary by more than ∼ 10 %.Therefore, it remains unclear whether the variations that we observe in Figure 5 d) are real physical effects, or experimental artifacts due to the difficulty to measure small pile angles during long times 7 .We believe that accurate measurement of S should be obtained from flow curves with an initial inclination angle θ S below the threshold angle θ C , to avoid any influence of the transition between the creep regime and the preceding fast avalanche regime.

V. CONCLUSION
In conclusion, we have measured the flow of dense colloidal suspensions, in microfluidic drums after an initial inclination angle, for different salt concentrations.The flowing curves show two regimes: a "fast avalanche" regime above a critical angle θ c , and a "slow creep" regime (logarithmic in time) below θ c .We observe that the flowing behavior is strongly modified by the ionic strength of the suspension.As the salt concentration increases, the initial time delay τ S before the fast avalanche regime increases, the speed of this regime ∆θ/∆t decreases, and the critical angle θ c increases.All those observations are well explained by the fact that ions added in solution screen the repulsive double-layer electrostatic forces between the colloidal particles, which increases the effective microscopic friction between the particles µ p and the initial packing fraction of the suspension Φ 0 .We have independently verified with AFM measurements that the particles roughness r is consistent with the critical salt concentration C c ∼ 2.3 × 10 −2 mol • L −1 at which we observe a transition from "very flowing curves" to "almost not flowing curves" (with particles agglomerates).We have also verified with direct confocal microscopy observations that the initial packing fraction of the sedimented suspension increases from Φ 0 ∼ 51 % in deionized water to Φ 0 ∼ 61 % in solution with ionic strength 1 × 10 −2 mol • L −1 .This explains why increasing dilatancy effects are observed when more salt is added to the suspension.
Finally, even though all our measurements seem to indicate that the microscopic friction between the particles µ p is increased by the salt concentration C, we cannot conclude on the physical origin of this effective friction increase.Indeed, this effective friction might come from direct contact friction (if the particles surfaces roughness touch each others), or from indirect hydrodynamic interactions (either long-range pore pressure effects, or shortrange lubrication effects).Numerical simulations tend to show that contact friction dominates over long-range hydrodynamics at high volume fraction [51] (Φ ≥ 40 %), and over both long-range and short-range hydrodynamics at low viscous number [45] (J ≤ 10 −1 ).However, only direct measurement of the normal and tangential forces between two colloidal particles (such as those obtained with quartz-tuning fork atomic force microscopy [52,53], or lateral force microscopy [54,55]), in different ionic strength suspensions, would be able to experimentally confirm this result in our system.

FIG. 1 .
FIG. 1. a) Schematic representation of the horizontal videomicroscopy apparatus, made to visualize the flow of the dense colloidal suspensions in vertically held microfluidic drums.b) Example of experimental image measured during the pile relaxation experiment: an array of 8 drums filled with ∼25 % of colloidal suspension made of 2.12 µm particles in 1 × 10 −3 mol • L −1 NaCl solution.The measured angle θ is shown in blue on the first drum.Note that the field of view has been cropped for readability, and that real images used during experiments have 20 drums in the field of view.

FIG. 3 .
FIG.3.Average flow curves for dense suspensions of 2.12 µm silica particles inclined with same starting angle θS = 30 • , in salty solutions with different concentrations.Values of the NaCl concentration is indicated in legend: the bottom curve (dark blue) corresponds to the solution with the lowest concentration (deionized water, ionic strength estimated ∼ 1 × 10 −5 mol • L −1 ), and the top curve (light blue) corresponds to the solution with the highest concentration (ionic strength = 5 × 10 −2 mol • L −1 ).Each curve is obtained by averaging pile angles θ(t) measured on the 20 drums visible in a single experiment.To help visualize the typical data dispersion, the 20 individual flow curves used to compute the average are showed for the NaCl concentration 1 × 10 −2 mol • L −1 (fine brown curves).

FIG. 5 .
FIG. 5. Measured values of a) the starting time τS, b) the critical angle θc, c) the avalanche speed ∆θ/∆t, and d) the slope of the creep regime in semilogarithmic time scale S, as a function of the ionic strength of the suspension.Color points are individual measurements, black stars are the averaged values.Gray dashed lines corresponds to a linear fit where each average point is weighted by the inverse of its standard deviation.

TABLE I .
Average best fitting parameters for surface potential ψ0 and Debye length λD obtained from AFM repulsive force measurements (errorbars are obtained from the fit accuracy and data dispersion).