Toward a structural understanding of turbulent drag reduction: nonlinear coherent states in viscoelastic shear flows

Coherent Structures in Plane Channel Flow of Dilute Polymer Solutions with Vanishing Inertia

When subjected to sufficiently strong velocity gradients, solutions of long, flexible polymers exhibit flow instabilities and chaotic motion, often referred to as elastic turbulence. Its mechanism differs from the familiar, inertia-driven turbulence in Newtonian fluids and is poorly understood. Here, we demonstrate that the dynamics of purely elastic pressure-driven channel flows of dilute polymer solutions are organized by exact coherent structures that take the form of two-dimensional traveling waves. Our results demonstrate that no linear instability is required to sustain such traveling wave solutions and that their origin is purely elastic in nature. We show that the associated stress profiles are characterized by thin, filamentlike arrangements of polymer stretch, which is sustained by a solitary pair of vortices. We discuss the implications of the traveling wave solutions for the transition to elastic turbulence in straight channels and propose ways for their detection in experiments.

Polymers and Plastrons in Parallel Yield Enhanced Turbulent Drag Reduction

Despite polymer additives and superhydrophobic walls being well known as stand-alone methods for frictional drag reduction in turbulent flows, the possibility of employing them simultaneously in an additive fashion has remained essentially unexplored. Through experimental friction measurements in turbulent Taylor–Couette flow, we show that the two techniques may indeed be combined favorably to generate enhanced levels of frictional drag reduction in wall-bounded turbulence. We further propose an additive expression in Prandtl–von Kármán variables that enables us to quantitatively estimate the magnitude of this cooperative drag reduction effect for small concentrations of dissolved polymer.

Viscoelastic Pipe Flow is Linearly Unstable

Newtonian pipe flow is known to be linearly stable at all Reynolds numbers. We report, for the first time, a linear instability of pressure-driven pipe flow of a viscoelastic fluid, obeying the Oldroyd-B constitutive equation commonly used to model dilute polymer solutions. The instability is shown to exist at Reynolds numbers significantly lower than those at which transition to turbulence is typically observed for Newtonian pipe flow. Our results qualitatively explain experimental observations of transition to turbulence in pipe flow of dilute polymer solutions at flow rates where Newtonian turbulence is absent. The instability discussed here should form the first stage in a hitherto unexplored dynamical pathway to turbulence in polymer solutions. An analogous instability exists for plane Poiseuille flow.

Simple framework for understanding the universality of the maximum drag reduction asymptote in turbulent flow of polymer solutions

Self-consistent direct numerical simulations of turbulent channel flows of dilute polymer solutions exhibiting friction drag reduction (DR) show that an effective Deborah number defined as the ratio of polymer relaxation time to the time scale of fluctuations in the vorticity in the mean flow direction remains O(1) from the onset of DR to the maximum drag reduction (MDR) asymptote. However, the ratio of the convective time scale associated with streamwise vorticity fluctuations to the vortex rotation time decreases with increasing DR, and the maximum drag reduction asymptote is achieved when these two time scales become nearly equal. Based on these observations, a simple framework is proposed that adequately describes the influence of polymer additives on the extent of DR from the onset of DR to MDR as well as the universality of the MDR in wall-bounded turbulent flows with polymer additives.

Inertio-elastic focusing of bioparticles in microchannels at high throughput

Controlled manipulation of particles from very large volumes of fluid at high throughput is critical for many biomedical, environmental and industrial applications. One promising approach is to use microfluidic technologies that rely on fluid inertia or elasticity to drive lateral migration of particles to stable equilibrium positions in a microchannel. Here, we report on a hydrodynamic approach that enables deterministic focusing of beads, mammalian cells and anisotropic hydrogel particles in a microchannel at extremely high flow rates. We show that on addition of micromolar concentrations of hyaluronic acid, the resulting fluid viscoelasticity can be used to control the focal position of particles at Reynolds numbers up to Re≈10,000 with corresponding flow rates and particle velocities up to 50 ml min(-1) and 130 m s(-1). This study explores a previously unattained regime of inertio-elastic fluid flow and demonstrates bioparticle focusing at flow rates that are the highest yet achieved.

