Polymer stretching by turbulence

Large is different: Nonmonotonic behavior of elastic range scaling in polymeric turbulence at large Reynolds and Deborah numbers

We use direct numerical simulations to study homogeneous and isotropic turbulent flows of dilute polymer solutions at high Reynolds and Deborah numbers. We find that for small wave numbers k, the kinetic energy spectrum shows Kolmogorov-like behavior that crosses over at a larger k to a novel, elastic scaling regime, E(k) ∼ k^-ξ, with ξ ≈ 2.3. We study the contribution of the polymers to the flux of kinetic energy through scales and find that it can be decomposed into two parts: one increase in effective viscous dissipation and a purely elastic contribution that dominates over the nonlinear flux in the range of k over which the elastic scaling is observed. The multiscale balance between the two fluxes determines the crossover wave number that depends nonmonotically on the Deborah number. Consistently, structure functions also show two scaling ranges, with intermittency present in both of them in equal measure.

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.

Tumbling dynamics of inertial inextensible chains in extensional flow

This paper investigates the effect of inertia on the dynamics of elongated chains to go beyond the overdamped case that is often used to study such systems. For that purpose, numerical simulations are performed considering the motion of freely jointed bead-rod chains in an extensional flow in the presence of thermal noise. The coil-stretch transition and the tumbling instability are characterized as a function of three parameters: the Péclet number, the Stokes number, and the chain length. Numerical results show that the coil-stretch transition remains when inertia is present and that it depends nonlinearly on the Stokes and Péclet numbers. Theoretical and numerical analyses also highlight the role of intermediate stable configurations in the dynamics of elongated chains: chains can indeed remain trapped for a certain time in these configurations, especially while undergoing a tumbling event.

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.

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.

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.

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.

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.

Nonlinear trans-resonant waves, vortices and patterns: From microresonators to the early Universe

Perturbed wave equations are considered. Approximate general solutions of these equations are constructed, which describe wave phenomena in different physical and chemical systems. Analogies between surface waves, nonlinear and atom optics, field theories and acoustics of the early Universe can be seen in the similarities between the general solutions that govern each system. With the help of the general solutions and boundary conditions and/or resonant conditions we have derived the basic highly nonlinear ordinary differential equation or the basic algebraic equation for traveling waves. Then, approximate analytic resonant solutions are constructed, which describe the trans-resonant transformation of harmonic waves into traveling shock-, jet-, or mushroom-like waves. The mushroom-like waves can evolve into cloud-like and vortex-like structures. The motion and oscillations of these waves and structures can be very complex. Under parametric excitation these waves can vary their velocity, stop, and change the direction of their motion. Different dynamic patterns are yielded by these resonant traveling waves in the x-t and x-y planes. They simulate many patterns observed in liquid layers, optical systems, superconductors, Bose-Einstein condensates, micro- and electron resonators. The harmonic excitation may be compressed and transformed inside the resonant band into traveling or standing particle-like waves. The area of application of these solutions and results may possibly vary from the generation of nuclear particles, acoustical turbulence, and catastrophic seismic waves to the formation of galaxies and the Universe. In particular, the formation of galaxies and galaxy clusters may be connected with nonlinear and resonant phenomena in the early Universe. (c) 2001 American Institute of Physics.

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.

Turbulent dynamics of polymer solutions

We study properties of dilute polymer solutions. The probability density function (PDF) of polymer end-to-end extensions R in turbulent flows is examined. We show that if the value of the Lyapunov exponent lambda is smaller than the inverse molecular relaxation time 1/tau then the PDF has a strong peak at the equilibrium size R0 and a power tail at R>>R0. This confirms and extends the results of J. L. Lumley [Symp. Math. 9, 315 (1972)]. There is no essential influence of polymers on the flow in this regime. At lambdatau>1 the majority of molecules is stretched to the linear size R(op)>>R0, which can be much smaller than the maximal length of the molecules due to their back reaction.