No Arabic abstract
We present a thermodynamically consistent constitutive model for fluid-saturated sediments, spanning dense to dilute regimes, developed from the basic balance laws for two phase-mixtures. The model can represent various limiting cases, such as pure fluid and dry grains. It is formulated to capture a number of key behaviors such as: (i) viscous inertial rheology of submerged wet grains under steady shearing flows, (ii) the critical state behavior of grains, which causes granular Reynolds dilation/contraction due to shear, (iii) the viscous thickening of the fluid response due to the presence of suspended grains, and (iv) the Darcy-like drag interaction observed in both dense and dilute mixtures, which gives rise to complex fluid-grain interactions under dilation and flow. The full constitutive model is combined with the basic equations of motion for each mixture phase and implemented in the material point method (MPM) to accurately model the coupled dynamics of the mixed system. Qualitative results show the breadth of problems which this model can address. Quantitative results demonstrate the accuracy of this model as compared with analytical limits and experimental observations of fluid and grain behaviors in inhomogeneous geometries.
Fine particle suspensions (such as cornstarch mixed with water) exhibit dramatic changes in viscosity when sheared, producing fascinating behaviors that captivate children and rheologists alike. Recent examination of these mixtures in simple flow geometries suggests inter-granular repulsion is central to this effect --- for mixtures at rest or shearing slowly, repulsion prevents frictional contacts from forming between particles, whereas, when sheared more forcefully, granular stresses overcome the repulsion allowing particles to interact frictionally and form microscopic structures that resist flow. Previous constitutive studies of these mixtures have focused on particular cases, typically limited to two-dimensional, steady, simple shearing flows. In this work, we introduce a predictive and general, three-dimensional continuum model for this material, using mixture theory to couple the fluid and particle phases. Playing a central role in the model, we introduce a micro-structural state variable, whose evolution is deduced from small-scale physical arguments and checked with existing data. Our space- and time-dependent model is implemented numerically in a variety of unsteady, non-uniform flow configurations where it is shown to accurately capture a variety of key behaviors: (i) the continuous shear thickening (CST) and discontinuous shear thickening (DST) behavior observed in steady flows, (ii) the time-dependent propagation of `shear jamming fronts, (iii) the time-dependent propagation of `impact activated jamming fronts, and (iv) the non-Newtonian, `running on oobleck effect wherein fast locomotors stay afloat while slow ones sink.
Discrete particle simulations are used to study the shear rheology of dense, stabilized, frictional particulate suspensions in a viscous liquid, toward development of a constitutive model for steady shear flows at arbitrary stress. These suspensions undergo increasingly strong continuous shear thickening (CST) as solid volume fraction $phi$ increases above a critical volume fraction, and discontinuous shear thickening (DST) is observed for a range of $phi$. When studied at controlled stress, the DST behavior is associated with non-monotonic flow curves of the steady-state stress as a function of shear rate. Recent studies have related shear thickening to a transition between mostly lubricated to predominantly frictional contacts with the increase in stress. In this study, the behavior is simulated over a wide range of the dimensionless parameters $(phi,tilde{sigma}$, and $mu)$, with $tilde{sigma} = sigma/sigma_0$ the dimensionless shear stress and $mu$ the coefficient of interparticle friction: the dimensional stress is $sigma$, and $sigma_0 propto F_0/ a^2$, where $F_0$ is the magnitude of repulsive force at contact and $a$ is the particle radius. The data have been used to populate the model of the lubricated-to-frictional rheology of Wyart and Cates [Phys. Rev. Lett.{bf 112}, 098302 (2014)], which is based on the concept of two viscosity divergences or textquotedblleft jammingtextquotedblright points at volume fraction $phi_{rm J}^0 = phi_{rm rcp}$ (random close packing) for the low-stress lubricated state, and at $phi_{rm J} (mu) < phi_{rm J}^0$ for any nonzero $mu$ in the frictional state; a generalization provides the normal stress response as well as the shear stress. A flow state map of this material is developed based on the simulation results.
We present a new strategy for introducing population balances into full-chain constitutive models of living polymers with linear chain architectures. We provide equations to describe a range of stress relaxation processes covering both unentangled systems (Rouse-like motion) and well entangled systems (reptation, contour length fluctuations, chain retraction, and constraint release). Special attention is given to the solutions that emerge when the breaking time of the chain becomes fast compared to various stress relaxation processes. In these fast breaking limits, we reproduce previously known results (with some corrections) and also present new results for nonlinear stress relaxation dynamics. Our analysis culminates with a fully developed constitutive model for the fast breaking regime in which stress relaxation is dominated by contour length fluctuations. Linear and nonlinear rheology predictions of the model are presented and discussed.
We investigate the structure of a dilute mixture of amphiphilic dimers and spherical particles, a model relevant to the problem of encapsulating globular guest molecules in a dispersion. Dimers and spheres are taken to be hard particles, with an additional attraction between spheres and the smaller monomers in a dimer. Using Monte Carlo simulation, we document the low-temperature formation of aggregates of guests (clusters) held together by dimers, whose typical size and shape depend on the guest concentration $chi$. For low $chi$ (less than $10%$), most guests are isolated and coated with a layer of dimers. As $chi$ progressively increases, clusters grow in size becoming more and more elongated and polydisperse; after reaching a shallow maximum for $chiapprox 50%$, the size of clusters again reduces upon increasing $chi$ further. In one case only ($chi=50%$ and moderately low temperature) the mixture relaxed to a fluid of lamellae, suggesting that in this case clusters are metastable with respect to crystal-vapor separation. On heating, clusters shrink until eventually the system becomes homogeneous on all scales. On the other hand, as the mixture is made denser and denser at low temperature, clusters get increasingly larger until a percolating network is formed.
We present a general theory for predicting the interaction potentials between DNA-coated colloids, and more broadly, any particles that interact via valence-limited ligand-receptor binding. Our theory correctly incorporates the configurational and combinatorial entropic factors that play a key role in valence-limited interactions. By rigorously enforcing self-consistency, it achieves near-quantitative accuracy with respect to detailed Monte Carlo calculations. With suitable approximations and in particular geometries, our theory reduces to previous successful treatments, which are now united in a common and extensible framework. We expect our tools to be useful to other researchers investigating ligand-mediated interactions. A complete and well-documented Python implementation is freely available at http://github.com/patvarilly/DNACC .