Turbulent dynamics of polymer solutions

Nonequilibrium interactions between multi-scale colloids regulate the suspension microstructure and rheology

Understanding nonequilibrium interactions of multi-component colloidal suspensions is critical for many dynamical settings such as self-assembly and material processing. A key question is how the nonequilibrium distributions of individual components influence the effective interparticle interactions and flow behavior. In this work, we develop a first-principle framework to study a bidisperse suspension of colloids and depletants using a Smoluchowski equation and corroborated by Brownian dynamics (BD) simulations. Using nonlinear microrheology as a case study, we demonstrate that effective depletion interactions between driven colloids are sensitive to particle timescales out of equilibrium and cannot be predicted by equilibrium-based pair potentials like Asakura-Oosawa. Furthermore, we show that the interplay between Brownian relaxation timescales of different species plays a critical role in governing the viscosity of multi-component suspensions. Our model highlights the limitations of using equilibrium pair potentials to approximate interparticle interactions in nonequilibrium processes such as hydrodynamic flows and presents a useful framework for studying the transport of driven, interacting suspensions.

The statistical geometry of material loops in turbulence

Material elements - which are lines, surfaces, or volumes behaving as passive, non-diffusive markers - provide an inherently geometric window into the intricate dynamics of chaotic flows. Their stretching and folding dynamics has immediate implications for mixing in the oceans or the atmosphere, as well as the emergence of self-sustained dynamos in astrophysical settings. Here, we uncover robust statistical properties of an ensemble of material loops in a turbulent environment. Our approach combines high-resolution direct numerical simulations of Navier-Stokes turbulence, stochastic models, and dynamical systems techniques to reveal predictable, universal features of these complex objects. We show that the loop curvature statistics become stationary through a dynamical formation process of high-curvature folds, leading to distributions with power-law tails whose exponents are determined by the large-deviations statistics of finite-time Lyapunov exponents of the flow. This prediction applies to advected material lines in a broad range of chaotic flows. To complement this dynamical picture, we confirm our theory in the analytically tractable Kraichnan model with an exact Fokker-Planck approach.

Continuous Pathway between the Elasto-Inertial and Elastic Turbulent States in Viscoelastic Channel Flow

Viscoelastic plane Poiseuille flow is shown to become linearly unstable in the absence of inertia, in the limit of high elasticities, for ultradilute polymer solutions. While inertialess elastic instabilities have been predicted for curvilinear shear flows, this is the first ever report of a purely elastic linear instability in a rectilinear shear flow. The novel instability continues up to a Reynolds number (Re) of O(1000), corresponding to the recently identified elasto-inertial turbulent state believed to underlie the maximum-drag-reduced regime. Thus, for highly elastic ultradilute polymer solutions, a single linearly unstable modal branch may underlie transition to elastic turbulence at zero Re and to elasto-inertial turbulence at moderate Re, implying the existence of continuous pathways connecting the turbulent states to each other and to the laminar base state.

Entropic characterization of the coil-stretch transition of polymers in random flows

Polymer molecules in a flow undergo a coil-stretch phase transition on an increase of the velocity gradients. Model-independent identification and characterization of the transition in a random flow has been lacking so far. Here we suggest to use the entropy of the extension statistics as a proper measure due to strong fluctuations around the transition. We measure experimentally the entropy as a function of the local Weisenberg number and show that it has a maximum, which identifies and quantifies the transition. We compare the new approach with the traditional one based on the theory using either linear Oldroyd-B or nonlinear finite extensible nonlinear elastic polymer models.

Effect of internal friction on the coil-stretch transition in turbulent flows

