Institute of Fundamental Technological Research

The long-time behavior of highly elastic fibers in a shear flow is investigated experimentally and numerically. Characteristic attractors of the dynamics are found. It is shown that for a small ratio of bending to hydrodynamic forces, most fibers form a spinning elongated double helix, performing an effective Jeffery orbit very close to the vorticity direction. Recognition of these oriented shapes, and how they form in time, may prove useful in the future for understanding the time history of complex microstructures in fluid flows and considering processing steps for their synthesis.

Stokes equations, shear flow, elastic fibers

Three-dimensional dynamics of flexible fibers in shear flow are studied numerically, with a qualitative comparison to experiments. Initially, the fibers are straight, with different orientations with respect to the flow. By changing the rotation speed of a shear rheometer, we change the ratio A of bending to shear forces. We observe fibers in the flow-vorticity plane, which gives insight into the motion out of the shear plane. The numerical simulations of moderately flexible fibers show that they rotate along effective Jeffery orbits, and therefore the fiber orientation rapidly becomes very close to the flow-vorticity plane, on average close to the flow direction, and the fiber remains in an almost straight configuration for a long time. This ``ordering'' of fibers is temporary since they alternately bend and straighten out while tumbling. We observe numerically and experimentally that if the fibers are initially in the compressional region of the shear flow, they can undergo a compressional buckling, with a pronounced deformation of shape along their whole length during a short time, which is in contrast to the typical local bending that originates over a long time from the fiber ends. We identify differences between local and compressional bending and discuss their competition, which depends on the initial orientation of the fiber and the bending stiffness ratio A. There are two main finding. First, the compressional buckling is limited to a certain small range of the initial orientations, excluding those from the flow-vorticity plane. Second, since fibers straighten out in the flow-vorticity plane while tumbling, the compressional buckling is transient - it does not appear for times longer than 1/4 of the Jeffery period. For larger times, bending of fibers is always driven by their ends.

Stokes flow, flexible fibers, bending, buckling, orientational order

We present a numerical study of the dynamics of an elastic fibre in a shear flow at low Reynolds number, and seek to understand several aspects of the fibre's motion using the equations for slender-body theory coupled to the elastica. The numerical simulations are performed in the bead-spring framework including hydrodynamic interactions in two theoretical schemes: the generalized Rotne-Prager-Yamakawa model and a multipole expansion corrected for lubrication forces. In general, the two schemes yield similar results, including for the dominant scaling features of the shape that we identify. In particular, we focus on the evolution of an initially straight fibre oriented in the flow direction and show that the time scales of fibre bending, curling and rotation, which depend on its length and stiffness, determine the overall motion and evolution of the shapes. We document several characteristic time scales and curvatures representative of the shape that vary as power laws of the bending stiffness and fibre length. The numerical results are further supported by an interpretation using an elastica model.

Stokesian dynamics, particle/fluid flow, slender-body theory

The three-dimensional dynamics of a single non-Brownian flexible fiber in shear flow is evaluated numerically, in the absence of inertia. A wide range of ratios A of bending to hydrodynamic forces and hundreds of initial configurations are considered. We demonstrate that flexible fibers in shear flow exhibit much more complicated evolution patterns than in the case of extensional flow, where transitions to higher-order modes of characteristic shapes are observed when A exceeds consecutive threshold values. In shear flow, we identify the existence of an attracting steady configuration and different attracting periodic motions that are approached by long-lasting rolling, tumbling, and meandering dynamical modes, respectively. We demonstrate that the final stages of the rolling and tumbling modes are effective Jeffery orbits, with the constant parameter C replaced by an exponential function that either decays or increases in time, respectively, corresponding to a systematic drift of the trajectories. In the limit of C→0, the fiber aligns with the vorticity direction and in the limit of C→∞, the fiber periodically tumbles within the shear plane. For moderate values of A, a three-dimensional meandering periodic motion exists, which corresponds to intermediate values of C. Transient, close to periodic oscillations are also detected in the early stages of the modes.

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.

Poiseuille flow, Stokes equations, vesicles, flexible fibers

