首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Flux-profile relationships based on surface-layer similarity theory are used to derive relationships between the Monin-Obukhov stability parameter = z/L and the bulk Richardson number Ri b . In contrast to previous studies, the roughness length for heat, z 0h ,is assumed unequal to the roughness length for momentum, z 0m .For the stable case, an analytic expression of in terms of Ri b can be derived and in the unstable case, the solution is obtained through a simple iterative process.Errors introduced from the simplification of z 0h = z 0m are evaluated and are shown to be very significant in most cases. Thus, this error in many practical applications may invalidate the intended solution.  相似文献   

2.
It is shown that for the purpose of trajectory simulation, the vertical velocityw L (t) of a fluid element, which is moving in a system (such as a forest canopy, or the unstably stratified atmospheric surface layer) whose turbulent velocity scale w is height-dependent, must be chosen from a frequency-distribution which is asymmetric aboutw L = 0. If the gradient w /z varies only slowly with height, correct trajectories may be obtained by adding a bias (where L is the length scale) to a fluctuating velocity chosen from a symmetric distribution with variance w 2(z).  相似文献   

3.
The standard deviation of temperature T is proposed as a temperature scale and as a velocity scale to describe the behaviour of turbulent flows in the Atmospheric Surface Layer (ASL), instead of * andu * of the Monin—Obukhov similarity theory, and ofT f andU f used for free convection stability conditions. On the basis of experimental evidence reported in the literature, it is shown that T T f andv * U f in the free convection region, and T * andv * U * in nearneutral and stable conditions. This implies that the proposed scales can be applied for all stabilities. Furthermore, a new length scale is proposed and its relation with Obukhov length is given. Also, a simple semi-empirical expression is presented with which T andv * can be evaluated in a rather simple way. Some examples of practical applications are given, e.g., a stability classification for unstable conditions.  相似文献   

4.
Summary This work deals with the Linke turbidity factor, based on total spectrum observations of the direct solar beam and aerosol turbidity parametersa a , , and based on observations in broad spectral bands. Diurnal and seasonal variations of these turbidity parameters were analyzed for the period 1975 to 1991.Annual variations of these parameters show low values in winter and high values in both spring and summer. The extinction coefficients decrease with increase of both wavelength and optical airmass. Trend analysis shows an increase in aerosol extinction coefficient below 0.63 m, and a slight decrease for longer wave-lengths.Linear regression relations are also constructed to estimate botha a and whenT L is available. The relations show thata a can be estimated with errors below 20%. The relation with the parameter, may give better results when it is estimated by assigning a fixed value of .Nomenclature AV Monthly and total average of the measured parameter - a Atmospheric optical thicknes - a a Aerosol optical thickness - a r Mean value of optical thickness of. Rayleigh atmosphere over all wavelengths - a o Ozone optical thickness - a o Ozone absorption coefficient - a w Water vapor optical thickness - COR Correlation coefficient of the linear relation in percentage - Ex1, Ex2, Ex3 Aerosol extinction coefficients in the bands .2–.53, .53–.63, .63–.695, respectively - I (o) Normal incident direct solar radiation, under clear sky condition - I o Extraterrestrial insolation at normal incidence - m r Relative (optical) air mass - NO Number of the observations used in either making the relation or the verification - RMSE Root mean square error of the linear relation - RMSE% Percentage value of the root mean square error relative to the average measured value AV - T L Linke turbidity factor - T Dry bulb temperature in °C - u o Ozone layer thickness, cm - Z Zenith angle - Ångström wavelength exponent - Ångström turbidity coefficient - Wavelength - Y The year number after 1975 With 5 Figures  相似文献   

5.
Equilibrium evaporation beneath a growing convective boundary layer   总被引：1，自引：1，他引：0  
Expressions for the equilibrium surface Bowen ratio ( s ) and equilibrium evaporation are derived for a growing convective boundary layer (CBL) in terms of the Bowen ratio at the top of the mixed layer i and the entrainment parameter A R . If AR is put equal to zero, the solution for s becomes-that previously obtained for the zero entrainment or closed box model. The Priestley-Taylor parameter is also calculated and plotted in terms ofA R and i . Realistic combinations of the atmospheric parameters give values of in the range 1.1 to 1.4.  相似文献   

