arxiv: v1 [physics.flu-dyn] 10 Dec 2018

Transcription

1 Folding of he frozen-in-fluid di-voriciy field in wo-dimensional hydrodynamic urbulence arxiv: v1 [physics.flu-dyn] 10 Dec 2018 E.A. Kuznesov a,b,c and E.V. Sereshchenko c,d,e b P.N. Lebedev Physical Insiue, RAS, Moscow, Russia c L.D. Landau Insiue for Theoreical Physics, RAS, Moscow, Russia a Novosibirsk Sae Universiy, Novosibirsk, Russia d A.S. Khrisianovich Insiue of Theoreical and Applied Mechanics, SB RAS, Novosibirsk, Russia e Far-Easern Federal Universiy, Vladivosok, Russia The voriciy roor field B = ro ω (di-voriciy) for freely decaying wo-dimensional hydrodynamic urbulence due o a endency o breaking is concenraed in he viciniy of he lines corresponding o he posiion of he voriciy quasi-shocks. The maximum value of he di-voriciy a he sage of quasi-shocks formaion increases exponenially in ime, while he hickness l() of he maximum area in he ransverse direcion o he vecor B decreases in ime also exponenially. I is numerically shown ha () depends on he hickness according o he power law: () l α (), where he exponen α 2/3. This behavior indicaes in favor of folding for he divergence-free vecor field of he di-voriciy. PACS: Cv, a, Ra I. INTRODUCTION In he wo-dimensional developed hydrodynamic urbulence a high Reynolds numbers, Re 1, as i was found ou in [1 3], he Kraichnan direc cascade [4] wih a consan ensrophy flux is formed due o he appearance of he voriciy quasi - shocks (or jumps), because of he compressibiliy of coninuously disribued lines of he field of he voriciy roor, B = roω, ofen called he di-voriciy afer Kida [5]. This compressibiliy propery follows direcly from he equaion for B, B = ro[v B], divv = 0, (1) which has he form of a frozenness equaion. From his equaion we have ha B changes only by virue of he velociy componen v n normal o he di-voriciy vecor (in consequence of he vecor produc). The velociy componen v n due o frozenness B deermines he change in he posiion of he force lines of his field by means of he Lagrangian rajecories, which are defined from he equaions dr d = v n(r,) = 0, r =0 = a. In he general siuaion, divv n 0, and herefore he mapping r = r(a,) as he soluion of hese equaions is urned ou o be compressible. This fac also follows from he Liouville equaion for he mapping Jacobian J = (x,y)/ (a x,a y ) (as a measure of he infiniesimally small area variaion): dj d = divv n J. Thus, Jacobian J can ake arbirary values, including zero. This is he reason for he compressibiliy of coninuously disribued divoriciy lines and, accordingly, he endency o breaking, ha resuls in he formaion of voriciy quasi-shocks. In he case of freely decaying urbulence, his process is dominan, leading o a srong anisoropy of he urbulence specrum because of he presence of jes generaed by quasishocks [1, 2]. This process urns ou o be he fases, as a resul he urbulence specrum of he direc cascade a he iniial sage forms a power dependence on he wave numberk wih he Kraichnanype behaviorforhe specrum: E k k 3 (see he original paper by Kraichnan [4]), even in he presence of pumping, as shown by numerical experimens[8]. A he same ime, he formaion of he voriciy quasi-shocks is exponenial; in accordance wih his, he regions of he maximum of di-voriciy are decreased in he direcion perpendicular o he lines of he consan voriciy. As demonsraed in he numerical experimens[1 3]), for

