Yves Dubief

On coherent-vortex identification in turbulence

The identification issue of coherent vortices is investigated on the basis of direct numerical simulation (DNS) and large-eddy simulations (LES) of turbulent flows. It is first shown that the pressure Laplacian is positive within a low-pressure tube of small cross section enclosed by convex isobaric surfaces in a uniform-density flow. Since this quantity is related to the second invariant Q of ∇ u , the Q criterion (region where Q is positive) is a necessary condition for the existence of such tubes. This eduction scheme is compared to other classical methods in incompressible simulations of isotropic turbulence: a mixing layer, a channel flow and a backward-facing step. Q-isosurfaces turn out to display very nice coherent vortices. This criterion is also used in combination with a conditional sampling method to discuss the characteristics of quasi-longitudinal vortices in a manipulated channel flow. The contribution of near-wall vortical structures to velocity and vorticity fluctuations is clearly isolat...

On the coherent drag-reducing and turbulence-enhancing behaviour of polymers in wall flows

Numerical simulations of turbulent polymer solutions using the FENE-P model are used to characterize the action of polymers on turbulence in drag-reduced flows. The energetics of turbulence is investigated by correlating the work done by polymers on the flow with turbulent structures. Polymers are found to store and to release energy to the flow in a well-organized manner. The storage of energy occurs around near-wall vortices as has been anticipated for a long time. Quite unexpectedly, coherent release of energy is observed in the very near-wall region. Large fluctuations of polymer work are shown to re-energize decaying streamwise velocity fluctuations in high-speed streaks just above the viscous sublayer. These distinct behaviours are used to propose an autonomous regeneration cycle of polymer wall turbulence, in the spirit of Jimenez & Pinelli (1999).

Direct numerical simulation of polymer-induced drag reduction in turbulent boundary layer flow

We describe a method for direct numerical simulation of polymer-induced friction drag reduction in turbulent boundary layers. The effect of the polymer additives that induce spatial variations of skin-friction drag is included in the momentum equation through a continuum constitutive model for the viscoelastic stress, which is based on the evolution of a parameter describing the fluid microstructure. We demonstrate that the turbulence structure and polymer microstructure evolve asynchronously as one moves in the streamwise direction. We observe an initial development length, which is followed by a quasisteady region where variations in drag reduction are weak. High drag reduction behavior can be present at short downstream distances from the inflow plane.

Numerical simulation of turbulent drag reduction using rigid fibres

We present a study of the drag reduction induced by rigid fibres in a turbulent channel flow using direct numerical simulation. The extra stresses due to the fibres are calculated with the well-known constitutive equation involving the moments of the orientation vector. Drag reductions of up to 26% are calculated, with the largest drag reductions observed using non-Brownian fibres and semi-dilute concentrations. These findings suggest that elasticity is not necessary to achieve turbulent drag reduction. Flow statistics show trends similar to those observed in simulation of polymeric drag reduction: Reynolds stresses are reduced, velocity fluctuations in the wall-normal and spanwise directions are reduced while streamwise fluctuations are increased, and streamwise vorticity is reduced. We observe strong correlations between the fibre stresses and inter-vortex extensional flow regions. Based on these correlations and instantaneous visualizations of the flow field, we propose a mechanism for turbulent drag reduction by rigid fibre additives.

Simulated polymer stretch in a turbulent flow using Brownian dynamics

We examine the phenomenon of polymer drag reduction in a turbulent flow through Brownian dynamics simulations. The dynamics of a large number of single polymer chains along their trajectories is investigated in a Newtonian turbulent channel flow. In particular, the FENE, FENE-P and multimode FENE models with realistic parameters are used to investigate the mechanisms of polymer stretching. A topological methodology is applied to characterize the ability of the flow to stretch the polymers. It is found using conditional statistics that at moderate Weissenberg number Wi the polymers, that are stretched to a large fraction of their maximum extensibility, have experienced a strong biaxial extensional flow. When Wi is increased other flow types can stretch the polymers but the few highly extended molecules again have, on average, experienced a biaxial extensional flow. Moreover, highly extended polymers are found in the near-wall regions around the quasi-streamwise vortices, essentially in regions of strong biaxial extensional flow.

Elasto-inertial turbulence