6.
A numerical model of airflow above changes in surface roughness and thermal conditions is extended to include cases with stable thermal stratification within the internal boundary-layer. The model uses a mixing-length approach with empirical forms for M and H.Results are presented for some basic cases and an attempt is then made to compare results given by the model with the experimental results of Rider, Philip and Bradley. Tolerable agreement is achieved. The importance of roughness change and thermal stability effects in the diffusion of heat and moisture near a leading edge is emphasised.Notation A Refers to Taylor (1970) - B Businger-Dyer constant (= 16.0) in forms for M and H - C Constant in form for in stable case - c p Specific heat at constant pressure - E Scaled absolute humidity - g Acceleration due to gravity - H Upward vertical heat flux - H 0, H 1 Surface heat fluxes for x <0, x0 - H E Upward latent heat flux - k Von Kármán's constant (= 0.4) - K H K W Eddy transfer coefficients for heat and water vapour - L Monin-Obukhov length - L H Latent heat of evaporation for water - m Ratio of roughness lengths ( = z 1/z 0) - RPB Refers to Rider et al. (1964) - RL* Non-dimensional parameter (see Equations (9), (20a), (22a), (24a)) - R* Net radiation less ground heat flux (see Equations (15), (16)) - T Scaled temperature - T 1 Downstream scaled surface temperature - u 0 u 1(x) Surface friction velocities for x <0, x0 - U, W Horizontal and vertical mean velocities - x, z Horizontal and vertical co-ordinates - Z i Local roughness length - z 0, z i Roughness lengths for x < 0, x 0 - Temperature - 0, 1 Surface temperatures for x<0, x0 - E Non-dimensional absolute humidity gradient - H Non-dimensional temperature gradient of heat flux - M Non-dimensional wind shear - = M = H = E an assumption used in stable conditions - Air density - Absolute humidity - w Density of water - Kinematic shear stress - Logarithmic height scale (= ln(z+z 1)/z 1)  相似文献   

7.
Many applied dispersion models require the knowledge of boundary-layer parameters such as sensible heat flux,Q H , friction velocity,u *, and turbulent energy components, w and v . Formulas are suggested for calculating these parameters over a wide variety of types of ground surfaces, based on simple observations of wind speed near the ground and fractional cloud cover, and specification of constants such as roughness length, albedo, and soil moisture availability. Observations ofu *,Q H , w , and v during field experiments in St. Louis and Indianapolis are used to test the formulas for urban sites. Relative errors of about ±20% in the predictions are seen to occur whenu *,Q H , w , and v are large. However, when these quantities are small (e.g.,u * < 0.2 m/s), the errors in the predictions are as large as the mean value of the quantity itself.In addition, it is concluded from studies of available field data and theories that the magnitude of w is not well-known at elevations above about 100m during the late afternoon and night. Some simple parameterizations for w . are suggested that are consistent with the observed steady decrease in ground-level concentration in the afternoon and the sudden increase in concentration that can occur a few hours after sunset due to wind shears associated with a low-level jet, for continuous plumes emitted from moderate to tall stacks.  相似文献   

