Cross-streamline migration of a semiflexible polymer in a pressure driven flow

Cross-streamline migration and near-wall depletion of elastic fibers in micro-channel flows

The complex dynamics of elastic fibers in viscous fluids are central to many biological and industrial systems. Fluid-structure interactions underlying these dynamics govern the shape and transport of flexible fibers, and understanding these interactions can help tune flow properties in applications such as microfluidic separation, printing and clogging. In this work, we use slender-body theory to study micromechanical dynamics that arise from the coupling between the elastic backbone of a fiber and the local straining flow that contributes to filament flipping and cross-streamline migration. The resulting transverse drift is unbiased in either direction in simple shear flow. However, a non-uniform shear rate results in bias towards regions of high shear, which we connect to the shape transitions during flips. We discover a depletion layer that forms near the boundaries of pressure-driven channel flow due to the competition between such a cross-streamline drift and steric exclusion from the walls. Finally, we develop scaling laws for the curvature of filaments during flip events, demonstrating the origin of the drift bias in non-uniform flows, and confirm this behavior from our simulations. Put together, these results shed light on the role of a local and dominant coupling between elasticity and viscous resistance in dictating long-term dynamics and transport of elastic fibers in confined flows.

Conformations and dynamic behaviors of confined wormlike chains in a pressure-driven flow

The conformations and dynamic behaviors of wormlike chains confined by a slit in a pressure-driven flow were investigated using dissipative particle dynamics method. The wormlike chains exhibit varying conformations due to the varying shear stresses across the slit. The wormlike chain solution can be well described by the power-law fluid, and the power-law index decreases with the increase in chain rigidity. We also presented that the wormlike chain undergoes tumbling motion in the vicinity of the wall in the presence of pressure-driven flow. We also found that the wormlike chains can migrate both away from the wall and slightly away from the slit center, and the migration away from the slit center increases as the chain rigidity is increased because of hydrodynamic interactions induced in a more rigid wormlike chain.

Stability and Deformation of Vesicles in a Cylindrical Flow

In this work, we used dissipative particle dynamics to study the stability, deformation, and rupture of polymer vesicles confined in cylindrical channels under the flow field. The morphological evolution, elongation, and rupture of vesicles and the corresponding mechanisms were intensively investigated. Bullet-like vesicles, leaking vesicles, spherical micelles, hamburger-like micelles, and bilayers were observed by changing the degree of confinement and dimensionless shear rate. We found that increasing the dimensionless shear rate and the degree of confinement can cause the deformation or rupture of polymeric vesicles. The asphericity parameter was utilized to describe the degree of elongation of vesicles deviating from the sphere in the direction of the flow. The results show that the aggregates are more likely to be spherical when the confinement is weak, while they become elongated bullet-like shapes when the confinement is strong. The investigation of dynamics reveals that the degree of confinement and the dimensionless shear rate can affect the chain stretching and reorganization during the process of vesicle elongation. Furthermore, the rupture time of the vesicle shows a nonlinear decrease with an increase in the dimensionless shear rate, and the confinement also contributes to the rupture. The results are very useful for guiding the application of vesicles in a flow environment.

Migration of active filaments under Poiseuille flow in a microcapillary tube

We present a comprehensive study of active filaments confined in a cylindrical channel under Poiseuille flow. The activity drives the filament towards the channel boundary, whereas external fluid flow migrates the filament away from the boundary. This migration further shifts towards the centre for higher flow strength. The migration behaviour of the filaments is presented in terms of the alignment order parameter that shows the alignment grows with shear and activity. Further, we have also addressed the role of length of filament on the migration behaviour, which suggests higher migration for larger filaments. Moreover, we discuss the polar ordering of filaments as a function of distance from the centre of channel that displays upstream motion near the boundary and downstream motion at the centre of the tube.

Understanding conformational and dynamical evolution of semiflexible polymers in shear flow

