Deformation of flow-stabilized solids

Nonaffine deformation under compression and decompression of a flow-stabilized solid

Carlos P. Ortiz, Robert Riehn and Karen E. Daniels Department of Physics, North Carolina State University, Raleigh, NC, 27695, USA

Understanding the particle-scale transition from elastic deformation to plastic flow is central to making predictions about the bulk material properties and response of disordered materials. To address this issue, we perform experiments on flow-stabilized solids composed of micron-scale spheres within a microfluidic channel, in a regime where particle inertia is negligible. Each solid heap exists within a stress gradient imposed by the flow, and we track the positions of particles in response to single impulses of fluid-driven compression or decompression. We find that the resulting deformation field is well-decomposed into an affine field, with a constant strain profile throughout the solid, and a non-affine field. The magnitude of this non-affine response decays with the distance from the free surface in the long-time limit, suggesting that the distance from jamming plays a significant role in controlling the length scale of plastic flow. Finally, we observe that compressive pulses create more rearrangements than decompressive pulses, an effect that we quantify using the statistic for non-affine motion. Unexpectedly, the time scale for the compression response is shorter than for decompression at the same strain (but unequal pressure), providing insight into the coupling between deformation and cage-breaking.

  • July 24, 2019

1 Introduction

Understanding how structural rearrangements in disordered solids differ from crystalline solids is central [1, 2, 3] to achieving control of material properties such as resistance to flow [4], sound propagation  [5], heat capacity [6] , and dielectric constants [7]. For large deformations, the microscopic response differs non-perturbatively from the predictions of linear elasticity [8]. Instead of linear deformations, phenomena such as shear banding [9], yielding and plastic rearrangements [10], and non-local effects [11] are present. Recent experiments have explored non-affine deformations in 3D sheared colloidal glasses [12], 3D emulsions [13], and 2D foams [14]. For sufficiently slow deformations, it is an open question whether the flow behavior [15] is controlled by the jamming transition, where moduli vanish as the packing approaches a critical packing fraction [16].

In this paper, we present experiments quantifying the particle-scale deformation of flow-stabilized solids: particle heaps formed under controlled hydrodynamic stress (see Fig. 1). These quasi-2D heaps are assembled via the slow accumulation of micron-scale particles against a barrier within a microchannel, and are found to be stable above Péclet number  [17]. For lower Péclet numbers (slower flows) the particles reversibly evaporate away from the solid. We have previously observed that the elastic modulus of the solid is proportional to the confining stress provided by the fluid flow [18]. However, the amount of deformation of the solid in response to a flow perturbation is dependent on the sign of the perturbation: for piles prepared under identical conditions, compressions result in smaller strains than decompressions [18]. At the bulk scale, this effect can be understood by considering an excluded volume equation of state, as in thermal systems, under the assumption of locally affine deformations. In this paper, we investigate how the particle-scale dynamics lead to deviations from the excluded volume argument at high deformations.

Figure 1: (a) Schematic of channel geometry, not to scale, illustrating barrier with overflow. Streamlines are scaled versions of a calculation solving Stokes’ equations on the true device dimensions. (b) Image of central region of a heap of bidisperse particles. Dark regions are dimly fluorescent 600 nm particles and bright regions are brightly fluorescent 710 nm particles. White box indicates the location of the region in which deformations are analyzed.

Our experiments use sterically- and electrostatically-stabilized Brownian microspheres with a short Debye length (3 nm), so that the net interparticle interaction is well-approximated by a hard-sphere potential except near contact. A bidisperse mixture of particle sizes (5:4 diameter ratio) suppresses the nucleation of crystal domains. Using fluorescence microscopy, we measure the particle-scale deformation fields and characterize the response of the heap under compressive and decompressive loads created by changing the hydrodynamic stress. We characterize the influence of cooperative motions by measuring the degree to which the the deformation field locally deviates from global affine deformations.

The affine (or homogeneous) component of the deformation field is the part that can be described by an affine transformation: rotation, shear, extension, or compression [19]. After identifying the affine component of the deformation, the residual (or inhomogeneous) component is the non-affine deformation. For linearly-elastic materials, only affine deformations are present, but non-affine deformations can arise due to either thermally-driven cage-breaking events [20] or local rearrangements [2]. Here, we quantify two effects: the total non-affine rearrangements, and the spatiotemporal dynamics of the response. We observe, as expected [21, 22], that non-affine deformation fields typically exhibit mesoscale correlations. Furthermore, the non-affine deformations are about twice as large for compressive deformations as compared with decompressive deformations of similar size, and happen over a shorter time scale. For both compressive and decompressive deformations, non-affine deformations continue after affine deformations have completed.

