Large-scale motions, important in turbulent shear flows, are frequently
attributed to the interaction of structures at smaller scale. Here we show
that, in a turbulent channel at Re_{\tau} \approx 550, large-scale motions can
self-sustain even when smaller-scale structures populating the near-wall and
logarithmic regions are artificially quenched. This large-scale self-sustained
mechanism is not active in periodic boxes of width smaller than Lz ~ 1.5h or
length shorter than Lx ~ 3h which correspond well to the most energetic large
scales observed in the turbulent channel

A. Townsend

Carlo Cossu

F. Waleffe

G. Pujals

J. Jiménez

T. R. Bewley

Yongyun Hwang

English

arXiv

International audienceLarge-scale motions, important in turbulent shear flows, are frequently attributed to the interaction of structures at smaller scales. Here we show that, in a turbulent channel at Reτ≈550, large-scale motions can self-sustain even when smaller-scale structures populating the near-wall and logarithmic regions are artificially quenched. This large-scale self-sustained mechanism is not active in periodic boxes of width smaller than Lz≈1.5h or length shorter than Lx≈3h which correspond well to the most energetic large scales observed in the turbulent channel

Hwang, Yongyun

Cossu, Carlo

HAL: Hyper Article en Ligne

Self-sustained process at large scales in turbulent channel flow

Large-scale motions, important in turbulent shear flows, are frequently attributed to the interaction of structures at smaller scales. Here we show that, in a turbulent channel at Reτ≈550, large-scale motions can self-sustain even when smaller-scale structures populating the near-wall and logarithmic regions are artificially quenched. This large-scale self-sustained mechanism is not active in periodic boxes of width smaller than Lz≈1.5h or length shorter than Lx≈3h which correspond well to the most energetic large scales observed in the turbulent channel

