Exact statistical theory ofisotropic turbulence1Ran zhengShanghai Institute of Applied Mathematics and Mechanics,Shanghai University, Shanghai 200072,P.R.ChinaAbstractSome physically significant consequences of recent advances in the theory ofself-preserved homogenous statistical solutions of the Navier-Stokes equationsare presented.Keywords: isotropic turbulence, Karman-Howarth equation, exact solution1 IntroductionHomogeneous isotropic turbulence is an idealized concept of turbulence assumed to be governed by a statistical law that is invariant under arbitrary translation (homogeneity), rotation or reflection (isotropic) of the coordinate system. This is an idealization of real turbulent motions, which are observed in nature or produced in a laboratory have much more complicated structures. This idealization was first introduced by Taylor (1935) to the theory of turbulence and used to reduce the formidable complexity of statistical expression of turbulence and thus make the subject feasible for theoretical treatment. Up to the present, a large amount of theoretical work has been devoted to this rather restricted kind of turbulence. Remarkable progress has been achieved so far in discovering the various nature of turbulence, but nevertheless our understanding of the fundamental mechanics of turbulence is still partial and unsatisfactory (Tatsumi, 1980)The assumption of self-preservation permits an analytical determination of the energy decay in isotropic turbulence. This idea of similarity and self-preservation has played an important the development of turbulence theory for more than a half-century. The traditional approach to search for similarity solutions in turbulence has been assume the existence of a single length and velocity scale, then ask whether and under what conditions the turbulent motions admit to such solutions. Excellent contributions had been given to this direction. Von Karman and Howarth (1938) firstly deduced the basic equation, also presented a particular set of solutions of Karman-Howarth equation for the final decaying turbulence. Later on, two Russian authors, Loitsiansky (1939) and Millionshtchikov (1941), have each discussed the solutions of Karman-Howarth equation which are obtained when the term describing the effect of the triple 1 National Natural Science Foundation of China (10272018, 10572083) velocity correlation is ignored. Their work is an extension of the “small Reynolds number” solution first put forward by von Karman and Howarth.H.Dryden (1943) gave a useful review on this subject. Great deep research on the solutions of Karman-Howarth equation was conducted by Sedov (1944). He showed that one could use the separability constraint to obtain the analytical solution of Karman-Howarth equation. Sedov’s solution could be expressed in terms of the confluent hypergeometric function. Batchelor (1948), who was also the first to carry out a detailed study of this problem, based his analysis on the assumption that the Loitsianskii integral is a dynamic invariant- a widely accepted hypothesis that was later discovered to be invalid. Batchelor conclude that the only complete self-preserving solution that was internally consistent existed at low turbulence Reynolds numbers where the turbulent kinetic energy consistent with the final period of turbulent decay. Batchelor (1948) also found a self-preserving solution to the Karman-Howarth equation in the limit of infinite Reynolds numbers for which Loitsiansky’s integral was an invariant. Objections were later raised to the use of Loitsiansky’s integral as a dynamic invariant: at high Reynolds numbers this integral can be shown to be a weak function of time (see Proudman & Reid 1954 and Batchelor & Proudman 1956). Saffman (1967) proposed an alternative dynamic invariant which yields another power law decayin the infinite Reynolds numbers limit (see Hinze 1975). While the results of Batchelor and Saffman formally constitute complete self-preserving solutions to the inviscid Karman-Howarth equation, it must be kept in mind that they only exhibit partial self-preservation with respect to the full viscous equation. Recently, George (1992) revived this issue concerning the existence of complete self-preserving solutions in isotropic turbulence. In an interesting paper he claimed to find a complete self-preserving solution, valid for all Reynolds numbers. George’s (1992) analysis is based on the dynamic equation for the energy spectrum rather than on the Karman-Howarth equation. Strictly speaking, the solution presented by George was an anlternative self-preserving solution to that of Karman-Howarth and Batchelor since he relaxed the constraint that the triple longitudinal velocity correlation be self-similarity in the classical sense. Speziale and Bernard (1992) reexamined this issue form a basic theoretical and computational standpoint. Several interesting conclusions had been drawn from their analysis.From the long history stories, we know that: Decaying homogeneous and isotropic turbulence is one of the most important and extensively explored problems in fluid dynamics and a central piece in the study of turbulence. Despite all the effort, a general theory describing the decay of turbulence based on the first principles has not yet been developed. (Skrbek and Steven R.Stalp, 2000). It thus appears that the theory of self-preservation in homogeneous turbulence has many interesting features that have not yet been fully understood and are worth of further study. (see Speziale and Bernard ,1992,pp.665).This paper offers a unified investigation of isotropic turbulence, based on the exact solutions of Karman-Howarth equation. Firstly, we will point out that new complete solution set may be exist if we adopt the Sedov (1944) method. Based on the above finds, new results could be obtained for the analysis on turbulencefeatures, such as the scaling behaviour, the spectrum, and also the large scale dynamics.2 Self-preservation solution under Sedov’s separability constraintFor complete self-preserving isotropic turbulence, the Karman-Howarth equation will have a solution if Reynolds number based on the Taylor microscale is constant as first noticed by Dryden (1943). However, this equation also has solutions where Reynolds number based on the Taylor microscale is time dependents when separability is invoked. The separability condition implies that each side of the equation is equal to zero individually, yielding differential equations from which explicit solution for the correlation functions may be determined depending on the choice of parameters. These solutions were first discovered by Sedov (1944) and later compared with experimental data by Korneyev & Sedov (1976). Here, we will discuss the possible new complete solutions under Sedov’s separability constraint.The two-point double longitudinal velocity correlations read as ( named Sedov equation)02242122=+⎟⎟⎠⎞⎜⎜⎝⎛++f a d df a d f d ξξξξ (1) with boundary conditions ()10=f ,()0=∞f .In the following analysis, we introduce alternative two parameters denoted by σ,1a ， here122a a =σ (2) The complete new set of the solution of the equation (1) with the boundary condition could be given as following ：The first kind of solution ：if 25=σ，thus ()214ξξa e f −= (3)The secondary kind of solution ：if 45−=σκ，thus ()⎟⎠⎞⎜⎝⎛−=−2144,25,2521ξσξξa F e f a (4) The third kind of solution ：if σκ−=45，thus ()⎟⎠⎞⎜⎝⎛=−2144,25,21ξσξξa F e f a (5) The forth kind of solution ：if 45=σ，thus ()⎟⎠⎞⎜⎝⎛=−2144,23,4321ξξξa F e f a (6) where (z F ,,)γαis the confluent hypergeometric function and the definition of the existing parameter κwill be given appendix. From the asymptotic expansions and the limiting forms of the confluent hypergeometric function, we could deduce the existence conditions of these solutions ：For all four kind of solutions ：,01>a For the secondary kind of soultions ：0>σ；For the third kind of solution ：250<<σ； The details could be seen in appendix .A simple comparison shows that the special solution found by Sedov (1944) is belonging to the secondary kind of our new set of solution.3 Some theoretical results based on the exact solutionsA unified investigation of isotropic turbulence, based on the above exact solutions of Karman-Howarth equation could be given. New results could be obtained for the analysis on turbulence features, such as the scaling behaviour, the spectrum, and also the large scale dynamics, some results could be seen in the following references [19,20,21,22] .References[1]Bareenblatt,G.J.& Garilov,A.A. 1974. Sov. Phys. J. Exp. Theor. Phys. 38,399-402.[2]Batchelor,G.K. 1948.Q. Appl. Maths. 6,97-116.[3]Batchelor, G.K. 1953 The Theory of Homogeneous Turbulence Turbulence. Cambrige University Press.[4]Batchelor,G.K. & Proudman,I. 1956. Phil. Trans. R. Soc. Lond. A.248,369-405. [5]Dryden,J.L. 1943. Q. Appl. Maths. 1,7-42.[6]George,W.K. 1992. Phys. Fluids A 4,1492-1509.[6]Hinze,J.O. 1975 Turbulence. McGraw-Hill.[7]Karman,T.Von & Howarth,L. 1938. Proc. R. Soc. Lond. A164,192-215.[8]Korneyev & Sedov, L.I. 1976. Fluid Mechanics-Soviet Research 5,37-48.[9]Lesieur,M. 1990 Turbulence in Fluids, 2nd Edn. Martinus Nijhoff.[10]Lin,C.C. 1948. Proc. Natl. Acad. Sci. 34,540-543.[11]Millionshtchikov, M. 1941. Dokl. Akad. Nauk SSSR 32,615-618.[12]Monin,A.S. & Yaglom,A.M. 1975 Statistical Fluid Mechanics: Mechanics of Turbulence, vol.2, MIT Press.[13]Proundman, I. & Reid, W. H. 1954. Philos. Trans. R. Soc. London, A 247,163-189.[14]Saffman,P.G. 1967. J. Fluid Mech. 27,581-594.[15]Sedov, L.I. 1944. Dokl.Akad.Nauk SSSR 42,116-119.[16]Speziale,C.G. & Bernard, P.S. 1992. J. Fluid Mech. 241,645-667.[17]Skbek,L. & Steven, R.S. 2000. Phys. Fluids 12,1997-2019.[18]Tatsumi, T. 1980. Advances in Applied Mechanics,39-133.[19]Ran zheng, Exact solutions of Karman-Howarth equation. -200510-311[20]Ran zheng, Scales and their interaction in isotropic turbulence. -200510-205.[21]Ran zheng, Dynamic of large scales in isotropic turbulence. -200509-230.[22]Ran zheng, On von Karman’s decaying turbulence theory. -200511-233.Appendix: Solutions of the correlation coefficientsThe ideas of similarity and self-preservation were firstly introduced by von Karman (1938). Following the methods adopted by Sedov (1944,1951), the two point double longitudinal velocity correlations satisfied02242122=+⎟⎟⎠⎞⎜⎜⎝⎛++f a d df a d f d ξξξξ （1） with the boundary condition()10=f()0=∞fThe complete solution are given in this paper, these are ：When 25=σ，()214ξξae f −= When 45−=σκ，()⎟⎠⎞⎜⎝⎛−=−2144,25,2521ξσξξa F e f a When σκ−=45，()⎟⎠⎞⎜⎝⎛=−2144,25,21ξσξξa F e f a When 45=σ，()⎟⎠⎞⎜⎝⎛=−2144,23,4321ξξξa F e f a The detailed calculation is given as following:A lot of useful partial differential equations can be reduced to confluent hypergeometric equations. Let ()ςm k P ,is the solution of Whittaker equation as that defined by Whittaker and Waston041412222=⎥⎦⎤⎢⎣⎡−++−+W m k d W d ςςς （2） where()()()()z h P e z z y m z f ,κβ= （3）After some reduction, the equation of ()z y reads()()02222=⋅+⎥⎦⎤⎢⎣⎡′++′′′−z y g dzdy z f z h h dz y d β （4） where()()12212g f z h h z z f f f g +⎟⎠⎞⎜⎝⎛′+′′′+++′+′′−′=ββββ ⎟⎟⎠⎞⎜⎜⎝⎛−+−⎟⎠⎞⎜⎝⎛′=4412221h h m h h g κ The solutions of above equation could be deduced in terms of Whittaker function. We discussed this equation in following special case:()λaz z f = （5）()λAz z h = （6） The equation under this condition reads as()0221122=+⎥⎦⎤⎢⎣⎡−−−+−z qy dzdy z z dz y d λλαβλ （7） where()()()2222222224124z m z A z A q −+++++⎟⎟⎠⎞⎜⎜⎝⎛−=−−λλββκλαβλαλλλThe solution of this equation is()()λκαβλAz P e z z y m z ,= （8）For isotropic turbulence, the corresponding parameters satisfied421=−−βλ （9） 11=−λ （10）221a =−λα （11） 04222=⎟⎟⎠⎞⎜⎜⎝⎛−A αλ （12） ()04122=⎟⎠⎞⎜⎝⎛−++m λλββ （13） ()222a A =+κλαβλ （14） Hence, we have2=λ （15）81a −=α （16） 25−=β （17） 43±=m （18） 41a A ±= （19） ⎭⎬⎫⎩⎨⎧−±=45212a a κ （20） From above analysis, we can introduce two parameters to classification turbulence, they are: 1212,a a a =σ。