Ryohei Seto

Discontinuous Shear Thickening of Frictional Hard-Sphere Suspensions

Discontinuous shear thickening (DST) observed in many dense athermal suspensions has proven difficult to understand and to reproduce by numerical simulation. By introducing a numerical scheme including both relevant hydrodynamic interactions and granularlike contacts, we show that contact friction is essential for having DST. Above a critical volume fraction, we observe the existence of two states: a low viscosity, contactless (hence, frictionless) state, and a high viscosity frictional shear jammed state. These two states are separated by a critical shear stress, associated with a critical shear rate where DST occurs. The shear jammed state is reminiscent of the jamming phase of granular matter. Continuous shear thickening is seen as a lower volume fraction vestige of the jamming transition.

Shear thickening, frictionless and frictional rheologies in non-Brownian suspensions

Particles suspended in a Newtonian fluid raise the viscosity and also generally give rise to a shear-rate dependent rheology. In particular, pronounced shear thickening may be observed at large solid volume fractions. In a recent article [R. Seto et al., Phys. Rev. Lett. 111, 218301 (2013)], we have considered the minimum set of components to reproduce the experimentally observed shear thickening behavior, including discontinuous shear thickening. We have found frictional contact forces to be essential and were able to reproduce the experimental behavior by a simulation including this physical ingredient along with viscous lubrication. In the present article, we thoroughly investigate the effect of friction and express it in the framework of the jamming transition. The viscosity divergence at the jamming transition has been a well known phenomenon in suspension rheology, as reflected in many empirical laws for the viscosity. Friction can affect this divergence, and in particular the jamming packing fraction is reduced if particles are frictional. Within the physical description proposed here, shear thickening is a direct consequence of this effect: As the shear rate increases, friction is increasingly incorporated as more contacts form, leading to a transition from a mostly frictionless to a mostly frictional rheology. This result is significant because it shifts the emphasis from lubrication hydrodynamics and detailed microscopic interactions to geometry and steric constraints close to the jamming transition.Particles suspended in a Newtonian fluid raise the viscosity and also generally give rise to a shear-rate dependent rheology. In particular, pronounced shear thickening may be observed at large solid volume fractions. In a recent article [R. Seto et al., Phys. Rev. Lett. 111, 218301 (2013)], we have considered the minimum set of components to reproduce the experimentally observed shear thickening behavior, including discontinuous shear thickening. We have found frictional contact forces to be essential and were able to reproduce the experimental behavior by a simulation including this physical ingredient along with viscous lubrication. In the present article, we thoroughly investigate the effect of friction and express it in the framework of the jamming transition. The viscosity divergence at the jamming transition has been a well known phenomenon in suspension rheology, as reflected in many e...

Discontinuous shear thickening in Brownian suspensions by dynamic simulation

Significance Colloidal shear thickening is a phenomenon by which a dense mixture (or suspension) of sub–micrometer-sized particles in a fluid becomes more viscous when strained. It has important (mostly negative) consequences as dense suspensions are found in many industrial processes. Its origin is a long-standing unsolved problem in soft-matter physics. We perform dynamic simulations that show that the shear thickening may be explained as a stress-induced transition from a flow of lubricated near-contacting particles to a flow of a frictionally contacting network of particles. Our numerical data quantitatively match previous experimental data. Dynamic particle-scale numerical simulations are used to show that the shear thickening observed in dense colloidal, or Brownian, suspensions is of a similar nature to that observed in noncolloidal suspensions, i.e., a stress-induced transition from a flow of lubricated near-contacting particles to a flow of a frictionally contacting network of particles. Abrupt (or discontinuous) shear thickening is found to be a geometric rather than hydrodynamic phenomenon; it stems from the strong sensitivity of the jamming volume fraction to the nature of contact forces between suspended particles. The thickening obtained in a colloidal suspension of purely hard frictional spheres is qualitatively similar to experimental observations. However, the agreement cannot be made quantitative with only hydrodynamics, frictional contacts, and Brownian forces. Therefore, the role of a short-range repulsive potential mimicking the stabilization of actual suspensions on the thickening is studied. The effects of Brownian and repulsive forces on the onset stress can be combined in an additive manner. The simulations including Brownian and stabilizing forces show excellent agreement with experimental data for the viscosity η and the second normal stress difference N2.

Restructuring of colloidal aggregates in shear flow Coupling interparticle contact models with Stokesian dynamics

A method to couple interparticle contact models with Stokesian dynamics (SD) is introduced to simulate colloidal aggregates under flow conditions. The contact model mimics both the elastic and plastic behavior of the cohesive connections between particles within clusters. Owing to this, clusters can maintain their structures under low stress while restructuring or even breakage may occur under sufficiently high stress conditions. SD is an efficient method to deal with the long-ranged and many-body nature of hydrodynamic interactions for low Reynolds number flows. By using such a coupled model, the restructuring of colloidal aggregates under shear flows with stepwise increasing shear rates was studied. Irreversible compaction occurs due to the increase of hydrodynamic stress on clusters. Results show that the greater part of the fractal clusters are compacted to rod-shaped packed structures, while the others show isotropic compaction.Graphical abstract

