arxiv:chao-dyn/ v1 11 Jun 1997

Transcription

1 arxiv:chao-dyn/ v1 11 Jun 1997 ON THE TRIPLE CORRELATIONS IN HELICAL TURBULENCE 1 O.G.Chkhetiani Space Research Institute, Russian Academy of Sciences, Profsoyusnaya 84/32, Moscow, , Russia e mail:ochkheti@mx.iki.rssi.ru Abstract The evolution of correlation characteristics in homogeneous helical turbulence is considered. Additional Kármán Howarth type equations, describing the evolution of the mixed correlation tensor of the velocity and vorticity are obtained. In the helical scaling region, the solution of obtained equation gives the exact relation between antisymmetric component of a rank-three tensor and the average dissipation of helicity. This relation is a helical analogue of Kolmogorov s known 4 5 law [1]. Helicity plays an important role in the evolution and stability of turbulent and laminar flows [2]. In [3] the conception of the helical cascade was introduced and the limiting cases of parallel energy and helicity flows along the spectrum E(k) ε 2/3 k 5/3, H(k) η ε 1/3 k 5/3, corresponds to Kolmogorov cascade, and a helicity flux with no energy flux. E(k) η 2/3 k 7/3, H(k) η 2/3 k 4/3 pure helical cascade ( ε, η - dissipation of energy and helicity). Many years pure helical cascade is considered as mainly theoretical exercise without any practical significance. However, today we have nature and laboratory results that is a direct evidence of minus seven-third law. In turbulent MHD flow in channel [4], tropical pretayphoon atmosphere [5] spectra of energy E(k) k 7/3, quantitatively and qualitatively corresponds to the helical cascade observe. In particular in tropical zone the helical cascade observes on scales l h > 6 7 km [5]). 1 Slightly extended version of the article On the third moments in helical turbulence, published in JETP Lett, 63(10), 802,

2 Such a spectrum arises a problem of helicity generation. In [6, 7] was shown, that for total helicity generation the simultaneous presence of rigid body rotation and temperature dencity stratification or other inhomogeneous factors be present simultaneously. In [7] also was emphasized, that in this case helicity is pumping in system on all turbulent scales. Helicity also generates in Eckman boundary layers [8], in convective liquids [9, 10, 11] and at aerodynamic flow [12]. The dissipative mechanism of helicity generation, first demonstrated in [13] should be noted also. Direct measurements have provided in turbulent boundary layer, mixing layer [14] and lattice turbulence [15, 16, 14, 17, 18] are demonstrated the nonzero values of mean helicity. Is noted also, that small helicity perturbations, made by an artificial source before lattice are picked up and amplify turbulence by a flow at lattice. Moreover, small, but the finite mean helicity values at lattice were received even at the absence of an obvious its source. Such phenomenon can be explained as influence of weak uncollinearity of lattice apertures on the mirror symmetry of turbulent flow [14]. In [11] helical cascades in a stratified and compressible turbulent media were studied in detail. It was noted that there exists an inner turbulent scale that separates the regions of Kolmogorov and helicity cascades: l h ε η. As rule, helicity cascade is observed on large scales (atmosphere) and for more small scales we observe a Kolmogorov minus five third law. However, if the helicity source lies on small scales [19] we have reverse situation: helical scaling region lies on small scales (k l h ) (see Figure 1). The scale invariance of characteristic functional allows to obtaine a form of helical scaling in stratified turbulent fluid [11] E(k) = η 2/3 k 7/3 Ψ ( (kl ) 1/2, kl h ), (1) E TT (k) = Nη 1/3 k 4/3 Ψ TT ( (kl ) 1/2, kl h ) L = l h (l h /L BO ) 4, L BO Bolgiano Obukhov scale, and in compressible fluid [11] E(k) = γ 2γ γ 1 (γ 1)(3γ 1) ρ 3γ 1 C 2 3γ 1 0 η 2γ 3γ 1 k 7γ 1 2 3γ 1 Ψ (klh, kl), (2)

