No Arabic abstract
Hydrodynamic interactions as modeled by Multi-Particle Collision Dynamics can dramatically influence the dynamics of fully flexible, ring-shaped polymers in ways not known for any other polymer architecture or topology. We show that steady shear leads to an inflation scenario exclusive to ring polymers, which depends not only on Weissenberg number but also on contour length of the ring. By analyzing velocity fields of the solvent around the polymer, we show the existence of a hydrodynamic pocket which allows the polymer to self-stabilize at a certain alignment angle to the flow axis. This self-induced stabilization is accompanied by transitioning of the ring to a non-Brownian particle and a cessation of tumbling. The ring swells significantly in the vorticity direction, and the horseshoe regions on the stretched and swollen ring are effectively locked in place relative to the rings center-of-mass. The observed effect is exclusive to ring polymers and stems from an interplay between hydrodynamic interactions and topology. Furthermore, knots tied onto such rings can serve as additional stabilization anchors. Under strong shear, the knotted section is pulled tight and remains well-localized while tank-treading from one horseshoe region to the opposite one in sudden bursts. We find knotted polymers of high contour length behave very similarly to unknotted rings of the same contour length, but small knotted rings feature a host of different configurations. We propose a filtering technique for rings and chains based on our observations and suggest that strong shear could be used to tighten knots on rings.
Molecular dynamics simulations confirm recent extensional flow experiments showing ring polymer melts exhibit strong extension-rate thickening of the viscosity at Weissenberg numbers $Wi<<1$. Thickening coincides with the extreme elongation of a minority population of rings that grows with $Wi$. The large susceptibility of some rings to extend is due to a flow-driven formation of topological links that connect multiple rings into supramolecular chains. Links form spontaneously with a longer delay at lower $Wi$ and are pulled tight and stabilized by the flow. Once linked, these composite objects experience larger drag forces than individual rings, driving their strong elongation. The fraction of linked rings generated by flow depends non-monotonically on $Wi$, increasing to a maximum when $Wisim1$ before rapidly decreasing when the strain rate approaches the relaxation rate of the smallest ring loops $sim 1/tau_e$.
Red blood cells (RBCs) -- erythrocytes -- suspended in plasma tend to aggregate and form rouleaux. During aggregation the first stage consists in the formation of RBC doublets [Blood cells, molecules, and diseases 25, 339 (1999)]. While aggregates are normally dissociated by moderate flow stresses, under some pathological conditions the aggregation becomes irreversible, which leads to high blood viscosity and vessel occlusion. We perform here two-dimensional simulations to study the doublet dynamics under shear flow in different conditions and its impact on rheology. We sum up our results on the dynamics of doublet in a rich phase diagram in the parameter space (flow strength, adhesion energy) showing four different types of doublet configurations and dynamics. We find that membrane tank-treading plays an important role in doublet disaggregation, in agreement with experiments on RBCs. A remarkable feature found here is that when a single cell performs tumbling (by increasing vesicle internal viscosity) the doublet formed due to adhesion (even very weak) remains stable even under a very strong shear rate. It is seen in this regime that an increase of shear rate induces an adaptation of the doublet conformation allowing the aggregate to resist cell-cell detachment. We show that the normalized effective viscosity of doublet suspension increases significantly with the adhesion energy, a fact which should affect blood perfusion in microcirculation.
We investigate the chain conformation of ring polymers confined to a cylindrical nanochannel using both theoretical analysis and three dimensional Langevin dynamics simulations. We predict that the longitudinal size of a ring polymer scales with the chain length and the diameter of the channel in the same manner as that for linear chains based on scaling analysis and Flory-type theory. Moreover, Flory-type theory also gives the ratio of the longitudinal sizes for a ring polymer and a linear chain with identical chain length. These theoretical predictions are confirmed by numerical simulations. Finally, our simulation results show that this ratio first decreases and then saturates with increasing the chain stiffness, which has interpreted the discrepancy in experiments. Our results have biological significance.
BiconeDrag is a software package that allows one to perform a flow field based data processing of dynamic interfacial rheology data pertaining to surfactant laden air-fluid interfaces obtained by means of a rotational bicone shear rheometer. MATLAB and Pyth
The tumbling dynamics of individual polymers in semidilute solution is studied by large-scale non-equilibrium mesoscale hydrodynamic simulations. We find that the tumbling time is equal to the non-equilibrium relaxation time of the polymer end-to-end distance along the flow direction and strongly depends on concentration. In addition, the normalized tumbling frequency as well as the widths of the alignment distribution functions for a given concentration-dependent Weissenberg number exhibit a weak concentration dependence in the cross-over regime from a dilute to a semidilute solution. For semidilute solutions a universal behavior is obtained. This is a consequence of screening of hydrodynamic interactions at polymer concentrations exceeding the overlap concentration.