Compressive consolidation of strongly aggregated particle gels

The compressive yield stress of particle gels shows a highly nonlinear dependence on the packing fraction. We have studied continuous compression processes and discussed the packing-fraction dependence with the particle-scale rearrangements. The two-dimensional simulation of uniaxial compression was applied to fractal networks, and the required compressive stresses were evaluated for a wide range of packing fractions that approached close packing. The compression acts to reduce the size of the characteristic structural entities (i.e., the correlation length of the structure). We observed three stages of compression: (I) Elastic-dominant regime; (II) single-mode plastic regime, where the network strengths are determined by the typical length scale and the rolling mode; and (III) multimode plastic regime, where sliding mode and connection breaks are important. We also investigated the way of losing the fractal correlation under compression. It turns out that both fractal dimension Df and correlation length ξ start to change from the early stage of compression, which is different from the usual assumption in theoretical models.

Microstructure and thickening of dense suspensions under extensional and shear flows

Dense suspensions are non-Newtonian fluids which exhibit strong shear thickening and normal stress differences. Using numerical simulation of extensional and shear flows, we investigate how rheological properties are determined by the microstructure which is built under flows and by the interactions between particles. By imposing extensional and shear flows, we can assess the degree of flow-type dependence in regimes below and above thickening. Even when the flow-type dependence is hindered, nondissipative responses, such as normal stress differences, are present and characterise the non-Newtonian behaviour of dense suspensions.

Novel Aspects of Evolution of the Stokes Parameters for an Electromagnetic Wave in Anisotropic Media

Polarization of a plane electromagnetic wave travelling through a medium is studied in the slowly-varying field envelope approximation. It is shown that the problem is identical to the 4-momentum evolution of a negatively-charged massless relativistic particle in an electromagnetic field. The approach is exemplified by the resonant oscillations of circular polarization in a medium embedded in a static magnetic field and a modulated electric field. The effect of dissipation in the medium is discussed. It is shown that the Rabi oscillations are stable below a threshold depending on the absorption coefficient. Above it, oscillations disappear.

Crystallization kinetics of binary colloidal monolayers

Experiments and simulations are used to study the kinetics of crystal growth in a mixture of magnetic and nonmagnetic particles suspended in ferrofluid. The growth process is quantified using both a bond order parameter and a mean domain size parameter. The largest single crystals obtained in experiments consist of approximately 1000 particles and form if the area fraction is held between 65-70% and the field strength is kept in the range of 8.5-10.5 Oe. Simulations indicate that much larger single crystals containing as many as 5000 particles can be obtained under impurity-free conditions within a few hours. If our simulations are modified to include impurity concentrations as small as 1-2%, then the results agree quantitatively with the experiments. These findings provide an important step toward developing strategies for growing single crystals that are large enough to enable follow-on investigations across many subdisciplines in condensed matter physics.

A theoretical framework for steady-state rheometry in generic flow conditions

We introduce a general decomposition of the stress tensor for incompressible fluids in terms of its components on a tensorial basis adapted to the local flow conditions, which include extensional flows, simple shear flows, and any type of mixed flows. Such a basis is determined solely by the symmetric part of the velocity gradient and allows for a straightforward interpretation of the non-Newtonian response in any local flow conditions. In steady homogeneous flows, the material functions that represent the components of the stress on the adapted basis generalize and complete the classical set of viscometric functions used to characterize the response in simple shear flows. Such a general decomposition of the stress is effective in coherently organizing and interpreting rheological data from laboratory measurements and computational studies in nonviscometric steady flows of great importance for practical applications. The decomposition of the stress in terms with clearly distinct roles is also useful in developing constitutive models.We introduce a general decomposition of the stress tensor for incompressible fluids in terms of its components on a tensorial basis adapted to the local flow conditions, which include extensional flows, simple shear flows, and any type of mixed flows. Such a basis is determined solely by the symmetric part of the velocity gradient and allows for a straightforward interpretation of the non-Newtonian response in any local flow conditions. In steady homogeneous flows, the material functions that represent the components of the stress on the adapted basis generalize and complete the classical set of viscometric functions used to characterize the response in simple shear flows. Such a general decomposition of the stress is effective in coherently organizing and interpreting rheological data from laboratory measurements and computational studies in nonviscometric steady flows of great importance for practical applications. The decomposition of the stress in terms with clearly distinct roles is also useful in de...

Viscosity of Rigid and Breakable Aggregate Suspensions Stokesian Dynamics for Rigid Aggregates

Suspensions of rigid aggregates have been investigated by Stokesian dynamics. In our recent work (Seto R, Botet R, Briesen H. Phys Rev E 84:041405, 2011), the motions of freely suspended aggregates in shear flows were determined by considering the force and torque balance, and forces and moments acting on the contact points within aggregates were also evaluated. Here, by comparing the obtained results with a bond strength between particles, the sustainable sizes of the aggregates under shear flows have been estimated, which leads to a viscosity-shear rate relation. Our method allows us to see not only the power-law shear thinning for fractal aggregates but also some deviations due to the finite-size effects.