8.
During spring and autumn, many lakes in temperate latitudes experience intensive convective mixing in the vertical, which leads to almost isothermal conditions with depth. Thus the regime of turbulence appears to be similar with that characteristic of convective boundary layers in the atmosphere. In the present paper a simple analytical approach, based on boundary-layer theory, is applied to convective conditions in lakes. The aims of the paper are firstly to analyze in detail the temperature distribution during these periods, and secondly to investigate the current system, created by the horizontal temperature gradient and wind action. For these purposes, simple analytical solutions for the current velocities are derived under the assumption of depth-constant temperatures. The density-induced current velocities are shown to be small, in the order of a few mm/sec. The analytical model of wind-driven currents is compared with field data. The solution is in good qualitative agreement with observed current velocities under the condition that the wind field is steady for a relatively long time and that residual effects from former wind events are negligible.The effect of the current system on an approximately depth-constant temperature distribution is then checked by using the obtained current velocity fields in the heat transfer equation and deriving an analytical solution for the corrected temperature field. These temperature corrections are shown to be small, which indicates that it is reasonable to describe the temperature distribution with vertical isotherms.Notation T temperature - t time - x, y, z cartesian coordinates - molecular viscosity - h , v horizontal and vertical turbulent viscosity - K h ,K v horizontal and vertical turbulent conductivity - Q heat flux through the water surface - D depth - u, v, w average current velocity components inx, y andz directions - f Coriolis parameter - p pressure - density - g gravity acceleration - a constant in the freshwater state equation - h s deviation from the average water surface elevation - L *,H * length and depth scale - U *,W * horizontal and vertical velocity scale - T temperature difference scale - bottom slope - u * friction velocity at the water surface - von Karman constant - L Monin-Obukhov length scale - buoyancy parameter - l turbulence length scale - C 1,C 2,C 3 dimensionless constants in the expressions for the vertical turbulent viscosity - , dimensionless vertical coordinate and dimensionless local depth - angle between surface stress direction andx-axis - T bx ,T by bottom stress components - c bottom drag coefficient  相似文献   

9.
It is shown that the ratio of standard deviation of lateral velocity to the friction velocity, /u *, and therefore wind direction fluctuations, are sensitive to mesoscale terrain properties. Under neutral conditions, /u * is almost 40% larger in rolling terrain than over a horizontal surface. In the lee of a low mountain, the fluctuations may be 2.5 times as strong as over horizontal terrain. In contrast, vertical velocity fluctuations are little influenced by mesoscale terrain features.Now with Air Weather Service, Offutt AFB, Omaha, Nebraska.  相似文献   

10.
Surface-layer features with different prevailing wind directions for two distinct seasons (Southwest Monsoon and Northeast Monsoon) on the west coast of India are studied using data obtained from tower-based sensors at a site located about 500 m from the coast. Only daytime runs have been used for the present analysis. The surface boundary-layer fluxes have been estimated using the eddy correlation method. The surface roughnessz 0 obtained using the stability-corrected wind profiles (Paulson, 1970) has been found to be low for the Southwest monsson season. For the other season,z 0 is relatively high. The drag coefficientC D varies with height in the NE monsoon season but not in the season with lowz 0. This aspect is reflected in the wind profiles for the two seasons and is discussed in detail. The scaling behaviour of friction velocityu * and the turbulence intensity of longitudinal, lateral and vertical winds u, v and w, respectively) are further examined to study their dependence on fetch. Our study shows that for the non-dimensional case, u/u* and v/u* do not show any surface roughness dependence in either season. On the other hand, for w/u* for the season with lowz 0, the values are seen to agree well with that of Panofskyet al. (1977) for homogeneous terrain whereas for the other season with highz 0, the results seem to conform more to the values observed by Smedman and Högström (1983) for coastal terrain. The results are discussed in the light of observations by other investigators.  相似文献   

11.
The standard deviation of vertical two-point longitudinal velocity fluctuation differences is analyzed experimentally with eleven sets of turbulence measurements obtained at the NASA 150-m ground-winds tower site at Cape Kennedy, Florida. It is concluded that /u *0 is proportional to (fz/u *0)0.22, where the coefficient of proportionality is a function of fz/u *0 and u *0/fL 0. The quantities f and L0 denote the Coriolis parameter and the surface Monin-Obukhov stability length, respectively; u *0 is the surface friction velocity; z is the vertical distance between the two points over which the velocity difference is calculated; and zz is the mean height of the mid-point of the interval z above natural grade. The results of the analysis are valid for 20<-u *0/fL 0<2000.  相似文献   

12.
A theory is offered for the drag and heat transfer relations in the statistically steady, horizontally homogeneous, diabatic, barotropic planetary boundary layer. The boundary layer is divided into three regionsR 1,R 2, andR 3, in which the heights are of the order of magnitude ofz 0,L, andh, respectively, wherez 0 is the roughness length for either momentum or temperature,L is the Obukhov length, andh is the height of the planetary boundary layer. A matching procedure is used in the overlap zones of regionsR 1 andR 2 and of regionsR 2 andR 3, assuming thatz 0 L h. The analysis yields the three similarity functionsA(),B(), andC() of the stability parameter, = u */fL, where is von Kármán's constant,u * is the friction velocity at the ground andf is the Coriolis parameter. The results are in agreement with those previously found by Zilitinkevich (1975) for the unstable case, and differ from his results only by the addition of a universal constant for the stable case. Some recent data from atmospheric measurements lend support to the theory and permit the approximate evaluation of universal constants.  相似文献   