Fluorescence anisotropy decay measurements and all atom molecular dynamics simulations are used to characterize the orientational motion and preferential interaction of a peptide, N-acetyl-tryptophan-amide (NATA) containing two peptide bonds, in aqueous, urea, guanidinium chloride (GdmCl), and proline solution. Anisotropy decay measurements as a function of temperature and concentration showed moderate slowing of reorientations in urea and GdmCl and very strong slowing in proline solution, relative to water. These effects deviate significantly from simple proportionality of peptide tumbling time to solvent viscosity, leading to the investigation of microscopic preferential interaction behavior through molecular dynamics simulations. Examination of the interactions of denaturants and osmolyte with the peptide backbone uncovers the presence of strongest interaction with urea, intermediate with proline, and weakest with GdmCl. In contrast, the strongest preferential solvation of the peptide side chain is by the nonpolar part of the proline zwitterion, followed by urea, and GdmCl. Interestingly, the local density of urea around the side chain is higher, but the GdmCl distribution is more organized. Thus, the computed preferential solvation of the side chain by the denaturants and osmolyte can account for the trend in reorientation rates. Analysis of water structure and its dynamics uncovered underlying differences between urea, GdmCl, and proline. Urea exerted the smallest perturbation of water behavior. GdmCl had a larger effect on water, slowing kinetics and stabilizing interactions. Proline had the largest overall interactions, exhibiting a strong stabilizing effect on both water–water and water–peptide hydrogen bonds. The results for this elementary peptide system demonstrate significant differences in microscopic behavior of the examined solvent environments. For the commonly used denaturants, urea tends to form disorganized local aggregates around the peptide groups and has little influence on water, while GdmCl only forms specific interactions with the side chain and tends to destabilize water structure. The protective osmolyte proline has the strongest and most specific interactions with the tryptophan side chain, and also stabilizes both water–water and water–peptide hydrogen bonds. Our results strongly suggest protein or peptide denaturation triggered by urea occurs by direct interaction, whereas GdmCl interacts favorably with side chains and destabilizes peptide–water hydrogen bonds. The stabilization of biopolymers by an osmolyte such as proline is governed by favorable preferential interaction with the side chains and stabilization of water.

molecular dynamics simulations, fluorescence anisotropy, peptides, orientational motion

We use numerical simulations of a bead–spring model chain to investigate the evolution of the conformations of long and flexible elastic fibers in a steady shear flow. In particular, for rather open initial configurations, and by varying a dimensionless elastic parameter, we identify two distinct conformational modes with different final size, shape, and orientation. Through further analysis we identify slipknots in the chain. Finally, we provide examples of initial configurations of an 'open' trefoil knot that the flow unknots and then knots again, sometimes repeating several times.

knots, low Reynolds number flows, multipole method

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.

Shear flows, Chaotic dynamics, Vortex dynamics, Numerical solutions, Periodic solutions

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.

Self-organisation, Supramolecular assemblies

Dynamics of single flexible non-Brownian fibers, tumbling in a Poiseuille flow between two parallel solid plane walls, is studied with the use of the HYDROMULTIPOLE numerical code, based on the multipole expansion of the Stokes equations, corrected for lubrication. Fibers, which are closer to a wall, more flexible (less stiff) or longer, deform more significantly and, for a wide range of the system parameters, they faster migrate towards the middle plane of the channel. For the considered systems, fiber velocity along the flow is only slightly smaller than (and can be well approximated by) the Poseuille flow velocity at the same position. In this way, the history of a fiber migration across the channel is sufficient to determine with a high accuracy its displacement along the flow.

Stokes equations, flexible fiber, Poiseuille flow, solid walls

