Victor Steinberg

Elastic turbulence in a polymer solution flow

Turbulence is a ubiquitous phenomenon that is not fully understood. It is known that the flow of a simple, newtonian fluid is likely to be turbulent when the Reynolds number is large (typically when the velocity is high, the viscosity is low and the size of the tank is large). In contrast, viscoelastic fluids such as solutions of flexible long-chain polymers have nonlinear mechanical properties and therefore may be expected to behave differently. Here we observe experimentally that the flow of a sufficiently elastic polymer solution can become irregular even at low velocity, high viscosity and in a small tank. The fluid motion is excited in a broad range of spatial and temporal scales, and we observe an increase in the flow resistance by a factor of about twenty. Although the Reynolds number may be arbitrarily low, the observed flow has all the main features of developed turbulence. A comparable state of turbulent flow for a newtonian fluid in a pipe would have a Reynolds number as high as 105 (refs 1, 2). The low Reynolds number or ‘elastic’ turbulence that we observe is accompanied by significant stretching of the polymer molecules, resulting in an increase in the elastic stresses of up to two orders of magnitude.

Efficient mixing at low Reynolds numbers using polymer additives

Mixing in fluids is a rapidly developing area in fluid mechanics, being an important industrial and environmental problem. The mixing of liquids at low Reynolds numbers is usually quite weak in simple flows, and it requires special devices to be efficient. Recently, the problem of mixing was solved analytically for a simple case of random flow, known as the Batchelor regime. Here we demonstrate experimentally that very viscous liquids containing a small amount of high-molecular-weight polymers can be mixed quite efficiently at very low Reynolds numbers, for a simple flow in a curved channel. A polymer concentration of only 0.001% suffices. The presence of the polymers leads to an elastic instability and to irregular flow, with velocity spectra corresponding to the Batchelor regime. Our detailed observations of the mixing in this regime enable us to confirm several important theoretical predictions: the probability distributions of the concentration exhibit exponential tails, moments of the distribution decay exponentially along the flow, and the spatial correlation function of concentration decays logarithmically.

Elastic turbulence in curvilinear flows of polymer solutions.

Following our first report (A Groisman and V Steinberg 2000 Nature 405 53), we present an extended account of experimental observations of elasticity-induced turbulence in three different systems: a swirling flow between two plates, a Couette–Taylor (CT) flow between two cylinders, and a flow in a curvilinear channel (Dean flow). All three set-ups had a high ratio of the width of the region available for flow to the radius of curvature of the streamlines. The experiments were carried out with dilute solutions of high-molecular-weight polyacrylamide in concentrated sugar syrups. High polymer relaxation time and solution viscosity ensured prevalence of non-linear elastic effects over inertial non-linearity, and development of purely elastic instabilities at low Reynolds number (Re) in all three flows. Above the elastic instability threshold, flows in all three systems exhibit features of developed turbulence. They include: (i) randomly fluctuating fluid motion excited in a broad range of spatial and temporal scales and (ii) significant increase in the rates of momentum and mass transfer (compared with those expected for a steady flow with a smooth velocity profile). Phenomenology, driving mechanisms and parameter dependence of the elastic turbulence are compared with those of the conventional high-Re hydrodynamic turbulence in Newtonian fluids. Some similarities as well as multiple principal differences were found. In two out of three systems (swirling flow between two plates and flow in the curvilinear channel), power spectra of velocity fluctuations decayed rather quickly, following power laws with exponents of about −3.5. It suggests that, being random in time, the flow is rather smooth in space, in the sense that the main contribution to deformation and mixing (and, possibly, elastic energy) is coming from flow at the largest scale of the system. This situation, random in time and smooth in space, is analogous to flows at small scales (below the Kolmogorov dissipation scale) in high-Re turbulence.

Continuous particle size separation and size sorting using ultrasound in a microchannel

Continuous separation and size sorting of particles and blood cells suspended in a microchannel flow due to an acoustic force are investigated both numerically and experimentally. Good agreement in the measured particle trajectories in a microchannel flow subjected to the acoustic force with those obtained by the numerical simulations up to the fitting parameter is found. High separation efficiency, particularly in a three-stage microdevice (up to 99.975%), for particles and blood cells leads us to believe that the device can be developed into commercially useful set-up. The novel particle size sorting microdevice provides an opportunity to replace rather expensive existing devices based on specific chemical bonding with an ultrasound cell size sorter that can be considerably improved by adding many stages for multistage size sorting.

Mechanism of elastic instability in Couette flow of polymer solutions: Experiment

Experiments on flow stability and pattern formation in Couette flow between two cylinders with highly elastic polymer solutions are reported. It is found that the flow instabilities are determined by the elastic Deborah number, De, and the polymer concentration only, while the Reynolds number becomes completely irrelevant. A mechanism of such “purely elastic” instability was suggested a few years ago by Larson, Shaqfeh, and Muller [J. Fluid Mech. 218, 573 (1990)], referred to as LMS. It is based on the Oldroyd-B rheological model and implies a certain functional relation between De at the instability threshold and the polymer contribution to the solution viscosity, ηp/η, that depends on the polymer concentration. The elastic force driving the instability arises when perturbative elongational flow in radial direction is coupled to the strong primary azimuthal shear. This force is provided by the “hoop stress” that develops due to stretching of the polymer molecules along the curved streamlines. It is found...

