No Arabic abstract
The tumbling of a rigid rod in a shear flow is analyzed in the high viscosity limit. Following Burgers, the Master Equation is derived for the probability distribution of the orientation of the rod. The equation contains one dimensionless number, the Weissenberg number, which is the ratio of the shear rate and the orientational diffusion constant. The equation is solved for the stationary state distribution for arbitrary Weissenberg numbers, in particular for the limit of high Weissenberg numbers. The stationary state gives an interesting flow pattern for the orientation of the rod, showing the interplay between flow due to the driving shear force and diffusion due to the random thermal forces of the fluid. The average tumbling time and tumbling frequency are calculated as a function of the Weissenberg number. A simple cross-over function is proposed which covers the whole regime from small to large Weissenberg numbers.
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.
The intrinsic viscosity of a dilute dispersion of rigid rods is studied using a recently developed direct numerical simulation (DNS) method for particle dispersions. A reentrant transition from shear-thinning to the 2nd Newtonian regime is successfully reproduced in the present DNS results around a Peclet number ${rm Pe}=150$, which is in good agreement with our theoretical prediction of ${rm Pe}=143$, at which the dynamical crossover from Brownian to non-Brownian behavior takes place in the rotational motion of the rotating rod. The viscosity undershoot is observed in our simulations before reaching the 2nd Newtonian regime. The physical mechanisms behind these behaviors are analyzed in detail.
As 2D materials such as graphene, transition metal dichalcogenides, and 2D polymers become more prevalent, solution processing and colloidal-state properties are being exploited to create advanced and functional materials. However, our understanding of the fundamental behavior of 2D sheets and membranes in fluid flow is still lacking. In this work, we perform numerical simulations of athermal semiflexible sheets with hydrodynamic interactions in shear flow. For sheets initially oriented in the flow-gradient plane, we find buckling instabilities of different mode numbers that vary with bending stiffness and can be understood with a quasi-static model of elasticity. For different initial orientations, chaotic tumbling trajectories are observed. Notably, we find that sheets fold or crumple before tumbling but do not stretch again upon applying greater shear.
We present numerical results for the dynamics of a single chain in steady shear flow. The chain is represented by a bead-spring model, and the smoothed profile method is used to accurately account for the effects of thermal fluctuations and hydrodynamic interactions acting on beads due to host fluids. It is observed that the chain undergoes tumbling motions and that its dimensionless frequency F depends only on the Peclet number Pe with a power law. The exponent of Pe clearly changes from 2/3 to 1 around the critical Peclet number, indicating that the crossover reflects the competition of thermal fluctuation and shear flow. The presented numerical results agree well with our theoretical analysis based on Jefferys work.
We study the linear response of interacting active Brownian particles in an external potential to simple shear flow. Using a path integral approach, we derive the linear response of any state observable to initiating shear in terms of correlation functions evaluated in the unperturbed system. For systems and observables which are symmetric under exchange of the $x$ and $y$ coordinates, the response formula can be drastically simplified to a form containing only state variables in the corresponding correlation functions (compared to the generic formula containing also time derivatives). In general, the shear couples to the particles by translational as well as rotational advection, but in the aforementioned case of $xy$ symmetry only translational advection is relevant in the linear regime. We apply the response formulas analytically in solvable cases and numerically in a specific setup. In particular, we investigate the effect of a shear flow on the morphology and the stress of $N$ confined active particles in interaction, where we find that the activity as well as additional alignment interactions generally increase the response.