Hydrodynamic properties of fibrinogen molecules were theoretically calculated. Their shape was approximated by the bead model, considering the presence of flexible side chains of various length and orientation relative to the main body of the molecule. Using the bead model, and the precise many-multipole method of solving the Stokes equations, the mobility coefficients for the fibrinogen molecule were calculated for arbitrary orientations of the arms whose length was varied between 12 and 18 nm. Orientation averaged hydrodynamic radii and intrinsic viscosities were also calculated by considering interactions between the side arms and the core of the fibrinogen molecule. Whereas the hydrodynamic radii changed little with the interaction magnitude, the intrinsic viscosity exhibited considerable variation from 30 to 60 for attractive and repulsive interactions, respectively. These theoretical results were used for the interpretation of experimental data derived from sedimentation and diffusion coefficient measurements as well as dynamic viscosity measurements. Optimum dimensions of the fibrinogen molecule derived in this way were the following: the contour length 84.7 nm, the side arm length 18 nm, and the total volume 470 nm3, which gives 16% hydration (by volume). Our calculations enabled one to distinguish various conformational states of the fibrinogen molecule, especially the expanded conformation, prevailing for pH < 4 and lower ionic strength, characterized by high intrinsic viscosity of 50 and the hydrodynamic radius of 10.6 nm. On the other hand, for the physiological condition, that is, pH = 7.4 and the ionic strength of 0.15 M NaCl, the semi-collapsed conformation dominates. It is characterized by the average angle equal to = 55, intrinsic viscosity of 35, and the hydrodynamic radius of 10 nm. Additionally, the interaction energy between the arms and the body of the molecule was predicted to be 4 kT units, confirming that they are oppositely charged than the central nodule. Results obtained in our work confirm an essential role of the side chains responsible for a highly anisotropic charge distribution in the fibrinogen molecule. These finding can be exploited to explain anomalous adsorption of fibrinogen on various surfaces.

The expansion of a plume generated during laser ablation is studied with the Direct Simulation Monte Carlo method. The plume is a mixture of four disparate molecular mass components and expands in vacuum or into ambient gas. The time dependence of deposition rate is studied and the transition from an initial vacuum-like to a diffusion-like regime of expansion in ambient gas is shown. The lack of stoichiometry increases with the ratio of molecular masses of ablated particles and at disparate masses the stoichiometry is seriously affected. Ambient gas worsens the stoichiometry unless it supplies particles compensating the backward and sideward flows of plume constituents.

laser deposition, plume expansion, DSMC

The paper presents the results of numerical simulation of processes aimed at production of nanostructures with the use of oil emulsions in water. The appropriate molecular models of water and oil, as well as the model of the substance which would sediment at the water – oil interface, are looked for. Such substance, after suitable solidification, would become the main component of the produced material. For the described simulations, the Molecular Dynamics method has been used throughout this paper.

thin layers, contact surface, nonmixing liquids

Współczesne technologie materiałowe są jedną z najszybciej rozwijających się dziedzin nauki i techniki. Szczególnie prężnie rozwijaj się nano- i mikrotechnologie. Jedna z takich nowoczesnych nanotechnologii, badana obecnie w kilku europejskich ośrodkach, wykorzystuje efekt gromadzenia się substancji na granicy faz emulsji. Użyta emulsja ma bardzo drobną strukturę nano-kropli oleju w wodzie. Trzecia substancja, dzięki odpowiednim właściwościom molekularnym, osiadając na powierzchni styku faz cieczy pokrywa powierzchnie kropel oleju. Po usunięciu emulsji substancja ta, zachowując strukturę, zostaje utwardzona i tworzy nano-materiał. W prezentowanej technologii emulsja spełnia rolę matrycy, na której powstaje struktura wytwarzanego materiału. Technik rozdrabniania emulsji jest wiele; stosuje się m.in. aparaty miksujące (homogenizatory), które w przepływie ścinającym rozrywają krople oleju na mniejsze lub mikrokanały o podobnym działaniu [2]. Nanomateriały o prezentowanej strukturze będą miały wiele interesujących właściwości, takich jak lekkość, elastyczność czy wytrzymałość mechaniczną, co zapewniają silne oddziaływania międzyatomowe w układzie oraz porowatość substancji. Tworzywa o takiej budowie mogą znaleźć zastosowanie w różnych dziedzinach np. medycynie czy aerodynamice. Celem naszych prac nad omawianą technologią wytwarzania nanomateriałów było stworzenie numerycznego modelu zjawiska powstawania nanostruktur w emulsjach. Ponieważ tak drobne układy wymagają modelowania na poziomie atomowym, do opisu procesów zachodzących w cieczach posłużono się metodą Dynamiki Molekularnej. Bazując na symulacjach numerycznych zbudowano modele molekularne cieczy tworzących emulsję oraz zaproponowano kilka typów substancji, które mogłyby wytworzyć pożądaną warstwę na granicy faz. Zaproponowane modele testowano numerycznie, poszukując kombinacji oddziaływań międzyatomowych zapewniającej powstawanie oczekiwanej nanostruktury