A clear description of the conformational and dynamical evolution of polymer chains in shear flow is the fundamental basis of microfluidic separations and macroscopic rheological behaviors. We employ graph theory analysis to analyze the local deformation and dynamics of linear polymer chains with different rigidities in shear flow based on the simulation trajectories that record the instantaneous conformations and dynamics. Our results show that all semiflexible chains experience quasi-periodic tumbling motions when the shear strain overwhelms the U-shape (or S-shape) deformation energy barrier. More interestingly, the contact map provides solid evidence for the asymmetric deformation in the whole tumbling motion. In the stretching process: at small and intermediate shear strength, flexible polymers show a quasi-affine deformation while semiflexible ones are initially unfolded from the center of the chains, then both of them follow the extension with half dumbbell- or dumbbell-like ends; at high shear strength, all polymer chains present only a dumbbell-like extension. In the collapse process, all chains prefer to initiate the folding from chain ends. This finding can facilitate our understanding on how semiflexible polymer chains relax and dissipate the stress in shear flow.

Polymer margination in uniform shear flows

We address the issue of polymer margination (migration towards surfaces) in uniform shear flows through extensive LBMD (lattice-Boltzmann molecular dynamics) simulations. In particular we consider the effect of monomer size, a on the chain's overall margination tendency for chains of length N = 16, 32 monomers in flows at multiple shear rates [small gamma, Greek, dot above]. We observed higher margination of chains with larger radii monomers in comparison to smaller radii monomer chains of the same length N. We quantify this effect by considering various measures such as the distribution of the maximum extent of the chain into the channel bulk, zm, distribution of its center of mass in the direction normal to the surface, zc and the distributions of the chain's radius of gyration in directions parallel and perpendicular to the surface i.e. Rx, Ry and Rz respectively.

Flow Behavior of Chain and Star Polymers and Their Mixtures

Star-shaped polymers show a continuous change of properties from flexible linear chains to soft colloids, as the number of arms is increased. To investigate the effect of macromolecular architecture on the flow properties, we employ computer simulations of single chain and star polymers as well as of their mixtures under Poiseuille flow. Hydrodynamic interactions are incorporated through the multi-particle collision dynamics (MPCD) technique, while a bead-spring model is used to describe the polymers. For the ultradilute systems at rest, the polymers are distributed uniformly in the slit channel, with a weak dependence on their number of arms. Once flow is applied, however, we find that the stars migrate much more strongly towards the channel center as the number of arms is increased. In the star-chain mixtures, we find a flow-induced separation between stars and chains, with the stars located in the channel center and the chains closer to the walls. In order to identify the origin of this flow-induced partitioning, we conduct additional simulations without hydrodynamic interactions, and find that the observed cross-stream migration originates from a combination of wall-induced hydrodynamic lift forces and viscoelastic effects. The results from our study give valuable insights for designing microfluidic devices for separating polymers based on their architecture.

Topology-Sensitive Microfluidic Filter for Polymers of Varying Stiffness

The separation of polymers based on their size, rigidity, and topology is an essential but also highly challenging task for nanoscience and engineering. Using hybrid molecular dynamics simulations that correctly take into account hydrodynamics, we have designed microfluidic channels for separating linear from ring polymers in dilute solutions. We establish that the transport velocity of the polymers is independent of their topology and rigidity when the channel walls are smooth and repulsive. However, when the walls are decorated with attractive spots arranged on lines parallel to the flow, ring polymers exhibit an order of magnitude higher transport velocity compared to linear chains. The spots induce a homeotropic-like reorientation of ring polymers close to walls leading to a tank treading motion along them, whereas linear chains are immobilized upon adsorption. This mechanism becomes more enhanced with increasing polymer rigidity. The presented technique holds thus promise for reliably separating nanoparticles based on their topology.

Lateral migration of electrospun hydrogel nanofilaments in an oscillatory flow

The recent progress in bioengineering has created great interest in the dynamics and manipulation of long, deformable macromolecules interacting with fluid flow. We report experimental data on the cross-flow migration, bending, and buckling of extremely deformable hydrogel nanofilaments conveyed by an oscillatory flow into a microchannel. The changes in migration velocity and filament orientation are related to the flow velocity and the filament's initial position, deformation, and length. The observed migration dynamics of hydrogel filaments qualitatively confirms the validity of the previously developed worm-like bead-chain hydrodynamic model. The experimental data collected may help to verify the role of hydrodynamic interactions in molecular simulations of long molecular chains dynamics.

Dynamics of flexible fibers and vesicles in Poiseuille flow at low Reynolds number