2 Experimental Setup

Our experiments begin by assembling a microsphere heap by flowing a dilute suspension against a barrier (see Fig. 1a). The microchannel is fabricated to have a height  nm, higher than the height of a barrier ( nm), so that the fluid overflow accumulates particles against the barrier of width . The heights are chosen to create a quasi-2D heap, shallow enough to suppress both stacked and non-stacked bilayer phases [23]. The suspension is pumped into the channel by compressing a reservoir at the inlet using a low pressure, piezoelectrically actuated, digital regulator (AirCom PRE1-UA1), at  kPa above atmospheric pressure. After two hours of equilibration, the heap is deep ( angle of repose) and contains approximately 40,000 particles, as shown in Fig. 1b. The coordinate system takes parallel to the barrier and perpendicular to the barrier, with the origin at center of the barrier; the fluid flow is in the direction.

The dilute, aqueous suspension is prepared at a concentration of fluorescent microspheres. The particles are a bidisperse mixture of equal concentrations of 600 nm and 710 nm polystyrene microspheres (6% polydispersity, elastic modulus  GPa from Bangs Laboratories). We use steric and electrostatic stabilization (sulfate functionalized surface with -potential  mV and coated with Triton X-100) to provide reversible inter-particle and channel-particle interactions. The suspending fluid is a density-matched aqueous solvent at pH 5.4, buffered by citric acid to prevent crystallization and with sucrose to provide density-matching. The later suppresses segregation and sedimentation effects, important both at the barrier and at the inlet reservoir. Because the total particle brightness scales approximately with the particle volume, and we are working near the diffraction limit, the 600 nm particles appear dimmer than the 710 nm particles. This effect aids in tracking the motion of the particles.

We quantify the affine and non-affine deformation due to a pulse of either compression () or decompression (). The decompression pulse results from a change from 10 kPa to 0.5 kPa over a time scale of 10 ms, and the compression pulse returns the pressure to 10 kPa. This corresponds to for the compression pulse and for the decompression pulse. Imaging occurs in two phases. A full view of the initial heap is accessible with a objective, while experiments quantifying the particle motions require visualizing a zoomed-in region using a objective with a beam expander. The image was recorded by a CCD camera with  m pixels and an exposure time 36 . As shown in Fig. 1b, the zoomed measurement region is of size and is located adjacent to the barrier. For each pulse, we first acquire an image of a region of interest at the center of the heap prior to the pressure change, allowing us to extract the initial configuration of particle positions. Additional images, taken at 27 Hz, characterize the particle-scale response of the heap to the change in pressure. After a wait of 100 seconds, long enough for particles to settle onto a new, equilibrated configuration, we repeat this process for the compression pulse.

3 Image Processing

Figure 2: Illustration of image analysis process. (a) Sample micrograph, with grayscale indicating raw data. (b) Edge detection between bright and dark regions.(c) Wiener Deconvolved Image.

For either compressive or decompressive pulses, we first compare the initial and final configurations (separated by 100 sec), and quantify both the total deformation and the non-affine deformation. Second, using the series of frames immediately following the pulse, we track individual particles to identify non-affine effects on the local scale. Below, these are referred to as “long-time” and “short-time” dynamics, and require slightly different image-processing. For the long-time dynamics, the total distance traveled is on the order of a few particle diameters. Therefore, we first subtract the total affine deformation before performing particle tracking using the Blair-Dufresne implementation [24] of the Grier-Crocker particle tracking algorithm [25].

Particle identification:

Fig. 2 summarizes how we obtain particle positions beginning from a raw image. We identify the location of each particle by performing a Wiener deconvolution on the raw image, using a Gaussian approximation to the point-spread function with full-width at half maximum of 540 nm. This value is found to maximize the contrast in the output image, as measured from the ratio of the standard deviation to the mean intensity, but is more effective at locating the large (bright) particles than the small (dim) particles. The resulting deconvolved image allows us to detect the centroid of each particle using Matlab’s local extended maxima method.

Total deformation:

We estimate the total affine deformation due to a single pulse by making a coarse-grained measurement of the particle displacements between an initial image and a final image. These two images are created by averaging 10 initial images and 10 final images . We divide into horizontal strips of width and compute a cross-correlation with to determine its displacement. We find that the cross-correlation is sharply-peaked function for strips of at least this width.

Due to the large total strains, we perform particle pair matching between initial and final configurations based on particle positions from which the total affine deformation has already been subtracted. After this adjustment, pair identification proceeds as in the one-step particle tracking [24], with the size of the search region selected to correspond to the estimate of the maximum non-affine displacement amplitude, plus an estimate of the error in the affine strain.

