arxiv: v2 [physics.class-ph] 4 Apr 2009

Transcription

1 arxiv: v2 [physics.class-ph] 4 Apr 2009 Geometric evolution of the Reynolds stress tensor in three-dimensional turbulence Sergey Gavrilyuk and Henri Gouin Abstract The dynamics of the Reynolds stress tensor is determined by an evolution equation coupling geometrical effects and turbulent source terms. The effects of the mean flow geometry are shown up when the source terms are neglected. Then, the Reynolds stress tensor is expressed as the sum of three tensor products of vector fields, which are governed by a distorted gyroscopic equation. Along the mean flow trajectories and in the directions of the vector fields, the fluctuations of velocity are determined by differential equations whose coefficients depend only on the mean flow deformation. 1 Introduction The governing equations of barotropic turbulent compressible fluids are see, for example [1], [2], [3] ρ t + ρ U i, i = 0, ρ U i t + ρ U i U j + p δ ij + ρ u i u j, j = 0, 1 ρ u i u j t + ρ u i u j U k, k + ρu i u j U j,k + ρ u i u j U i,k = S ij. Here brackets mean the averaging, coma means the derivation with respect to the Eulerian coordinates x = {x i }, i {1, 2, 3}, index t means the partial derivative with respect to time, ρ is the fluid density, U = {U i }, i {1, 2, 3} is the mass average velocity, p is the pressure, u = {u i }, i {1, 2, 3} is the velocity fluctuation verifying ρu = 0. Repeated indices mean summation. The term S = {S ij } represents turbulent sources. We introduce the Reynolds stress tensor R = ρ u u, R ij = ρu i u j. Université d Aix-Marseille & C.N.R.S. U.M.R. 6595, IUSTI, Project SMASH, 5 rue E. Fermi, Marseille Cedex 13 France, sergey.gavrilyuk@polytech.univ-mrs.fr Université d Aix-Marseille & C.N.R.S. U.M.R. 6181, M2P2, Case 322, Av. Escadrille Normandie-Niemen, Marseille Cedex 20 France, henri.gouin@univ-cezanne.fr 1

2 The system 1 can also be rewritten in the tensorial form ρ + div ρ U = 0, t ρ du + p + divrt = 0, 2 dr T U U + R divu + R + R = S, where d/ means the material derivative with respect to the mean motion d = t + UT. The superscript T means the transposition. By using the mass conservation law, the equation for the volumic Reynolds tensor R can be rewritten as the equation for the specific or per unit mass Reynolds tensor where dp + U T U P + P = S ρ, 3 P = R ρ. The aim of this note is to understand the structure of homogeneous Reynolds equation 3 corresponding to S = 0. The reason is twofold: - first, in the numerical study of compressible turbulent flows, this is a natural step in applying the splitting-up technique see for example [4], - second, the system 2 also appears as an exact asymptotic model of weakly shearing flows of long waves over a flat bottom [5] h + div hu = 0, t h du + gh divr T = 0, 4 dr T U U + R divu + R + R = 0. In Eq. 4, h is the fluid depth, the average pressure is given by p = gh 2 /2, g is the gravity acceleration, and R = h 0 h Ũ U Ũ Udz, hu = Ũdz, 2 0

3 where Ũ is an instantaneous velocity. Equations are written for three-dimensional long waves and the production term is zero in the limit of weakly shearing flows. Also, Eqs. 4 are hyperbolic. Finally, we will focus on the equation of the Reynolds tensor per unit mass dp + U T U P + P = 0. 5 The particular case rotu = 0 was investigated in [6]. In such a case U/ T = U/ and Eq. 5 corresponds to a two-covariant tensor convected by the mean flow. This means that P has a zero Lie derivative d L with respect to the velocity field U and the tensor P 0, image of P in Lagrange coordinates t,x, only depends on X = {X i }, i {1, 2, 3} d L P dp + U P + P U = 0, P = F T 1 P0 XF 1, where F = / X is the deformation gradient of the mean motion [7]. The aim of the paper is the study of the homogeneous Reynolds tensor equations structure 5 in the case rotu 0. 2 Geometrical properties of the Reynolds tensor evolution The Reynolds stress tensor P is symmetric, semi-positive definite and consequently can be rewritten in a local basis of orthonormal eigenvectors in the form P = λ 2 α e α e α λ 2 α e α e T α. The eigenvalues λ 2 α, α {1, 2, 3} are non-negative; the case λ2 α > 0 is a generic one. For the two-dimensional case, λ Let us denote Then, From Eq. 6, we deduce a α = λ α e α, λ α > 0 α {1, 2, 3}. P = dp = a α a α da α By using Eq. 7, Eq. 5 can be written daα a α a T α. 6 a T α + a α T daα. 7 + U [ a α a T α + daα + U T a α aα] T =