The dynamics of flexible fibers and vesicles in unbounded planar Poiseuille flow at low Reynolds number is shown to exhibit similar basic features, when their equilibrium (moderate) aspect ratio is the same and vesicle viscosity contrast is relatively high. Tumbling, lateral migration, accumulation and shape evolution of these two types of flexible objects are analyzed numerically. The linear dependence of the accumulation position on relative bending rigidity, and other universal scalings are derived from the local shear flow approximation.

Dynamics of flexible fibers in shear flow

Dynamics of flexible non-Brownian fibers in shear flow at low-Reynolds-number are analyzed numerically for a wide range of the ratios A of the fiber bending force to the viscous drag force. Initially, the fibers are aligned with the flow, and later they move in the plane perpendicular to the flow vorticity. A surprisingly rich spectrum of different modes is observed when the value of A is systematically changed, with sharp transitions between coiled and straightening out modes, period-doubling bifurcations from periodic to migrating solutions, irregular dynamics, and chaos.

Hydrodynamic segregation in a bidisperse colloidal suspension in microchannel flow: A theoretical study

Colloids in suspension exhibit shear-induced migration towards regions of low viscous shear. In dense bidisperse colloidal suspensions under pressure driven flow large particles can segregate in the center of a microchannel and the suspension partially demixes. To develop a theoretical understanding of these effects, we formulate a phenomenological model for the particle currents based on the work of Phillips et al. [Phys. Fluids 4, 30 (1992)]. We also simulate hard spheres under pressure-driven flow in two and three dimensions using the mesoscale simulation technique of multi-particle collision dynamics. Using a single fit parameter for the intrinsic diffusivity, our theory accurately reproduces the simulated density profiles across the channel. We present a detailed parameter study on how a monodisperse suspension enriches the channel center and quantitatively confirm the experimental observation that a binary colloidal mixture partially segregates into its two species. In particular, we always find a strong accumulation of large particles in the center. Qualitative differences between two and three dimensions reveal that collective diffusion is more relevant in two dimensions.

Lateral migration of flexible fibers in Poiseuille flow between two parallel planar solid walls

Dynamics of non-Brownian flexible fibers in Poiseuille flow between two parallel planar solid walls is evaluated from the Stokes equations which are solved numerically by the multipole method. Fibers migrate towards a critical distance from the wall zc, which depends significantly on the fiber length N and bending stiffness A. This effect can be used to sort fibers. Three types of accumulation are found, depending on a shear-to-bending parameter Γ. In the first type, stiff fibers deform only a little and accumulate close to the wall, where their tendency to drift away from the channel is balanced by the repulsive hydrodynamic interaction with the wall. In the second type, flexible fibers deform significantly and accumulate far from the wall. In both types, the fiber shapes at the accumulation positions are repeatable, while in the third type, they are very compact and non-repeatable. The difference between the second and third accumulation types is a special case of the difference between the regular and irregular modes for the dynamics of migrating fibers. At the regular mode, far from walls, the fiber tumbling frequency satisfies Jeffery's expression, with the local shear rate and the aspect ratio close to N.

Mobility Gradient Induces Cross-Streamline Migration of Semiflexible Polymers

Many aspects of modern material science and biology rely on the strategic manipulation and understanding of polymer dynamics in confining micro- and nanoflow. We directly observe and analyze nonequilibrium structural and dynamic properties of individual semiflexible actin filaments in pressure-driven microfluidic channel flow using fluorescence microscopy. Different conformational shapes, such as filaments fluctuating in an elongated manner, parabolically bent, as well as tumbling, are identified. With increasing flow velocity, a strong center-of-mass migration toward the channel walls is observed. This significant migration effect can be explained by a shear rate dependent spatial diffusivity due to a gradient in chain mobility of the semiflexible polymers.

Shape diagram of vesicles in Poiseuille flow

Soft bodies flowing in a channel often exhibit parachutelike shapes usually attributed to an increase of hydrodynamic constraint (viscous stress and/or confinement). We show that the presence of a fluid membrane leads to the reverse phenomenon and build a phase diagram of shapes-which are classified as bullet, croissant, and parachute-in channels of varying aspect ratio. Unexpectedly, shapes are relatively wider in the narrowest direction of the channel. We highlight the role of flow patterns on the membrane in this response to the asymmetry of stress distribution.