Rotating Rayleigh–Bénard convection: asymmetric modes and vortex states

We present optical shadowgraph flow visualization and heat transport measurements of Rayleigh–Benard convection with rotation about a vertical axis. The fluid, water with Prandtl number 6.4, is confined in a cylindrical convection cell with radius-to-height ratio Γ = 1. For dimensionless rotation rates 150 R c (Ω) much less than those predicted by linear stability analysis for a laterally infinite system and qualitatively consistent with finite-aspect-ratio, linear-stability calculations of Buell & Catton (1983). As in the calculations, the forward bifurcation at onset is to states of localized flow near the lateral walls with azimuthal periodicity of 3 m et al . (1992), with a frequency that is finite at onset but goes to zero as Ω goes to zero. At Ω = 2145 we find primary and secondary stability boundaries for states with m = 4, 5, 6, and 7. Further, we show that at higher Rayleigh number, there is a transition to a vortex state where the vortices form with the symmetry of the existing azimuthal periodicity of the sidewall state. Aperiodic, time-dependent heat transport begins for Rayleigh numbers at or slightly above the first appearance of vortices. Visualization of the formation and interactions of thermal vortices is presented, and the behaviour of the Nusselt number at high Rayleigh numbers is discussed.

Concentration dependence of the longest relaxation times of dilute and semi-dilute polymer solutions

The longest relaxation times of polymer solutions of semi-flexible T4 DNA and flexible 18 M molar mass polyacrylamide (PAAm) in dilute and semi-dilute concentration range are studied by the polymer extension relaxation of stretched single DNA molecules and by the stress relaxation of PAAm solutions measurements. For both polymer solutions, the longest relaxation time normalized by the value at infinite dilution with the same solvent viscosity τ/τ0 increases with increasing concentration. In the dilute regime, the longest relaxation time increases just slightly with increasing concentration as τ/τ0=[1+cA−2(cA)1.5+2(cA)2] as well as the empirical relation of τ/τ0=exp(cA) up to c∼3c∗ with A≈0.5[η], where c∗ is the overlap concentration, in accord with the theory and previous experiments. For the semi-dilute solutions, the scaling of τ/τ0 with concentration shows two different exponents in two concentration regions, corresponding to the unentangled and entangled regimes. The exponents are consistent with thos...

Dynamics of a vesicle in general flow

An approach to quantitatively study vesicle dynamics as well as biologically-related micro-objects in a fluid flow, which is based on the combination of a dynamical trap and a control parameter, the ratio of the vorticity to the strain rate, is suggested. The flow is continuously varied between rotational, shearing, and elongational in a microfluidic 4-roll mill device, the dynamical trap, that allows scanning of the entire phase diagram of motions, i.e., tank-treading (TT), tumbling (TU), and trembling (TR), using a single vesicle even at λ = ηin/ηout = 1, where ηin and ηout are the viscosities of the inner and outer fluids. This cannot be achieved in pure shear flow, where the transition between TT and either TU or TR is attained only at λ>1. As a result, it is found that the vesicle dynamical states in a general are presented by the phase diagram in a space of only 2 dimensionless control parameters. The findings are in semiquantitative accord with the recent theory made for a quasi-spherical vesicle, although vesicles with large deviations from spherical shape were studied experimentally. The physics of TR is also uncovered.

Vertical pairing of identical particles suspended in the plasma sheath.

It is shown experimentally that vertical pairing of two identical microspheres suspended in the sheath of a radio-frequency (rf) discharge at low gas pressures (a few Pa) appears at a well-defined instability threshold of the rf power. The transition is reversible, but with significant hysteresis on the second stage. A simple model which uses measured microsphere resonance frequencies and takes into account, in addition to the Coulomb interaction between negatively charged microspheres, their interaction with positive-ion-wake charges, seems to explain the instability threshold quite well.

Fluid vesicles in flow.

We review the dynamical behavior of giant fluid vesicles in various types of external hydrodynamic flow. The interplay between stresses arising from membrane elasticity, hydrodynamic flows, and the ever present thermal fluctuations leads to a rich phenomenology. In linear flows with both rotational and elongational components, the properties of the tank-treading and tumbling motions are now well described by theoretical and numerical models. At the transition between these two regimes, strong shape deformations and amplification of thermal fluctuations generate a new regime called trembling. In this regime, the vesicle orientation oscillates quasi-periodically around the flow direction while asymmetric deformations occur. For strong enough flows, small-wavelength deformations like wrinkles are observed, similar to what happens in a suddenly reversed elongational flow. In steady elongational flow, vesicles with large excess areas deform into dumbbells at large flow rates and pearling occurs for even stronger flows. In capillary flows with parabolic flow profile, single vesicles migrate towards the center of the channel, where they adopt symmetric shapes, for two reasons. First, walls exert a hydrodynamic lift force which pushes them away. Second, shear stresses are minimal at the tip of the flow. However, symmetry is broken for vesicles with large excess areas, which flow off-center and deform asymmetrically. In suspensions, hydrodynamic interactions between vesicles add up to these two effects, making it challenging to deduce rheological properties from the dynamics of individual vesicles. Further investigations of vesicles and similar objects and their suspensions in steady or time-dependent flow will shed light on phenomena such as blood flow.