13.
The Langevin equation is used to derive the Markov equation for the vertical velocity of a fluid particle moving in turbulent flow. It is shown that if the Eulerian velocity variance wE is not constant with height, there is an associated vertical pressure gradient which appears as a force-like term in the Markov equation. The correct form of the Markov equation is: w(t + t) = aw(t) + b wE + (1 – a)T L ( wE 2)/z, where w(t) is the vertical velocity at time t, a random number from a Gaussian distribution with zero mean and unit variance, T L the Lagrangian integral time scale for vertical velocity, a = exp(–t/T L), and b = (1 – a 2)1/2. This equation can be used for inhomogeneous turbulence in which the mean wind speed, wE and T L vary with height. A two-dimensional numerical simulation shows that when this equation is used, an initially uniform distribution of tracer remains uniform.  相似文献   

14.
The maintenance of an elevated inversion in steady flow above a cold, rotating surface is shown to be possible for a certain range of the buoyancy number bfV g, where b is the buoyant acceleration appropriate to the density deficiency of the fluid above the inversion, f is Coriolis parameter and V gis geostrophic velocity (so that fV gis also horizontal pressure gradient in kinematic units). The height of the inversion lid is determined by a balance of surface stress and buoyancy, in a way which may be deduced from laboratory experiments. With the aid of such empirical evidence a theory is constructed for the layer below the inversion lid. The cross-isobar angle of ground-level stress is found to increase with the buoyancy number, to a limiting value of 90, by which time the inversion descends to the ground. Under typical conditions, a temperature difference of order 10C is necessary to eliminate the possibility of an equilibrium, elevated inversion lid and reduce ground level wind stress to a vanishingly small value.Woods Hole Oceanographic Institution Contribution #3011On leave from the University of Waterloo, Ontario  相似文献   

15.
A Random Displacement Model (RDM) and a Langevin Equation Model (LEM) are used to simulate point releases in a complex flow around a building. The flow field is generated by a three-dimensional finite element model that uses the standardk- model to parameterize the turbulence. The RDM- and LEM-calculated concentration fields are compared, with particular emphasis on the structure in regions with high turbulence and/or recirculation. RDM and LEM results are similar qualitatively, but RDM tends to predict lower concentration levels. In part this is due to the higher early-time diffusion. However, the expected convergence at later times is prevented by the interaction of the diffusion with the strongly inhomogeneous mean flow.Notation a i coefficient in the Langevin equation - b ij coefficient in the Langevin equation - C 0 the universal constant associated with the Lagrangian structure function - H building height (22.5 m) - K eddy viscosity - K k eddy viscosity used in the definition of the off-diagonal Reynolds stresses - k turbulent kinetic energy - LEM Langevin Equation Model - p 1 local unit vector in thexy-plane, orthogonal tos - p 2 local unit vector, orthogonal to boths andp 1 - RDM Random Displacement Model - s local unit vector in the streamline direction - T local decorrelation time (Lagrangian time scale) - U magnitude of the local Eulerian mean wind velocity - u s total velocity in the streamline direction - u 1 velocity component in thexy-plane, orthogonal to the streamline direction - u 2 velocity component orthogonal to bothu s andu 1 - i mean Eulerian wind velocity - W i stochastic vector-valued Wiener process - x unit vector inx-direction - y unit vector iny-direction - z unit vector inz-direction - angle between thexy-plane and the mean wind streamline - angle between the projection in thexy-plane of the streamline and thex-axis - ij the Kronecker delta function - rate of turbulence dissipation - i/ga the part ofa i that contains mean wind and turbulence gradients - ij inverse of a Reynolds stress tensor component - ij shorthand for a quantity that defines a part of i/ga - i shorthand for a quantity that defines a part of i/ga - ij Reynolds stress tensor component  相似文献   