4 The vector da α /+ U/a α can be developped in the local basis {a α }, α {1, 2, 3} of eigenvectors; one obtains da α + U a α = A βα a β, α {1, 2, 3} 9 β=1 where A βα, α, β {1, 2, 3} are the scalar components to be determined. By using Eq. 9, Eq. 8 leads to A αα g αα + α β=1 A αβ + A βα g αβ = 0, where g αα = 2a α a T α and g αβ = g βα = a α a T β + a β a T α, α, β {1, 2, 3}, are six independent symmetric tensors. Consequently, A αα = 0 and A αβ + A βα = 0, α, β {1, 2, 3}. Equation 8 is equivalent to da α + U a α = Λ i πe α with π =A 32 e 1 + A 13 e 2 + A 21 e 3, α {1, 2, 3}, 10 where diagonal matrix Λ and antisymmetric matrix i π are determined in the basis e α, α {1, 2, 3} as Λ = λ λ λ 3, i π = 0 A 21 A 13 A 21 0 A 32 A 13 A 32 0 The vectors a α, α {1, 2, 3} are orthogonal, a T α a β = 0, α β. If we assume the orthogonality at initial instant, this is equivalent to. a T α da β + a T β da α = So, Eqs yield α β {1, 2, 3}, a T α Λ i πe β U a β + a T β Λ i πe α U a α = 0, or 2λ α λ β e T α D e β = a T α Λ i πe β + a T β Λ i πe α λ 2 αe T α i πe β + λ 2 βe T β i πe α 12 4

5 where D = 1 2 U + T U is the rate of deformation tensor corresponding to the mean flow. We denote the mixed product of three vectors {a,b,c} as a,b,c a T b c. Hence, Eq. 12 is 2λ α λ β e T α D e β = λ 2 α λ2 β eα, π,e β = λ 2 β λ2 α π,eα,e β = λ 2 β λ2 α π T e γ, where {α, β, γ} is a cyclic permutation of the triplet {1, 2, 3}. Finally, we get Equation 10 can be written Since π T e γ = 2λ αλ β e T α D e β λ 2 β. 13 λ2 α λ α de α + dλ α e α + λ α U e α Λ i πe α = 0, α {1, 2, 3}. 14 e T de α α = 0, e T α Λ i πe α = λ α e T α i πe α = 0, by multiplying the left side of Eq. 14 with e T α, we get dλ α + e T α D e α λα = 0. By multiplying the left side of Eq. 14 with the projector I e α e T α, we get de α λ α + λ α I eα e T U α e α Λ i πe α = 0. Due to the fact there exists a vector Π such that vector Π verifies the condition de α = Π e α, λ α i Πe α + λ α I eα e T U α e α Λ i πe α = By multiplying Eq. 15 with e T β where β α,we get λ α Π T e γ = λ β π T e γ λ α e T β 5 U e α. 16

6 By replacing Rel. 13 into Eq. 16 we get T λ 2 U U β et α + e β Π T e γ = λ 2 β λ 2 α = e T β λ 2 β We can now formulate the following result: e T β U e α T U + λ 2 U α e α λ 2 β λ2 α Theorem 1. The Reynolds tensor can be written in the form R = ρ λ 2 α e α e α.. The eigenvectors e α and the eigenvalues λ α verify the equations: de α = Π e α, α {1, 2, 3} d Ln λ 2 α = 2 µ α, 17 where µ α = e T α D e α, Π T e γ = e T β λ 2 β T U + λ 2 U α e α λ 2 β λ 2. α The triplet {α, β, γ} corresponds to a cyclic permutation of the triplet {1, 2, 3}. Equation 17 1 is similar to the equations of a rigid body [8]. The vectors e α form a natural moving frame G = {e α } 3 depending on the mean deformation. The eigenvalues λ 2 α of the Reynolds tensor are determined by the evolution equation Note that if λ α are initially positive, they will be positive at any time. Hence, the tensor P is always positive definite. Due to the mass conservation law 2 1 and Eq. 17 2, we obtain d 3 ρ 2 λ 2 α = 0. 6