Short-time particle tracking:

In order to obtain particle trajectories during the full duration of the dynamics, we make several assumptions about the nature of valid trajectories. We limit the displacement per frame to  ; this value is consistent with the total affine deformation rate determined above. In addition, we consider a particle’s identified size (brightness) in order to either split incorrect trajectories or or reconnect broken trajectories.

4 Results

In previous experiments [18], we observed that flow-stabilized solids exhibit a nonlinear stress-strain relationship in which the magnitude of the deformation of the surface of the flow-stabilized solid is well-described by


The success of this description is somewhat surprising, as Eqn. 1 does not contain any information about the distribution of stresses or strains throughout the flow-stabilized solid. The stress field within the solid is anticipated to be similar to that in a sedimentation experiment where particles “on top” of the sedimented material apply some stress on lower layers (in the limit of shallow sediments without side walls). We believe that the success in describing our experiments is due to the universality of the van-der-Waals thermal argument. However, that argument breaks down if non-affine motions occur, and we anticipate that the anticipated lower stress at the upstream (“top”) surface of the flow-stabilized solid is not fully characterized by the van-der-Waals argument. In the following, we first identify the distributions of particle displacements in the asymptotic long-term limit, before following individual trajectories through compression and decompression. Our particular interest is in the associated particle-scale non-affine motions and their dependence on the sign of .

4.1 Total deformation

Figure 3: Determination of the total (long-time) affine field from particle tracks for  kPa. (a) Image difference of blobs with size corresponding to of the position uncertainty in initial (blue) and final (yellow) configuration. (b) The associated affine deformation field. The dashed line is a linear regression with %.

We find displacements of individual particles in the heap immediately before and 100 s after a compressive/decompressive pulse through a two-step analysis. Following the homogeneous strain field assumption from our prior work [18], we first use the image cross-correlation analysis of images before and after deformation to obtain a global estimate of the affine strain field. We then use the affine transformation identified by the cross-correlation analysis as a scaffold for the matching of particles in the images before and after deformation. In Fig. 3a, we show an example of the particle locations after Wiener deconvolution and prior to finding the centroids, for both (red, before compression) and (white, after compression). By tracking each centroid, we can plot the local displacement as a function of -position within the heap.

As shown in In Fig. 3b, the mean behavior is linear, confirming that the overall assumption of an affine deformation was sufficiently accurate. The best fit line to these points provides a measure of the strain: with . The precision of these measurements is insufficient to estimate the expected higher-order (quadratic) term, although we expect one to be present due to a depth-dependent stress field. The observed linear behavior, combined with Eqn. 1, suggests that the packing fraction is close to invariant along the -direction.

Figure 4: Deformation field for a compression given by a pressure change at the inlet of  kPa. The axes in both figures are the -position and the -position in units of particle diameters. (a) Full deformation field. (b) Affine deformation field. (b) Non-affine deformation field, magnified by a factor of two.

Figure 5: Deformation field for a decompression given by a pressure change at the inlet of  kPa. The axes in both figures are the -position and the -position in units of particle diameters. (a) Full deformation field. (b) Affine deformation field. (b) Non-affine deformation field, magnified by a factor of two.

To obtain the total non-affine deformation field, we subtract the local affine motion from each displacement vector, as in [2]. Figs. 4 and 5 (both compression and decompression) show the total, affine, and non-affine displacement fields, for comparison. Importantly, we observe bands of correlated motions, as expected from [21, 22]. Because the total deformation field is not robust in tracking individual bead pairs over long times, we next examine the short-time dynamics.

4.2 Dynamics of individual particle tracks

Figure 6: Particle Tracks and Packing Fraction. (a) Particle tracks during compression. Tracks are colored to introduce contrast. (b) Particle tracks during decompression.

Using the estimate of strain provided by Fig. 3b, we track the fast dynamics arising from a compressive/decompressive pulse. For each frame, we first subtract the estimated affine deformation, based on the fraction of total strain which should have accumulated so far (see §3). This temporary adjustment allows for the correct association of particle centroids from frame to frame. Once the particle trajectories have been assembled based on these adjusted trajectories, we return to using the original positions detected for each particle centroid. The analyses that follow are based on the non-adjusted deformations that result from that tracking.

Sample trajectories are shown in Fig. 6. While the total deformation field is of the similar magnitude under compression and decompression, we find a more pronounced scrambling of the particle trajectories under compression, as compared to decompression. Below, we quantify both the affine and non-affine contributions to these trajectories.

Affine deformations:

Data was binned within strips along the -direction, providing ensembles of particles-dynamics sampled as a function of depth. Figs. 7a and 8a show the mean deformation field as a function of -position for compression and decompression, respectively. For both deformation directions, we find an exponential-like asymptotic approach to the final displacement magnitude. The depth-dependence of the asymptotic value of (Fig. 7b and 8b) demonstrates the same linear relationship originally shown in Fig. 3b. The resulting slope () quantifies the dynamics of affine reorganization. We find a marked difference between decompression (Fig. 7c) and compression (Fig. 8c) in that decompressions are far slower than compressions, and that the strain curves for decompression collapse better onto a single dynamic curve.

To quantify the difference, we make the Ansatz of a single-exponential approach to the asymptotic deformation


where is the particle displacement after a time interval , is a characteristic time scale of particle rearrangements, and is the asymptotic strain. Note that the value here is a fitting parameter; we find its value to be consistent with the estimate from the long-time dynamics. As shown in both panels (d), this exponential form is a good fit for the decompression pulses with  s. For compressive deformations, a single-exponential form is less consistent with the observed dynamics. Instead, there appears to be a two-step process of compression in which the viscous stress increase acts nearly instantaneously throughout the solid, while stresses due to particle-particle contacts propagate at a distinct speed of sound from the immobile barrier on which the solid is formed. Given the two-step nature of the process under compression, we establish an upper bound on the relaxation time scale of  s.

Figure 7: Affine Deformation Time-Series Analysis During Decompression. (a) Magnitude of correlated displacement field as a function of time at varying distances upstream of the barrier. (b) Long-time displacement amplitude , as a function of distance upstream of the barrier, both in units of particle diameters. (c) Correlated displacement amplitude normalized by the long-time displacement amplitude. (d) Log-linear plot of growth curves in (c).

Figure 8: Affine Deformation Time-Series Analysis During Compression. (a) Magnitude of correlated displacement field as a function of time at varying -positions. (b) Long-time displacement amplitude , as a function of distance upstream of the barrier, both in units of particle diameters. (c) Correlated displacement amplitude normalized by the long-time displacement amplitude. (d) Log-linear plot of growth curves in (c).

Figure 9: Magnitude of non-affine deformation field as a function of time and depth for (a) compression and (b) decompression.

Non-affine deformations:

We identify the non-affine contribution to the deformation field by subtracting the affine portion associated with the best-fit instantaneous value of the strain, which we designate


To quantify the resulting non-affinity, we use the measure [2] defined by


Here, is the set of local displacement vectors connection nearest neighbors, and is the time immediately before the pressure step was applied. Fig. 9 shows the time-evolution of the non-affine displacement as a function of -position during compression and decompression, respectively. In both graphs, grows and ultimately saturates. Interestingly, the magnitude of the non-affine field scales linearly with depth as demonstrated by the collapse of data series shown in Fig. 10. This a surprising finding in light of the assumed constant strain throughout the flow-stabilized solid. Furthermore, the magnitude of non-affine deformations is approximately twice as large under compression than under decompression at near identical asymptotic strain .

Figure 10: Plots of as a function of time, with the -axis scaled either by the size of the large particles or the -position at which the non-affine deformation is being probed. The legend in panel b applies to all panels. (a) Compression. scaled by . (b) Compression. scaled by . (c) Decompression. scaled by . (d) Decompression. scaled by .

We observe that the growth of with time is smooth. We are able to determine a characteristic time for the approach to the asymptotic value of by fitting a single exponential approach, as we did for the affine deformation field. In doing so, we neglect the low background value of in steady-state flow-stabilized solids arising from Brownian motion. We find  s, and  s. Therefore, the non-affine field significantly lags the affine field for decompression ( s). For compression, where a single time scale is less well defined, and an upper bound on the affine time scale is  s, the non-affine field also lags the affine deformation.

5 Discussion

We have observed particle-scale non-affine motions within flow-stabilized solids, and examined how their spatiotemporal dynamics depend on whether the deformation is compressive or decompressive. We observed the typical swirling regions often associated with non-affine deformations, arising through cooperatively rearranging regions. The magnitude of these effects is nearly twice as large under compression than under decompression, in spite of very similar total strains.

We observe that compressive pulses (large ) generate more non-affine deformation, which is able to dissipate the effect of the pulse more quickly. Because the non-affine fields for both compressive/decompressive deformations occur after similar delays with respect to the affine deformations, is suggests that they are triggered by the affine deformations. In the context of caging behavior, this suggests that the affine deformation distorts the cages provided by the neighboring particles and thereby makes Brownian cage-breaking (non-affine deformation) more likely. Remarkably, this is the case even though the strain is approximately the same for decompression and compression.