2 2 ypical iniial condiions an increase of he divoriciy modulus consiss orders of magniude, and he ransverse size of he maximal area B decreases significanly. The explanaion of his growh is relaed o he possibiliy of parial inegraion of he equaion (1) in erms of mapping r = r(a,) [1]: B(r,) = (B 0(a) a )r(a,), (2) J whereb 0 (a)isheiniialb, whichhasameaning of he Cauchy invarian analogue. A similar formula for he hree-dimensional Euler equaions is basic in he so-called vorex line represenaion [9]. The key poin here for undersanding is he compressibiliy of he divoriciy field and he possibiliy of J o vanish. As is known, breaking in he gas dynamics occurs due o he compressibiliy of he gas. The formaion of quasi-wo-dimensional causics occurs when approaching he breaking poin(see, e.g. [10]). The formaion of he voriciy quasi-shocks happens similarly. I is necessary o menion ha probably firs ime appearance of he voriciy quasi-shocks was observed in he numerical simulaions by Lilly in 1971 [6]. I was one of he main moivaion for Saffman o sugges he specrum E(k) k 4 [7] differen from ha suggesed by Kraichnan [4]. The firs explanaion of hese facs was given in he papers [1] based on he represenaion (2) in In his paper, we invesigae how he maximum value of he di-voriciy varies depending on he hickness of he maximum area in order o find ou wheher his process can be considered as a fold formaion (wha is a fold - see, for example, [11]). As a resul of numerical simulaion on he grid 16384x16384, we found ha beween he maximum value of and he hickness of l, a he sage of exponenial growh, here is a power law dependence: l α, where he exponen α is close o 2/3. I should be noed ha such a dependence was found while formaion of he pancake-ype vorex srucures arising for inviscid hree-dimensional flows [12]. This resul indicaes ha he formaion of quasi-shocks can be considered as a process of folding for a divergen-free vecor field - he di-voriciy field. If for he hree-dimensional Euler equaion, he appearance of a power dependence beween he maximum voriciy ω max and he pancake hickness l of he form ω max l 2/3 could be aribued o he Kolmogorov ype relaion, hen he dependence l 2/3 indicaes ha in he wo-dimensional Euler hydrodynamics we are dealing wih folding ha is no relaed o he Kolmogorov behavior.the frozenness of boh he voriciy field in he hree-dimensional Euler equaion and he divoriciy field combines boh of hese cases - he formaion of hree-dimensional vorex srucures of pancake ype and quasi-shocks of voriciy for wo-dimensional flows. Despie he incompressibiliy of he velociy field, boh fields - hree dimensional voriciy field and divoriciy are compressible (see [12]). Apparenly, due o his propery of frozen-in-fluid fields, we can expec ha he 2/3 law is universal for any fields of his ype. II. NUMERICAL RESULTS In his paper for he numerical inegraion of he wo-dimensional Euler equaions wrien in erms of voriciy (Helmholz equaion), ω +(v )ω = 0, (3) we use algorihm described in deail in our previous papers [2, 3, 8]. The inegraion domain for he equaion (3) was a square box wih sizes 2π 2π, he boundary condiions were periodic a boh coordinaes. Velociy v and voriciy ω were found hrough he sream funcion ψ using he sandard formulas. For a given value of ω he sream funcion was found by reversing he Laplace operaor wih he help of he fas Fourier ransform and hen Bmax FIG. 1: Disribuion of di-voriciy B a = 12. Bmax

3 3 7,0 1 6,0 5,0 4,0 3,0 2,0 y=1.65exp(0.16x) y=1.65exp(0.16x) 1,0 1,0 FIG. 2: The dependence of he maximum di-voriciy on ime (on he lef is he usual scale, on he righ is he logarihmic). The poins correspond o he numerical resuls, and he line o exponenial fiing. 0,30 1,000 0,25 0,20 y=0.67exp(-0.25x) y=0.67exp(-0.25x) 0,15 0,100 0, FIG. 3: Dependence of hickness on ime (on he lef is he usual scale, on he righ is he logarihmic). The poins correspond o he numerical resuls, and he line o exponenial fiing. he velociy was deermined. As in he our previous works, wo ses of Gaussian vorices wih posiive and negaive voriciy wih zero oal voriciy were used as iniial condiions. The size of each pair was random in he range of , he locaion of vorices was also random. Unlike previous works [2, 3], we limied he number of vorices o 8 (4 posiive and 4 negaive) o more accuraely deermine he required dependences - he B field and is geomeric characerisics: maximum posiions, longiudinal and ransverse quasi-shock sizes, ec. A firs, he spaioemporal dependences of voriciy were numerically found and hen he emporal evoluion of he field of di-voriciy was deermined. Analysis of he disribuion of he di-voriciy field showed ha is main mass is concenraed in he small viciniy of he lines wih a maximum value of B in he form of narrowing wo-dimensional ribbons ha form a complex ne and a large imes becomes urbulen. Fig. 1 shows he srucure of B a = 12. Zoom shows ha beween he maximum lines of di-voriciy values of B is significanly less han he maximum. For voriciy his corresponds o a errace wih seps of variable heigh. Each of hese seps is a voriciy quasi-shock. As already noed, he maximum ampliude of he di-voriciy a he iniial sage increases exponenially and hen afer reaching is maximum performs small oscillaions near his maximum. Fig. 2 shows he dependence of on ime a he iniial sage. As can be seen from hese figure, exponenial growh is observed a imes from = 4 o = 12 wih he growh rae To find he ransverse size in he viciniy of he maximum, he Hessian marix i j B was calculaed a he maximum poin and hen by is eigenvalues he ransverse (relaive o he ribbon) size ( l) and longiudinal l 2 were defined as l i = (2 λ i / ) 1/2, where λ i are he eigenvalues of he i j B marix. Fig. 3 shows he dependence of he hickness on ime: a firs pracically does no change and hen saring wih = 4 here is an exponenial decrease o = 14 wih he nega-