nanostruktury, emulsje, nanomateriały

The very serious problem connected with long distance transport of gase-
ous fuels is connected with the fact that detonation may occur inside the duct if some
air leaks into it. Detonation is particularly dangerous for compressors which “push”
the gas through pipelines. It is therefore necessary to use some “detonation dampers”
to protect these compressors.
The commonly used detonation damper has a form of a matrix of narrow chan-
nels, placed across the pipe transporting gas. Detonation wave is supposed to be ex-
tinguished due to cooling by cold walls of these channels. To achieve efficient
damping the channels should be very narrow. Our earlier simulations [1, 2, 3] indi-
cate, that preferable widths should be of the order of 0.005mm. This is not acceptable
for practical reasons – in real applications it is close to 0.5mm. In such channels cool-
ing the gas by heat transfer is inefficient and cannot extinguish the flame (cooling ef-
fect is proportional to perimeter of the cross-section and amount of gas to be cooled
to cross-section surface).
Rarefaction wave generated behind sharp increase of channel cross-section is
an alternative phenomenon, which may be used to cool the burning gas [1]. In our
earlier papers we tested several channel shapes with increase of cross-section [2, 3].
Here we compare channels 0.5mm wide, without and with sharp increases and de-
creases of cross-section. Such channels were found to be the best for practical appli-
cations.

detonation waves, detonation damping, narrow channels

One of the important contemporary technological problems is connected withnecessity of extinguishing detonations, which may occur inpipelines transporting gaseousfuels. To achieve this goal usually a matrix of narrow channels is placed across the flowinside the pipeline. In our recent papers [1], [2] we have shown, that channels with sharpchanges of cross-section should be more efficient in this respect than traditionally usedstraight channels with constant cross-section area. In the present paper we demonstratehow detonation behaves in channels with changes of cross-section under realistic conditions– if the channel cross-section is of dimensions acceptable technologically. At the same timewe take into account the fact, that if friction and heat exchange at the walls are present,gas flowing through the channels accelerates and its densitydecreases considerably. Theresult of our considerations is a selection of, possibly, optimum shape of the channels ofa detonation damper.

One of the important contemporary technological problems is connected with necessity of extinguishing detonation, which may occur in pipelines transporting gaseous fuels. To achieve this goal usually a matrix of very narrow channels is placed inside the pipeline, perpendicularly to its axis. In our recent paper (Walenta and Slowicka (2016)) we have shown, that channels with sharp changes of cross-section should be more efficient in this respect than traditionally used straight channels with constant cross-section area. In this paper we demonstrate how detonation behaves in the channels, in which gas flows under realistic conditions – when friction and heat exchange are present. We take into account the fact, that gas flowing through such channels accelerates and its density decreases considerably.

detonation waves, detonation damping, narrow channels

In the present paper we show the dependence of the shock structure in a dense, noble gas on each of the three non-dimensional parameters: non-dimensional initial density, non-dimensional initial temperature and non-dimensional shock velocity. It will also be demonstrated, that the length scale, most suitable for measuring the thickness of the shock wave in a dense gas, is the sum of the mean free path (calculated the same way as for a dilute gas) and the diameter of a single gas molecule.