3 E(k) 7 3 E(k) l 1 h pretyphoon atmosphere k l 1 h laboratory MHD flow k Figure 1: C 0 sound velocity. In the shock wave limit γ 1 we have k 3 spectrum instead Kadomtsev Petviashvili k 2 E(k) = eηc 1 0 k 3 Ψ (kl h, kl) The presence of helicity in turbulent system naturally puts a problem of its influence to correlation properties of flow, in particular on high correlations (and accordingly about possibility of its measurement). Not doing of hypotheses about nature of the cascade we shall consider evolution of the correlation characteristics in homogeneous helical turbulence. At the initial moment of time t = 0 exist nonzero average helicity. Double correlation v 1i v 2j has a form v 1i v 2j = A(r)δ ij + B(r)r i r j + Cε ikl r l, r = r 2 r 1. (3) Obviously v 2 = 3A(0), vrotv = 6C(0). We shall consider correlation of the third rank tensor b ik,l = v 1i v 1k v 2l = v 2i v 2k v 1l General form of similar tensor following [20, 21] b ik,l = S 1 (r)δ ik n l + S 2 (r)δ il n k + S 3 (r)δ nl n i + S 4 (r)n i n k n l 3

4 + S 5 (r)ε ikl + S 6 (r)ε ikt n t n l + S 7 (r)ε ilt n t n k + S 8 (r)ε klt n t n i, (4) Where n = r/ r. The antisymmetric part tensor of the third rank (proportional to ε... ) usually was not considered, as it differ from 0 only at helicity presence. Tensor b ik,j satisfies to the following conditions of symmetry [20, 21] And incompressibility condition Whence we receive, that Thus b ik,j = b ki,j, b i,kj = b i,jk, b i,kj = b kj,i (5) b ik,j r j = 0, (6) S 5 = S 6 = 0, S 7 = S 8. b ik,l = D(r)δ ik n l + E(r)(δ il n k + δ nl n i ) + F(r)n i n k n l +S(r) (ε ilt n t n k + ε klt n t n i ), (7) (The symmetric part of tensor is in detail considered in [22, 20, 21] and, as we shall see below, does not influence on received results). We shall show, that the antisymmetric part, proportional S(r), is directly connected with flow helicity. We shall consider mixed point-to-point correlation tensor, component determined by product of velocity and vorticity v 1i w 2j, and having in homogeneous isotropic case the following kind: ( v 1i w 2j = 2C(r) + r C(r) ) δ ij + 1 C(r) r i r j r r r ( 5B(r) + r B(r) ) ε ijl r l. (8) r It is easy to obtain the equations for its components t v 1iw 2j = v 1l v 1i w 2j v 1i v 2l w 2j + v 1i w 2l v 2j x 1l x 1l x 1l x 1i p 1 w 2j + ν 1 v 1i w 2j + ν 2 v 1i w 2j (9) 4

5 These equations, are similar to Kármán Howarth equations and are additional to them at the helicity presence in turbulent system. Owing to solenoidality and regularity of correlation functions at r = 0 p 1 w 2j = 0 In the case of homogeneous turbulence we can extract derivative from averaging braces t h ij = ε jst 2 r l r s [b li,t b lt,i ] + 2ν r h ij (10) After simple transformations (the symmetric components tensor b li,t give the zero contribution) we receive t h ij = 2 {[ 2 S r + 2 S 2 r r And for i = j we receive 2ν C(r) = t r 4 r ] 6S n r 2 i n j (r 4 r C(r) ) [ 2 S r r ] } S δ ij + ν h ij, (11) r + 2 ( r 3 S(r) ). (12) r 4 r Assuming, that at r the order of functions C(r)/ r and S(r) higher, than 1/r 4 and 1/r 3 accordingly, and that νr 4 C(r)/ r+r 3 S(r) tends to zero at r = 0, we shall find a helical analogue of Loytsansky invariant, which is conserved on a final stage of free turbulence decay [22] Λ h = 0 C(r)r 4 dr = const. (13) This invariant concerns with the conservation of the mean product of the vortex momentum and angular momentum P M. Having presented C(r) in form C(r) = 1 6 ( vrotv + Ĉ(r) ), We receive η + 2ν = (r 4 rĉ(r) ) tĉ(r) r 4 r ( r 3 S(r) ), (14) r 4 r