4 4 0,7 0,6 l 2 1,000 l 2 0,5 0,4 y=0.27exp(0.06x) y=0.27exp(0.06x) 0,3 0,2 0,1 0,100 FIG. 4: Dependence of he longiudinal size l 2 on ime (on he lef is he usual scale, on he righ is he logarihmic). The poins correspond o he numerical resuls, and he line o exponenial fiing. ive growh rae A = 14 he number of grid poins on he hickness of maximal divoriciy line was abou 20. For longer imes, his number decreases. So a = 20 he hickness reaches a value of , which is comparable o he grid size 2π/16384 = , i.e. he furher analysis becomes already incorrec. I should be noed ha for > 14 exponenialgrowhofb max sopsand sauraion akes place accompanying by small oscillaions (see ref. [1, 2]). As for he longiudinal scale l 2, i grows exponenially up o = 14 wih he growh rae 0.06 (see Fig. 4). A large imes l 2 growh sops, and hen here is even an exponenial fall, which in our opinion is associaed wih he discreeness of he compuaion grid and herefore is incorrec. I is imporan o noe ha up o = 14 he oal value ofhe growh raes(for boh and l 2 ) urns ou o be negaive: = y=1.27x ,10 0,15 0,20 0,25 FIG. 5: The dependence of he maximum divoriciy on he hickness l. The poins correspond o numerical resuls, and he line corresponds o he power dependence l 2/ , which indicaes he compressibiliy of he region of he maximal value. Similar behavior is observed in hree-dimensional hydrodynamics when forming pancake-ype vorex srucures [13]. The obained dependences for and hickness l show ha a he exponenial sage beween hese values here is a power dependence = Cl α wih α = 0.16/( 0.25) = /3, C is a consan (see Fig. 5). I is worh noing ha his dependence of on l in he form of he 2/3 law was also verified for anoher iniial condiions (recall ha he posiions of vorices and heir sizes were random). This allows ones o believe ha his relaion can be considered as universal. III. CONCLUSION The main conclusion of his work is ha he formaion of he power dependence of on he hickness l - he 2/3 law - can be considered as folding, we emphasize, for he divergence-free vecor field of he di-voriciy. As already noed in he inroducion, he same scaling was found for he formaion of he pancake-ype vorex srucures for hreedimensional Euler hydrodynamics [12, 13]. In his case, he scaling 2/3 was se beween he maximal voriciy and he pancake hickness: ω max l 2/3. Boh hese phenomena are unied by he frozenness propery of he voriciy field for he hree-dimensional flows and he di-voriciy for wo-dimensional hydrodynamics. Is his law universal for any frozenin-fluid fields for incompressible flows? This quesion remains sill open; especially i is in-

5 5 eresing for magneic hydrodynamics where, in he non-dissipaive limi, he magneic field is also frozen. This work was carried ou wih he suppor of he Russian Foundaion of Basic Research (gran no ). [1] E.A. Kuznesov, V. Naulin, A.H. Nielsen, and J.J. Rasmussen, Phys. Fluids 19, (2007); Theor. Compu. Fluid Dyn. 24, (2010). [2] A.N. Kudryavsev, E.A. Kuznesov, E.V. Sereshchenko, JETP Leers, 96, (2013). [3] E.A. Kuznesov, E.V. Sereshchenko, JETP Leers 102, (2015) [4] R. Kraichnan, Phys. Fluids, 11, 1417 (1967). [5] S. Kida, J. Phys. Soc. Jpn. 54, 2840 (1985). [6] D.K. Lilly, J. Fluid Mech. 45, 395 (1971). [7] P.G. Saffman, Sud. Appl. Mah. 50, 377 (1971). [8] E.A. Kuznesov, E.V. Sereshchenko, JETP Leer 105, 83-88, (2017). [9] E.A. Kuznesov, V.P. Ruban, JETP Leers, 67, (1998); E.A. Kuznesov, V.P. Ruban, Phys. Rev. E, 61, 831 (2000). [10] S.F. Shandarin, Ya.B. Zeldovich, Rev. Mod. Phys. 61, 185 (1989). [11] V.I. Arnold, Caasrophe Theory (3rd ed.), Berlin: Springer-Verlag, [12] D.S. Agafonsev, E.A. Kuznesov and A.A. Mailybaev, Phys. Fluids 30, (2018). [13] D.S. Agafonsev, E.A. Kuznesov and A.A. Mailybaev, Phys. Fluids 27, (2015).