7 Consequently, system 17 admits a scalar invariant along the trajectories of mean flow. This invariant was earlier obtained in [6] in a different form. Let us introduce the turbulent specific energy e T = 1 2 trp = 1 2 λ 2 α. In the incompressible isochoric case, we have d ρ / = 0; the turbulent energy is minimal in the isotropic case when the three eigenvalues λ 2 α are equal λ 2 1 = λ2 2 = λ2 3 = λ2. In this case, the orthonormal eigenvectors e α, α {1, 2, 3} of the Reynolds stress tensor P are also the orthonormal eigenvectors of the mean rate of deformation tensor D. The eigenvalues of D are µ α, α {1, 2, 3}. In the compressible isotropic case e T = ρ 2/3 κ, κ = 3λ 2 /2 ρ 2/3, and κ is a classical invariant of isotropic turbulence. In presence of shock waves the quantity ρ 2 3 λ 2 α is not conserved through shocks; it increases like the classical entropy in compressible fluid dynamics. The estimation of the jump of turbulence entropy in isotropic case was given in [9]. The governing equations 2 admit the energy conservation law t ρ 1 2 U 2 + e i + e T + div ρ U 1 2 U 2 + e i where the internal specific energy e i is defined by 1 de i = p d ρ + p I + RU = 0, and the mean pressure p is supposed to be a given function of ρ. Indeed, using 17 2 we immediately obtain 1 1 ρ t 2 U 2 + e i + e T + div ρ U 2 U 2 + e i + p I + RU = ρ 2 3 Conclusion ρ de T +trrd = d λ 2 α + ρ tr λ 2 αµ α = 0. The equations of fluid turbulent motions take three equations into account: the equation of the mass balance 1 1, the balance equation of the average momentum 1 2, and the Reynolds stress tensor equation of evolution 1 3 ; this last 7

8 equation has been the object of our study. If the turbulent sources are neglected, the turbulent fluid motion is a superposition of the mean motion and turbulent fluctuations. The eigenvectors of the Reynolds tensor carry the fluctuations associated with the mean flow deformation. The amplitude of turbulent deformations is defined by the eigenvalues of the Reynolds stress tensor. Locally, the equations for the directions of turbulent fluctuations describe a small solid whose rotation is given by a gyroscopic type equation Equation The amplitude evolution of turbulent deformations is determined by the diagonal values µ α of the mean deformation tensor D expressed in the eigenvector basis of the Reynolds stress tensor. The turbulence increases with the time when µ α < 0, and decreases when µ α > 0. In the particular case of incompressible fluid motions we have tr D = 0, and hence there always exists a direction in which the turbulence is increasing while in other directions it is decreasing. These mathematical deductions are confirmed by experiments [10]. Acknowledgement: This article is dedicated to the memory of Professor Vladimir M. Teshukov who drew our attention to his model of weakly sheared flows during his visit in Marseille in the fall of year References [1] S. B. Pope, Turbulent flows, Cambridge University Press, [2] D. Wilcox, Turbulence Modeling for CFD, DCW Industries, [3] B. Mohammadi, O. Pironneau, Analysis of the K-epsilon turbulence model, Research in Applied Mathematics, John Wiley & Sons, New York, [4] C. Berthon, F. Coquel, J.M. Hérard, M. Uhlmann, An approximate solution of the Riemann problem for a realisable second-moment turbulent closure, Shock Waves, 11, [5] V.M. Teshukov, Gas dynamic analogy for vortex free-boundary flows, Journal of Applied Mechanics and Technical Physics, 48, [6] J.-F. Debiève, H. Gouin, J. Gaviglio, Evolution of the Reynolds stress tensor in a shock wave-turbulence interaction, Indian Journal of Technology 20 March [7] J. Serrin, Mathematical principles of classical fluid mechanics, in: Encyclopedia of Physics VIII/1, Springer, Berlin, 1959, pp [8] J. E. Marsden, T.S. Ratiu, Introduction to mechanics and symmetry, Series in Applied Mathematics, 17, Springer, Berlin, [9] S. Gavrilyuk, R. Saurel, Estimation of the turbulence energy production across a shock wave, The Journal of Fluid Mechanics, 549,

9 [10] J.-F. Debiève, H. Gouin, J. Gaviglio, Momemtum and temperature fluxes in a shock wave-turbulence interaction, in: Proceedings of the ICHMT/IUTAM Symposium on the Structure of Turbulence and Heat and Mass Transfer, Z. Zarić, Ed., p , Hemisphere Publishing Corporation, London