16.
We propose a simple model for esimating the average number of occurrences per unit time (c o) that a threshold concentration c o is exceeded. It is based on the joint probability density of the observed concentration c(t) and its time derivative (t) under the assumption that c(t) is a stationary time series; this assumption leads to the hypothesis that c(t) and (t) are statistically independent. Adopting plausible forms of the frequency distributions of c and , we apply the model to diffusion from an infinite area source and from an elevated point source, both in the neutral boundary layer, and obtain simple results for (c o) and the average duration of one excursion above c o as functions of c o, the mean and the standard deviation of the concentration, and surface-layer variables.The National Center for Atmospheric Research is sponsored by the National Science Foundation.  相似文献   

17.
The I-atom sensitised decomposition of ozone in air at 1 atm pressure and ambient temperature has been investigated. Iodine atoms were produced by photolysis of I2 using visible light or of CH3I using ultraviolet light. In both cases, the quantum yield for O3 decomposition was 1.25 (±0.11) per I atom. An important role is proposed for the self-reaction of IO radicals leading to higher oxides of iodine, IO+IO(+M)I2O2(+M)higher oxides, which predominated over the bimolecular reaction leading to regeneration of I atoms, IO+IO2I+O2, with k 2a/k 2b4. Simple computer modelling calculations indicate that reaction (2a) may be important in determining the fate of photolabile iodine species in the atmosphere. The consequences for the behaviour of radioiodine releases are also discussed.  相似文献   

18.
Summary A zonally averaged global energy balance model with feedback mechanisms was constructed to simulate (i) the poleward limits of ITCZ over the continent and over the ocean and (ii) a simple monsoon system as a result of differential heating between the continent and the ocean. Three numerical experiments were performed with lower boundary as (1) global continent, (2) global ocean and (3) continent-ocean, with freezing latitudes near the poles. Over the continent, midlatitude deserts were found and the ITCZ migrates 25° north and south with seasons. Over a global swamp ocean results do not show migration of ITCZ with time but once the ocean currents are introduced the ITCZ migrates 5° north and south with seasons. It was found that the seasonal migration of ITCZ strongly depends on the meridional distribution of the surface temperature. It was also found that continent influences the location of the oceanic ITCZ. In the tropics northward progression of quasi-periodic oscillations called events are found during the pre- and post-monsoon periods with a period of 8 to 15 days. This result is consistent with the observed quasi-periodic oscillations in the tropical region. Northward propagation of the surface temperature perturbation appears to cause changes in the sensible heat flux which in turn causes perturbations in vertical velocity and latent heat flux fields.List of Symbols vertical average - 0 zonal average - vertical mean of the zonal average - 0s zonal average at the surface - 0a zonal average at 500 mb level - latitude We now define the various symbols used in the model rate of atmospheric heating due to convective cloud formation (K/sec) - dp/dt (N/m2/sec) - density - potential temperature (K) - rate of rotation of the earth (rad/sec) - empirical constant - humidity mixing ratio - * saturated humidity mixing ratio - opacity of the atmosphere - 1,2 factors for downward and upward effective black body long wave radiation from the atmosphere - Stefan-Boltzmann constant - emissivity of the surface - D subsurface temperature (K) - a specific volume - 0xs ,0ys eastward and northward components of surface frictional stress - * vertical velocity at the top of the boundary layer (N/m2/sec) - P Thickness of the boundary layer (mb) - nondimensional function of pressure - P pressure - P a pressure of the model atmosphere (N/m2) - P s pressure at the surface (N/m2) - t time (sec) - U eastward wind speed (m/sec) - V northward wind speed (m/sec) - surface water availability - T absolute temperature (K) - heat addition due to water phase changes - g acceleration due to gravity (m2/sec) - a radius of the earth (m) - R gas constant for dry air (J/Kg/K) - C p specific heat of air at constant pressure (J/Kg/K) - k R/C p - L latent heat of condensation (J/Kg) - f coriolis parameter (rad/sec) - H s H 0s (1) +H 0s (2) +H 0s (3) +H 0s (4) +H 0s (5) (J/m2/Sec)=sum of the rates of vertical heat fluxes per unit surface area, directed toward the surface - H a H 0a (1) +H 0a (2) +H 0a (3) +H 0a (4) (J/m2/Sec)=sum of the rates of heat additions to the atmospheric column per unit horizontal area by all processes - H 0s (1) ,H 0a (1) heat flux due to short wave radiation - H 0s (2) ,H 0a (2) heat flux due to long wave radiation - H 0s (3) ,H 0a (3) heat flux due to small scale convection - H 0s (4) heat flux due to evaporation - H 0a (4) heat flux due to condensation - H 0s (5) heat flux due to subsurface conduction and convection - e * saturation vapor pressure - R solar constant (W/m2) - r a albedo of the atmosphere - r s albedo of the surface - b 2 empirical constant (J/m2/sec) - c 2 empirical constant (J/m2/sec) - e 2 nondimensional empirical constant - f 2 empirical constant (J/m2/sec) - factor proportional to the conductive capacity of the surface medium - a s constant used in Sellers model - b s positive constant of proportionality used in the Sellers model (kg m2/J/sec2) - K HT coefficient for eddy diffusivity of heat (m2/sec) - K HE exchange coefficient for water vapor (m2/sec) - h depth of the water column (m) - z height (m) - V 0ws meridional component of surface current (m/sec) - n cloud amount - G 0,n long wave radiation form the atmosphere for cloud amount n (W/m2) - B 0 long wave radiation from the surface (W/m2) - S 0,n short wave radiation from the atmosphere for cloud amount n (W/m2) - A n albedo factor for a cloud amount n - R f1 large scale rainfall (mm/day) - R f2 small scale rainfall (mm/day) With 22 Figures  相似文献   