A polymer in a turbulent flow undergoes the coil-stretch transition when the Weissenberg number, i.e. the product of the Lyapunov exponent of the flow and the relaxation time of the polymer, surpasses a critical value. The effect of internal friction on the transition is studied by means of Brownian dynamics simulations of the elastic dumbbell model in a homogeneous and isotropic, incompressible, turbulent flow and analytical calculations for a stochastic velocity gradient. The results are explained by adapting the large deviations theory of Balkovsky et al. [Phys. Rev. Lett., 2000, 84, 4765] to an elastic dumbbell with internal viscosity. In turbulent flows, a distinctive feature of the probability distribution of polymer extensions is its power-law behaviour for extensions greater than the equilibrium length and smaller than the contour length. It is shown that although internal friction does not modify the critical Weissenberg number for the coil-stretch transition, it makes the slope of the probability distribution of the extension steeper, thus rendering the transition sharper. Internal friction therefore provides a possible explanation for the steepness of the distribution of polymer extensions observed in experiments at large Weissenberg numbers.

Elastic turbulence in entangled semi-dilute DNA solutions measured with optical coherence tomography velocimetry

The flow instabilities of solutions of high molecular weight DNA in the entangled semi-dilute concentration regime were investigated using optical coherence tomography velocimetry, a technique that provides high spatial (probe volumes of 3.4 pL) and temporal resolution (sub μs) information on the flow behaviour of complex fluids in a rheometer. The velocity profiles of the opaque DNA solutions (high and low salt) were measured as a function of the distance across the gap of a parallel plate rheometer, and their evolution over time was measured. At lower DNA concentrations and low shear rates, the velocity fluctuations were well described by Gaussian functions and the velocity gradient was uniform across the rheometer gap, which is expected for Newtonian flows. As the DNA concentration and shear rate were increased there was a stable wall slip regime followed by an evolving wall slip regime, which is finally followed by the onset of elastic turbulence. Strain localization (shear banding) is observed on the boundaries of the flows at intermediate shear rates, but decreases in the high shear elastic turbulence regime, where bulk strain localization occurs. A dynamic phase diagram for non-linear flow was created to describe the different behaviours.

Dynamics of comb-of-comb-network polymers in random layered flows

We analyze the dynamics of comb-of-comb-network polymers in the presence of external random flows. The dynamics of such structures is evaluated through relevant physical quantities, viz., average square displacement (ASD) and the velocity autocorrelation function (VACF). We focus on comparing the dynamics of the comb-of-comb network with the linear polymer. The present work displays an anomalous diffusive behavior of this flexible network in the random layered flows. The effect of the polymer topology on the dynamics is analyzed by varying the number of generations and branch lengths in these networks. In addition, we investigate the influence of external flow on the dynamics by varying flow parameters, like the flow exponent α and flow strength W_{α}. Our analysis highlights two anomalous power-law regimes, viz., subdiffusive (intermediate-time polymer stretching and flow-induced diffusion) and superdiffusive (long-time flow-induced diffusion). The anomalous long-time dynamics is governed by the temporal exponent ν of ASD, viz., ν=2-α/2. Compared to a linear polymer, the comb-of-comb network shows a shorter crossover time (from the subdiffusive to superdiffusive regime) but a reduced magnitude of ASD. Our theory displays an anomalous VACF in the random layered flows that scales as t^{-α/2}. We show that the network with greater total mass moves faster.

Polymer stretching in the inertial range of turbulence

We study the deformation of flexible polymers whose contour length lies in the inertial range of a homogeneous and isotropic turbulent flow. By using the elastic dumbbell model and a stochastic velocity field with nonsmooth spatial correlations, we obtain the probability density function of the extension as a function of the Weissenberg number and of the scaling exponent of the velocity structure functions. In a spatially rough flow, as in the inertial range of turbulence, the statistics of polymer stretching differs from that observed in laminar flows or in smooth chaotic flows. In particular, the probability distribution of polymer extensions decays as a stretched exponential, and the most probable extension grows as a power law of the Weissenberg number. Furthermore, the ability of the flow to stretch polymers weakens as the flow becomes rougher in space.

Large-deviation statistics of vorticity stretching in isotropic turbulence

