Biomedical subjects
Publications and source records attributed to Victor S L'vov.
Drag reduction by bubbles in stationary turbulent flows is sensitive to the compressibility of the bubbles. Without this dynamical effect the bubbles only renormalize the fluid density and viscosity, an effect that by itself can only lead to a small percentage of drag reduction. We show in this paper that the dynamics of bubbles and their effect on the compressibility of the mixture can lead to a much higher drag reduction.
We construct a simple analytic model for wall-bounded turbulence, containing only four adjustable parameters. Two of these parameters are responsible for the viscous dissipation of the components of the Reynolds stress tensor. The other two parameters control the nonlinear relaxation of these objects. The model offers an analytic description of the profiles of the mean velocity and the correlation functions of velocity fluctuations in the entire boundary region, from the viscous sublayer, through the buffer layer, and further into the log-law turbulent region. In particular, the model predicts a very simple distribution of the turbulent kinetic energy in the log-law region between the velocity components: the streamwise component contains a half of the total energy whereas the wall-normal and cross-stream components contain a quarter each. In addition, the model predicts a very simple relation between the von Kármán slope k and the turbulent velocity in the log-law region v+ (in wall units): v+=6k. These predictions are in excellent agreement with direct numerical simulation data and with recent laboratory experiments.
Drag reduction by polymers in wall turbulence is bounded from above by a universal maximal drag reduction (MDR) velocity profile that is a log law, estimated experimentally by Virk as V+(y+) approximately 11.7logy+ - 17. Here V+(y+) and y+ are the mean streamwise velocity and the distance from the wall in "wall" units. In this Letter we propose that this MDR profile is an edge solution of the Navier-Stokes equations (with an effective viscosity profile) beyond which no turbulent solutions exist. This insight rationalizes the universality of the MDR and provides a maximum principle which allows an ab initio calculation of the parameters in this law without any viscoelastic experimental input.
We address the additive equivalence discovered by Virk and co-workers: drag reduction affected by flexible and rigid rodlike polymers added to turbulent wall-bounded flows is limited from above by a very similar maximum drag reduction (MDR) asymptote. Considering the equations of motion of rodlike polymers in wall-bounded turbulent ensembles, we show that although the microscopic mechanism of attaining the MDR is very different, the macroscopic theory is isomorphic, rationalizing the interesting experimental observations.
Drag reduction by microbubbles is a promising engineering method for improving ship performance. A fundamental theory of the phenomenon is lacking, however, making actual design quite haphazard. We offer here a theory of drag reduction by microbubbles in the limit of very small bubbles, when the effect of the bubbles is mainly to normalize the density and the viscosity of the carrier fluid. The theory culminates with a prediction of the degree of drag reduction given the concentration profile of the bubbles. Comparisons with experiments are discussed and the road ahead is sketched.
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.
Drag reduction by polymers in turbulent flows raises an apparent contradiction: the stretching of the polymers must increase the viscosity, so why is the drag reduced? A recent theory proposed that drag reduction, in agreement with experiments, is consistent with the effective viscosity growing linearly with the distance from the wall. With this self-consistent solution the reduction in the Reynolds stress overwhelms the increase in viscous drag. In this Rapid Communication we show, using direct numerical simulations, that a linear viscosity profile indeed reduces the drag in agreement with the theory and in close correspondence with direct simulations of the FENE-P model at the same flow conditions.
We elucidate the mechanism of drag reduction by polymers in turbulent wall-bounded flows: while momentum is produced at a fixed rate by the forcing, polymer stretching results in the suppression of momentum flux to the wall. On the basis of the equations of fluid mechanics we develop the phenomenology of the "maximum drag reduction asymptote" which is the maximum drag reduction attained by polymers. Based on Newtonian information only we demonstrate the existence of drag reduction, and with one experimental parameter we reach agreement with the experimental measurements.
We suggested a multizone shell (MZS) model for wall-bounded flows accounting for the space inhomogeneity in a piecewise approximation, in which the cross-sectional area of the flow, S, is subdivided into j zones. The area of the first zone, responsible for the core of the flow, S1 approximately S/2, and the areas of the next j zones, S(j), decrease toward the wall like S(j) proportional, variant 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 proportional to 2(-n). The MZS model includes a set of complex variables V(j)(t), j=1,2, em leader, infinity, describing the amplitudes of the near-wall coherent structures of the scale s(j) approximately 2(-j) and responsible for the mean velocity profile. The suggested MZS equations of motion for u(nj)(t) and V(j)(t) preserve the actual conservation laws (energy, mechanical, and angular momenta), respect the existing symmetries (including Galilean and scale invariance), and account for the type of nonlinearity 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 toward the universal logarithmic profile near the flat-plane boundary layer as Re--> infinity.
The weak version of universality in turbulence refers to the independence of the scaling exponents of the nth order structure functions from the statistics of the forcing. The strong version includes universality of the coefficients of the structure functions in the isotropic sector, once normalized by the mean energy flux. We demonstrate that shell models of turbulence exhibit strong universality for both forced and decaying turbulence. The exponents and the normalized coefficients are time independent in decaying turbulence, forcing independent in forced turbulence, and equal for decaying and forced turbulence. We conjecture that this is also the case for Navier-Stokes turbulence.
We address the phenomenon of drag reduction by a 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-Loéve 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.
We propose a one-fluid analytical model for a turbulently flowing dilute suspension, based on a modified Navier-Stokes equation with a k-dependent effective density of suspension rho(eff)(k) and an additional damping term proportional, variant gamma(p)(k), representing the fluid-particle friction (described by Stokes law). The statistical description of turbulence within the model is simplified by a modification of the usual closure procedure based on the Richardson-Kolmogorov picture of turbulence with a differential approximation for the energy transfer term. The resulting ordinary differential equation for the energy budget is solved analytically for various important limiting cases and numerically in the general case. In the inertial interval of scales, we describe analytically two competing effects: the energy suppression due to the fluid-particle friction and the energy enhancement during the cascade process due to decrease of the effective density of the small-scale motions. An additional suppression or enhancement of the energy density may occur in the viscous subrange, caused by the variation of the extent of the inertial interval due to the combined effect of the fluid-particle friction and the decrease of the kinematic viscosity of the suspensions. The analytical description of the complicated interplay of these effects supported by numerical calculations is presented. Our findings allow one to rationalize the qualitative picture of the isotropic homogeneous turbulence of dilute suspensions as observed in direct numerical simulations.
In anisotropic turbulence, the correlation functions are decomposed in the irreducible representations of the SO(3) symmetry group (with different "angular momenta" l). For different values of l, the second-order correlation function is characterized by different scaling exponents zeta(2)(l). In this paper, we compute these scaling exponents in a closure approximation. By linearizing the closure equations in small anisotropy we set up a linear operator and find its zero modes in the inertial interval of scales. Thus the scaling exponents in each l sector follow from solvability condition, and are not determined by dimensional analysis. The main result of our calculation is that the scaling exponents zeta(2)(l) form a strictly increasing spectrum at least until l=6, guaranteeing that the effects of anisotropy decay as power laws when the scale of observation diminishes. The results of our calculations are compared to available experiments and simulations.
Motivated by the large effect of turbulent drag reduction by minute concentrations of polymers, we study the effects of a weakly space-dependent viscosity on the stability of hydrodynamic flows. In a recent paper [Phys. Rev. Lett. 87, 174501, (2001)], we exposed the crucial role played by a localized region where the energy of fluctuations is produced by interactions with the mean flow (the "critical layer"). We showed that a layer of a weakly space-dependent viscosity placed near the critical layer can have a very large stabilizing effect on hydrodynamic fluctuations, retarding significantly the onset of turbulence. In this paper we extend these observations in two directions: first we show that the strong stabilization of the primary instability is also obtained when the viscosity profile is realistic (inferred from simulations of turbulent flows with a small concentration of polymers). Second, we analyze the secondary instability (around the time-dependent primary instability) and find similar strong stabilization. Since the secondary instability develops around a time-dependent solution and is three dimensional, this brings us closer to the turbulent case. We reiterate that the large effect is not due to a modified dissipation (as is assumed in some theories of drag reduction), but due to reduced energy intake from the mean flow to the fluctuations. We propose that similar physics act in turbulent drag reduction.
A theory of clustering of inertial particles advected by a turbulent velocity field caused by an instability of their spatial distribution is suggested. The reason for the clustering instability is a combined effect of the particles inertia and a finite correlation time of the velocity field. The crucial parameter for the clustering instability is the size of the particles. The critical size is estimated for a strong clustering (with a finite fraction of particles in clusters) associated with the growth of the mean absolute value of the particles number density and for a weak clustering associated with the growth of the second and higher moments. A new concept of compressibility of the turbulent diffusion tensor caused by a finite correlation time of an incompressible velocity field is introduced. In this model of the velocity field, the field of Lagrangian trajectories is not divergence free. A mechanism of saturation of the clustering instability associated with the particles collisions in the clusters is suggested. Applications of the analyzed effects to the dynamics of droplets in the turbulent atmosphere are discussed. An estimated nonlinear level of the saturation of the droplets number density in clouds exceeds by the orders of magnitude their mean number density. The critical size of cloud droplets required for cluster formation is more than 20 microm.
We consider shell models that display an inverse energy cascade similar to two-dimensional turbulence (together with a direct cascade of an enstrophylike invariant). Previous attempts to construct such models ended negatively, stating that shell models give rise to a "quasiequilibrium" situation with equipartition of the energy among the shells. We show analytically that the quasiequilibrium state predicts its own disappearance upon changing the model parameters in favor of the establishment of an inverse cascade regime with Kolmogorov scaling. The latter regime is found where predicted, offering a useful model to study inverse cascades.
The statistics of two-dimensional turbulence exhibit a riddle: the scaling exponents in the regime of inverse energy cascade agree with the K41 theory of turbulence far from equilibrium, but the probability distribution functions are close to Gaussian-like in equilibrium. The skewness S identical with S3(R)/S(3/2)(2)(R) was measured as S (exp) approximately 0.03. This contradiction is lifted by understanding that two-dimensional turbulence is not far from a situation with equipartition of enstrophy, which exists as true thermodynamic equilibrium with K41 exponents in space dimension of d= 4 / 3. We evaluate the skewness S( d) for 4 / 3 < or =d< or =2, showing that S(d)=0 at d= 4 / 3, and that it remains as small as S (exp) in two dimensions.
A variation principle is suggested to find self-similar solitary solutions (referred to as solitons) of shell model of turbulence. For the Sabra shell model the shape of the solitons is approximated by rational trial functions with relative accuracy of O(10(-3)). It is found how the soliton shape, propagation time t(n) (from a shell n to shells with n --> infinity), and the dynamical exponent z(0) (which governs the time rescaling of the solitons in different shells) depend on parameters of the model. For a finite interval of z the author discovered quasisolitons which approximate with high accuracy corresponding self-similar equations for an interval of times from -infinity to some time in the vicinity of the peak maximum or even after it. The conjecture is that the trajectories in the vicinity of the quasisolitons (with continuous spectra of z) provide an essential contribution to the multiscaling statistics of high-order correlation functions, referred to in the paper as an asymptotic multiscaling. This contribution may be even more important than that of the trajectories in the vicinity of the exact soliton with a fixed value z(0). Moreover there are no solitons in some regions of the parameters where quasisolitons provide a dominant contribution to the asymptotic multiscaling.