CC BY
62
УДК 532.526.4
ТЕОРИЯ ТУРБУЛЕНТНОСТИ И МОДЕЛИРОВАНИЕ ТУРБУЛЕНТНОГО ПЕРЕНОСА В АТМОСФЕРЕ ЧАСТИ 1, 2
Трунев Александр Петрович к. ф.-м. н., Ph.D.
Директор, A&E Trounev IT Consulting, Торонто,
Канада
В работе представлена полностью замкнутая модель турбулентного пограничного слоя, полученная из уравнения Навье-Стокса. Фундаментальные константы пристенной турбулентности, включая постоянную Кармана, определены из теории. Эта модель была развита для ускоренного и неизотермического пограничного слоя над шероховатой поверхностью. Исследованы численные решения системы уравнений турбулентного преноса примесей в приземном слое атмосферы для больших масштабов
Ключевые слова: АТМОСФЕРНАЯ ТУРБУЛЕНТНОСТЬ, ТУРБУЛЕНТНЫЙ ПЕРЕНОС, УСКОРЕННЫЕ ТЕЧЕНИЯ,
ПОГРАНИЧНЫЙ СЛОЙ, ШЕРОХОВАТАЯ ПОВЕРХНОСТЬ, ПРИЗЕМНЫЙ СЛОЙ АТМОСФЕРЫ, ТУРБУЛЕНТНЫЙ ПЕРЕНОС АЭРОЗОЛЕЙ
CONTENTS
1 Introduction
2 Theory of turbulence
2.1 Formal principles of theory
2.2 Transformation of hydrodynamic equations
2.2.1 Input equations
2.2.2 Transformed equations
2.3 Turbulent boundary layer equations
2.4 Theory of turbulent incompressible flows
2.4.1 Turbulent boundary layer in zero pressure gradient
2.4.2 Non-linear theory
2.4.3 Mean velocity profile in turbulent flow over smooth surface
2.4.4 Streamwise turbulent intensity profiles
2.4.5 Normal and transversal turbulent intensity profiles
2.5 Dynamic roughness surface equation
3 Model of turbulent flows over rough surface
3.1 Empirical models of turbulent flow over rough surfaces
3.2 Model of wall roughness effect on turbulent flow
3.3 Modelling of roughness density effect. 3D roughness elements
3.4 Modelling of roughness density effect. 2D roughness elements
3.5 Model of total length of frontal and leeward re-circulation zones
3.6 About efficiency of air pollutants removal
4 Theory of turbulent flows in temperature and pressure gradients
4.1 Thermal boundary layer
4.2 Turbulent flow in the stratified surface layer
4.3 Turbulent boundary layer in pressure gradient
4.4 Planetary boundary layer
5 Turbulent transports in the atmospheric surface layer
UDC 532.526.4
THEORY OF TURBULENCE AND SIMULATION OF TURBULENT TRANSPORT IN THE ATMOSPHERE PARTS 1, 2
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 over rough surface. Numerical solutions of equations system of turbulent transport of admixtures in a surface layer of the atmosphere for a large scale have been studied
Keywords: ACCELERATED FLOW, AEROSOL TURBULENT TRANSPORT, ATMOSPHERIC STRATIFIED FLOW, ATMOSPHERIC TURBULENCE, ATMOSPHERIC SURFASE LAYER, BOUNDARY LAYER, ROUGH SURFACE, TURBULENT TRANSPORT
5.1 Model of turbulent transport in the stratified surface layer
5.2 Numerical modelling of air pollutants turbulent transport in region with several roads
5.2.1 Turbulent transport model description
5.2.2 Turbulent transport over the region with several parallel roads
5.2.3 Air pollutants turbulent transport in Central Sochi
5.2.4 Turbulent transport of lead aerosols in Central Sochi
5.2.5 Air pollution impact on human health in the Sochi region
5.3 Model of aerosols turbulent transport
6 Dynamics of boundary layer
6.1 Boundary layer structure
6.2 Laminar boundary layer
6.3 Transition to turbulence
6.3.1 Continuous transition to turbulence
6.3.2 3D Transition to turbulence
6.3.3 2D Transition to turbulence
6.4 Spectral characteristics of turbulent flows
1. Introduction
Lately an interest for the problems of estimation of atmosphere air quality in the cities has increased. Firstly it should be connected with the policy of environment protection which is carried out by developed countries. In 1992, II UN Conference on Environment protection and development adopted the Declaration which enunciated principles of stable development. In particular a strategy has been worked out on limitation and reduction of carbon dioxide emission which can effected catastrophically on the world climate. In 1993 in Russia a national plan on environment protection and stable development was worked out. The main principles of this plan became the basic in the policy of management in the regions [2]. In these declarations and plans atmospheric air is considered a resource, qualitative consumption of which should be guaranteed not just for the present generation of human beings but for future generations as well.
The problem of air pollution is tightly connected with the development of industry, transport and energetic. Continuous process of coal, natural gas and organic fuel burning aimed at obtaining electrical energy and heat, wide spread of automobile transport, waste of chemical plants and metallurgical works - all this leads to accumulation of different chemical compounds in the atmosphere, which affect the atmosphere composition in a planetary scale. Waste of nitric oxide, sulphur oxide and carbon oxide badly affect different components of biosphere. Health of the urban population in big industrial cities is aggravating because of air pollution. It has been found that air pollution results in building and monuments destruction. Among natural factors of air pollution there are volcanic activity and wind erosion, the latter is partially connected with agricultural development: extension of land under crops and soil destruction because of intensive exploitation.
Practically all industrial countries exercise environment control, in the sense, that waste of harmful substances should be limited, including the atmosphere. For the effective management a special control of air quality is needed, which includes a certain waste volume limitation, taking into consideration the
current composition of atmosphere. In its turn air composition is determined by continuous and periodical measurements in the local points, as well as at mobile stations i.e. by monitoring. Even a well-organised network of the observation stations cannot secure true information of the air quality in the whole physical volume in a given region. It is necessary to have a suitable model in order to apply data of observations to the areas where monitoring stations are not available. Thus, the problem arises of mathematical modelling of the quality of free air, tightly connected with a solution of a problem of diffusion of admixtures in the atmosphere under a given emission [3].
Thanks to acknowledging by the community and government, the importance of solution of the problem of ensuring permanent monitoring of air pollution, such powerful research trends as EUROTRAC-program which includes over 250 research groups in 24 European countries [4] - have been supported. It should be noted that the main aim of EUROTRAC - is co-ordination of programs of scientific researches of transboundary transport and chemical transformation of admixtures in the troposphere over Europe.
Mathematical modelling of the air quality is getting more and more effective instrument in the analysis of atmosphere condition, due to the rapid development of electronic computers and decrease in their cost, and also perfection of mathematical models of transport of gaseous fluid- and solid dispersed components of pollution. Systems of free air quality modelling have been created in big cities such as Paris [5], Lisbon [6], Budapest [7], on a planetary scale [8], and over huge regions as Western Europe or northern latitudes of Eurasia [9-10], as well as in small towns such as Oxford or Cambridge [11] and even in the central street of London [12].
The mane peculiarity of the above mentioned models [5-12] and others is that modelling of admixtures transport is carried out on the base of systems of diffusion equations with coefficients depended on parameters of atmospheric turbulence. It should be noted that modelling of air flows in the boundary layer of atmosphere is a complicated problem, solution of which depends on the theoretical ideas about turbulence. Thus, a simplified parameterization of planetary boundary layer [see 3, 8, 10], which allows completely eliminate calculations of complicated atmospheric flows, using data of meteorological parameters - is a popular method in description of global transport of admixture in the atmosphere.
But this approach cannot be realized in solution the problem of transport of admixture in the lower layers of atmosphere, which is very important for the urban air quality control.
Indeed, transboundary transport of admixtures influences mainly a formation of background value of air pollution level in the regions, while to determine the local level of pollution it is necessary, first, to take into consideration local
centres of emission. But among all sources of pollution, automobile transport in most regions is an outright winner in the way of intensive pollution of the atmosphere. For example, traffic contribution into air pollution in Sochi (Russia) makes up over 80% [2, 13-14] that can be compared with Athens (75% from total mass of wastes NOx) and London (63% NOx) [15].
So, a problem of transport of admixtures into the areas having one or more highways is one of the priorities. Typical space scale in such task should be several hundreds meters from axis of the road and dozens meters in height. Thus the process of turbulent diffusion of traffic wastes is localized in the surface layer of atmosphere and greatly depends on parameters of turbulence as well as on the condition of stratification. Besides heterogeneity of ground surface, including artificial roughness in a form of buildings, trees, etc [16] is a very important factor in this problem.
Turbulent transport of moment, heat and mass in a surface layer greatly depends on heat convection [17]. Due to a great scale of atmospheric flows, convection in the surface layer is observed in a kind of turbulent motions of complex structure [18-20]. Nevertheless a turbulent flow in a thermal stratified surface layer is considered as a spatial homogeneous steady flow with well determined average parameters, which depend only on a coordinate across a boundary layer. Above mentioned approach to the description of turbulent flows in a stratified flow has been developed by Monin & Obukhov [21, 22]. Monin and Obukhov’s similarity theory has been further developed and has got experimental confirmation by numerous researches [23-27].
Modern models of turbulent transport of admixtures in the lower part of atmosphere are based, mainly, on different modifications of k -e model [28-29], which in its turn represents expansion of mean according to Reynolds equations of Navier-Stokes for the purpose of their closing. Some closures for turbulent stratified flow and with due regard for the planet rotation have been presented in [30-35]. A review of the main methods of atmospheric boundary layer modelling is given by authors [36-38].
As it is known, the main idea in the description of mean turbulent motions of a fluid consists in representations of a vector of flow velocity as a sum of vectors of average, by time, velocity of flow and velocity of pulsate motion: u = (u) + u', where by definition
The averaging time interval in (1.1) is suggested quite enough to exclude the chaotic pulsate movement of a fluid. Using the averaging procedure (1.1) in the
(1.1)
Navier-Stokes model one can derive the momentum equation for turbulent incompressible flow as follows
dUt dUiU]. l dP V2TT 1 dtv (1 2)
—L +--------L +--= vV Ui +--------------------------------------L (1.2)
dt dxj p dxi p dx]
where U = ( u}, P = (p) is the mean pressure, ti] =-p(uiuis the Reynolds stress tensor.
This approach was developed over 100 years ago by O. Reynolds [39], who presented equations of turbulent flows (1.2), widely used nowadays in engineer application as, for example, models [28-29].
It should be noted that according to Reynolds vectors (u), u ' are not
supposed to be of any solutions of the Navier-Stokes equations. So, resultant equations (1.2) in common case are unclosed and thus additional physical ideas are needed to close the model. In the simplest cases model closures for incompressible flows in a boundary layer is carried out the base of the Boussinesq's conception of eddy viscosity [40]. Thus for a turbulent flow in OX direction we have
ï’xcz =-p{uxu\) = vt U (13)
where z is the coordinate normal to the wall
Prandtl [41-42] offered a handy formula of turbulent viscosity, using a parameter of the mixing length, by analogy with the length of molecules run, as follows
= pl2
ЭЦ.
dz
(1.4)
where l is the mixing length.
The initial Prandtl’s model (1.4) is notable for its extraordinary simplicity, because the parameter of a mixing length is estimated out of independent hypotheses based on the similarity theory. Prandtl’s ideas turned out to be fruitful and constructive, so most of the models of turbulent flows in the boundary layer, in any case, are based on the hypotheses of eddy viscosity and a mixing length.
Kolmogorov [43] suggested more general expression of eddy viscosity, considering it dependent on an average kinetic energy of turbulence, which can be estimated on the base of hypothetical equations, in the form
¿Ex
dz
Кz = - p(uu) = PnTlkk1/2 (1.5)
where Cm is the constant, lk is the turbulent length scale, k = (ukuk)/2 is the turbulent kinetic energy.
Launder et al [28] suggested another form of expression (1.5) introduced the
turbulent dissipation rate eT = k3 2 /lk, thus with this parameter the eddy viscosity (1.5) can be rewritten as
The various modifications of the k - e model are based mainly on the eddy viscosity model (1.6). The standard k -e model is included the continuity equation for the mean velocity, the momentum equation (1.2), the turbulent transport equations for 'substances' k, eT, and a generalised form of the first equation (1.5) for Reynolds stress tensor:
The main problem, peculiar to the approach mentioned above, is difficulty to control various physical effects which are the key in the process of turbulent transport. As far as, according to Reynolds point of view, functions (u), u' are not the solutions of initial equations, their physical sense is not obvious enough to connect them with the forces functioning in the system. So, model of turbulence [28-35] as well as some others are notable for numerous of different constants and parameters introduced for co-ordination of estimated values with experimental data.
So, in a case of atmospheric stratified flows, following Rodi [29], k -e model can be written as
where p0 is the reference density, gt is the gravitational acceleration, T,T are the mean temperature and temperature fluctuations respectively, Prt is the turbulent Prandtl number.
The closures for model (1.8) can be derived by the common method explained by Rodi [29], as follows
vT = Cm k2 / eT
(1.6)
(1.7)
(1.8)
/
2
3
kd„, -(u,r)=
Pr Эх,.
vT ЭТ
Эк Эк Э
— + U,----------------= —
dt 1 Эх, Эх,.
vT Эк sк Эх,
V к 1 У
+ П + Г — £т
(1.9)
Эет + U д£т
Эt 1 Эх,.
_Э_
Эх,
^ v т Эет ^
se Эх,
Ve 1 У
+с„ек- (П+Ce Г)—сг,к кк
where n, r are stress and buoyancy production of the turbulent kinetic energy, respectively. The last terms can be expressed as
П = v T
'Эи Эи,Л
____l_ +___J_
Эх, Эх,
v J
Эи г = pgt v'т- (i.io)
Эх, 1 Prt Эх,
here b = -p ~l(dp / d T) is the coefficient of expansion, b = 1/ T for the perfect gases.
Model (1.8-1.10) depends on the six parameters C1e,C2e,C3e, Prf,aE,sk. Two of them, Cse and the turbulent Prandtl number, Prt, are not constants for arbitrary stratification, but its dependent on the stability parameter.
As it has been shown by Apsley & Castro [44] the k - e model requires modification in the atmospheric boundary layer. For instance, in a case of a neutral stratification, the standard k - e model (1.8) - (1.10) yields a very deep boundary layer and large friction velocity. They supposed the new form of the k - e model extended for the case of a stable stratification.
Thus, in general case the problem of closing the models of turbulent transport has not got a comprehensive solution, because the parameters of atmospheric turbulence depend on roughness of ground surface, local heat flux, pres-sure-gradient and velocity of the planet rotation (or Coriolis forces). So, a turbulent flow in the boundary layer of the atmosphere is a three-dimensional and not stationary and can be described in particular cases only, for example, in surface layer, in which the conditions of homogeneity and quasi stationary for the mean turbulent fluxes of the heat and momentum are carried out (these hypotheses are the bases of the Monin-Obukhov’s similarity theory [21-22]).
The aim of this work is development of the theory of turbulent diffusion of admixtures in the lower part of atmosphere, which can be applied to describe the flows with arbitrary stratification. In the work process it became vivid that hypotheses about influence of buoyancy forces and roughness of ground surface on turbulent boundary layer should have been revised. It became possible mainly due to simplicity of description of external forces within the theory of turbulence [45-48].
In this paper the basic principles of the developing theory of turbulence are given and the idea of mechanisms of influence of static and dynamic roughness on turbulent flows parameters in the boundary layer is substantiated. The averaging method of turbulent flow parameters in the boundary layer has been suggested, application of which allows distinguishing among all solutions of Na-vier-Stokes equations the functions, similar, according to the characteristics, to those which are to be observed in experiments. Formulas are given of hydrodynamic equations transformation to curvilinear nonstationary system of coordinates, connected with chosen surface, which is modelled of the dynamic roughness. The formal transformation of hydrodynamic equations in the surface layer is described. On the base of the equations of dynamics of viscous, heat-conducting fluid and an equation of diffusion, by application of transformations of the type mentioned, the system of dynamic equations for random amplitudes of flow parameters has been obtained: admixture velocity, pressure, temperature and concentration. For realistic modelling of atmospheric flows, the buoyancy forces are taking into consideration, produced by thermal expansion of the air, for the description of which Boussinesq approximation is adopted. The analysis of similarity of the system of equations in the case of steady, non-isothermal flow in longitudinal pressure gradient has been accomplished.
2 Theory of turbulence
2.1 Formal principles of theory
Modern models of turbulent flows are mainly based on the Navie-Stokes equations averaged according to Reynolds as it was mentioned in Introduction. It is clear, that any model, based on presentation of the velocity field as a sum of mean velocity and velocity of pulsate motion: u = (u) + u', in which vectors (u), u' are not the solution of initial equations, should not be closed. And this is the main problem of turbulent flows modelling - mean equations are not adequate to initial Navie-Stokes model, and therefore are not complete. Further locking of equations obtained, is practically based on experimental data interpretations, so models of turbulent flows are notable for excess adjusting parameters which are introduced to register an effect of real forces, functioning in the system. In spite of notable success in this sphere, Cantwell [52] passed a remark that with the exclusion of some results obtained from reasons of dimensionality the simplest problems for turbulent flow with the simplest out of possible boundary conditions, cannot be solved up to now. Therefore the methods of direct numerical simulation (DNS) of turbulent flows have been developed for the last decade [53], but their application to the atmospheric boundary layer is still a problem due to a lack of resources of modern computers [54].
It should be noted that if vector (u) is somehow a solution of initial equations, then the model which has been obtained by Reynolds’ method, will be closed without any other additional hypothesis. In this case the main target is a
search for a suitable solution, u , the characteristics of which should be conformed to the experimentally estimated profile of the mean velocity. Algorithm of finding such type of solutions for the flows in a turbulent boundary later was suggested in the papers [45-48]. It consists in transformation of Navier-Stokes initial equations so, that they would contain eddy viscosity, as in the Prandtl theory of mixing length. Let’s consider this algorithm in detail.
Note that the field of velocity in a turbulent flow is not a regular function, but it is not a random function as well, because it doesn’t depend on random parameters. The main idea is to introduce the random parameters into the hydrodynamic equations, the flow velocity evidently depends. This is possible in a particular case of flows in a boundary layer, for which a surface layer transformation is performed, i.e. presentation of a flow velocity vector, u = (u,v,w) , as u = u( x, y, z/h( x, y. t), t), where z is the normal to the wall variable, h = h( x, y, t) - is the dynamic roughness surface adjoining to the wall but not coinciding with it. Physical interpretation and the model of the dynamic roughness surface will be given below.
The dynamic roughness surface in a turbulent flow is advised to characterise by the set of the random continuous parameters h, ht, hx, hy , which have a meaning of the height, inclination and velocity of transport of surface elements, with a well-known function of distribution fs = fs (h, hx, hy, ht).
We suppose that h = z / h = const and 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 - Figure 2.1. Let us consider subregion dVs , lying in the representative region of the flow, dV , in which random parameters
h, ht, hx, hy are changed in the intervals
(h;h + dh) ,(ht; ht + dht), (hx; hx + dhx) ,(hy;hy + dhy) accordingly. In common case the subregion dVs is a multiply connected domain. The volume of this subregion is given by
dVs = dVfs (h, hx, hy, ht )dhdhxdhydht.
The random amplitude of 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,hy,ht) = lim — \ u(x,y,h,t)dxdydz (2.1)
SV®dVs OV
OV
where dV 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 func-
tions, u = u(h, t, h, hx, hy, ht), immediately follow from the hydrodynamic equations of viscous fluid and its are adduced below in subsection 2.3. Statistical moment of an order m of the random function u(h, t, h, hx, hy, ht) are determined as
follows
u m (Z t) = J u m (h t, К hx, hy, ht )fs (h, hx, hy, ht )dhdhxdhydht
(2.2)
Thus, in the given theory the mean velocity (and any other mean value) is determined in two steps. On a first steps accordingly to equation (2.1) the random functions u(h, t, h, hx, hy, ht) are calculated, on the second step accordingly to
the equation (2.2) the mean values u m (z, t) are calculated. The similar algorithm has been proposed by Trunev & Fomin [55], and Trunev [56] for the problem of impingement erosion.
The transformation (2.1) can also be developed for the turbulent flows over rough surfaces. One can assume, that the rough surface can be described by a function z = r(x,y) and that the dynamic roughness surface presented as h( x, y, t) = r (x, y)+h( x, y, t), where h( x, y, t) is the thickness of viscous sublayer. Then, in a case of smooth surface we have at z = 0: u = 0, and in a case of a rough surface at z = r: u = 0. Therefore the factor of static roughness enters in the model as a random parameter. The basic formulas of transformation of flow parameters in the surface layer over a rough surface are given by
~ 1 г
u(h, t, r,h, hx, h , ht) = lim — J u(x,y, h, t)dxdydz
7 SV®dVs OV
(2.3)
SV
um (z, t) = J uгт (h, t, r, h, hx, hy, ht) fs (r, h, hx, hy, ht )drdhdhxdhydht
Figure 2.1: Selection of coordinate system at the description of turbulent flows
Let us consider two examples. If u = u* ln(z / r), then the mean velocity, calculated on second equation (2.3) can be written as follows
u(z) = |u(z /r)fs(r)dr = u* ln(z /r*) ,
where r* = exp(J ln(r) fs (r)dr). Assuming, that u = u* ln(z / r) is some solution of the Navier-Stokes equations one can discover, that (u) = ~(z)is also the solution of these equations. If the random amplitude of velocity is u = u*(z/r)b and this function is the solution of Navier-Stokes equations, then the mean velocity profile is also the solution (u) = ~(z) = u*(z / r*)b at fixed value of roughness parameter: r = r* = (J r ~bfs (r )dr) .
The main problem on this way is how to estimate the multiple density of a probability distribution function fs = fs (r, h, hx, hy, ht)? Nevertheless, for the solutions presented by the logarithmic function one can suppose that
u = J ¿/(hi, t, r, h, hx, hy, ht) fs (r, h, hx, hy, ht )drdhdhxdhydht =
= ~(z1/ h*, r*, h*, h^, h*)
where the parameters with stars can be estimated from the experimental data or calculated from the theory of turbulence. Practically the roughness parameter r* should be given as an input value and all another parameters can be calculated from the similarity theory considered in sections 2.6-2.7.
2.2. Transformation of hydrodynamic equations
2.2.1. Input equations
We shall consider an air flow containing a scalar impurity. Air is assumed as a viscous, heat-conducting, incompressible gas in a rather slow turbulent motion. It is well known fact that the surface layer comprises one-tenth of the planetary boundary layer, in which the earth’s rotation effect can be neglected (see Arya [18]). Thus, the model of the turbulent flow in the atmospheric surface layer can be written as follows:
V. u = 0 (2.4)
+ (u.V)u + Vp = vV2u + (p - p 0)
at: p 0 p0
— + (u.V)T = — VT at v ' Pr
f + (u.V)f = SV- V2f at Sc
where p is the air density, u = (u,v, w) is the flow velocity vector, v is the kinematics viscosity, p is the pressure with the exception of the hydrostatic atmospheric pressure, g is the gravity acceleration vector, p0 is the equilibrium
density, T is the temperature, Pr is the Prandtl number, f is the mass concentration of an impurity, Sc = v / D is the Schmidt number, D is the molecular diffusion coefficient. The hydrostatic equation and the standard Boussinesq approximation for the density fluctuations are given by
VP0 = gP0(P0. U p - p0 =-p0b(T - T0), (2.5)
where b = -p -\dp / 0T)p is the coefficient of expansion, b = 1/ T for the perfect gases.
The coordinates system should be determined in such a way that the Z -axis is directed opposite to the vector of the gravity acceleration. The relief of a ground surface is given by the equation z = r (x, y) - see Figure 2.1.
Boundary conditions for the flow parameters are set on the ground surface and on the top of boundary layer are set as follows:
z = r (x, y): u = 0, T = Tg, f = f g (2.6)
z = H: u = ^0(1,0,0), T = T0, f = f 0-
where Tg is the surface temperature, f g is the impurity concentration over the ground, H is the boundary layer height, U0 is the wind velocity at the height z = H, T0,f 0 are the temperature and the impurity concentration at the height z = H respectively.
The boundary conditions (2.6) are spatial homogeneous on horizontal coordinates. This simplification is accepted only at a derivation of the main equations of turbulent flows. Hereinafter the boundary conditions will vary as required, to result them in conformity with a type of soluble problems.
2.2.2. Transformed equations
The transformation (2.3) can be applied to equations (2.4) in the form S(h, t, r, h, hx, hy, ht) = lim -1 J S(x, y, h, t)dxdydz
7 SV®dVs dv J
UV SV
where S = (u,p,T,<~~),S = (u,p,T,f). The equations for random functions S = (u,T ,T,<T )can be obtained directly from system (2.4) recorded in curvilinear non-stationary coordinates (x, y ,h, t), where h = z / h( x, y, t). The transition to curvilinear non-stationary coordinates system in hydrodynamic equations is explicitly described by [59-60] and other. Following Pulliam&Steger [59] equations
(2.4) can be presented in the form:
^ + |- (E - E,.) + 0^ (F - F,) + J (G - G,.) + J [h,Q+K (E - E,.) + hy (F - F, )] = B (2.7)
at ax ay ah
where J is the Jacobian of transformation, J = h 1 ф 0, J 1 = h ф 0,
■O 0 ^ P 0 U ' v 0 P ^ P oW ^
о a P 0u 2 + p P 0u v P 0uW- hhxp
P 0v O 0 w , E = P 0vu P 0 wu , F = P 0v2 + p P 0w v , G = P 0v W - h hyp P 0wW + p
T uT v T WT
v f 0 v uf 0 v v f 0 v Wf 0
E =
œ 0 " œ 0 " œ 0 " œ 0 "
Txv h & 0
tx F = ’ v ty Tzy , Gv =1 i ht 1 i yi h , B = 0 - g(P - P0)
n Pr-1Tx n Pr-1Ty nPr 'h T 0
1 x J Df y v Dh f i 0 v 0 0
Here tkl is the tensor of viscous stress, Tkl = m cosity, k, l = 1,2,3 ; h x = -JhK, h y = -Jh hy, h z
du, du
+ -
dx, dx,
m is the dynamic vis-
J, W = w - h(ht + hxu + hv v), i = x, y, z
In the curvilinear coordinates system it is necessary to execute replacements in terms with gradients as follows:
a
d dh d
--------®-------1----------
dx, dx, dxj dh
-P • 1 о d 1 d
for j = 1,2;-------------------®---------
dz h dh
The equations (2.7) in a common case look rather cumbersome and here are not considered. The transformation of the Navier-Stokes model to the non-stationary curvilinear coordinate system is only the convenient method to allocate some thin effect of rough surface in the viscous flow. The turbulent flow over a smooth surface can be considered as the limiting case of flow over a rough surface, when the influence of viscous sublayer exceeds effect of the static roughness.
At first sight it seems that the transition to the new coordinate system (x, y,h, t) is connected with a selection of dynamic roughness surface, and thus it can't be defined as the univalent transformation. Nevertheless, as it will be shown hereinafter, the requirements imposed at calculation of the mean velocity profile allow us to define the mean parameters h*,h*x,h*,h* accurate to an undefined factor, which connects to the scale of the turbulent boundary layer.
2.3 Turbulent boundary layer equations
Let us consider the special type of solution of transformed equations (2.7) which depends on time and normal variable h as it often supposed in the turbulent boundary layer theory. Thus put
d d
— (E - Ev) = — (F - Fv) = 0 d x dy
in the left part of (2.7). In this case eq. (2.7) can be presented in the form
+ J dh (G - Gv) + j[htQ+hx (E - Ev)+hy (F - Fv )]=B (2.8)
dt dh
To derive the turbulent boundary layer model one can apply the averaging operator in the form (2.3) with an arbitrary averaging volume 5V to equation
(2.8) to conserve the commutative properties of the averaging operator with the space and time differential operators. Then one can consider the limit of all terms of the averaged equation at 5V ® dVs . At this step the theorem about two limits of the continuous function can be used (since the differential operators can be considered as some limits). Finally we have
-dQ- + J — (G - G v ) + J ht Q + hx (E - E v ) + hy (F - Fv )] = B (2 9)
where parameters with tilde are defined similar to S = (U,p,~,C). If
S = (u,p,T,C) is the solution of the transformed Navier-Stokes equations (2.8) in
any sense, then we have the turbulence model closures automatically as follows:
1 f uu
lim I U (x,yh, tK (x,y, h, t)dxdydz = utuk + 05lk
5V ®dVs 5V
SV
1 Г ~ ~
lim — I u(x,y,h,t)T(x,y,h,t)dxdydz = uT
SV ®dVs SV
SV
1 f ~ ~
lim — I u(x,y, h, t)f (x,y, h, t)dxdydz = uf
SV®dVs OV
SV
Here 3d /2 is the kinetic energy of turbulent fluctuations in the small volume dVs, S lk is the Kronecker delta: S lk = 0 for i * k, S lk = 1 for i = k.
Therefore (2.8) and (2.9) are the identical equations. The dissipative terms in
(2.9) can be written as follows
E =
f 0 " f 0 " f 0 ^
2mhxuv mhx ~+mhy~h
mhyUh + phx ~ mhxw -mT /r F =-h v h 2h t ßyW -m T/r Gv = (1+rhf)1 dh dr m T ¡Ж
nPr~XTn vPr-1hT7h vPr-1~
Dhx~ \ ' / v DhT , d~ V r /
where uh = du /dh , ....
Substituting in (2.9) the expressions of all terms finally we have the dynamic equations for random flow parameters as follows:
Note, that the parameters of a dynamic roughness in equations (2.10), are not already the functions of space variables or time. Really, in virtue of transformation (10), the values of these parameters are fixed in intervals from r up to r + dr , from h up to h + dh , from ht up to ht + dht , from hx up to hx + dhx , from hy up to hy + dhy . These values, thus, are considered as the random parameters, and the law of their distribution in specific intervals is described by a known function fs = fs (r, h, hx, hy, ht).
As we can see from the derived equations (2.10) there are the factors in the second derivatives terms, which depend on a distance up to a rigid surface. It should be noted also, that the equation (2.9) is not in the strong conservation form, as, for example, it is given by Pulliam & Steger [59]. Therefore the number of terms in a square brackets, breaking conservation of this system are chosen in the left part of equations (2.8) and (2.9). Such allocation of non-divergent terms is stipulated by the purposes of modelling of the eddy viscosity, which, in our opinion, arises in a boundary layer from the transformation of a tensor of viscous stresses near the dynamic roughness surface. It is obvious in the case of viscous flow over a rigid rough surface and is connected with an adhesion of a viscous flow to a rigid surface of any configuration. In the turbulent flow over a smooth surface the eddy viscosity is simulated by analogy to a more widespread type of turbulent flows, as in a special case, when r ® 0. Thus the eddy viscosity is connected (mathematically) with the transformation of a tensor of viscous stresses to coordinate mapping which brings rigid surface onto coordinate surface.
Cu W Эи N C
-------I------------I---------;
Ct h Ch p0h с
(2.10)
where W = w - hF , F = ht + hxu + hy v, P = p + в, n = ^hl + h2y , N = (-Ц hx,-h hy ,1).
In this model the Reynolds stress can be calculated as follows
ïk (z> 0 = J p[(~ - ~ )(Sk - )+^k 1/, (r K hx, hy, K )drdhdlhdhydht
Therefore the random function u = u(h, t, r,...) gives the main contribution in the non-diagonal components of the Reynolds stress. Now we take it as granted because we haven't any contradictions. Hence, the first assumption of this theory is that the turbulence interaction between the hydrodynamic fields can be described with the solutions S = (u,~,~, C) as well as with the solutions S = (u, p, T, C). The second assumption is that it's possible to neglect longitudinal and transversal gradients of flow parameters in a comparison with gradients across a boundary layer, at least for steady turbulent flow.
For the diffusion equation it is possible to derive the boundary layer model by the simplified way. Let us suppose that in the last equation (2.4)
C = C (h ( x, y, z, t ), t ), then we have
— + (u. V)f - — V2f = — + (ht + u. Vh) — - D(Vh)2 - DV2h — =0
at V Sc y dt yit Udh dh2 dh
In partial case when h = z / K( x, y, t ) thus
Vh = h _1(-hhx ,-hhy ,1); V 2h = -hh _1V2 h + 2h(h _1VK)2,
And therefore the last equation can be written as
a— w a— = d (i 2 2)a — 2Dn2h a— dvV2k a— at k ah k2 h ah2 k2 ah k ah
This equation can be transformed to the form of the last equation (2.10). According to definition
f (h, t, r,h, hx, hy, ht ) = ^lim —— J f (h, t)dxdydz,
S— ®d—sS—
S—
Using the identity h V2h = V(h 'Vh) + h 2(Vh)2, and averaging all terms, finally we have
= ^2(1 + w2h2)—- lim — fdzDh f \V(h~lVh)dxdy
J'l 1 J 3n2 h2 5V®dv. AT/J ^
dt h dh h2 dh2 h2 dh S— ®d—s S—r~ 'dh
AS
Where DS = LxLy. But the last term is annulled if region DS = LxLy is large enough (the divergence theorem). Therefore we have an equation
dC W dC D 2 2ч d2C Dn2h dC n d „ 2 2ч dC nn2h dC
----+---------------------------------------------------------------------= —г(1 + П h )-^ ^-=-Г (1 + П h )-2
dt h dh h dh h dh Sch dh dh Sch dh
which is identical to the last equation (2.10).
The hydrodynamic part of the system (2.10) can be transformed to a form convenient for integration. For this purpose one can multiple the second equation (2.10) by scalar way on vectors
N, N = (hx, hy ,0), n2 = (hy ,-hx ,0)
Expressing the pressure gradient through other parameters one can derive
dP _ 2vrLdF-gh (u - r ) - (211)
dh h dh n2 (r p0) n2 dt ( . )
dF W dF n2h dP n d „ 2 2x dF 2nn2h dF
---+----------------_——(1 + n h )-------------^-----
dt h dh 0h dh h2 dh dh h2 dh
dY W dY n d „ 2 ^ dY vn2h dY
----+-------------_—r—(1 + n h )----^----
dt h dh h2 dh dh h2 dh
where N2 _ 1 + n2h2, Y _ hyu - hxv .
The system (2.11) can be closed using the continuity equation, which is represented for this purpose as follows
+ f_ 0 (2.11, а)
dh
Let us consider some similarity properties of the system (2.11-2.11^) in case of the steady turbulent flow. Put
dF/dt _dY/dt _d#/dt _ 0
in (2.11). Then the dimensionless components of a flow velocity and their combination can be written as follows
u + _ ~ / u*,v+ _ v/ u*, w + _ w / u*, y _ Y / nu* _ u + sin a - v + cosa, j _ F / nu* _ u + cosa + v + sin a + w0 +,
Z _ j - w +,
where u* _ ux _^tj~p is the friction velocity, tw is the wall shear stress, X _ z /1, 1 _ h / n,a _ arctan(hy / hx), w+ _ ht / nu*, w0 _ ht / n is the second scale of velocity in the turbulent boundary layer.
Substituting a gradient of pressure from the first equation (2.11) in second one, expressing the buoyancy force through the temperature according to approximation (2.5) and using dimensionless variables finally we have:
I_j (2-12)
/1 е2\ d2j ,.+ Xl + BT4
(1 + X2)-т2 + (l Ж + 2X)-%т +ь
dX2 dX 1+x2
(1+X2) dy+(1+Z+X) % _ 0
dX2 dX
where 1 +_ 1v / u*, T+_ (Tg - ~) / T* is the dimensionless temperature, 5 _ -gb T*v / u*3, T* _ qH / (p cpu*) is the turbulent scale of temperature, qH is the heat flux from the ground to the air, cp is the specific heat at the constant pressure of the gas.
It's obvious, that the solutions of the equations system (2.12) are not evidently dependent on the choice of the scale h, and depend on the combinations of random parameters
1 _ h / n, a _ arctan(hy / hx), w+ _ ht / nu*.
Thus, the turbulent steady flow equations (2.12) are applicable to the arbitrary scales flows, including the atmospheric flows. This property is also applied to the unsteady flows, if the dimensionless time is determined as follows t _ u*t /1. For example, the equation for the contravariant component of velocity in common case can be written as
d 2w + + W +9W - ± A (1+X 2) _ _L_ dW+ + hX (p/ p 0 -1) (213)
dXdt dX2 1 + dX X ) dX2 1 + X2 dt 1 + X2 )
where W+_ J~/u*, b0 _ 1 g/ u*2.
It should be noted that the hydrodynamic part of the system (2.10), as it follows from (2.11), describes flow, which has all components of velocity. That is the essential difference of this model from the standard models [28-29], which are based on the Navier-Stokes equations averaged accordingly to Reynolds. Using this feature it is possible to simulate the turbulent intensity profiles in the turbulent boundary layer as well as the mean velocity, temperature and impurity concentration.
2.4 Theory of turbulent incompressible flows
2.4.1 Turbulent boundary layer in zero pressure gradient
The turbulent boundary layer over a smooth surface is the well investigated flow [51, 61-62]. We shall consider the steady turbulent flow of incompressible fluid in a boundary layer in zero pressure gradients. Let's assume in equations (2.11)
9Y _0
dt dt dt
Then the normal pressure gradient can be written as follows
dP _ 2vp dF dh h dh
Using this expression, we can transform the equations (2.11) to a form convenient for integration:
+ F_ 0 (2.14)
dh
W dW V d „ 2 2x d2W
(1 + nh )^Y~ h dh 2 h2 dh dh 2
Wu dY v d 2 2 dY vn2h dY
-------_—r (1 + n h )--------------^-----
h dh h2 dh dh h2 dh
Here the second equation is obtained by a substitution O from the continuity equation (1.11, a) in the second equation (1.11). The obtained equations can be investigated in general case. Integrating these equations one time, we have
dF Aj exp[-I(h)] dY A2 exp[-I(h)]
' dh
dh
1 + n 2h2
a/1 + n 2h2
d 2W _ Aj exp[ -1 (h )] dh2
1 + n 2h2
h çv Wdh
2 2
vJ 1+n v
The velocity components can be written as
where A are some constants, I _ - — H-‘ ’ vJ 1
(2.15)
[(F - ht )hx + Yhy , v _ n 2 [(F - ht )hy - Yhx , _ W + h F
Then using (2.15), one can derive
du n 2hvA,e 1 n 2hyA2e
+ (2.16)
dh 1 + n2h2 ^1 + n 2h
dV _ n 2hyA1e 1 n~2hxA2e"1 _ A1he
dV 1 + nV -JÏ+nV ’ dv 1 + n2v2
The further analysis of turbulent incompressible flow in a boundary layer will be based on equations (2.16), which will be used for the computing of the mean velocity and turbulent intensity profiles.
I
2.4.2. Non-linear theory
The second equation in system (2.14) is non-linear. We'll consider the numerical method of its solution. In case of the flow over a smooth surface the boundary conditions can be written as follows
h _ 0: J~(0) _ 0, dW/dh _-ht, d2W/dh2 _-A15 (2.17)
where A1 is the free parameter required to obtain the limited value of the inh rv PWdh
tegral I(h) _ - v J0 1 + n 2 2 for h ® ¥ . This condition is used to obtain the logarithmic asymptotic of the mean velocity profile. The parameter A1 has a clear physical sense. because this parameter is directly proportional to the normal pressure gradient on the wall:
_ 2v^_3*, thus ^ ^A.
dh h dh dh h_0 h
0.01 0.1 1 10 100 1000 if 0.01 0.1 1 10 100 R
Figure 2.2: a) The normalised function I(X)/Rt computed for various Rt _-0.83;0.026;0.83;3.32;26.56 - the solid lines 1-5 respectively; b) The normalised integral I0(Rt)/ Rt depending on the dynamic roughness parameter for Rt > 0
First and second boundary conditions (2.17) can be derived from the definition of Wu and from the boundary conditions for the viscous flow velocity on a rigid surface.
To minimise the number of the independent random parameters the general solution of the second equation (2.14) can be written as W _-(ht /n)zi(nh,Rt),
where the universal function z 1 depends on the composition of random parameters (the dynamic roughness Reynolds number):
hht
R = „2 , 2ч (2.18)
' V (И2х + И2)
and satisfies to the equation
(1 + Х 2) dX 3 + (RС 1 +2Х) dX2 0 (2Л9)
with boundary conditions at
X _ 0: Z 1(0) _ 0, dz 1/dX _ 1, d2Z1/dX2 _ a (2.20)
where X _ nh , a is also the free parameter required to obtain the limited value of the integral I(X) at X . Note, that (2.19) can be derived from the second equation (2.14) and that a _ 1+ pX (0)/2 pu*2. Consequently the integral I (X) depends on the composition of random parameters Rt and can be calculated as
I (X, Rt) R' xf+ 'R) dX (2.21)
The integral (2.21) has been computed in the range - 2.5 < Rt < 700 together with problem (2.19-2.20)). Fourth-order scaled Runge-Kutta algorithms and the shooting method (see [63]) have been used to get the numerical solution. Note that for Rt = 0 this problem has the analytical solution:
Z1(X) _ j0X (1 +a arctan X)dX.
In this case I (X) _ 0, hence one can suggest that function
I (X) _ limR, ® I(X)/ R, _ |0{ Z—dr
has a limited value for X . It is possible if only a _ -2 / p, thus
Z1(X) _ j0X (1 - p2arctan X )dX
This solution has been used as the initial position in the shooting method.
The normalized function I(X)/Rt is shown in Figure 2.2,a for various R't _-0.83;0.026;0.83;3.32;26.56 - solid lines 1-5 respectively. As it is shown the function I(X) / Rt is simple and smooth function such arctan(X). The calculated limited value I0(Rt) _ limX®¥ I(X,Rt) is shown in Figure 2.2,b by the symbols together with the approximated line. To simplifier the numerical modelling of the mean velocity profile over a rough surface the function I( X, R )has been approximated as
2
I (X, R) = - I0(R )arctan[(0.4 + 0.02Rt 3/4)£ ] (2.22)
p
where I0( R) is given by
10 (Rt) / Rt = 1.38 - 1.13 arctan[0.4 ln(1 + Rtq)],
J 1, 0 < Rt < 100
4 = 11 - (1.5Rt - 150)10-4, 100 < Rt < 700
Finally note that for the negative value of the parameter Rt in the range Rt < -2 5 the numerical procedure becomes unstable one. In this case the value I0(Rt) _ limX®¥ I(X, Rt) increases considerably with the small decreasing of the dynamic roughness parameter Rt. Since this branch of the integral I0(Rt) will not be used in the analysis, therefore data for the negative value Rt < 0 is not presented in Figure 2.2, b and has been neglected in approximation formula
(2.22).
2.4.3. Mean velocity profile in turbulent flow over smooth surface
The turbulent boundary layer over a smooth surface is the best example for the theoretical consideration and modelling according to the model (2.16). In this case the streamwise velocity gradient can be written in the standard form using the inner layer variables z + _ zux / v, u + _ u / ux, and boundary condition for the mean velocity gradient on the smooth wall:
at z —— 0: du / dz —— 1.
Besides one can require, that at the great distance from the wall the profile of mean velocity is described by the logarithmic function, i.e.:
at z + —¥ du + / dz + — 1 / k z +,
where k is the Karman constant (we use the parameters with stars instead of the random parameters according to subsection 2.1, thus it's the reason why the Karman constant has been taken). Finally we have got for the streamwise velocity gradient
du+ Ae 1 e10 1____________________________________ (2 23)
dz+ 1 + (z+ /1+ )2 + kl 1 + (z+ /1+ )2 ‘
where A _ 1 - eh / kl +, 1 + _ hur / nv .
The first term in the right part (2.23) has the essential value mainly close to the wall (if A ^ 0) and the second one gives the main contribution in the logarithmic layer. To derive the mean velocity profile we should firstly define the parameter A _ 1 -eh / kl +. Note, from first equation (2.16) and (2.23) it follows
that
. 2 du+ (0) . d v+ (0)
A = cos a----------^ + cosa sin a-------------+-L
dz dz
Our suggestion about the dynamical roughness structure 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, Cantwell, Coles & Dimotakis [64]) where the visualisation of the coherent structure in the turbulent boundary layer is presented, and subsection 2.5 here).
Thus for the mean flow A _ 1 - e10/ kl +_ 0, then the length scale 1 + _ hur / nv can be found as the solution of the next equation
k _ < exp[I0 (Rt)] / Rt (2.24)
where Rt _ 1+ w0 / ux, w0 _ ht / n is the second scale of the turbulent velocity. For an arbitrary value w0 the equation (2.24) has two roots or hasn't any roots and only if dk / dRt _ 0 this equation has one root. Hence for the uniqueness of the mean velocity profile should be done
1 dk eh dI0 e1
w+ dRt Rt dRt R,
0 (2.25)
The numerical solution of the equation (2.25) with I0(Rt) determined from
(1.22) gives R, _ R* »1.22 and therefore the predicted values of the turbulent theory constants are given by (for k _ 0.41)
w*+ _ kR* exp(-I0) _ 0.14 , 1 0_ R* / w*+ _ 8.71 (2.26)
The fundamental parameter of length for the turbulent boundary layer is defined from here: 10_ 1 +v/ u* » 8.71v / u*, that almost coincides with the peak of turbulence production, obtained by Klebanoff [49] and Laufer [50]. The fundamental scale of length in the turbulent boundary layer, determined as 10_ R* / w*+ _ 8.71 also coincides with the peak of the frequency diagram of instantaneous thickness of a viscous sublayer [62] - see Fig. 2.3.
If A _ 0 then 1 +_ 1 + in (2.23) and this equation can be presented in the form:
du + eJ°-I
dz + kl+V1 + (z + /1+)2 where I(z* / Л,*,R/) = |0{ R xf+f:)dX.
(2.27)
20
£ 10
15
0
5
♦
0
5
10
15
20 z+
Figure 2.3: Frequency diagram of instantaneous thickness of a viscous sublayer in the turbulent boundary layer [62]
Integrated the first equation (2.27) we have:
Therefore, with the given constant k another constant of the mean velocity logarithmic profile can be calculated from the second equation (2.28). It gives c0 = 5.015 for k = 0.41. The velocity profile calculated with (2.27) for
k = 0.41; l += 8.71 is shown in Figure 2.4, a by the solid line (1). For this profile
the shooting parameter in (2.20) is estimated as a = -2/p - 0.27 » 0.9066. The predicted profile (1) has been compared with the mean velocity profile computed on the model of the transitional layer proposed by Van Driest [65] - the solid line (2). As explained by Cebeci & Bradshaw [51] the Van Driest's model can be written in the form
The standard logarithmic profile can be derived from here at z + >> l +:
ln z + + c.
0 5
(2.28)
du (-< u' w '>)12 du 2
— = —---------------------------------------------------------------, -< u'w'>+v— = uT2 (2.29)
dz k z[1 — exp( - z / ld )]’ dz y ’
Here ld is the damping length, ld = 26n / ut. The model (2.29) is based on the
Prandtl theory of mixing length. Thus, this model depends on two parameters:
the Karman constant and the damping length.
Figure 2.4: a) Mean velocity profiles in the turbulent boundary layer for inner region calculated on eq. (2.27) - (1), and with Van Driest model (2.29) -(2); b) comparison of computed profiles (1-2) with DNS data by Kuroda et al [66] - (3), and with experimental data by Nagano et al [68] and Smith [69] - (5); c) mean velocity defect low in the outer region: solid line 1 is calculated on eq. (2.34) for v* = 0.4 and experimental data (2,3) by Nagano et al [68] obtained at x = 0.525m;1.125m accordingly; d) mean velocity defect low in the outer region: solid line 1 is calculated on eq. (2.35) for v* = 0.27; e = 0.8 with experimental data (2,3) by Nagano et al [68]
The profile computed on (1.29) coincides with the predicted profile 1 in the viscous sublayer and in the logarithmic layer but differs a bit in the transitional layer (see Figure 2.4, a). This difference can be explained by the pressure gradient effect. Figure 2.4,b demonstrates the comparison of both profiles (1,2) with several data bases: 3 - the direct numerical simulation of the turbulent flow in the two-dimensional channel (Re = 2980) by Kuroda et al [66]; 4 - the turbulent boundary layer in zero pressure gradient (the Reynolds number based on the
momentum thickness Ree = 1040) by Nagano et al [67,68] and 5 - the turbulent
boundary layer mean velocity profile (Re = 13052) presented by Smith [69]. Note that both profiles (1,2) well correlated with computed and experimental data.
In the upper layer for z + > H + = Hu*Jv the mean velocity profile should be constant in contrast to the logarithmic profile which diverges at z + ® ¥. As it is well known in the outer region of the turbulent boundary layer the mean velocity profile can be described by the defect low:
U+ =—---------------= — ln — + FH
ut k H H
H0
where the universal function FH (g) has the properties:
(2.30)
dF 1
lim F = 0; lim FH (g) = — — ln g (2.31)
g®0 dg g®¥ k b 7
Here the first condition means that the universal function does not change boundary conditions on a rigid wall. The second condition can be used to obtain the limiting value of the mean velocity u = U0.
One can suggest that the mixed layer turbulence is generated in the same way as the wall turbulence. Then the new dynamic roughness surface can be introduced and the equation system which is similar to (2.16) can be derived. In the case of the mixing layer we can put I = I0 and l = g*H, where g* is the parameter. Using the new characteristic scale and equation (2.30) the velocity gradient can be written as
d-f = ± — _L (1 + z- 2 r + _L S (2.32)
dg kg kg*K ' kg* (1 + z )
where z = (z — z0)/g*H = (g — g0)/g*, z0 = H/2 is the middle position of the mixing layer.
Obviously, that the first term in the right (2.32) part corresponds to the logarithmic profile in the inner layer. Therefore the combined mean velocity profile can be written as follows
-2
U = u+n (z+ ) - — (Arsh( z ) + Arsh(z0))+^-— (arctan z + arctan z0) (2.33)
k k
where u+n (z+)is the mean velocity profile in the inner layer, which is given
by equation (2.27), Arsh(z) = ln(z W1 + z2), z0 = 1/2^*.
Using the asymptotic formula Arsh(z) = ln(2z) for z ® ¥ one can derive from (2.33):
Uо+= со - -1 k
/
ln-----------+ Arsh(z0)
ç* Re* v 07
+
л/1
+ z,
k
P _
—+ arctan z0 2 0
where Re* = Hut / v is the dynamic Reynolds number. Then, subtracting the equation (2.33) from both parts of this expression we have got the mean velocity defect low in the upper layer
ln(2v / V*) Arsh(z) yfl-
+ z,
P _
— - arctan z
(2.34)
Thus, the defect low (2.34) depends on the parameter ç*, which can precisely be determined, for instance, in case of turbulent flow in a flat channel from the additional condition u+ (1) = 0. However in the boundary layer this condition seems artificial as well as concept of the external boundary. It can be understood if the mean velocity defect low (2.34) to compare with the experimental data - see Figure 2.4, c. Figure 2.4,c shows that the computed profile (solid line 1) comes to zero when the parameter z / H > 3, whereas the experimental data concentrates near zero for z / H = 1. It should be noted, that the parameter of the profile 1 is calculated by minimisation of root-mean-square deviation of experimental points from the computational curve. That is reached for ç* = 0.4 .
To obtain the best correlation with the experimental data, the first condition (2.31) has been changed and. Then the mean velocity defect low can be rewritten as
ln(2ç / ç*) Arsh(z) ^д/l
+ z,
P _
— - arctan z 2
(2.35)
k kk where e is the parameter.
The mean velocity defect low computed on (2.35) (the solid line 1) and the experimental data by Nagano et al (1992) are shown in Figure 2.4,d. As it has been established ç* = 0.27; e = 0.8. After this correction the mean velocity defect low in the form (2.35) is in a good agreement with the experimental data - see Figure 2.4.d.
It is necessary to take into account the contribution of the upper layer universal function gradient to the mean velocity gradient on the wall, because the first condition (2.31) is broken. For this purpose the expression (2.33) can be modified as follows
(z+ ) - — (Arsh(z) + Arsh(z0 )) + ^ 0 (arctan z + arctan z0 )
kk
Differentiating both parts of this expression and calculating the derivative near the wall, we have
du_ = dul-------------------------------------------------1 e = 1 (2.36)
dz dz kç* Re* -^1 + Z02
Therefore the contribution of the upper layer universal function gradient to the mean velocity gradient on the smooth wall decreases with the growth of the dynamic Reynolds number. Using (2.36) the mean velocity profile in the turbulent boundary layer can be written as
u + = e0u+(z + ) —-(Arsh(z) + Arsh(z0)) + —--— (arctanz + arctanz0) (2.37)
k k
where e0 = 1 + (1 - e)/ kg* Re* yfl + zf » 1 + 0.9/ Re*.
The constants of the theory of turbulence should be re-normalized due to
(2.36) as follows
w+ = e0kR* exp(-10) = 0.14e0, l += R* / w+ = 8.71/e0 (2.38)
It is obvious, that the uncertainty of constants in the turbulent boundary layer on a flat plate is stipulated by the fact, that this flow is two-dimensional one, because it develops from the laminar boundary layer through the transition layer to the developed turbulent layer, down to the verge of separation (see chapter 6).
2.4.4 Streamwise turbulent intensity profiles
The variations of the velocity gradient around the mean value given by the first equation (2.27) are the production of two terms. One of them depends on the variations of the parameter A, which can estimated as
5A = 5(1 - eh / kl +) = -5w0+ eh / kRt* = -5w0 + / w+ ; another one depends on the fluctuating part of the wall shear stress tw or boundary condition on the wall -see Figure 2.5. Thus, in common case the velocity gradient can be written as
du + Ae"1 + (E - A)e - (2 39)
dz
1 + (z + / l +)2 д/1 + (z + /1 + )2
Boundary condition on the wall: z + ® 0: du + /dz + ® E where E = 1 for the mean flow.
Therefore the fluctuating streamwise velocity gradient can be written as
ddu + 5Ae-1 5A(5E / 5A - 1)e(2 40)
+ r, , + . „ 2 (2.40)
dz 1 +
(Z + / 1 +)2 д/1 + (z+ / 1 0 + )2
Let us suggest that for the considered turbulent flow 5E / 5A = const, then we have:
(du+)2 = (SA)2 f 2(z+, 8E / SA) (2.41)
Here f (z +,5E / 5A) is the solution of the value boundary problem:
df = e-1 + (5E / 5A - 1)e-1 (242)
dz
where f (0) = f (50) = 0 in case of the turbulent boundary layer on the flat plate, and f (0) = 0, f (50) = f0 for the turbulent flow in the two-dimensional channel.
dU/dz
Figure. 2.5: Frequency diagram of the wall shear stress in the turbulent boundary layer [62]. Mean velocity gradient dU / dz = 2000 c_1
One constant of the mean-squared value of the fluctuating streamwise velocity profile can be found out from the boundary value problem (2.42) solution and another one can be determined from the experiment. The predicted stream-
wise turbulent intensity profile, 8U + = ^{(ôA)2) f ( z+ )and the profile by Kuroda
et al [66] computed by the direct numerical simulation of the turbulent flow in the two-dimensional channel for Re = 2980, are shown in Figure 2.6,a. The constants determined for this case are given by ^((ôA)2) = 0925, ôE/ ôA = 0444. Figure
2.6,b illustrates the comparison of the computed streamwise turbulent intensity profile (solid line) with the experimental data by Nagano et al [68] (the line with symbols). In this case the best correlation of the computed and experimental data is for *J((ôA)2) = 0.735, ôE / ôA = 0.56.
Estimating the contribution of the outer region of the turbulent boundary layer to the streamwise turbulent intensity, one can assume, that the velocity fluctuations in the upper layer depend on the fluctuations of the velocity gradient on the wall. Then the generalized form of the expression (2.37) can be written as
u+ = e0u+ (z+) -—A)e°(Arsh(z) + Arsh(z0))+ (---) ^ ~° (arctanz + arctanz0) (2.43)
k k
where u+n(z+) is given by the equation (2.39). Varying parameters A, E in this profile, we have
Su+ = e08u+n (z + ) - —dAe° (Arsh(z) + Arshz0)) + )-) ^ ~0 (arctanz + arctanz0) (2.44)
k k
where 8u+n (z+) is calculated from (2.40). Hence, in this case it is possible to
present the streamwise turbulent intensity profile as follows 8U+ =^ ((8A)2) f (z),
where
f(z) = e0f (z+) -(dE/ dA—(Arsh(z) + Arsh(z0)) +)/ 8A—)g^ ~0 (arctanz + arctanz0) (2.45)
k k
SIT
A, a)
/ Л
\
/ N
/
(4 Ч- \
СГ
\и
J
10
10
Figure 2.6: a) Comparison of computed streamwise turbulent intensity profile with DNS data by Kuroda et al [66], b) experimental data by Nagano et al [67, 68]
The streamwise turbulent intensity profile computed on equation (2.45) and the experimental data by Nagano et al [68] are shown in Figure 2.7,a. Comparing the computed profiles in Figures 2.6, b and 2.7,a one can conclude that the boundary conditions on the external border for the turbulent intensity profile is differ for this two models. The profile computed on the model (2.45) has zero value together with the first derivative and then at z > H saves zero value, whereas the profile (2.42) approaching to zero has the nonzero decline.
The parameters of the profile (2.45) are given by
-J((8A)2) = 0.7, 8E / 8A = 0.568, thus its practically have the same values as for the
profile shown in Figure 2.6,b. The fluctuation of the velocity profile declination on the wall are connected with the fluctuation of the shear stress, 8E = 8tW / tW. This value can be defined from the data shown in Figures 2.6, 2.7. For the turbulent flow in the flat channel this value is given by 8tW / tW = 8A(8E / 8A) = 0.41.
The same value approximately we have for the turbulent boundary layer on the flat plate.
SU+ a) u+ b)
ЗУ
f\ 20
J V 10 /
n у
1 1(1 1(1(1 Z+ 0.1 1 10 100 z+
Fig. 2.7: a) The streamwise turbulent intensity profile in the turbulent boundary layer calculated on (2.45) - the solid line, and the experimental data [68]; b) the mean velocity profile in the turbulent boundary layer, calculated on
(2.37) for g* = 0.27; e = 0.79, and the experimental data [68] obtained in the crosssection x = 0.525 m
2.4.5. Normal and transversal turbulent intensity profiles
The flow in the turbulent boundary layer is not characterised by significant value of the normal and transversal (parallel to the wall) mean velocity components. Therefore the turbulent fluctuation of velocity in the specified directions can be considered as the random deviations of flow velocity vector from the mean value. The main parameter of these deviations is direct proportional to the dimensionless pressure normal gradient on the wall. Assuming that in the equation system (2.16) a = p /2, let us rewrite the last two equations as follows
W ■ < % - T+T (2.46)
where v+ = aw+ is the turbulent velocity scale characterised the flow velocity pulsation in the normal and transversal directions. Suggesting that Rt = R*, then the turbulent velocity scale can be written through the characteristic length scale, thus v+ = aR* /1+. At fixed turbulent theory parameters Rt = Rf, 1+ = 1+ the solutions of equations (2.46) depend on the parameter a, which is proportional to the dimensionless normal to the wall gradient of pressure, a = 1+Px (0)/2 pu*2.
Integrating equations (2.46) the velocity components in the inner layer can be written as follows
v + /v+=J (e 1 ~e,20)^ + e-J0 arctanX (2.47)
0 j i + x2
/+ X (e -1 - e -, e -1’ , x 2)
/V0 -J 1 + x 2 + 2 ln(1 + X )
-+x
It is necessary to add the compensatory functions, which depend on the upper layer variable z = (z — z0)/ g*H, to these expressions, thus we have
f (e ~1 — e ~10 )df
v + /v+ = 1--------~2--+ e ~10 arctan f + e, (arctan z + arctan z0)
0 1 + f
\ (e~1 — e~10)fdf e~h 2 e 1 + z2
' ' '
, + f (e - e )XdX , e , „ .£=24.^,
/ v0 — I + ln(1 + X ) + ln
0 I 1 + X2 2 v ъ 7 2
+ X2 ' 2 7 ' 2 “1 + z02
where parameter e1 is determined from the boundary conditions:
lim ® v + = v + ; lim ® w + = w+
Z®¥ ¥ 5 Z®¥ ¥
Therefore we have e1 = -e-Io and two additional equations for calculating of the parameters of this problem:
¥(eI0-1 - 1)dX . z v 1e0 (2 ,0)
J--------7—X2--- arctan zo =------— (2-48)
t 1 + X v0
f(+ _ln(1 + z 02) +ln-J0 1 + +2 2 K 0) r v0
If all parameters of the inner and outer layers are fixed, i.e.,
Rt = R*,w+ » 0.14,1+ = 1+, then the normal and transversal velocity fluctuations depend only on the parameter a = 1+ p (0)/2pu*2. In this case the turbulent intensity of the transversal velocity can be written as follows:
S v + = da[v + (z+) - v +ut(z)] (2.49)
where v +n is the function of the inner layer calculated according to the first
equation (2.47) at v+ = 1. The function of the outer layer is determined as follows
v +ut = e1 (arctan z + arctan z0).
The transversal turbulent intensity profile computed on (2.49) is shown in the left part of Figure 2.8 together with DNS data by Kuroda et al (1989). The
normal pressure gradient fluctuation parameter has been estimated from these data as 8a »1.1 / w+ » 7.87 and consequently v+ = w0+ 8a » 1.1.
v +
0.8 0.4
0^ ■ ■ ,
0.1 1 10 100 I+ 0.1 1 10 100 Z +
Figure 2.8: Turbulent intensity profiles of transversal velocity components in the turbulent boundary layer in the two-dimensional channel computed on (1.49) - (1.50) - solid lines, and DNS data by Kuroda et al [66]
The similar expression for the turbulent intensity of the normal velocity component is given by
8w + = w+&[< (z+) - w+out (Z )] (2.50)
where the function of the inner layer can be calculated on second equation (2.47) at v+ = 1. The outer layer function is determined in this case as follows
w+ut = eh [ln(1 + z2) - ln(1 + z02) + arctan z + arctan z0 ].
The normal turbulent intensity profile calculated on (2.50) is shown in the right part of Figure 2.8 with DNS data by Kuroda et al [66]. As it has been estimated the turbulent normal velocity scale v+ » 0.75, i.e., it is not equal to the transversal velocity scale. This is obviously connected to the influence of dynamic roughness parameters, because in common case a ^ p / 2. However it is difficult to calculate this influence, so far as we have the non-linear problem.
2.5. Dynamic roughness surface equation
The dynamic roughness surface is connected with the coherent structures in the turbulent boundary layer reported by Kline et al [70]. The review of the coherent structure problem, including the evaluation of the basic structure scales has been submitted by Cantwell [52]. To derive the dynamic roughness surface equation, the common method developed by Trunev et al [71-73] can be used. This method was proposed for the simulation of rough surfaces which are formed on the rigid body in the impingement erosion process, including the sputtering of a rigid body by ionic bombardment.
The Reynolds number calculated on the dynamic roughness parameters is given by
=X+w+ = hth
Rt = X0 w0 =
0 0 /7*2 y*2.
V (hx + hy )
This statistical relationship between the dynamic roughness parameters corresponds to the special type of the Nave-Stocks equation solution transformation
(2.4). Let’s consider small variation of the second velocity scale around the mean value w+ = kR*exp[-/0(R*)]» 0.14. In this case A * 0 in the first equation (2.16) and this parameter can be written as follows
A = Aj cosa = w0"a cosa
where a = arctan(hy / hx )is the parameter characterising the dynamic roughness structure. If A = 0, then a = p/2. This case corresponds to the special type of the dynamic roughness composed of furrows elongated along the streamlines of the mean flow. If a * p /2, then tan a = hy / hx » X x/ X y. It can be considered as
the relationship between the scales of the dynamic roughness elements in the X-and Y-directions. In this case a fluctuates around the mean value, a = p/2, and produces the velocity and pressure fluctuation. The velocity gradient on the wall fluctuates with a as follows
du + / dz + = w+a cosa + A2 sin a
For a = p /2 we have du + / dz + = A2 = w 0/ w0*. Therefore in this case the velocity gradient on the wall is given by
du + / dz + = (w 0/w0)(w+a cosa + sin a) (2.51)
as
The pressure gradient on the wall also depends on the second velocity scale
(0) = (2pu*2 /X+)aw+ (2.52)
The velocity components near the smooth wall for z ® 0 can be written as follows (see Kutateladze [62])
u @ K1 z, v @ K2 z, w @ K3z2,
where K = K (x, y, t) are the viscous sublayer functions. Substituting the velocity approximation formulas in the x-component of momentum equation and supposing that p = p(h) one can derive
Э K1 (0)
z——1 - z . „ cosa - zn
Эл
pl2
^Э 2K1 Э 2K1 Л
------1 н-------1
Эх2 Эу2
: z 2о(1)
where o(1) is the restricted value for z = 0. Using equations (2.51)-(2.52) finally we have the dynamic roughness surface equation in case of the steady turbulent boundary layer over a smooth surface
ЭК 2u* w+ a
cosa = n
Эх2 Эу2
(2.53)
dt v (X+)3
\ /
where K = (w0/ w0)(w+acosa + sin a).
The steady turbulent flow dynamic roughness is realised for w 0 = w *. In this case equation (2.53) can be transformed into the quasi-linear differential equation
П Ka
^Э2a Э2a^ Эх2 + Эу2
+ n Kaa (ax + ^ )=■
2u2w0a Эa
-----^cosa+Ka —
nl+ a Эt
(2.54)
where Ka = -w+a sin a + cosa, K„„ = -w+a cosa - sin a .
The point in which Ka = 0 is the singular point of the equation (2.54). In this point tana* = —(w+a)-1. As it has been estimated in the numerical experiments w + a » -0.127 for the mean flow, therefore a* » p/2 - w+|a|. For the stationary case,
i.e. for a » p / 2, the equation (2.54) can be written in the quasi-elliptical form
(1 - a)
Эх2 Эу2
V У
were a = (a - p / 2)/ \w +a|, k02 = 2u*2 / n21
{cx2x + a2
)-k0
a,
(2.55)
Using the function a = a( x, y)one can calculate the dynamic roughness surface parameters as follows
w+ a
ahy = o, ht = wUh2x + ку » w¡hy
(2.56)
In the special case when a << 1, equation (2.55) has the periodical solution
a = a 0cos()
where a 0 is the amplitude, ky = V2( A+ )-32 w* / v is the wave number in the y-
direction. Therefore the transversal length scale of coherent structures can be estimated as
1 * = 2p / ky
+ 3/2n / u* » 114n / u*
The predicted length scale is in a good agreement with the experimental value, Xy» 100v /u*, obtained by Kline et al. [70]. This type of coherent structures corresponds to the furrows considered above.
In the special case when |1 - a l<< 1, the periodical solution of the equation (2.55) is given by
a = 1 + a 0cos(kyy)
where ky = V2( X+)-32 u*/ a0v . Therefore in this case X y depends on the amplitude a0. The periodical solution of the first equation (2.56) can be written as
h(x, y,t) = h(kxx - kyyt - b(kyyt^ (2.57)
b (s) = ^ a n | cosn sds
n=1
1 + 7 *
kx = |w 0 a|ky, yt = y + w 0t .
The transversal phase velocity of the dynamic roughness surface disturbances can be determined as cy = w0*. The dynamic roughness length scale X x
depends on the amplitude a0 as follows
X x = X *ya0 /|w + a| » 898a0v / u*.
For a 0»1 the estimated streamwise length scale of coherent structures
agrees with the experimental value Xx» 1000v /u* obtained by Blackwelder &
Eckelmann [74], and discussed by Cantwell [52].
In this paper the problems of non-linear theory of turbulent boundary layer have been studied. The algorithm of numerical solution of the problem has been considered. Equation is deduced, connecting constants of non-linear theory.
A fundamental parameter of the turbulent boundary layer length has been determined, which coincides with the position of velocity peak of turbulence energy generation according to Klebanoff [49] and Laufer’s [50] data. It is shown that the profile of an average velocity in the boundary layer can be described satisfactory, using only one constant. The Karman constant can be used for this purpose. The second constant of the logarithmic profile can be estimated within this theory. Velocity profile, calculated according to the model suggested, conforms well to the data of direct numerical modelling, to experimental data and models of other authors. Results of velocity intensity pulsation modelling have been presented, as well as their compliance with the results of direct numerical modelling and experimental data. A model of dynamic roughness in turbulent boundary layer has been suggested. It has also been shown that in a stationary case there are two types of periodic solutions. One solution corresponds to dynamic roughness in a kind of furrows, stretched along the main flow. The vis-
cous flow over the structures is a physical mechanism of formation of logarithmic profile of velocity. The second solution corresponds to the perturbations of limited amplitude which have a limited length in the direction of the mean flow. It is shown that the parameters of dynamic roughness, having been calculated on the base of this model, coincide with the data of experiments.
(To be continued)
References
[1] Liepmann, H.W., The Rise and Fall of Ideas in Turbulence, American Scientist, 67, pp. 221-228, 1979.
[2] Amirkhanov, M.M, Lukashina, N.S. & Trunev, A. P., Natural recreation resources, state of environment and economical and legal status of coastal resorts, Publishing House "Economics", Moscow, 207 p., 1997 (in Russian).
[3] Marchuk, G. I., Mathematical Modelling in the Environmental Problem, "Nauka", Moscow, 1982 (in Russian).
[4] Borrell P.M., Borrell P., Cvitas T. & Seiler W., Transport and transformation of pollutants in the troposphere. Proc. EUROTRAC Symp., SPB Academic Publishing, Hague, 1994.
[5] Jaecker-Voirol A., Lipphardt M., Martin B., Quandalle, Ph., Salles, J., Carissimo, B., Dupont, E., Musson-Genon, L., Riboud, P.M., Aumont, B., Bergametti, G., Bey, I., Toupance, G., A 3D regional scale photochemical air quality model - application to a 3 day summertime episode over Paris, Air Pollution IV. Monitoring, Simulation and Control, eds. B. Caussade, H. Power & C.A. Brebbia, Comp. Mech. Pub., Southampton, pp. 175-194, 1996.
[6] Borrego, C., Coutinho, M., Carvalho, A.C.& Lemos, S., A modelling package for air quality management in Lisbon, Air Pollution V. Modelling, Monitoring and Management, eds. H. Power, T. Tirabassi & C.A. Brebbia, CMP, Southampton-Boston, pp. 35-44, 1997.
[7] Bozo, L.& Baranka, G., Air quality modelling over Budapest, Air Pollution IV. Monitoring, Simulation and Control, eds. B. Caussade, H. Power & C.A. Brebbia, Comp. Mech. Pub., Southampton, pp. 31-36, 1996.
[8] Marchuk, G.I. & Aloyan, A.E., Global Admixture Transport in the Atmosphere, Proc. Rus. Acad. Sci., Phys. Atmosphere and Ocean, 31, pp. 597-606, 1995.
[9] Moussiopoulos N., Air pollution models as tools to integrate scientific results in environmental policy, Air Pollution III, Vol.1. Theory and Simulation, eds. H. Power, N. Moussiopoulos & C.A. Brebbia, Comp. Mech. Publ., Southampton, pp. 11-18, 1995.
[10] Pekar M., Regional models LPMOD and ASIMD. Algorithms, parametrization and results of application to Pb and Cd in Europe scale for 1990, EMEP/MSC-E Report 9/96, Aug, 78 p., 1996.
[11] Carruthers, D.J, Edmunds, H.A., McHugh, C.A., Riches, P.J. & Singles, R.J., ADMS Urban - an integrated air quality modelling system for local government, Air Pollution V. Modelling, Monitoring and Management, eds. H. Power, T. Tirabassi & C.A. Brebbia, CMP, SouthamptonBoston, pp. 45-58, 1997.
[12] Ni Riain, C., Fisher, B., Martin, C. J. & Littler J., Flow field and Pollution Dispersion in a Central London Street, Proc. of the 1st Int. Conf on Urban Air Quality: Monitoring and Modelling, ed. R. S. Sokhi, Kluwer Academic Publishers, pp. 299-314, 1998.
[13] Lukashina, N.S. & Trunev, A. P., Principles of Recreation Ecology and Natural Economics, Russian Academy of Sciences, Sochi, 273 p., 1999 (in Russian).
[14] Lukashina, N.S., Amirkhanov, M.M, Anisimov, V.I. & Trunev, A.P., Tourism and environmental degradation in Sochi, Russia, Annals of Tourism Research, 23, pp. 654-665, 1996.
[15] Lenhart, L. & Friedrich, R. European emission data with high temporal and spatial resolution, Air Pollution III Vol.2: Air Pollution Engineering and management, eds. H. Power, N. Moussiopoulos & C.A. Brebbia. Comp. Mech. Pub., Southampton, pp.285-292, 1995.
[16] Oke, T.R., Street design and urban canopy layer climate, Energy and Buildings, 11, pp. 103-111, 1988.
[17] Zilitinkevich, S. Non-local turbulent transport: pollution dispersion aspects of coherent structure of convective flows, Air Pollution III, Vol.1. Air Pollution Theory and Simulation, eds. H. Power, N. Moussiopoulos & C.A.Brebbia, Comp. Mech. Publ., Southampton, pp.-53-60, 1995.
[18] Arya, S. P., Introduction to Micrometeorology, Academic Press, San Diego, 307 p., 1988.
[19] Stull, R. B, An Introduction to Boundary Layer Meteorology, Kluwer Academic Publishers, Dordrecht, 666 p., 1988.
[20] Kaimal, J. C.& Finnigan, J. J., Atmospheric Boundary Layer Flows: Their Structure andMeas-urements, Oxford University Press, 289 p., 1994.
[21] Monin, A.S. & Obukhov, A.M., Basic Laws of Turbulent Mixing in the Atmospheric surface layer, Trudy Geofiz. Inst. Akad. Nauk SSSR 24 (151), pp. 163-187, 1954.
[22] Monin, A. S., The Atmospheric Boundary Layer, Ann. Rev. Fluid Mech., 22, 1970.
[23] Businger, J.A., Wyangaard, J.C., Izumi, Y. & Bradley, E.F. Flux Profile Relationships in the Atmospheric Surface Layer, J. Atmos. Sciences, 28. pp. 181-189, 1971.
[24] Businger, J. A., A Note on the Businger-Dyer Profile, Boundary-Layer Meteorol., 42, pp. 145— 151, 1988.
[25] Yaglom, A.M., Data on Turbulence Characteristics in the Atmospheric Surface Layer, Izv. Acad. Sci. USSR, Phys. Atmosphere and Ocean, 10, pp. 566-586, 1974.
[26] Dyer, A. J., A Review of Flux-Profile Relationships, Boundary-Layer Meteorol., 7, pp. 363-372, 1974.
[27] Van Ulden, A. & Holtslag, A. A. M., Estimation of Atmospheric Boundary Layer Parameters for Diffusion Applications, J. Clim. Appl. Meteorol,. 24, pp. 1196-1207, 1985.
[28] Hanjalic, K.& Launder, B. E., A Reynolds Stress Model of Turbulence and its Application to Thin Shear Flows, J. Fluid Mech, 52, pp. 609-638, 1972.
[29] Rodi, W. Calculation of Stably Stratified Shear-layer Flows with a Buoyancy-extended k - e Turbulence Model, Turbulence and Diffusion in Stable Environments, ed. J. C. R. Hunt, Clarendon Press, Oxford, pp. 111-143, 1985.
[30] Mellor, G. L. & Yamada, T. A, Hierarchy of Turbulence Closure Models for Planetary Boundary Layers, J. Atmos. Sci., 31, pp. 1792-1806, 1974.
[31] Mellor, G. L. & Yamada, T., Development of a Turbulence Closure Model for Geophysical Fluid Problems, Rev. Geophys. Space Phys., 20, pp.851-875, 1982.
[32] Wyngaard, J. C. & Cote, O. R., The Evolution of a Convective Planetary Boundary Layer -a Higher-order-closure Model Study, Boundary-Layer Meteorol., 7, pp. 289-308, 1974.
[33] Zeman, O. & Lumley, J. L., Modelling Buoyancy Driven Mixed Layers, J. Atmos. Sci., 33, pp.1974-1988, 1976.
[34] Deardorff, J. W. & Willis, G. E., Further Results from a Laboratory Model of the Convective Boundary Layer, Boundary-Layer Meteorol, 32, pp. 205-236, 1985.
[35] Enger, L., A Higher Order Closure Model Applied to Dispersion in a Convective PBL, Atmos. Environ., 20, pp. 879-894, 1986.
[36] Holt, T. & Raman, S., A Review and Comparative Evaluation of Multilevel Boundary Layer Parameterisations for First Order and Turbulent Kinetic Energy Closure Schemes, Rev. Geophys. Space Phys., 26, pp. 761-780, 1988.
[37] Danilov, S.D., Koprov, B. M. & Sazonov, L. A., Atmospheric Boundary Layer and the Problem of Its Description (Review), Proc. Rus. Acad. Sci., Phys. Atmosphere and Ocean, 31, pp. 187204, 1995.
[38] Hurley, P. J., An Evaluation of Several Turbulence Schemes for the Prediction of Mean and Turbulent Fields in Complex Terrain, Boundary-Layer Meteorol., 83, pp. 43-73, 1997.
[39] Reynolds, O., On the dynamically theory of incompressible viscous fluids and the determination of the criterion, Philos. Trans. R. Soc. London, A 186, 123, 1895.
[40] Boussinesq, J., Theorie de l'ecoulement tourbillat, Mem. Pres. Acad. Sci., 23, p. 46, 1877.
[41] Prandtl, L., Bericht uber untersuchungen zur ausgebildeten turbulenz, Z. Angew. Math. Mech., 5, pp.136-139, 1925.
[42] Prandtl, L., Neuere Ergebnisse der Turbulenzforschung, VDI -Ztschr., 77, 5, p.105, 1933.
[43] Kolmogorov, A.N., The Equations of Turbulent Motion in an Incompressible Fluid, Izv. Acad. Sci. USSR, Phys., 6, pp. 56-58, 1942.
[44] Apsley, D. D. & Castro, I. P., A Limited-Length-Scale-Model for the Neutral and Stable-Stratified Atmospheric Boundary Layer, Boundary Layer Meteorol., 83, pp. 75-98, 1997.
[45] Trunev, A. P., Diffuse processed in turbulent boundary layer over rough surface, Air Pollution III, Vol.1. Theory and Simulation, eds. H. Power, N. Moussiopoulos & C.A. Brebbia, Comp. Mech. Publ., Southampton, pp. 69-76, 1995.
[46] Trunev, A. P., Similarity theory and model of turbulent dusty gas flow over large-scale roughness, Abstr. of Int. Conf On Urban Air Quality: Monitoring and Modelling, University of Hertfordshire, Institute of Physics, London, p. 3.8, 1996.
[47] Trunev, A. P., Similarity theory for turbulent flow over natural rough surface in pressure and temperature gradients, Air Pollution IV. Monitoring, Simulation and Control, eds. B. Caussade, H. Power & C.A. Brebbia, Comp. Mech. Pub., Southampton, pp. 275-286, 1996.
[48] Trunev, A. P., Similarity theory and model of diffusion in turbulent atmosphere at large scales, Air Pollution V. Modelling, Monitoring and Management, eds. H. Power, T. Tirabassi & C.A. Brebbia, CMP, Southampton-Boston, pp. 109-118, 1997.
[49] Klebanoff, P. S., Characteristics of turbulence in a boundary layer with zero pressure gradient, NACA Tech. Note, 3178, 1954.
[50] Laufer, J., The structure of turbulence in fully developed pipe flow, NACA Tech. Note, 2954, 1954.
[51] Cebeci, T. & Bradshaw, P., Physical and Computational Aspects of Convective Heat Transfer, Springer-Verlag, NY, 1984.
[52] Cantwell, Brian J., Organized motion in turbulent flow, Ann. Rev. Fluid Mech., 13, pp. 457-515, 1981.
[53] Kuroda, A., Direct Numerical Simulation of Couette-Poiseuille Flows, Dr. Eng. Thesis, the University of Tokyo, Tokyo, 1990.
[54] Coleman, G.N., Ferziger, J. R. & Spalart, P. R., A numerical study of the turbulent Ekman layer, J. Fluid Mech, 213. pp.313-348, 1990.
[55] Trunev, A. P. & Fomin, V. M., Continual model of impingement erosion, J. Applied Mech. Tech. Phys, 6, pp. 113-120, 1985.
[56] Trunev, A. P., Research of bodies erosion distraction in gas flows with admixture particles, Ph.D. Thesis, Inst. Theoretical and Appl. Mech., Novosibirsk, 1986.
[57] Nikolaevskii, V.N., The space averaging in the turbulence theory, Vortexes and Waves, ed. V.N. Nikolaevskii, Mir, Moscow, pp. 266-335, 1984 (in Russian).
[58] Landau, L.D. & Lifshitz, E. M., Hydrodynamics, 3rd ed., Nauka, Moscow, 1986 (in Russian).
[59] Pulliam, T. H. & Steger, J. L., Implicit Finite-Difference Simulations of three-dimensional Compressible Flow, AIAA Journal, 18, p. 159, 1980.
[60] Hirschel, E.H. & Kordulla, W., Shear Flow in Surface-Oriented Coordinates, Friedr. Vieweg & Sohn, Wiesbaden, 1986.
[61] Schlichting, H., Boundary Layer Theory, McGraw-Hill, NY, 1960.
[62] Kutateladze, S.S., The Wall Turbulence, Nauka, Novosibirsk, 1973 (in Russian).
[63] Hairer, E., Norsett, S.P. & Wanner, G., Solving Ordinary Differential Equations 1. Nonstiff Problems, Springer-Verlag, Berlin, 1987.
[64] 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, pp. 641-672, 1978.
[65] Van Driest, E.R., On turbulent flow near a wall, J. Aero. Sci., 23. p.1007, 1956.
[66] Kuroda, A., Kasagi, N. & Hirata, M., A Direct Numerical Simulation of the Fully Developed Turbulent Channel Flow, Proc. Int. Symp. on Computational Fluid Dynamics, Nagoya, pp. 1174-1179, 1989.
[67] Nagano, Y., Tagawa, M. & Tsuji, T., Effects of Adverse Pressure Gradients on Mean Flows and Turbulence Statistics in a Boundary Layer, Proc. 8th Symposium on Turbulent Shear Flows, 1992.
[68] 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.
[69] Smith, R.W., Effect of Reynolds Number on the Structure of Turbulent Boundary Layers, Ph.D. Thesis, Princeton University, Princeton, NJ, 1994.
[70] Kline, S.J , Reynolds, W.C., Schraub, F.A. & Runstadler P.W., The structure of turbulent boundary layers, J. Fluid Mech., 30, pp. 741-773, 1967.
[71] Kriklivy, V.V., Trunev, A.P. & Fomin, V.M., Investigation of two-phase flow in channel with damaging wall, J. Applied Mech. Tech. Phys., 1, pp. 82-87, 1985.
[72] Trunev, A. P. & Fomin, V.M., Surface instability during erosion in the gas-particles stream, J. Applied Mech. Tech. Phys., 3. pp. 78-84, 1986.
[73] Trunev, A. P., Evolution of the surface relief at sputtering by ionic bombardment, Interaction of nuclear particles with a rigid body, Moscow, Vol.1, Part 1, pp. 83-85, 1989.
[74] Blackwelder, R. F. & Eckelmann, H., Streamwise vortices associated with the bursting phenomena, J. Fluid Mech., 94, pp. 577-594, 1979.
