Poiseuille flow of soft glasses in narrow channels: from quiescence to steady state

Poiseuille Flow of Soft Polycrystals in 2D Rough Channels

Polycrystals are partially ordered solids where crystalline order extends over mesoscopic length scales, namely, the grain size. We study the Poisuielle flow of such materials in a rough channel. In general, similar to yield stress fluids, three distinct dynamical states, namely, flowing, stick-slip, and jammed can be observed, with a yield threshold dependent on channel width. Importantly, the interplay between the finite channel width, and the intrinsic ordering scale (the grain size) leads to a new type of spatiotemporal heterogeneity. In wide channels, although the average flow profile remains pluglike, at the underlying granular level, there is vigorous grain remodeling activity resulting from the velocity heterogeneity among the grains. As the channel width approaches typical grain size, the flowing polycrystalline state breaks up into a spatially heterogeneous mixture of flowing liquid like patches and chunks of nearly static grains. Despite these static grains, the average velocity still shows a parabolic profile, dominated by the moving liquidlike patches. However, the solid-liquid front moves at nearly constant speed in the opposite direction of the external drive.

On the existence of thermodynamically stable rigid solids

Customarily, crystalline solids are defined to be rigid since they resist changes of shape determined by their boundaries. However, rigid solids cannot exist in the thermodynamic limit where boundaries become irrelevant. Particles in the solid may rearrange to adjust to shape changes eliminating stress without destroying crystalline order. Rigidity is therefore valid only in the metastable state that emerges because these particle rearrangements in response to a deformation, or strain, are associated with slow collective processes. Here, we show that a thermodynamic collective variable may be used to quantify particle rearrangements that occur as a solid is deformed at zero strain rate. Advanced Monte Carlo simulation techniques are then used to obtain the equilibrium free energy as a function of this variable. Our results lead to a unique view on rigidity: While at zero strain a rigid crystal coexists with one that responds to infinitesimal strain by rearranging particles and expelling stress, at finite strain the rigid crystal is metastable, associated with a free energy barrier that decreases with increasing strain. The rigid phase becomes thermodynamically stable when an external field, which penalizes particle rearrangements, is switched on. This produces a line of first-order phase transitions in the field-strain plane that intersects the origin. Failure of a solid once strained beyond its elastic limit is associated with kinetic decay processes of the metastable rigid crystal deformed with a finite strain rate. These processes can be understood in quantitative detail using our computed phase diagram as reference.

Attraction-induced jamming in the flow of foam through a channel

We study the flow of a pressure-driven foam through a straight channel using numerical simulations, and examine the effects of a tuneable attractive potential between bubbles. We show that the effect of an attractive potential is to introduce a regime of jamming and stick-slip flow in a channel, and report on the behaviour resulting from varying the strength of the attraction. We find that there is a force threshold below which the flow jams, and upon further increasing the driving force, a crossover from intermittent (stick-slip) to smooth flow is observed. This threshold force below which the foam jams increases linearly with the strength of the attractive potential. By examining the spectra of energy fluctuations, we show that stick-slip flow is characterized by low frequency rearrangements and strongly local behaviour, whereas steady flow shows a broad spectrum of energy drop events and collective behaviour. Our work suggests that the stick-slip and the jamming regimes occur due to the increased stabilization of contact networks by the attractive potential - as the strength of attraction is increased, bubbles are increasingly trapped within networks, and there is a decrease in the number of contact changes.

Yielding of glass under shear: A directed percolation transition precedes shear-band formation

Under external mechanical loading, glassy materials, ranging from soft matter systems to metallic alloys, often respond via formation of inhomogeneous flow patterns, during yielding. These inhomogeneities can be precursors to catastrophic failure, implying that a better understanding of their underlying mechanisms could lead to the design of smarter materials. Here, extensive molecular dynamics simulations are used to reveal the emergence of heterogeneous dynamics in a binary Lennard-Jones glass, subjected to a constant strain rate. At a critical strain, this system exhibits for all considered strain rates a transition towards the formation of a percolating cluster of mobile regions. We give evidence that this transition belongs to the universality class of directed percolation. Only at low shear rates, the percolating cluster subsequently evolves into a transient (but long-lived) shear band with a diffusive growth of its width. Finally, the steady state with a homogeneous flow pattern is reached. In the steady state, percolation transitions also do occur constantly, albeit over smaller strain intervals, to maintain the stationary plastic flow in the system.