Turbulence is ubiquitous in nature, yet even for the case of ordinary Newtonian fluids like water, our understanding of this phenomenon is limited. Many liquids of practical importance are more complicated (e.g., blood, polymer melts, paints), however; they exhibit elastic as well as viscous characteristics, and the relation between stress and strain is nonlinear. We demonstrate here for a model system of such complex fluids that at high shear rates, turbulence is not simply modified as previously believed but is suppressed and replaced by a different type of disordered motion, elasto-inertial turbulence. Elasto-inertial turbulence is found to occur at much lower Reynolds numbers than Newtonian turbulence, and the dynamical properties differ significantly. The friction scaling observed coincides with the so-called “maximum drag reduction” asymptote, which is exhibited by a wide range of viscoelastic fluids.

Direct numerical simulation of polymer-induced drag reduction in turbulent boundary layer flow of inhomogeneous polymer solutions

Skin-friction drag reduction in turbulent boundary layer flow of inhomogeneous polymer solutions is investigated using direct numerical simulations. A continuum constitutive model (FENE-P) accounting for the effects of polymer microstructure and concentration is used to describe the effect of viscoelasticity. The evolution of wall friction along the streamwise direction is a function of the dynamics of the polymer distribution in the boundary layer. It is observed that polymer transport decreases drag reduction downstream compared to the homogeneous case. The fluctuations of polymer concentration are anti-correlated with those of the streamwise velocity. Concentration is largest in the low-speed streaks. The physical process creating this effect is primarily that of dilution of the high-speed streaks, where due to the local turbulence structure the dispersion of polymer is strongest. Thus, the polymer-induced drag reduction phenomenon is sustained primarily in the vicinity of the low-speed streaks where the injected polymer additive is most effective.

On the mechanism of elasto-inertial turbulence

Elasto-inertial turbulence (EIT) is a new state of turbulence found in inertial flows with polymer additives. The dynamics of turbulence generated and controlled by such additives is investigated from the perspective of the coupling between polymer dynamics and flow structures. Direct numerical simulations of channel flow with Reynolds numbers ranging from 1000 to 6000 (based on the bulk and the channel height) are used to study the formation and dynamics of elastic instabilities and their effects on the flow. The flow topology of EIT is found to differ significantly from Newtonian wall-turbulence. Structures identified by positive (rotational flow topology) and negative (extensional/compressional flow topology) second invariant Qa isosurfaces of the velocity gradient are cylindrical and aligned in the spanwise direction. Polymers are significantly stretched in sheet-like regions that extend in the streamwise direction with a small upward tilt. The Qa cylindrical structures emerge from the sheets of high polymer extension, in a mechanism of energy transfer from the fluctuations of the polymer stress work to the turbulent kinetic energy. At subcritical Reynolds numbers, EIT is observed at modest Weissenberg number (Wi, ratio polymer relaxation time to viscous time scale). For supercritical Reynolds numbers, flows approach EIT at large Wi. EIT provides new insights on the nature of the asymptotic state of polymer drag reduction (maximum drag reduction), and explains the phenomenon of early turbulence, or onset of turbulence at lower Reynolds numbers than for Newtonian flows observed in some polymeric flows.

Dilutional control of prothrombin activation at physiologically relevant shear rates.

The generation of proteolyzed prothrombin species by preassembled prothrombinase in phospholipid-coated glass capillaries was studied at physiologic shear rates (100-1000 s(-1)). The concentration of active thrombin species (α-thrombin and meizothrombin) reaches a steady state, which varies inversely with shear rate. When corrected for shear rate, steady-state levels of active thrombin species exhibit no variation and a Michaelis-Menten analysis reveals that chemistry of this reaction is invariant between open and closed systems; collectively, these data imply that variations with shear rate arise from dilutional effects. Significantly, the major products observed include nonreactive species arising from the loss of prothrombins phospholipid binding domain (des F1 species). A numerical model developed to investigate the spatial and temporal distribution of active thrombin species within the capillary reasonably approximates the observed output of total thrombin species at different shears; it also predicts concentrations of active thrombin species in the wall region sufficient to account for observed levels of des FI species. The predominant feedback formation of nonreactive species and high levels of the primarily anticoagulant intermediate meizothrombin (∼40% of total active thrombin species) may provide a mechanism to prevent thrombus propagation downstream of a site of thrombosis or hemorrhage.

Re-examining the logarithmic dependence of the mean velocity distribution in polymer drag reduced wall-bounded flow

A re-examination of the logarithmic dependence of the mean velocity distribution in polymer drag reduced flows shows that drag reducing polymers modify the von Karman coefficient and, in channel flow, eradicate the log-layer at high drag reductions. It is also found that the “ultimate profile,” corresponding to the state of maximum drag reduction is not logarithmic.