The necessity of extinguishing detonation, which may occur in pipelines transporting gaseous fuels, is nowadays a very important technological problem. The standard devices used for it consist of matrices of very narrow channels. Cooling the gas by cold walls of such channels may extinguish the flame and stop detonation. Detonation may also be extinguished if the cross-section of the channel transporting gas increases abruptly at some place. The desired effect may be achieved if the generated rarefaction waves decrease sufficiently the temperature of the flame (Teodorczyk et al. 1988, Cai et al. 2002, Dremin 1999, Walenta et al. 2004). It might be expected, that simultaneous use of both methods – using narrow channels with variable cross-section – should give even better results. Additional profit might come from the fact, that the flow in narrow channels is usually laminar; the abrupt increase of the cross-section would introduce some turbulence and this way enhance cooling by the walls. However, if the cross-section of the channel increases and decreases, the unwanted heating of the gas may occur. To estimate the net esult, which is not obvious, it is necessary to perform suitable simulations and experiments. The present paper is devoted to numerical simulation of the phenomenon.

extinguishing detonation, DSMC

The present paper is a continuation of our earlier work on the
structure of shock waves in dense media. Simple, monoatomic gas, argon, was con
sidered first [1]. It was found, that the shock thickness in dense argon, when re
lated to the mean distance between the molecules (as suggested by P.W. Bridgm
an [2]) decreased with increasing density. To illustrate the dependence of the thickness of the
shock wave on density of the medium, the results for liquid sulfur hexafluoride were
compared with experimental, shock tube results in gas phase at low density.

molecular dynamics simulations, complex liquids, shock waves

In our earlier paper [3] we reported the investigation of the shock structure in dense, monatomic gas - argon. Here we extend our work to dense molecular gases and to liquids. We investigate, in particular, the influence of the electric charges (electric dipoles, quadrupoles etc.) of the molecules on the shock wave structure.

Shock structure, Dense media, Molecular Dynamics

laser ablation, plume expansion, DSMC

shock waves, dense fluids, molecular dynamics simulations

Stokes equations, Euler-Bernoulli beam, elastica, elastic filament, shear flow, buckling

The research, reported in the present paper, was aimed at optimization of devices
extinguishing detonation, which may occur in pipelines transporting gaseous fuels. Such
device, “detonation damper”, is a matrix of narrow channels placed across the pipe. Detonation
is extinguished by heat exchange between the flame and cold walls of the channels.
Improvement of efficiency of the damper was achieved by modification of the channel shape,
so that rarefaction waves, cooling additionally the burning gas, appear in the flow.

Flow Control and Optimisation, Micro- and Nano- flows, Detonation

fluid dynamics, Stokes flows, particulate flows, hydrodynamic interactions, bending

Dynamics of a single non-Brownian flexible fiber in shear flow at low Reynolds number is investigated numerically. Initially, the fiber is straight and at the equilibrium. For different initial orientations and values of bending stiffness, three generic scenarios are observed: the fiber tends to: align along the vorticity direction, tumble within the plane perpendicular to vorticity, or perform a periodic motion superposed with translation along the flow.

Stokes equations, flexible fibers, shear flow

Dynamics of flexible fibers and vesicles in unbounded planar Poiseuille flow at the low-Reynolds-number are shown to exhibit similar basic features, when their equilibrium (moderate) aspect ratio is the same and vesicle viscosity contrast is relatively high. The lateral migration and accumulation of these two types of flexible objects are analyzed numerically.

Stokes equations, vesicles, flexible fibers

Modes of the dynamics of flexible fibers in shear flow in a plane perpendicular to vorticity are analyzed numerically.

Dynamics of flexible fibers in shear flow are analyzed numerically for a wide range of the ratios A of the fiber bending force to the hydrodynamic force. The Reynolds number is much smaller that unity, and the Péclet number much greater than one. A fiber is modeled as a chain of solid beads as in Ref. [1]. The centers of the consecutive beads are linked by springs. The equilibrium length of each spring is such that the consecutive beads almost touch each other. The spring constant is large, and the fiber’s length practically does not change with time. At the equilibrium, the fiber is straight. At a deformed configuration, there appear a bending force exerted on each bead, proportional to the bending parameter A. The details of the model and the numerical method can be found in Ref. [2]. Velocities of the beads at given positions are evaluated with the use of the HYDROMULTIPOLE numerical code, based on solving the Stokes equations by the multipole expansion [3]. The time-dependent positions of the beads are determined by the adaptive fourth order Runge-Kutta method.