Two-dimensional magnetic colloids under shear

Complex rheological properties of soft disordered solids, such as colloidal gels or glasses, inspire a range of novel applications. However, the microscopic mechanisms of their response to mechanical loading are not well understood. Here, we elucidate some aspects of these mechanisms by studying a versatile model system, i.e. two-dimensional superparamagnetic colloids in a precessing magnetic field, whose structure can be tuned from a hexagonal crystal to a disordered gel network by varying the external field opening angle θ. We perform Langevin dynamics simulations subjecting these structures to a constant shear rate and observe three qualitatively different types of material response. In hexagonal crystals (θ = 0°), at a sufficiently low shear rate, plastic flow occurs via successive stress drops at which the stress releases due to the formation of dislocation defects. The gel network at θ = 48°, on the contrary, via bond rearrangement and transient shear banding evolves into a homogeneously stretched network at large strains. The latter structure remains metastable after switching off of the shear. At θ = 50°, the external shear makes the system unstable against phase separation and causes a failure of the network structure leading to the formation of hexagonal close packed clusters interconnected by particle chains. At a microcopic level, our simulations provide insight into some of the mechanisms by which strain localization as well as material failure occur in a simple gel-like network. Furthermore, we demonstrate that new stretched network structures can be generated by the application of shear.

Length scales in glass-forming liquids and related systems: a review

The central problem in the study of glass-forming liquids and other glassy systems is the understanding of the complex structural relaxation and rapid growth of relaxation times seen on approaching the glass transition. A central conceptual question is whether one can identify one or more growing length scale(s) associated with this behavior. Given the diversity of molecular glass-formers and a vast body of experimental, computational and theoretical work addressing glassy behavior, a number of ideas and observations pertaining to growing length scales have been presented over the past few decades, but there is as yet no consensus view on this question. In this review, we will summarize the salient results and the state of our understanding of length scales associated with dynamical slow down. After a review of slow dynamics and the glass transition, pertinent theories of the glass transition will be summarized and a survey of ideas relating to length scales in glassy systems will be presented. A number of studies have focused on the emergence of preferred packing arrangements and discussed their role in glassy dynamics. More recently, a central object of attention has been the study of spatially correlated, heterogeneous dynamics and the associated length scale, studied in computer simulations and theoretical analysis such as inhomogeneous mode coupling theory. A number of static length scales have been proposed and studied recently, such as the mosaic length scale discussed in the random first-order transition theory and the related point-to-set correlation length. We will discuss these, elaborating on key results, along with a critical appraisal of the state of the art. Finally we will discuss length scales in driven soft matter, granular fluids and amorphous solids, and give a brief description of length scales in aging systems. Possible relations of these length scales with those in glass-forming liquids will be discussed.

Multi-scale lattice Boltzmann and mode-coupling theory calculations of the flow of a glass-forming liquid

We present a hybrid-lattice Boltzmann (LB) algorithm for calculating the flow of glass-forming fluids that are governed by integral constitutive equations with pronounced nonlinear, non-Markovian dependence of the stresses on the flow history. The LB simulation for the macroscopic flow fields is combined with the mode-coupling theory (MCT) of the glass transition as a microscopic theory, in the framework of the integration-through transients formalism. Using the combined LB-MCT algorithm, pressure-driven planar channel flow is studied for a schematic MCT model neglecting spatial correlations in the microscopic dynamics. The cessation dynamics after removal of the driving pressure gradient shows strong signatures of oscillatory flow both in the macroscopic fields and the microscopic correlation functions.