6 vi w i Where η = 2ν is the mean helicity dissipation [11]. x j x j The helicity flow on spectrum is constant as in Kolmogorov, and in helical interval [3]. Therefore it is possible with sufficient accuracy to neglect the change in time of Ĉ(r) in comparison with η. Product on r4 and integrate on r we receive the following expression for S(r) S(r) = η 60 r2 νr (C(r)). (15) r It is obvious, that in depth of an inertial interval the viscous member is small and we have simply S(r) = η 60 r2. (16) It s easy to obtaine the first terms of the Taylor expansion of C(r) for small r: C(r) = 1 ( v curl v + η ) 6 20ν r2 (17) This form of C(r) turn to zero the expression for S(r). Evidently, its concerned with the more high order to r dependence of S(r) ( r 4 ) in small scale viscous region. Received exact relation is similar to known 4/5 law, connecting thriple longitudinal correlation of velocity with dissipation of energy [1]. We shall note, that at its derivation we didn t use any hypotheses about character of scaling. If in turbulent to system on which or reasons arise helicity, occur additional nonzero components 2 pointed correlation tensor of the third rank, proportional to the mean helicity dissipation. In particular for correlation v i (x)v j (x)w j (x + r) we receive v i (x)v j (x)w j (x + r) = η 10 r i. (18) Having chosen axis z parallel to a vector r we receive b xy,z = b y,xz = 0, b zy,x = b yz,x = S(r), b xz,y = b z,xy = S(r), b xy,z = b y,xz = 0, b zy,x = b yz,x = S(r), b xz,y = b zx,y = S(r). The observation of such fine structure of high correlations of velocity should be the direct and obvious evidence of helicity presence in turbulent system. 6

7 It should be noted, that in [11] expression for point-to-point semiinvariant any rank in the helical scaling region is obtained v i1 (x)v i2 (x)...v in (x + r) = const n (ηr 2) n 3 Θ i1,...,i n (n), (19) Where Θ i1,...,i n (n) - an angular part of spectral tensor. (In nonhelical case accordingly (εr) n 3 ). In summary it would be necessary to express gratitude S.S.Moiseev and A.V.Belyan for stimulating discussions. The work is carried out at support RFFI (grant ), Russian Ministry of a Science, INTAS (grant INTAS ). References [1] A.N.Kolmogorov, Doklady Akademii Nauk (1941), 32 (1), 19 (in russian). [2] H.K.Moffat, A.Tsinober, Ann. Rev. Fluid Mech., 24, 281 (1992). [3] A.Brissaud A. Et all, Phys. Fluids, 16, 1363 (1973). [4] H.Branover, A.Eydelman, M.Nagorny, M.Kireev, in Progress in Turbulence Research Progress in Astronautics and Aeronautics, 162, 64 (1994). [5] I.N.Klepikov, I.V.Pokrovskaya, E.A.Sharkov, Issledovaniya Zemli iz Kosmosa, N 3, 13 (1995). [6] M.Steenbeck, F.Krause, and K.-H.Rädler, Z.Naturforsch A, 21, 364 (1966) [7] S.I.Vainshtein, Ya.B.Zeldovich, A.A.Ruzmaikin, Turbulent dynamo in Astrophysics, Moscow: Nauka Publ (in russian). [8] R.Hide, Geophys.Asrophys.Fluid Dyn., 48, 69 (1989). [9] D.K.Lilly, J.Atm.Sci., 43, 126 (1986). [10] M.V.Kurgansky, Fizika Atmosfery i Okeana 29, 464 (1993). 7

8 [11] S.S.Moiseev, O.G.Chkhetiani, JETP, 83(1), 192 (1996). [12] J.C.R.Hunt, F.Hussain, J.Fluid Mech., 229, 569 (1991) [13] J.C.André, M.Lesieur, J.Fluid Mech., bf 81, 181 (1977). [14] J.M.Wallace, J.L.Balint, in Topological Fluid Mechanics. Cambridge, Cambridge Univ.Press 1990, p.585. [15] E.A.Kit et all, Phys.Fluids, 30, 3323 (1987). [16] E.A.Kit et all, Fluid Dyn. Res., 3, 289 (1988). [17] T.Dracos et all, in Topological Fluid Mechanics. Cambridge, Cambridge Univ.Press 1990, p.564. [18] M.Kholmyansky et all, Fluid. Dyn. Res., 7, 65 (1991). [19] O.G.Chkhetiani, E.I.Golbraikh, S.S.Moiseev ( unpublished results) [20] I.O.Hince, Turbulence, Moscow: Science Publ [21] S.Chandraseckhar in Selected papers. v.3, Chickago Press 1989, p.498. [22] L.D.Landau, E.M.Livshits, Hydrodynamics, Moscow: Science Publ