A key feature of three-dimensional fluid turbulence is the stretching and realignment of vorticity by the action of the strain rate. It is shown in this paper, using the cumulant-generating function, that the cumulative vorticity stretching along a Lagrangian path in isotropic turbulence obeys a large deviation principle. As a result, the relevant statistics can be described by the vorticity stretching Cramér function. This function is computed from a direct numerical simulation data set at a Taylor-scale Reynolds number of Re(λ)=433 and compared to those of the finite-time Lyapunov exponents (FTLE) for material deformation. As expected, the mean cumulative vorticity stretching is slightly less than that of the most-stretched material line (largest FTLE), due to the vorticity's preferential alignment with the second-largest eigenvalue of strain rate and the material line's preferential alignment with the largest eigenvalue. However, the vorticity stretching tends to be significantly larger than the second-largest FTLE, and the Cramér functions reveal that the statistics of vorticity stretching fluctuations are more similar to those of the largest FTLE. In an attempt to relate the vorticity stretching statistics to the vorticity magnitude probability density function in statistically stationary conditions, a model Kramers-Moyal equation is constructed using the statistics encoded in the Cramér function. The model predicts a stretched-exponential tail for the vorticity magnitude probability density function, with good agreement for the exponent but significant difference (35%) in the prefactor.

Impact of the Peterlin approximation on polymer dynamics in turbulent flows

We study the impact of the Peterlin approximation on the statistics of the end-to-end separation of polymers in a turbulent flow. The finitely extensible nonlinear elastic (FENE) model and the FENE model with the Peterlin approximation (FENE-P) are numerically integrated along a large number of Lagrangian trajectories resulting from a direct numerical simulation of three-dimensional homogeneous isotropic turbulence. Although the FENE-P model yields results in qualitative agreement with those of the FENE model, quantitative differences emerge. The steady-state probability of large extensions is overestimated by the FENE-P model. The alignment of polymers with the eigenvectors of the rate-of-strain tensor and with the direction of vorticity is weaker when the Peterlin approximation is used. At large Weissenberg numbers, the correlation times of both the extension and of the orientation of polymers are underestimated by the FENE-P model.

Dynamics of generalized Gaussian polymeric structures in random layered flows

We develop a formalism for the dynamics of a flexible branched polymer with arbitrary topology in the presence of random flows. This is achieved by employing the generalized Gaussian structure (GGS) approach and the Matheron-de Marsily model for the random layered flow. The expression for the average square displacement (ASD) of the center of mass of the GGS is obtained in such flow. The averaging is done over both the thermal noise and the external random flow. Although the formalism is valid for branched polymers with various complex topologies, we mainly focus here on the dynamics of the flexible star and dendrimer. We analyze the effect of the topology (the number and length of branches for stars and the number of generations for dendrimers) on the dynamics under the influence of external flow, which is characterized by their root-mean-square velocity, persistence flow length, and flow exponent α. Our analysis shows two anomalous power-law regimes, viz., subdiffusive (intermediate-time polymer stretching and flow-induced diffusion) and superdiffusive (long-time flow-induced diffusion). The influence of the topology of the GGS is unraveled in the intermediate-time regime, while the long-time regime is only weakly dependent on the topology of the polymer. With the decrease in the value of α, the magnitude of the ASD decreases, while the temporal exponent of the ASD increases in both the time regimes. Also there is an increase in both the magnitude of the ASD and the crossover time (from the subdiffusive to the superdiffusive regime) with an increase in the total mass of the polymeric structure.

Mixing and clustering in compressible chaotic stirred flows

The effect of compressibility on the mixing of Lagrangian tracers is analyzed in chaotic stirred flows. Mixing is studied in terms of the finite-time Lyapunov exponents. Mixing and clustering of passive tracers surrounded by Lagrangian coherent structures is observed to increase with compressibility intensity. The role of the stirring rate and compressibility on mixing and clustering has been analyzed.

Statistics of polymer extensions in turbulent channel flow

