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.Comment: 5 pages, 3 figures, published version, Phys. Rev. Lett. 89, 208301 (2002

Drag Reduction by Polymers in Turbulent Channel Flows: Energy Redistribution Between Invariant Empirical Modes

We address the phenomenon of drag reduction by dilute polymeric additive to turbulent flows, using Direct Numerical Simulations (DNS) of the FENE-P model of viscoelastic flows. It had been amply demonstrated that these model equations reproduce the phenomenon, but the results of DNS were not analyzed so far with the goal of interpreting the phenomenon. In order to construct a useful framework for the understanding of drag reduction we initiate in this paper an investigation of the most important modes that are sustained in the viscoelastic and Newtonian turbulent flows respectively. The modes are obtained empirically using the Karhunen-Loeve decomposition, allowing us to compare the most energetic modes in the viscoelastic and Newtonian flows. The main finding of the present study is that the spatial profile of the most energetic modes is hardly changed between the two flows. What changes is the energy associated with these modes, and their relative ordering in the decreasing order from the most energetic to the least. Modes that are highly excited in one flow can be strongly suppressed in the other, and vice versa. This dramatic energy redistribution is an important clue to the mechanism of drag reduction as is proposed in this paper. In particular there is an enhancement of the energy containing modes in the viscoelastic flow compared to the Newtonian one; drag reduction is seen in the energy containing modes rather than the dissipative modes as proposed in some previous theories.Comment: 11 pages, 13 figures, included, PRE, submitted, REVTeX

Multi-Zone Shell Model for Turbulent Wall Bounded Flows

We suggested a \emph{Multi-Zone Shell} (MZS) model for wall-bounded flows accounting for the space inhomogeneity in a "piecewise approximation", in which cross-section area of the flow, $S$, is subdivided into "$j$-zones". The area of the first zone, responsible for the core of the flow, $S_1\simeq S/2$, and areas of the next $j$-zones, $S_j$, decrease towards the wall like $S_j\propto 2^{-j}$. In each $j$-zone the statistics of turbulence is assumed to be space homogeneous and is described by the set of "shell velocities" $u_{nj}(t)$ for turbulent fluctuations of the scale $\propto 2^{-n}$. The MZS-model includes a new set of complex variables, $V_j(t)$, $j=1,2,... \infty$, describing the amplitudes of the near wall coherent structures of the scale $s_j\sim 2^{-j}$ and responsible for the mean velocity profile. Suggested MZS-equations of motion for $u_{nj}(t)$ and $V_j(t)$ preserve the actual conservations laws (energy, mechanical and angular momenta), respect the existing symmetries (including Galilean and scale invariance) and account for the type of the non-linearity in the Navier-Stokes equation, dimensional reasoning, etc. The MZS-model qualitatively describes important characteristics of the wall bounded turbulence, e.g., evolution of the mean velocity profile with increasing Reynolds number, \RE, from the laminar profile towards the universal logarithmic profile near the flat-plane boundary layer as \RE\to \infty.Comment: 27 pages, 17 figs, included, PRE, submitte

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 two orders of magnitude higher than in a laminar shear flow. The results indicate that the polymer molecules get strongly stretched by the random motion of the fluid.Comment: 4 pages, 3 figure

Effects of heme on the structure of the denatured state and folding kinetics of cytochrome b562

Heme-linked proteins, such as cytochromes, are popular subjects for protein folding studies. There is the underlying question of whether the heme affects the structure of the denatured state by cross-linking it and forming other interactions, which would perturb the folding pathway. We have studied wild-type and mutant cytochrome b562 from Escherichia coli, a 106 residue four-α-helical bundle. The holo protein apparently refolds with a half-life of 4 μs in its ferrous state. We have analysed the folding of the apo protein using continuous-flow fluorescence as well as stopped-flow fluorescence and CD. The apo protein folded much more slowly with a half-life of 270 μs that was unaffected by the presence of exogenous heme. We examined the nature of the denatured states of both holo and apo proteins by NMR methods over a range of concentrations of guanidine hydrochloride. The starting point for folding of the holo protein in concentrations of denaturant around the denaturation transition was a highly ordered native-like species with heme bound. Fully denatured holo protein at higher concentrations of denaturant consisted of denatured apo protein and free heme. Our results suggest that the very fast folding species of denatured holo protein is in a compact state, whereas the normal folding pathway from fully denatured holo protein consists of the slower folding of the apo protein followed by the binding of heme. These data should be considered in the analysis of folding of heme proteins.This work was supported by the BBSRC under the SBDA initiative, the EU TMR Life Sciences programme (contract ERBFMRX-CT-98-0230), the MRC and the NIH (grant GM056250 to H.R.). P.D.B. thanks the BBSRC for an Advanced Research Fellowshi

Crude incidence in two-phase designs in the presence of competing risks.

BackgroundIn many studies, some information might not be available for the whole cohort, some covariates, or even the outcome, might be ascertained in selected subsamples. These studies are part of a broad category termed two-phase studies. Common examples include the nested case-control and the case-cohort designs. For two-phase studies, appropriate weighted survival estimates have been derived; however, no estimator of cumulative incidence accounting for competing events has been proposed. This is relevant in the presence of multiple types of events, where estimation of event type specific quantities are needed for evaluating outcome.MethodsWe develop a non parametric estimator of the cumulative incidence function of events accounting for possible competing events. It handles a general sampling design by weights derived from the sampling probabilities. The variance is derived from the influence function of the subdistribution hazard.ResultsThe proposed method shows good performance in simulations. It is applied to estimate the crude incidence of relapse in childhood acute lymphoblastic leukemia in groups defined by a genotype not available for everyone in a cohort of nearly 2000 patients, where death due to toxicity acted as a competing event. In a second example the aim was to estimate engagement in care of a cohort of HIV patients in resource limited setting, where for some patients the outcome itself was missing due to lost to follow-up. A sampling based approach was used to identify outcome in a subsample of lost patients and to obtain a valid estimate of connection to care.ConclusionsA valid estimator for cumulative incidence of events accounting for competing risks under a general sampling design from an infinite target population is derived

Polymer transport in random flow

The dynamics of polymers in a random smooth flow is investigated in the framework of the Hookean dumbbell model. The analytical expression of the time-dependent probability density function of polymer elongation is derived explicitly for a Gaussian, rapidly changing flow. When polymers are in the coiled state the pdf reaches a stationary state characterized by power-law tails both for small and large arguments compared to the equilibrium length. The characteristic relaxation time is computed as a function of the Weissenberg number. In the stretched state the pdf is unstationary and exhibits multiscaling. Numerical simulations for the two-dimensional Navier-Stokes flow confirm the relevance of theoretical results obtained for the delta-correlated model.Comment: 28 pages, 6 figure