Turbulent plane Poiseuille-Couette flow as a model for fluid slip over superhydrophobic surfaces

In this study, plane Poiseuille-Couette flow is simulated as a model for specified streamwise slip on one of the channel walls. The relative velocity between the two walls is set to be 1, 2, and 4 in viscous wall units. This is equivalent to the presence of a superhydrophobic surface at one of the channel walls that causes fluid to slip on the boundary. The results show that the streamwise slip forces turbulence in the near-wall region to tend towards a limiting one-component state. This leads to the suppression of small scale turbulence and laminarization close to the wall and then to drag reduction. The selective weakening of the streamwise vorticity close the wall and the observed decrease of turbulence kinetic energy production can then be considered as a consequence of this effect. Changes in the coherent structures, including a decrease of sweep events and increase of ejection events close to the wall where slip occurs, are also observed.

Polymer-induced drag enhancement in turbulent Taylor-Couette flows: direct numerical simulations and mechanistic insight

We report for the first time the polymer-induced breakdown of large-scale Taylor vortex structures leading to drag enhancement in viscoelastic turbulent Taylor-Couette flows. Specifically, we demonstrate that upon the addition of trace amounts of soluble high molecular weight macromolecules the Newtonian large-scale Taylor vortices are replaced by small-scale vortices in the inner and outer cylinder wall regions. This flow transition and a commensurate drag increase of up to 62% are facilitated by the presence of large polymeric normal stresses in a narrow region immediately close to the outer wall. A simple mechanism for this striking flow transition is proposed with the aim of paving the way for a mechanistic understanding of polymer-induced structure and drag modifications in high-Re turbulent curvilinear flows.

Spatiotemporal evolution of hairpin eddies, Reynolds stress, and polymer torque in polymer drag-reduced turbulent channel flows

To study the influence of dynamic interactions between turbulent vortical structures and polymer stress on turbulent friction drag reduction, a series of simulations of channel flow is performed. We obtain self-consistent evolution of an initial eddy in the presence of polymer stresses by utilizing the finitely extensible nonlinear elastic-Peterlin (FENE-P) model. The initial eddy is extracted by the conditional averages for the second quadrant event from fully turbulent Newtonian flow, and the initial polymer conformation fields are given by the solutions of the FENE-P model equations corresponding to the mean shear flow in the Newtonian case. At a relatively low Weissenberg number We(τ) (=50), defined as the ratio of the polymer relaxation time to the wall time scale, the generation of new vortices is inhibited by polymer-induced countertorques. Thus fewer vortices are generated in the buffer layer. However, the head of the primary hairpin is unaffected by the polymer stress. At larger We(τ) values (≥100), the hairpin head becomes weaker and vortex autogeneration and Reynolds stress growth are almost entirely suppressed. The temporal evolution of the vortex strength and polymer torque magnitude reveals that polymer extension by the vortical motion results in a polymer torque that increases in magnitude with time until a maximum value is reached over a time scale comparable to the polymer relaxation time. The polymer torque retards the vortical motion and Reynolds stress production, which in turn weakens flow-induced chain extension and torque itself. An analysis of the vortex time scales reveals that with increasing We(τ), vortical motions associated with a broader range of time scales are affected by the polymer stress. This is qualitatively consistent with Lumley's time criterion for the onset of drag reduction.

Numerical study of laminar-turbulent transition in particle-laden channel flow