19.
This paper is written to report observations of the structure of the atmospheric surface layer over a coastal industrialized equatorial area. The observations were recorded at Prai Industrial Park, Penang (5° 22′ N, 100° 23′ E) a relatively simple terrain area during the south-west monsoon season in the period of three months using slow response systems. The limitations of the instruments used and its effects on the results are discussed. Wind turbulence and temperature were measured on a 10 m tower and analyzed using eddy correlation method and Monin–Obukhov similarity relations to obtain the normalized standard deviation of longitudinal (σu/u), lateral (σv/u) and vertical wind velocity fluctuations (σw/u) with respect to stability parameter z/L. From the results of the analysis, we found that most of turbulence is generated by shear or mechanical force. It was found that the average neutral value of σu/u is 2.35, 1.98 for σv/u and 1.47 for σw/u with a significantly lower than the proportionality to the power of 1/3 during unstable atmospheric conditions, and thus do not obey Monin–Obukhov similarity theory. It was observed that σu/u and σv/u values increase linearly in the range of 0 < z/L < 2 and fairly well correlated while σw/u does not.  相似文献   

20.
Turbulence mechanisms at an agricultural site   总被引：8，自引：0，他引：8  
An extensive set of turbulence data from the 3- and 12-m heights taken over an agricultural site (Marsta, Sweden) are analyzed and compared with data from ideal sites.In unstable air, Monin-Obukhov similarity is found to be valid for the non-dimensional gradients of wind, m , temperature, h , and humidity, e , for (only a few data), for T /|T *|,/ E /|E *| and for the non-dimensionalized inertial subrange spectra of temperature and humidity. Where comparison is possible, the unstable data also agree with those found in the Kansas study, with one remarkable exception, the inertial subrange constant of the temperature spectrum, 1, being only 0.39, compared to the value 0.80 found at the Kansas site.On the stable side, most similarity predictions break down, with most of the data differing systematically from the corresponding Kansas results, the only exception being . The inertial subrange constants for temperature, 1, and for humidity, 1 are found to have the same values, 0.39 and 0.30, respectively, as they do on the unstable side. Remarkable similarity is found for the shape of the stable u- and - and e-spectra. In addition, this shape is found to be identical with that found in Kansas. The peak wavelength of the stable u-, and -spectra is found to be about four times larger than it is for the corresponding Kansas spectra. This is interpreted to be a result of the increased macro-roughness at the Marsta site as compared with that at the Kansas site. A possible explanation for the low 1-value is discussed, suggesting that 1 is not a universal constant, but instead dependent on the turbulent structure.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号