We present direct numerical simulations of turbulent channel flow with passive Lagrangian polymers. To understand the polymer behavior we investigate the behavior of infinitesimal line elements and calculate the probability distribution function (PDF) of finite-time Lyapunov exponents and from them the corresponding Cramer's function for the channel flow. We study the statistics of polymer elongation for both the Oldroyd-B model (for Weissenberg number Wi<1) and the FENE model. We use the location of the minima of the Cramer's function to define the Weissenberg number precisely such that we observe coil-stretch transition at Wi ≈1. We find agreement with earlier analytical predictions for PDF of polymer extensions made by Balkovsky, Fouxon, and Lebedev [Phys. Rev. Lett. 84, 4765 (2000)] for linear polymers (Oldroyd-B model) with Wi <1 and by Chertkov [Phys. Rev. Lett. 84, 4761 (2000)] for nonlinear FENE-P model of polymers. For Wi >1 (FENE model) the polymer are significantly more stretched near the wall than at the center of the channel where the flow is closer to homogenous isotropic turbulence. Furthermore near the wall the polymers show a strong tendency to orient along the streamwise direction of the flow, but near the center line the statistics of orientation of the polymers is consistent with analogous results obtained recently in homogeneous and isotropic flows.

Control of particle clustering in turbulence by polymer additives

We study the clustering properties of inertial particles in a turbulent viscoelastic fluid. The investigation is carried out by means of direct numerical simulations of turbulence in the Oldroyd-B model. The effects of polymers on the small-scale properties of homogeneous turbulence are considered in relation with their consequences on clustering of particles, both lighter and heavier than the carrying fluid. We show that, depending on particle and flow parameters, polymers can either increase or decrease clustering.

Elastic turbulence in a curvilinear channel flow

We report detailed quantitative studies of elastic turbulence in a curvilinear channel flow in a dilute polymer solution of high molecular weight polyacrylamide in a high viscosity water-sugar solvent. Detailed studies of the average and rms velocity and velocity gradients profiles reveal the emergence of a boundary layer associated with the nonuniform distribution of the elastic stresses across the channel. The characteristic boundary width is independent of the Weissenberg number Wi and proportional to the channel width, which is consistent with the findings our early investigations of the boundary layer in elastic turbulence in different flow geometries. The nonuniform distribution of the elastic stresses across the channel and appearance of the characteristic spatial scales of the order of the boundary layer width of both velocity and velocity gradient in the correlation functions of the velocity and velocity gradient fields in a bulk flow may suggest that excessive elastic stresses, concentrated in the boundary layer, are ejected into the bulk flow similar to jets observed in passive scalar mixing in elastic turbulence observed recently. Finally, the experimental results show that one of the main predictions of the theory of elastic turbulence, namely, the saturation of the normalized rms velocity gradient in the bulk flow of elastic turbulence contradicts the experimental observations both qualitatively and quantitatively in spite of the fact that the theory explains well the observed sharp power-law decay of the velocity power spectrum. The experimental findings call for further development of theory of elastic turbulence in a bounded container, similar to what was done for a passive scalar problem.

Coil-stretch transition in an ensemble of polymers in isotropic turbulence

We study the statistical properties of ensembles of polymers in isotropic turbulence numerically in the one-way coupling regime. A linear polymer chain passively convected by turbulence is modeled by a line of beads, each of which is connected by a finitely extensible nonlinear elastic force and subject to Brownian motion. We find that when the Weissenberg number Wi(η)<1, the polymer chain has a coiled configuration, while for Wi(η)>10, it remains stretched for a much longer time than the typical time scale of the fluctuating turbulent velocity gradient. Various statistical quantities characterizing the ensemble of polymers, such as the mean, variance, autocorrelation time, and probability density function of the end-to-end vector distance, indicate that the coil-stretch transition occurs at Wi(η)=3-4. We also find that this trend is insensitive to the number of beads N(b) ( N(b)=20 or N(b)=2), provided that the parameters in the model with a small number of beads are properly generated from the one with a large number of beads (i.e., using the formula of Jin and Collins). Finally, the Wi(η) effects on the alignment of the end-to-end vector versus the principal axis of the rate of strain tensor and on the polymer elongation are examined from the viewpoint of local flow topology.

Low Reynolds number turbulence in nonlinear Maxwell-model fluids

A generalized nonlinear Maxwell model which had previously been analyzed for plane Couette geometry is here applied to a lid-driven cavity flow. The full three-dimensional hydrodynamical problem is treated numerically. Depending on the relevant model parameters, both smooth laminar and low Reynolds number turbulent flows are found, strikingly similar to the experimentally observed elastic turbulence phenomena in polymer solutions. Representative results of the calculated flow patterns, as well as measures for the turbulent nature of the flow are presented graphically.