We present direct numerical simulations of subcritical transition to turbulence in a particle-laden channel flow, with particles assumed rigid, spherical, and heavier than the fluid. The equations describing the fluid flow are solved with an Eulerian mesh, whereas those describing the particle dynamics are solved by Lagrangian tracking. Two-way coupling between fluid and particles is modeled with Stokes drag. The numerical code is first validated against previous results from linear stability: the nonmodal growth of streamwise vortices resulting in streamwise streaks is still the most efficient mechanism for linear disturbance amplification at subcritical conditions as for the case of a single phase fluid. To analyze the full nonlinear transition, we examine two scenarios well studied in the literature: (1) transition initiated by streamwise independent counter-rotating streamwise vortices and one three-dimensional mode and (2) oblique transition, initiated by the nonlinear interaction of two symmetric oblique waves. The threshold energy for transition is computed, and it is demonstrated that for both scenarios the transition may be facilitated by the presence of particles at low number density. This is due to the fact that particles may introduce in the system detrimental disturbances of length scales not initially present. At higher concentrations, conversely, we note an increase of the disturbance energy needed for transition. The threshold energy for the oblique scenario shows a more significant increase in the presence of particles, by a factor about four. Interestingly, for the streamwise-vortex scenario the time at which transition occurs increases with the particle volume fraction when considering disturbances of equal initial energy. These results are explained by considering the reduced amplification of oblique modes in the two-phase flow. The results from these two classical scenarios indicate that, although linear stability analysis shows hardly any effect on optimal growth, particles do influence secondary instabilities and streak breakdown. These effects can be responsible of the reduced drag observed in turbulent channel flow laden with heavy particles.

Dynamics of hairpin vortices and polymer-induced turbulent drag reduction

It has been known for over six decades that the dissolution of minute amounts of high molecular weight polymers in wall-bounded turbulent flows results in a dramatic reduction in turbulent skin friction by up to 70%. First principles simulations of turbulent flow of model polymer solutions can predict the drag reduction (DR) phenomenon. However, the essential dynamical interactions between the coherent structures present in turbulent flows and polymer conformation field that lead to DR are poorly understood. We examine this connection via dynamical simulations that track the evolution of hairpin vortices, i.e., counter-rotating pairs of quasistreamwise vortices whose nonlinear autogeneration and growth, decay and breakup are centrally important to turbulence stress production. The results show that the autogeneration of new vortices is suppressed by the polymer stresses, thereby decreasing the turbulent drag.

Mechanism of polymer drag reduction using a low-dimensional model

Using a retarded-motion expansion to describe the polymer stress, we derive a low-dimensional model to understand the effects of polymer elasticity on the self-sustaining process that maintains the coherent wavy streamwise vortical structures underlying wall-bounded turbulence. Our analysis shows that at small Weissenberg numbers, Wi, elasticity enhances the coherent structures. At higher Wi, however, polymer stresses suppress the streamwise vortices (rolls) by calming down the instability of the streaks that regenerates the rolls. We show that this behavior can be attributed to the nonmonotonic dependence of the biaxial extensional viscosity on Wi, and identify it as the key rheological property controlling drag reduction.

Pattern formation in Taylor-Couette flow of dilute polymer solutions: dynamical simulations and mechanism

We report spatiotemporal pattern formation in Taylor-Couette flow (i.e., flow between rotating cylinders) of viscoelastic dilute polymer solutions obtained for the first time from first-principles dynamical simulations. Solution structures with varying spatial and temporal symmetries, such as rotating standing waves, flames, disordered oscillations, and solitary vortex solutions which include diwhirls (stationary and axisymmetric) and oscillatory strips (axisymmetric or nonaxisymmetric), are observed, depending on the ratio of fluid relaxation time to the time period of inner cylinder rotation. The flow-microstructure coupling mechanisms underlying the pattern formation process are also discussed.

Subcritical finite-amplitude solutions for plane Couette flow of viscoelastic fluids

Plane Couette flow of viscoelastic fluids is shown to exhibit a purely elastic subcritical instability at a very small-Reynolds number in spite of being linearly stable. The mechanism of this instability is proposed and the nonlinear stability analysis of plane Couette flow of the Upper-Convected Maxwell fluid is presented. Above a critical Weissenberg number, a small finite-size perturbation is sufficient to create a secondary flow, and the threshold value for the amplitude of the perturbation decreases as the Weissenberg number increases. The results suggest a scenario for weakly turbulent viscoelastic flow which is similar to the one for Newtonian fluids as a function of Reynolds number.

Traveling waves in pipe flow

A family of three-dimensional traveling waves for flow through a pipe of circular cross section is identified. The traveling waves are dominated by pairs of downstream vortices and streaks. They originate in saddle-node bifurcations at Reynolds numbers as low as 1250. All states are immediately unstable. Their dynamical significance is that they provide a skeleton for the formation of a chaotic saddle that can explain the intermittent transition to turbulence and the sensitive dependence on initial conditions in this shear flow.