Initially, the fiber is aligned with the flow, and the springs are at the equilibrium. Owing to symmetry, the centers of all the beads move in the plane perpendicular to the vorticity direction. The fiber end-to-end vector tumbles, in a similar way as a rigid elongated body [4]. While the fiber turns, its shape evolves accordingly, and the center-of-mass oscillates across the flow. A surprisingly rich spectrum of different modes is observed when the value of A is systematically changed, with regular and chaotic trajectories. (For the details, see Ref. [5].)

For some ranges of small and large values of A, the center-of-mass trajectories are periodic with a single tumbling time τ, and there is no migration across the flow. For a certain range of small values of A, evolutions during every second tumbling time are the same, but consecutive tumbling times differ from each other, what leads to regular migrating trajectories of the fiber center-of-mass. For moderate values of A, a chaotic behavior is observed - a large sensitivity of the dynamics to a small change of A, with (typically) many irregular, erratic trajectories or (exceptionally) some regular migrating trajectories. At a moderate value of A, a transition is observed between a ‘straightening out’ mode of more stiff fibers to the coiled mode of more flexible fibers. In the straightening out mode, the fiber significantly changes its shape while tumbling - from almost straight and aligned with the flow to S-shaped. In the coiled mode, the fiber is always compact, it never aligns with the flow.

References
[1] E. Gauger, H. Stark, Numerical study of a microscopic artificial swimmer, Phys. Rev. E 74, 021907 (2006).
[2] A. M. Słowicka, E. Wajnryb, M. L. Ekiel-Jeżewska, Lateral migration of flexible fibers in Poiseuille flow between two parallel planar solid walls, Eur. Phys. J. E 36, 1-12 (2013).
[3] B. Cichocki, M. L. Ekiel-Jeżewska, E. Wajnryb, Lubrication corrections for three-particle contribution to short-time self-diffusion coefficients in colloidal dispersions, J. Chem. Phys. 111, 3265 (1999).
[4] G. Jeffery, The Motion of Ellipsoidal Particles Immersed in a Viscous Fluid, Proc. R. Soc. Lond. A. 102, 161-179 (1922).
[5] A. M. Słowicka, E. Wajnryb, M. L. Ekiel-Jeżewska Dynamics of flexible fibers in shear flow, arXiv [cond-mat.soft] (2015).

Stokes equations, flexible fibers, shear flow

Pulsed laser deposition is a method frequently used for creating thin films of various materials on
solid substrates. High energy laser pulse causes evaporation of the target material, forming a
plume which subsequently
expands and moves with high speed from the target. Thin film of the
evaporated material is deposited on the substrate placed at some distance in front of the target.
The behavior of the plume influences both the stoichiometry and homogeneity of the deposit
ed
layer
–
the final product of the process. Better understanding of the process of expansion of the
plume, variation of its structure as well as deposition of the material itself is therefore very
important and should give us opportunity for better contro
l of formation of the deposited layer.

laser ablation, plume expansion, DSMC

shock waves, dense fluids, molecular dynamics simulations

shock waves, dense media, molecular dynamics symulations

In this work, we consider a single non-Brownian mobile and flexible fiber immersed in Poiseuille flow in a channel consisting of two parallel infinite walls. The dynamics of the fiber is evaluated numerically from the Stokes equations by a multipole code HYDROMULTIPOLE. Investigating the fiber dynamics we found out that fibers migrate to a critical position across the channel. The distance between the wall and a limiting position depends on the fiber elongation and flexibility. For more stiff fibers the critical position results from the interplay between their tendency to drift away from the channel and the repulsive hydrodynamic interaction with the wall. For less stiff fibers the limiting position is not influenced by the presence of the wall. Differences between the critical position for different fibers can be used in the process of microfibers separation by the flow.

Stokes equations, hydrodynamic interactions, Poiseuille flow, flecible fiber

shock waves, dense fluids, molecular dynamics simulations

Stokes equations, Poiseuille flow, flexible fibers, migration