Statistics of entropy production in linearized stochastic systems

We consider a wide class of linear stochastic problems driven off the equilibrium by a multiplicative asymmetric force. The force breaks detailed balance, maintained otherwise, thus producing entropy. The large deviation function of the entropy production in the system is calculated explicitly. The general result is illustrated using an example of a polymer immersed in a gradient flow and subject to thermal fluctuations.

Universal scaling for polymer chain scission in turbulence

We report that previous polymer chain scission experiments in strong flows, long analyzed according to accepted laminar flow scission theories, were in fact affected by turbulence. We reconcile existing anomalies between theory and experiment with the hypothesis that the local stress at the Kolmogorov scale generates the molecular tension leading to polymer covalent bond breakage. The hypothesis yields a universal scaling for polymer scission in turbulent flows. This surprising reassessment of over 40 years of experimental data simplifies the theoretical picture of polymer dynamics leading to scission and allows control of scission in commercial polymers and genomic DNA.

Dynamical slowdown of polymers in laminar and random flows

The influence of an external flow on the relaxation dynamics of a single polymer is investigated theoretically and numerically. We show that a pronounced dynamical slowdown occurs in the vicinity of the coil-stretch transition, especially when the dependence on polymer conformation of the drag is accounted for. For the elongational flow, relaxation times are exceedingly larger than the Zimm relaxation time, resulting in the observation of conformation hysteresis. For random smooth flows, hysteresis is not present. Yet, relaxation dynamics is significantly slowed down because of the large variety of accessible polymer configurations. The implications of these results for the modeling of dilute polymer solutions in turbulent flows are addressed.

Role of elastic stress in statistical and scaling properties of elastic turbulence

The role of elastic stress in statistical and scaling properties of elastic turbulence in a polymer solution flow between two disks is discussed. The analogy with a small-scale magnetodynamics and a passive scalar turbulent advection in the Batchelor regime is used to explain the experimentally observed statistical properties, the flow structure, and the scaling of elastic turbulence. The emergence of a new length scale, namely, the boundary layer thickness, is observed and studied.

Statistics of tumbling of a single polymer molecule in shear flow

We present experimental results on statistics of polymer orientation angles relative to the shear plane and tumbling times in shear flow with thermal noise. The strong deviation of the probability distribution functions (PDFs) of the orientation angles from Gaussian PDFs was observed in good accord with theory. A universal exponential PDF tail for the tumbling times and its predicted scaling with Wi (that is, the dimensionless shear rate normalized by the polymer relaxation time) are also tested experimentally against numerics. The scaling relations of PDF widths for both angles as a function of Wi are verified and compared with numerics.

Polymer stress tensor in turbulent shear flows

The interaction of polymers with turbulent shear flows is examined. We focus on the structure of the elastic stress tensor, which is proportional to the polymer conformation tensor. We examine this object in turbulent flows of increasing complexity. First is isotropic turbulence, then anisotropic (but homogenous) shear turbulence, and finally wall bounded turbulence. The main result of this paper is that for all these flows the polymer stress tensor attains a universal structure in the limit of large Deborah number De >> 1. We present analytic results for the suppression of the coil-stretch transition at large Deborah numbers. Above the transition the turbulent velocity fluctuations are strongly correlated with the polymer's elongation: there appear high-quality "hydroelastic" waves in which turbulent kinetic energy turns into polymer potential energy and vice versa. These waves determine the trace of the elastic stress tensor but practically do not modify its universal structure. We demonstrate that the influence of the polymers on the balance of energy and momentum can be accurately described by an effective polymer viscosity that is proportional to the cross-stream component of the elastic stress tensor. This component is smaller than the streamwise component by a factor proportional to De2. Finally we tie our results to wall bounded turbulence and clarify some puzzling facts observed in the problem of drag reduction by polymers.

Polymers suppress the inverse transfers of energy and the enstrophy flux fluctuations in two-dimensional turbulence