In probing the spatial dependence to the magnitude of the non-affine deformations (Fig. 9), we observe that the degree of non-affinity increases with distance from the barrier. This effect can be rescaled by the position to indicate a universal behavior. The form of this dependence suggests , for pressure and modulus [18]. One interpretation is that the surface of the heap is less rigid (smaller ), and therefore more prone to undergoing non-affine deformations (higher ). Similar effects have been observed in numerical simulations [22], where increasingly non-affine displacements are present in proximity to unjamming.

The significance of the above conclusions to soft-matter particle assemblies is to reinforce the centrality of understanding non-affine rearrangements to link bulk properties of the material, such as its modulus and global stability, to local properties about the typical particle geometry and rearrangement timescales. Based on these results, this experimental setup opens the possibility to explore this connection, by studying multiple orders of magnitude of heap sizes, under dynamically tunable interaction potentials and heap geometry, maintaining the ability relate particle-scale rearrangement dynamics to bulk properties. By doing so, it should be possible to determine length and time scales at which localized and collective rearrangements have the greatest impact on bulk properties, and shed light on the general mechanisms by which it is feasible to control the bulk properties of soft matter systems.


We are grateful for support from the National Science Foundation through an NSF Graduate Fellowship, grants DMR-0644743, DMS-0968258, DMR-1121107, MRSEC/DMR-112090, and INSPIRE/EAR-1344280. Research was also supported by US Army Research Office–Division of Earth Materials and Processes grant 64455EV. This work was performed in part at the Cornell NanoScale Facility, a member of the National Nanotechnology Infrastructure Network, which is supported by the National Science Foundation (Grant ECCS-0335765). This work was also performed in part at North Carolina State University facilities: Nanofabrication Facility, Advanced Instrumentation Facility, and Education and Research Laboratory.



  • [1] Schall P, Weitz D A and Spaepen F 2007 Science 318 1895–9
  • [2] Falk M L and Langer J S 1998 Phys. Rev. E 57 7192–7205
  • [3] Falk M L and Langer J 2011 Annu. Rev. Condens. Matter Phys. 2 353–373
  • [4] Brady J F 1993 J. Chem. Phys. 99 567
  • [5] Kriegs H, Petekidis G, Fytas G, Penciu R S, Economou E N and Schofield A B 2004 J. Chem. Phys. 121 7849–54
  • [6] Lubchenko V and Wolynes P G 2007 Annu. Rev. Phys. Chem. 58 235–266
  • [7] Bradshaw-Hajek B H, Miklavcic S J and White L R 2009 Langmuir 25 1961–9
  • [8] Bocquet L, Colin A and Ajdari A 2009 Phys. Rev. Lett. 103 036001
  • [9] Hays C C, Kim C P and Johnson W L 2000 Phys. Rev. Lett. 84 2901–2904
  • [10] Hébraud P, Lequeux F, Munch J and Pine D 1997 Phys. Rev. Lett. 78 4657–4660
  • [11] Lu C Y D, Olmsted P D and Ball R C 2000 Phys. Rev. Lett. 84 642–645
  • [12] Chikkadi V and Schall P 2012 Phys. Rev. E 85 031402
  • [13] Knowlton E D, Pine D J and Cipelletti L 2014 Soft Matter 10 6931–6940
  • [14] Twardos M and Dennin M 2005 Phys. Rev. E 71 061401
  • [15] Lerner E, Düring G and Wyart M 2012 Proceedings of the National Academy of Sciences of the United States of America 109 4798–803
  • [16] Liu A J and Nagel S R 2010 Annual Review of Condensed Matter Physics 1 347–369
  • [17] Ortiz C P, Riehn R and Daniels K E 2013 Soft Matter 9 543
  • [18] Ortiz C P, Daniels K E and Riehn R 2014 Physical Review E 90 022304
  • [19] Wu H C 2005 Continuum Mechanics and Plasticity (Chapman and Hall/CRC Press)
  • [20] Weeks E and Weitz D 2002 Phys. Rev. Lett. 89 095704
  • [21] Léonforte F, Tanguy A, Wittmer J and Barrat J L 2006 Phys. Rev. Lett. 97 055501
  • [22] Ellenbroek W G, Hecke M V and Saarloos W V 2009 Physical Review E 80 61307
  • [23] Van Winkle D and Murray C 1986 Phys. Rev. A 34 562–573
  • [24] Blair D and Dufresne E The matlab particle tracking code repository URL
  • [25] Crocker J C and Grier D G 1996 J. Colloid Interface Sci. 179 298–310
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description