CC BY
78
УДК 532.526.4
ТЕОРИЯ И КОНСТАНТЫ ПРИСТЕННОЙ ТУРБУЛЕНТНОСТИ
Трунев Александр Петрович к. ф.-м. н., Ph.D., директор
A&E Trounev IT Consulting, Торонто, Канада
В работе представлена полностью замкнутая модель турбулентного пограничного слоя, полученная из уравнения Навье-Стокса. Фундаментальные константы пристенной турбулентности, включая постоянную Кармана, определены из теории. Эта модель была развита для ускоренного и неизотермического пограничного слоя. Профили средней скорости и температуры, вычисленные на основе модели находятся в согласии с экспериментальными данными
Ключевые слова: ТУРБУЛЕНТНЫЙ ПОГРАНИЧНЫЙ СЛОЙ, ЛОГАРИФМИЧЕСКИЙ ПРОФИЛЬ
UDC 532.526.4
THEORY AND CONSTANTS OF WALL TURBULENCE
Alexander Trunev Ph.D., Director
A&E Trounev IT Consulting, Toronto, Canada
The completely closed model of wall turbulence was derived directly from the Navier-Stokes equation. The fundamental constants of wall turbulence including the Karman constant have been calculated within a theory. This model has been developed also for the accelerated and non-isothermal turbulent boundary layer flows.
The estimated mean velocity, temperature and impurity concentration profiles as well as the spectral characteristics of the streamwise velocity component are to be shown in a good agreement with the experimental data
Keywords: TURBULENT BOUNDARY LAYER, LOGARITHMIC PROFILE
1. Introduction
Nikuradse [1] found that the turbulent boundary layer structure includes the region where the mean velocity increases with the distance from the wall as a logarithmic function
U 1 , u z ..s
— = -ln-^ + ^ (1)
uT - n
where ut is the friction velocity, ut =^Jt I p, t is the wall shear stress, p is the fluid density, z is the distance from the wall, n is the kinematic viscosity. Nikuradse [1] established that - = 0.4, c0 = 5.5 for turbulent flows in smooth pipes.
The logarithmic law (1) has been confirmed by many and many researchers of turbulence [2-7, 13-14,17]. Recently it was shown [5] that the log law can be
excellent agreed with the experimental pip flow data for
600 < zut In < 0.07Rut /n, where R is the pip radius. For this range the log law
constants were shown to be - = 0.436, c0 = 6.15. Also the log law was confirmed for a zero pressure gradient turbulent boundary layer flow in the precise wind-tunnel experiment [6] and the constants were obtained - = 0.38, c0 = 4.1. To estimate the Karman constant the normalised local slope of the mean velocity profile has been used as follows X = dln z +1dU + where z + = zut In is the inner layer variable - see Figure 1a.
So this universal law can be considered as a fundamental fact of the theory of turbulence. Nevertheless it may be wonder that this law never been derived from "first principles" excepting the similarity theory (see the typical arguments in ref. [2-8] and other).
To derive this law we can begin with simple assumption that the logarithmic profile of velocity is a fragment in solution of hydrodynamics equations of viscous fluid in a form of Navier-Stokes. Thus the way a question is put can be different: how the initial equations should be transformed so that in their solution boundary layer should contain logarithmic profile? This type of transformation was developed recently [9-11]. The idea of transformation is rather simple and as a matter, was driven at bringing into initial system equations of some set of random parameters, describing fluctuations of conventional boundary of viscous sublayer.
To understand this idea we should note that the logarithmic distribution (1) is a typical mean velocity profile measured in the turbulent flow bounded by a smooth wall. In a case of the turbulent flow over a rough wall Nikuradse [12] established (for sand-roughened pipes) that if the roughness height significantly exceeds the viscous sublayer thickness, then the mean velocity profile can be described by the logarithmic function as follows:
U 1 z
— = —ln— + cs (2)
ut - ks
where ks is the characteristic scale of the sand roughness, -, cs are empirical values (Nikuradse found that - = 0.4, cs = 8.5 for the completely rough regime).
The sentence "the roughness height significantly exceeds the viscous sublayer thickness" means only that ks >> n Iut or in another form k+ = ksux In >> 1. Obviously that k+ is the roughness Reynolds number and n I ut is the characteristic thickness of viscous sublayer. Compared equations (1) and (2) one can conclude that there is an universal model described the turbulent boundary layer flows over smooth and rough surfaces. As it was shown [9-11] this model can be derived directly from the Navier-Stokes equation and that the logarithmic distribution is an asymptotic solution of this model.
2 Principles of theory
Let r = r(x, y)is a function describing the roughness geometry and = h~(x, y, t)is the thickness of viscous sublayer. The dynamic roughness surface can be defined as follows h( x, y, t) = r (x, y)+~(x, y, t). Thus the mean dynamic roughness height is given by h = ks + 1+n Iut, where 1+ is the dimensionless thickness of viscous sublayer. It should be noted that in the case of a smooth plate for
r (x, y) = 0 the dynamic roughness surface equals to the viscous sublyer conventional surface i.e. h( x, y, t) = ~(x, y, t).
In Figure 1b the frequency diagram [13] of instantaneous thickness of a viscous sublayer in the turbulent boundary layer over a smooth surface is shown. The maximum of frequency diagram is established for z + = z + » 9. The first maximum of the inverse normalised local slope of the mean velocity profile A = dU + / d ln z + also locates approximately at z + » 9- see Figure 1a. Thus one can put 1+ » 9 .
For the flow over a smooth surface we have h = 1+v / ut and in this case the logarithmic velocity distribution in the form (1) is realised. When ks >> 1+v / uT and therefore h » ks then the mean velocity profile can be described by eq. (2). This two cases can be described automatically if the flow velocity vector, u = (u, v, w), can be written as follows
u = u( x, y, z/h( x, y. t), t) (3)
A a) b)
Figure 1: a) The inverse normalised local slope of the mean velocity profile A = dU + / d ln z + estimated from the van Driest model [14]- the solid line 6 and from the experimental data by Osterlund [6] (1-4). Horizontal solid line (5) corresponds to the estimated quantity A = 2.54; b) the frequency diagram of instantaneous thickness of viscous sublayer in the turbulent boundary layer over a smooth wall [13].
If the streamwise velocity has a logarithmic asymptotic, u ^ ln(z / h), then also the mean velocity has the logarithmic asymptotic, U = (u) ^ ln(z / h).
Unfortunately this very simple idea can’t be used directly by substituting expression (3) in the Navier-Stokes equation, because the dynamic roughness surface can’t be described by a regular function in a common case.
Therefore the dynamic roughness surface is advised to characterise by the set of the random continuous parameters h, ht, hx, hy , with a well-known function of
distribution fs = fs (h, hx,hy,ht). The reason why only the first odder derivatives of the function h = h( x, y, t) play role in this theory is explained by the volume averaging used for a filtration of random functions.
We suppose that h = z / h = const and let us consider a representative region of flow by a volume dV = LxLydz, where Lx, Ly are typical scales of flow in the x, y directions accordingly. Then let us consider subregion dVs lying in the representative region of the flow, dV, in which the random parameters h, ht, hx, hy are changed in the intervals (h; h + dh) , (ht; ht + dht), (hx; hx + dhx) , (hy; hy + dhy) accordingly. The volume of this subregion is given by
dVs = dVfs (h, hx, hy, ht )dhdhxdhydht.
The random function describing the flow velocity can be determined by the toting expression u = u( x, y, z/h( x, y. t), t) in the volume dVs:
~ 1 f
u(h, t, h,hx,h , ht) = lim — I u(x,y, h, t)dxdydz
7 SV®dVs OV
SV
where SV is an arbitrary volume enclosed in dV = LxLydz and containing dVs as a whole.
Obviously, u(h, t, h, hx, hy, ht) is the random function, because it depends on the
random parameters. The equations describing dynamics of random functions u = u(h, t, h, hx, hy, ht) immediately follow from the hydrodynamic equations of
viscous fluid [9-11]. Statistical moment of an order m of the random function u(h, t, h, hx, hy, ht) are determined as follows
~ m (z, t) = I ~ m (h, t, K K , hy , ht )fs (h K , hy , ht )dhdKdhydht .
This algorithm is a very effective in the wall turbulence problem because in this case the random functions ~(h, t,h,hx,hy, ht)with fixed variables h,hx ,hy, ht have
the properties which are similar to the mean flow parameters. Otherwise for the solutions presented for instance by the logarithmic function we can suppose that
~ = | ~(h, t, h, hx, hy, ht) fs (h, hx, hy, ht )dhdhxdhydht = u( z / h*, h*, h*y, h*),
where the parameters with stars can be estimated from the comparison of solutions with experimental data or calculated from a theory considered below.
A zero pressure gradient turbulent boundary layer over a smooth surface can be described by the completely closed equation system derived directly from the
Navier-Stokes equation (see [10-11] for details). Utilised the inner layer variables this model can be written in the form
+
dz
+
f 2 • 2 Л
cos a sin a
(4)
Rc dll + d(1 + x2)d2| = 0 , I = R, \Щ
tA dX2 dX dX2 J01 + x2
where u + = u/ ux, X = z + /1+, A = huT /ViJ h2x + h2y is the characteristic dimen-sionless scale of the viscous sublyer, a = arctan(hy / hx), v += v/ ux, w+ = w/uT, ^ = [(u + cosa + v + sin a + w0+)X - w+ ]/w+, w+ = ht / u*^h2x + hy; is the second scale of velocity in the wall turbulence, Rt = hht/v(hx2 + h2y) is the Reynolds
number calculated on the dynamical roughness parameters (it should be noticed that the random parameters Rt, w0+ ,a, A+ are not functions of time or space variables!).
The boundary conditions for Eq. (4) are set as follows
u+(0) = 0, ^ (0) = 0, ^' (0) = 1, ^ "(0) = a (5)
where a is a free (shooting) parameter.
It should be noted that model (4) is related to the Blasius equation in the particular case when R = 1/2, a = 0 and X << 1.
3 Constants of wall turbulence
Let us suppose that for z + ® ¥ the streamwise mean velocity profile has a logarithmic asymptotic, i.e.
du + / dz + ® 1/ k z + .
Surmising that limx®¥I(X,R) = I0(R) we have from the first equation (4):
du + / dz + = AW z +, and therefore 1+ = eI 0/k. The last equation gives the continuous spectrum of the turbulent scales A = A+(R). Used the dynamic roughness Reynolds number in the form Rt = w+A+ the equation A+ = e10 / k can be rewritten as follows
w0+ / k = Rt exp[-I0(Rt)]
For an uniqueness of the mean velocity profile we can suppose that for k = const the second turbulent velocity scale has a stable value at small variations of the parameter R, i.e. Sw+ /OR = 0. It gives Rt = R* »1.2193 . The turbulent boundary
layer fundamental scale can be defined as A+ = eh{Rt) / k . This quantity approximately equals to the experimental value A+ » 9.16 for k = 0.39.
The function 1/ 1+ = ke-10 can be considered as a spectral density. The inverse length scale versus the Reynolds number is shown in Figure 2a. This type of a spectral density is similar to the spectral density of the streamwise velocity fluctuations in the turbulent boundary layers. The function w+ = w0+ (R) is shown in Figure 2b. This type of spectrum is similar to the spectral density of the transversal velocity pulsation [15]. Both functions represent the constructive model of the hydrodynamic chaos in this theory of turbulence.
Figure 2: a) The inverse length scale 1/A+ versus the Reynolds number of dynamic roughness in double logarithmic scale. This type of a spectral density is similar to the spectral density of the streamwise velocity fluctuations in the turbulent boundary layers; b) the normalised turbulent velocity scale w+ /k versus the Reynolds number of dynamic roughness. This type of spectrum is similar to the spectral density of the normal to the wall velocity pulsation.
Our suggestion about the dynamical roughness structure in the wall region is that the parameter a fluctuates around the mean value a = p /2. This structure looks like furrows elongated along of the mean flow stream lines in the viscous sublayer (see, for instance, [16] where the visualisation of the coherent structure in the turbulent boundary layer is presented). In this case the first Eq. (4) can be integrated once and the standard logarithmic profile can be derived for
z + >> 10 +:
+ 1 , + 1 ¥ el0-I -1 x 1 , 1 +
u =— ln z + c0, c0 = — I —j= dX -— ln-------------------------.
k °’ 0 k{fi+X2 h k 2
Therefore, with the given constant k another constant of the mean velocity logarithmic profile can be calculated from the last equation.
4 Spectral characteristics of wall turbulence
The spectral density of the streamwise velocity pulsation can be defined as follows
1 T
IF(k)dk =|u'2 (w)dw = lim t Ju'2 (t)dt =
u12
where w is the characteristic radian frequency, k = w / U is the flow wave number (Taylor's frozen turbulence hypothesis). To compare the spectral density with experimental data let us suppose that ht ^ wh, Therefore the Reynolds number calculated on the dynamical roughness parameters depends on the frequency,
Rt = ah2 /vn2 = wl+2v/u2 = kHl+ 2u + /Re*
where n2 = h2 + h2, 1+ is the typical turbulent length scale of the streamwise velocity pulsation, H is the boundary layer height, Re* = HuT /v. Due to this relationship the spectral characteristic of the turbulent flow is related to the eigen spectrum of the value problem (4-5). In general case it can be like the power series
F (k) = X cl)"'
Practically we can test one first term of this series. Then the spectral density can be proposed in the form
F (k) = = *c,Hu]e -'0[ Rt (k)] (6)
г =1
Figure 3:The estimated spectral density of the streamwise velocity pulsation in the turbulent channel flow (solid lines) for z + = 4.9 (a), and z + = 11.7(b), and the experimental data [17]
where cK is the normalising factor. Suggesting that the flow wave number depends on the dynamical roughness Reynolds number as the linear function, i.e.
k ^ Rt, we have cK = 2.131+ 2u+(u'2)/ut Re*.
The spectral density given by Eq. 6 is shown in Figure 3 (solid lines) together with the experimental data [17] obtained in the turbulent channel flow. As it was established both spectral density parameters slowly depend on the distance from the wall. The best correlation was found in the viscose sublayer for z += 4.9 1+ = 6.76 cK = 1.01(Figure 3a) and for z += 11.7,1+ = 5.45, cK = 1.82 (Figure 3b).
5 Thermal turbulent boundary layer
The model of turbulent flow (4) was developed to estimate the turbulent transport of heat and impurities in the boundary layer. The mean temperature gradient can be written as follows (see [10-11] for details)
dT+ = Pr exp[- °-5R C '(x )ln(1 + 11 )1 (7)
dz' V1+x .2
where T + = (Tg - ~)/ T*, Tg is the surface temperature, T is the turbulent flow temperature, T* = qH / (p cpuT) is the turbulent scale of temperature, qH is the heat flux from the rigid surface to the flow, cp is the specific heat at constant pressure of the gas, X i= z + / 1T+, 1T + = (kh Pr)-1 is the main turbulent scale of the ther-
mal layer, Pt = w°+1+ Pr is the Peclet number.
The boundary conditions for equation (12) on a smooth wall and for the long distance from the wall are given by
z + = 0: T + = 0; dT+ / dz + = Pr (8)
z + ® ¥ dT+ / dz + ® 1 / k hz + , where kh is a constant which approximately equals to the Karman constant.
iji+ iji+
Figure 4. Mean temperature profiles in the turbulent boundary layer computed for the small (left) and for the large Prandtl numbers (right). The square symbols -the model of Sebeci [18]; the solid lines - the model (7-9).
In special case when Pr ® 0 a general solution of the boundary value problem (7)-(8) can be written as follows:
T+= -^Arsh(k hPr z+) (9)
k h
The mean temperature profiles computed on equations (7), (9) are shown in Figure 4 (right) by solid lines together with the profiles computed on the model of Sebeci [18] - the square symbols. For the best correlation the function Pt(Pr)was numerically approximated as follows Pt =-161/(1 + 0.1lnPr) in the range of the Prandtl number 2 < Pr < 103. The agreement between two models in general is good.
It should be noted that equation (7) is a typical two-scale model describing the interaction between the transition layer with l»10 and thermal boundary layer
with the main scale 1T + = (kh Pr)-1. Considering the mean flow as a "substance" with the Prandtl number Pr = 1 we can utilized equation (7) also to describe the mean velocity profile.
6 Multi-scale model
Analysing data shown in Figure 1a one can found the several layers of different scale:
• The viscous syblayer z + < 1/ k , A » z +;
• The transition layer which corresponds to the first maximum of A for
z + » 1+;
• The logarithmic layer where A »1/ k ;
• The mixing layer or the "wake region" [19] which corresponds to the second maximum of A for z » H/2;
• The free stream region where A ® 0 for z > H.
With this scales the model of turbulence (4) has been developed for the case of turbulent boundary layer flow in pressure gradient as follows
du + = £0exp(I0 -1) _ A+p+e-- | eJd%____________£0 + £(P+ ^1 + zo where
dz + k1+^1 + X2 k-J 1 + X2 0^1 + X2 kg*Re*V1 + z2 Kg*Re*(1 + z2)
£ 0 = 1 + (1 - £ )/(kg*Re* Jl + zf), g* = 0.27, p + = -n- — is the dimensionless pres-
p u* dx
sure gradient, z = (z - z0)/ g*H, z0 = H /2 is the middle position of the mixing layer.
Figure 5: The estimated mean velocity profiles in the turbulent boundary layer in adverse pressure gradients computed for p+=0.023; 0.0251 - the solid lines 1-2 accordingly, and experimental data (symbols) from the data base by Nagano et al [20].
The function of pressure gradient parameter has been estimated from the experimental data by Nagano et al [20]
3 2
£(p+) @ 118099p + - 2386.8p + + 36.579p ++ 0.7881 for adverse pressure gradients. The computed mean velocity profiles in the turbulent boundary layer in adverse pressure gradients are shown in Figure 5 together with the experimental data [20].
7 Final closure and conclusion
Generally speaking the Karman constant also should be defined as a function of this parameter because the two-scale model (7) depends on the length scale 1+ = 1/ k which can be related to the fundamental length scale 1+. Used a symmetry of the multiple density of a probability distribution function fS (h, hx, hy ,-ht) = fS (h, hx, hy, ht) one can defined the length scale corresponding to
a negative value Rt = -R* as follows 1+n = 1+ exp[/0(-R*)] and consequently the scale of velocity w+ = -R*/ 1+n. Both scales depend on the main scales w0+, 1+ due to a symmetry so that 1+n = w+ ,w+ =-1++. Therefore the Karman constant should be defined as follows
1
k = ,-------exp
10(R*) +10(-R*) 2
(10)
The formula (10) gives k = 0.3931 for R* = 1.2193. This result is presented in Figure 1a by the horizontal line 5.
We have demonstrated the completely closed theory of wall turbulence. The main parameter of this theory was established as the Reynolds number of dynamic roughness Rt = w+1+ = hht /(h2x + h2y) . This parameter has been calculated
within a theory based on the mean velocity profile asymptotic behaviour as R* = 1.2193. The Karman constant and the turbulent boundary layer fundamental length scale have been defined as functions of this parameter. The estimated parameters k = 0.3931, 1+ » 9.16 and mean velocity and temperature profiles are in a good agreement with the available experimental data.
References
1. Nikuradze, J. Gesetzmaessigkkeiten der turbulenten Stroemung in glatten Rohren,
Forsch. Arb. /ng. Wes., 356, 1932.
2. Schlichting, H. Boundary Layer Theory, McGraw-Hill, NY, 1979.
3. Purtell, L. P., Klebanoff, P. S. & Buckley, F. T. Turbulent Boundary Layers at Low Reynolds Number. Phys. Fluids, 24(5), 802-811, 1981.
4. Smith R. W. Effect of Reynolds Number on the Structure of Turbulent Boundary Layers, Ph.D. Thesis, Princeton University, 1994.
5. Zagarola, M. V. & Smits, A. J. A New Mean Velocity Scaling For Turbulent Boundary Layers, Proceedings of FEDSM’98, 1998 ASME Fluids Engineering Division Summer Meeting, June 21-25, 1998, Washington DC, 1998.
6. Osterlund, J.M. Experimental studies of zero pressure gradient turbulent boundary layer flow. Doctoral thesis. Stockholm, 1999.
7. Fernholz, H. H. & Finley, P .J. The incompressible zero-pressure-gradient turbulent boundary layer: An assessment of the data. Prog. Aerospace Sci. 32 , pp. 245 -311,1996.
8. Millikan, C. M. A Critical Discussion of Turbulent Flows in Channels and Circular Tubes, Proc. 5th /nt. Congr. Appl. Mech. (Wiley, NY) 386-392,1938.
9. Trunev A. P. Similarity theory for turbulent flow over natural rough surface in pressure and temperature gradients. Air Pollution /V. Monitoring, Simulation and Control ed. B. Caus-sade, H. Power and C. A. Brebbia (Comp. Mech. Pub., Southampton) 275-286. 1996.
10. Trunev A. P. Theory of Turbulence and Model of Turbulent Transport in the Atmospheric Surface Layer, Russian Academy of Sciences, Sochi, 160 p., 1999.
11. Trunev A. P. Theory of Turbulence and Turbulent Transport in the Atmosphere. 180 p., 2001.
12. Nikuradse J. Stromungsgesetze in Rauhen Rohren ForschHft. Ver. Dt. /ng., 361, 1933.
13. Kutateladze S. S. Wall Turbulence (Nauka, Novosibirsk) p.102, 1973.
14. Van Driest, E. R. On turbulent flow near a wall, J. Aero. Sci. 23,1007,1956.
15. Tennekes H. & Lumley J. L. A First Course in Turbulence, MIT Press, Cambridge, Massachusetts, 1972.
16. Cantwell B. J., Coles D. E. & Dimotakis P. E. Structure and entrainment in the plane of symmetry of a turbulent spot J. Fluid Mech, 87, 641-672, 1978
17. Hussain A. K. M. F. & Reynolds W. C. Measurements in fully development turbulent channel flow, J. Fluid /ng. 97, 568-80, 1975.
18. Cebeci T. A model for eddy conductivity and turbulent Prandtl number, J. Heat Transfer,
95, 227, 1973.
19. Coles, D. E. The law of the wake in the turbulent boundary layer. J.Fluid Mech.1 , pp.191 -226, 1956.
20. Nagano, Y., Kasagi, N., Ota, T., Fujita, H., Yoshida, H. & Kumada, M., Data-Base on Turbulent Heat Transfer, Department of Mechanical Engineering, Nagoya Institute of Technology, Nagoya, DATA No. FW BL004, 1992.