The addition of minute amounts of a flexible polymer to two-dimensional turbulence produced in fast-flowing soap films affects large scales and small scales differently. For large scales, the inverse transfers of energy are suppressed. For small scales, where mean quantities are barely affected, the enstrophy flux fluctuations are significantly reduced, making the flow less chaotic.

Passive scalar evolution in peripheral regions

We consider the evolution of a passive scalar (concentration of pollutants or temperature) in a chaotic (turbulent) flow. A universal asymptotic behavior of the passive scalar decay (homogenization) related to peripheral regions (near walls) is established. The passive scalar moments and their pair correlation function in the peripheral region are analyzed. A special case investigated in our paper is the passive scalar decay along a pipe.

Two-dimensional turbulence of dilute polymer solutions

We investigate theoretically and numerically the effect of polymer additives on two-dimensional turbulence by means of a viscoelastic model. We provide compelling evidence that, at vanishingly small concentrations, such that the polymers are passively transported, the probability distribution of polymer elongation has a power law tail: Its slope is related to the statistics of finite-time Lyapunov exponents of the flow, in quantitative agreement with theoretical predictions. We show that at finite concentrations and sufficiently large elasticity the polymers react on the flow with manifold consequences: Velocity fluctuations are drastically depleted, as observed in soap film experiments; the velocity statistics becomes strongly intermittent; the distribution of finite-time Lyapunov exponents shifts to lower values, signaling the reduction of Lagrangian chaos.

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

Nontrivial steady flows have recently been found that capture the main structures of the turbulent buffer layer. We study the effects of polymer addition on these "exact coherent states" (ECS) in plane Couette flow. Despite the simplicity of the ECS flows, these effects closely mirror those observed experimentally: Structures shift to larger length scales, wall-normal fluctuations are suppressed while streamwise ones are enhanced, and drag is reduced. The mechanism underlying these effects is elucidated. These results suggest that the ECS are closely related to buffer layer turbulence.

Polymers in 2D turbulence: suppression of large scale fluctuations

Small quantities of a long chain molecule or polymer affect two-dimensional turbulence in unexpected ways. Their presence inhibits the transfers of energy to large scales causing their suppression in the energy density spectrum. This also leads to the change of the spectral properties of a passive scalar which turns out to be highly sensitive to the presence of energy transfers.

Turbulent drag reduction and degradation of DNA

Turbulent drag reduction induced by lambda-DNA is studied. The double-stranded DNA is found to be a good drag reducer when compared with the other normal linear polymers. However, this drag reducing power disappears when the DNA denatures to form two single-strand molecules. Mechanical degradation of DNA is also different from that of the normal linear-chain polymers: DNA is always cut in half by the turbulence. Our results suggest that the mechanism for turbulent degradation of DNA is different from that of the normal flexible long-chain polymers.

Turbulence of polymer solutions

We investigate high-Reynolds-number turbulence in dilute polymer solutions. We show the existence of a critical value of the Reynolds number, which separates two different regimes. In the first regime, below the transition, the influence of the polymer molecules on the flow is negligible, so they can be regarded as passively embedded in the flow. This case admits a detailed investigation of the statistics of the polymer elongations. The second state is realized when the Reynolds number is larger than the critical value. This regime is characterized by the strong back reaction of polymers on the flow. We establish some properties of the statistics of the stress and velocity in this regime and discuss its relation to the drag reduction phenomenon.

Stretching of polymers in a random three-dimensional flow

Behavior of a dilute polymer solution in a random three-dimensional flow with an average shear is studied experimentally. Polymer contribution to the shear stress is found to be more than 2 orders of magnitude higher than in a laminar shear flow. The results indicate that the polymer molecules are strongly stretched by the random motion of the fluid.

Polymer stretching by turbulence

The stretching of a polymer chain by a large-scale chaotic flow is considered. The steady state which emerges as a balance of the turbulent stretching and anharmonic resistance of the chain is quantitatively described, i.e., the dependency on the flow parameters (Lyapunov exponent statistics) and the chain characteristics (the number of beads and the interbead elastic potential) is made explicit.