Open Archive Toulouse Archive Ouverte

   Open Archive TOULOUSE Archive Ouverte (OATAO)  OATAO is an open access repository that collects the work of Toulouse researchers and makes it freely available over the web where possible.  This is an author-deposited version published in : http://oatao.univ-toulouse.fr/ Eprints ID : 10493 To link to this article : DOI:10.1103/PhysRevLett.105.044505 URL : http://dx.doi.org/10.1103/PhysRevLett.105.044505 To cite this version : Hwang, Yongyun and Cossu, Carlo Self-sustained process at large scales in turbulent channel flow. (2010) Physical Review Letters, vol. 105 (n° 4). pp. 1-4. ISSN 0031-9007 Any correspondance concerning this service should be sent to the repository administrator: staff-oatao@listes-diff.inp-toulouse.fr Self-Sustained Process at Large Scales in Turbulent Channel FlowYongyun Hwang*LadHyX, CNRS-E´cole Polytechnique, 91128 Palaiseau, FranceCarlo Cossu†IMFT-CNRS, Alle´e du Professeur Camille Soula, 31400 Toulouse, Franceand De´partement de Me´canique, E´cole Polytechnique, 91128 Palaiseau, France(Received 16 February 2010; published 23 July 2010)Large-scale motions, important in turbulent shear flows, are frequently attributed to the interaction ofstructures at smaller scales. Here we show that, in a turbulent channel at Re  550, large-scale motionscan self-sustain even when smaller-scale structures populating the near-wall and logarithmic regions areartificially quenched. This large-scale self-sustained mechanism is not active in periodic boxes of widthsmaller than Lz  1:5h or length shorter than Lx  3h which correspond well to the most energetic largescales observed in the turbulent channel.DOI: 10.1103/PhysRevLett.105.044505Understanding the mechanisms by which turbulenceself-sustains in wall-bounded turbulent shear flows is stillan open challenge that is relevant for applications rangingfrom atmospheric dynamics to the design of low dragvehicles. One of the most robust features of these flowsis the presence of streamwise streaks, i.e., narrow stream-wise regions where the streamwise velocity is larger orsmaller than the mean velocity at the same distance fromthe wall. The existence of streaks in the buffer layer with acharacteristic average spanwise scale þz  100, in wallunits, has been revealed by early experimental observa-tions [1,2]. Numerical experiments furthermore reveal thatturbulence is not sustained in periodic domains smallerthan this spacing [3]. There is now a large consensus thatthese near-wall structures are generated via a self-sustainedprocess [4,5] based on the amplification of streaks from thevortices via the lift-up effect [6], followed by the break-down of the streaks [7,8] and the regeneration of thevortices. This buffer-layer process survives when turbulentstructures are artificially removed from the outer layers [9].Coherent structures, however, exist at larger scales. Ithas long been known that streamwise velocity correlationsare important up to lengths of the order of the outer lengthscale h of the flow (e.g. the channel half-width, the piperadius, or the boundary layer thickness), and have beenrelated to the presence of ‘‘bulges’’ or ‘‘large-scale mo-tions’’ [10]. Recently, it has also been realized that coher-ent streaks exist at even larger scales, with typical spanwiseand streamwise scales z  OðhÞ and x  Oð10hÞ, re-spectively [11–15]. These very-large-scale streaks, alsoreferred to as very-large-scale motions or global modesor superstructures, are important because they carry a verysignificant fraction of the turbulent kinetic energy and ofthe turbulent Reynolds stress [12,13,15,16], contradictingthe early view that the motions at very large scales areessentially inactive [17].The fundamental question in which we are interestedhere is as follows: What is the origin of the streaky motionsat very large scales? It has been conjectured that thesemotions could result from the concatenation of large-scalemotions [13,16]. These large-scale motions are, in turn,made of the aggregation of a huge number of hairpinvortices of smaller scales [12,18]. Even from a differentperspective, further numerical experiments [19] supportthe idea that very-large-scale motions can be directlyforced and can interact with near-wall structures, at muchsmaller scales, via a cosupporting cycle. In this mecha-nism, while the near-wall cycle would continue to existeven in the absence of large-scale structures, the reversewould not happen.The current wisdom is therefore that large- and very-large-scale motions would not exist in the absence ofsmaller-scale structures. An alternative way of thinkingis, however, emerging from recent findings that extendrecent hydrodynamic stability methodologies to fully de-veloped turbulent flows. It has recently been theoreticallypredicted that very-large-scale streamwise streaks can beamplified by a coherent lift-up effect which is able toextract energy from the mean flow at very large scales,without the mediation of coherent structures at smallerscales [20–24]. The existence of the coherent lift-up hasbeen confirmed in experiments where very-large-scale co-herent streaks have been artificially induced by forcingcoherent streamwise vortices [25–27]. It is therefore tempt-ing to conjecture that the very-large-scale motions couldresult from a self-sustained mechanism bearing some simi-larity to the processes observed in the buffer layer and intransitional flows (see, e.g., the discussions in [13,22,28]).The goal of the present study is to verify whether or not thisconjecture is true.We proceed by designing a ‘‘conceptual’’ numericalexperiment where we gradually remove the smaller scalesfrom the flow while still taking into account the associateddissipation. The ideal tool to perform such an experiment isthe large eddy simulation (LES), which also has the ad-vantage of keeping the computations in very-large domainsmanageable at high Reynolds numbers. In LES the very-small-scale motions are not resolved by the grid, and theresidual stress tensor associated with these small-scalefiltered motions is modeled with an eddy viscosity ij kkij=3 ¼ 2T Sij, where Sij is the symmetric part ofthe filtered velocity gradient. From the adopted perspectiveit is important to avoid any backward spectral energytransfer from the smaller to the larger scales. Therefore,we consider the ‘‘static’’ Smagorinsky model [29], wherethe eddy viscosity is simply modeled as T ¼ DðCS Þ2 S,with S ¼ ð2 Sij SijÞ1=2 being a measure of the local shearand CS  the local mixing length, given by the product ofthe filtering scale  to the Smagorinsky constant CS. Thevan Driest damping function D ¼ 1 eðyþ=AþÞ3 is alsoused to enforce a physical behavior near the wall (forfurther details on properties of the Smagorinsky model,the reader may refer, e.g., to [30,31]). DynamicSmagorinsky models, which are known to perform betterfor engineering applications, are not used here because wewant to avoid backscatter. The simulations are performedusing the public domain code DIABLO [32] which imple-ments the fractional-step method based on a semi-implicittime integration scheme and a mixed finite-difference andFourier discretization in space. The computational domainis chosen large enough (Lz ¼ 6h) to contain three to fourlarge-scale streaks in the spanwise direction and longenough (Lx ¼ 56h) to contain the longest observed stream-wise wavelengths of the very-large-scale motions. The gridsize of 384 65 96 points in, respectively, the stream-wise (x), wall-normal (y), and spanwise (z) directions is setup so as to resolve the near-wall cycle properly in thereference simulations [33]. A grid stretching is applied inthe wall-normal direction with the first grid point at yþ ¼1:8 and a maximum spacing þy ¼ 40 at the channelcenter.The simulations are performed by enforcing a constantmass flux corresponding to the constant Rem ¼ 20 133,where Rem ¼ Um2h= is the Reynolds number based onthe bulk velocity Um ¼ ð1=2hÞRhh Udy, the kinematicviscosity , and the channel half-height h. The referencecase is first considered by choosing the Smagorinsky con-stant CS ¼ 0:05 which was shown to provide the bestaccordance with statistics of direct numerical simulations[30]. The reference solution exhibits an average friction-based Reynolds number Re ¼ 550, where Re ¼ hu=and u is the friction velocity. The mean streamwise ve-locity profile and the turbulence intensities agree well withthose from direct numerical simulations at the same Re[34] as reported in Fig. 1. A snapshot of the streamwisevelocity fluctuation in the reference case is reported inFig. 2(a). The structure of this field is very similar to theone obtained from direct numerical simulations, withlarge-scale streaks embedded in clouds of smaller-scalestructures (see, e.g., [28]). The premultiplied spanwiseand streamwise velocity spectra kzEuuðkzÞ and kxEuuðkxÞcorresponding to the reference case, reported in Figs. 3(a)and 3(c), respectively, are also in qualitatively good accor-dance with those from direct numerical simulation at simi-(a) (b)+y +y+++rmsrmsrmswvu , ,)(++yU101 102051015200 100 200 300 400 50000.511.522.53FIG. 1. Comparison of mean statistics of the reference case(CS ¼ 0:05) with data from the direct numerical simulations inRef. [34] at the same Re ¼ 550: (a) mean streamwise velocityprofile UðyþÞ; (b) turbulence intensity profiles urmsðyþÞ,vrmsðyþÞ, and wrmsðyþÞ.FIG. 2 (color online). Isosurface (uþ ¼ 2) of the instanta-neous streamwise velocity fluctuation for the reference simula-tion with CS ¼ 0:05 (a) and for simulations performed withincreasing artificial dissipation of the small scales withCS ¼ 0:1 (b), CS ¼ 0:2 (c), and CS ¼ 0:3 (d).lar Re [11]. In particular, in the premultiplied spectra, thescales of the buffer-layer cycle þz  100 and þx  600are observed for yþ < 30, while in the outer region thepeaks corresponding to the large-scale motions are ob-served at z=h  1:5–2 and x=h  3–4 with long tailsin the range of x=h > 3–4 that represent the very-large-scale motions [16,35,36].To test if the large-scale structures can survive even inthe absence of smaller scales, we progressively damp thesmall-scale structures by simply repeating the simulationswith larger values of the Smagorinsky constant CS, whichonly affects the diffusion term and not the nonlinear ad-vective terms. In these additional simulations the Reynoldsnumber is kept at its reference value Rem ¼ 20 133. Thefriction-based Reynolds number is not strongly affected bythe change in CS: Re ¼ 494 is found for CS ¼ 0:1, Re ¼518 for CS ¼ 0:2, and Re ¼ 560 for CS ¼ 0:3.For increasing CS values, the small-scale coherent mo-tions are increasingly damped, as can be seen from theinstantaneous streamwise velocity fluctuation fields re-ported in Fig. 2. Only large- and very-large-scale streakymotions survive for CS  0:3 [Fig. 2(d)]; also, these mo-tions are quenched for CS ¼ 0:4. The effects observed onthe instantaneous velocity fields are confirmed by theanalysis of the premultiplied spectra reported in Fig. 3.For CS ¼ 0:1, the buffer-layer peak þz  100 in the span-wise premultiplied spectrum kzEuuðzÞ has been quenched,meaning that the near-wall cycle has been suppressed.However, structures with almost constant z=y survive inthe log layer for CS ¼ 0:1, as seen in Fig. 3(c). For largerCS, these structures are also quenched. Only one peaksurvives with z  1:5h for CS ¼ 0:2 and z ¼ 2h forCS ¼ 0:3, with z corresponding well to the large- andvery-large-scale motions ½z  ð1:5–2Þh apparent inFig. 3(a) and observed in numerical simulations and ex-periments [11,15].The analysis of the streamwise premultiplied spectrakxEuuðxÞ confirms this scenario. For CS ¼ 0:1 the þx 600 peak has been replaced by much longer structures(x  7h). In the channel center, structures with x  3hdominate, like in the reference simulation. The typicalstreamwise scales in the inner layer, no longer masked bythe near-wall cycle, are now longer than in the outer layer.These scales grow longer when CS is further increased [seeFigs. 3(d), 3(f), and 3(h)]. The streamwise scales survivingfor CS ¼ 0:3 are larger by a factor of 2 when compared tothe ‘‘natural’’ large-scale and very-large-scale motions,which can probably be attributed to the very importantartificial dissipation enforced at this CS. The analysis ofthe wall-normal structure of the modes associated withthese scales reveals that the modes associated with longerx, when compared to shorter x ones, (a) have turbulenceintensity peaks nearer to the wall, (b) penetrate more to thewall with a larger turbulence intensity, and (c) have asmaller ratio of the wall-normal to the streamwise turbu-lence intensity. These structures remain essentially outer-layer structures even for the longest x. We have verifiedthat the results described above do not depend on thedetails of the damping of the small-scale motions byrepeating the simulations with a different form of artificialextra dissipation.These findings, especially for the case CS ¼ 0:1, are inaccordance with recent numerical experiments, where itwas shown that large-scale motions in the outer layers arealmost independent of the buffer-layer processes [37,38].Even more importantly, these results prove that large- andvery-large-scale motions in the turbulent channel do notnecessarily rely on the existence of smaller-scale structureslike buffer-layer streaks and vortices or hairpin vortices.This is a strong indication that an independent self-FIG. 3. Spanwise premultiplied power spectrum kzEuuðzÞ[(a), (c), (e), and (g)] and streamwise premultiplied powerspectrum kxEuuðxÞ [(b), (d), (f), and (h)] for, respectively, thereference simulation with CS ¼ 0:05 [(a) and (b)] and for thecases CS ¼ 0:1 [(c) and (d)], CS ¼ 0:2 [(e) and (f)], and CS ¼0:3 [(g) and (h)]. The premultiplied spectra are extracted in theinner layer (yþ ¼ 16, 30, 70, 108; dashed lines) and in the outerlayer (y=h ¼ 0:38, 0.65, 1; solid lines).sustained process at large scales exists in the turbulentchannel and probably in the other turbulent canonical shearflows such as the pipe and the boundary layer, as waspreviously conjectured [13,20,22].The self-sustained process at large scales is associatedwith coherent large-scale streaks that undergo sinuousoscillations and break down in an apparently randomway. This phenomenology is similar to what is observedin minimal box realizations of the buffer-layer self-sustained process and in the subcritical laminar-turbulenttransition in the channel [3,4,7,39]. In these situations, theprocess can self-sustain only in periodic boxes sufficientlylarge for the streaks to grow by the lift-up effect andsufficiently long to allow the sinuous streak oscillationsto develop. We have therefore repeated the simulations forCS ¼ 0:3 in smaller horizontally periodic boxes. A con-stant spatial resolution (x, z) is used so as to maintainthe LES filter width, unchanged. A detailed analysis ofnumerous box combinations reveals that the large-scaleprocess is sustained only if the streamwise and spanwisebox sizes are larger than the minimal values Lx;min  3hand Lz;min  1:5h, respectively. These values are in verygood agreement with the spatial length scales of the mostenergetic outer-layer motions [11,15,20,35]. This similar-ity of the minimal box sizes to the characteristic scales ofthe most energetic structures observed in the ‘‘natural’’turbulent flows is also observed for buffer-layer structures[3], The minimal spanwise size Lz;min is also well in therange of scales where streaks are strongly amplified by thecoherent lift-up effect in turbulent channels [20,21].The analysis of the power input and dissipation, IðtÞ andDðtÞ, in the minimal box Lx;min  Lz;min for CS ¼ 0:3reveals pseudoperiodic oscillations with time period T ð13–15Þh=Ub associated with sinuous oscillations of thelow speed streaks in the whole channel. These oscillationsare typically followed, in an intermittent way, by a largeexcursion in the drag and dissipation that is very similar tothe ‘‘bursting’’ observed in the buffer-layer process [3].From the minimal box results we also conclude that theexistence of motions at very large scales (x > 3h) is notnecessary for the survival of the motions at large scales(x  3h). Whether the motions with x > 3h are active orare passively forced by the process at x  3h is still anopen question that is currently under scrutiny.We gratefully acknowledge the use of the DIABLO code,mainly developed by T. R. Bewley and J. Taylor, andpartial support from DGA, E´cole Polytechnique, andNORDITA, where part of this work was done.*yongyun@ladhyx.polytechnique.fr†carlo.cossu@imft.fr[1] S. J. Kline, W.C. Reynolds, F. A. Schraub, and P.W.Runstadler, J. Fluid Mech. 30, 741 (1967).[2] J. R. Smith and S. P. Metzler, J. Fluid Mech. 129, 27(1983).[3] J. Jime´nez and P. Moin, J. Fluid Mech. 225, 213 (1991).[4] J. Hamilton, J. Kim, and F. Waleffe, J. Fluid Mech. 287,317 (1995).[5] F. Waleffe, Phys. Fluids 9, 883 (1997).[6] M. T. Landahl, J. Fluid Mech. 212, 593 (1990).[7] F. Waleffe, Stud. Appl. Math. 95, 319 (1995).[8] W. Schoppa and F. Hussain, J. Fluid Mech. 453, 57 (2002).[9] J. Jime´nez and A. Pinelli, J. Fluid Mech. 389, 335(1999).[10] L. S. G. Kovasznay, V. Kibens, and R. F. Blackwelder,J. Fluid Mech. 41, 283 (1970).[11] J. Jime´nez, Annual Research Briefs, Center for TurbulenceResearch, Stanford University (1998).[12] C. D. Tomkins and R. J. Adrian, J. Fluid Mech. 490, 37(2003).[13] M. Guala, S. E. Hommema, and R. J. Adrian, J. FluidMech. 554, 521 (2006).[14] N. Hutchins and I. Marusic, J. Fluid Mech. 579, 1 (2007).[15] B. J. Balakumar and R. J. Adrian, Phil. Trans. R. Soc. A365, 665 (2007).[16] K. C. Kim and R. Adrian, Phys. Fluids 11, 417 (1999).[17] A. Townsend, The Structure of Turbulent Shear Flow(Cambridge University Press, Cambridge, England, 1976).[18] X. Wu and P. Moin, J. Fluid Mech. 630, 5 (2009).[19] S. Toh and T. Itano, J. Fluid Mech. 524, 249 (2005).[20] J. C. del A´lamo and J. Jime´nez, J. Fluid Mech. 559, 205(2006).[21] G. Pujals, M. Garcı´a-Villalba, C. Cossu, and S. Depardon,Phys. Fluids 21, 015109 (2009).[22] C. Cossu, G. Pujals, and S. Depardon, J. Fluid Mech. 619,79 (2009).[23] Y. Hwang and C. Cossu, J. Fluid Mech. 643, 333 (2010).[24] Y. Hwang and C. Cossu, J. Fluid Mech. (to be published).[25] O. Kitoh and M. Umeki, Phys. Fluids 20, 025107(2008).[26] G. Pujals, C. Cossu, and S. Depardon, in Sixth Symposiumon Turbulence and Shear Flow Phenomena (SeoulNational University, Seoul, Korea, 2009).[27] G. Pujals, C. Cossu, and S. Depardon, J. Turbul. (to bepublished).[28] J. Jime´nez, Rev. R. Acad. Cien. Serie A Mat. 101, 187(2007).[29] J. Smagorinsky, Mon. Weather Rev. 91, 99 (1963).[30] C. Ha¨rtel and L. Kleiser, J. Fluid Mech. 356, 327 (1998).[31] C. Meneveau and J. Katz, Annu. Rev. Fluid Mech. 32, 1(2000).[32] T. R. Bewley, P. Moin, and R. Temam, J. Fluid Mech. 447,179 (2001).[33] T.A. Zang, Phil. Trans. R. Soc. A 336, 95 (1991).[34] J. del A´lamo and J. Jime´nez, Phys. Fluids 15, L41 (2003).[35] J. C. del A´lamo, J. Jime´nez, P. Zandonade, and R.D.Moser, J. Fluid Mech. 500, 135 (2004).[36] J. P. Monty, N. Hutchins, H. C. H. Ng, I. Marusic, andM. S. Chong, J. Fluid Mech. 632, 431 (2009).[37] O. Flores and J. Jime´nez, J. Fluid Mech. 566, 357 (2006).[38] O. Flores, J. Jime´nez, and J. del A´lamo, J. Fluid Mech.591, 145 (2007).[39] S. C. Reddy, P. J. Schmid, J. S. Baggett, and D. S.Henningson, J. Fluid Mech. 365, 269 (1998).

On a self-sustained process at large scale in the turbulent channel flow

Abstract

Similar works

Full text

Available Versions

HAL: Hyper Article en Ligne

Open Archive Toulouse Archive Ouverte

Crossref

HAL-Polytechnique