首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 380 毫秒
1.
A wind-tunnel experiment has been used to investigate momentum absorption by rough surfaces with sparse random and clustered distributions of roughness elements. An unusual (though longstanding) method was used to measure the boundary-layer depth δ and friction velocity u * and thence to infer the functional relationship z 0/h = f(λ) between the normalised roughness length z 0/ h and the roughness density λ (where z 0 is the roughness length and h the mean height of the roughness elements). The method for finding u * is based on fitting the velocity defect in the outer layer to a functional form for the dimensionless velocity-defect profile in a canonical zero-pressure-gradient boundary layer. For the conditions investigated here, involving boundary layers over sparse roughness with strong local heterogeneity, this velocity-defect-law method is found to be more robust than several alternative methods for finding u * (uw covariance, momentum integral and slope of the logarithmic velocity profile).The experimental results show that, (1) there is general agreement in the relationship z 0/h = f(λ) between the present experiment with random arrays and other wind-tunnel experiments with regular arrays; (2) the main effect of clustering is to increase the scatter in the z 0/h = f(λ) relationship, through increased local horizontal heterogeneity; (3) this scatter obscures any trend in the z 0/h = f(λ) relationship in response to clustering; and (4) the agreement between the body of wind-tunnel data (taken as a whole) and field data is good, though with scatter for which it is likely that a major contribution stems from local horizontal heterogeneity in the field.  相似文献   

2.
A numerical model of airflow in the lowest 50–100 m of the atmosphere above changes in surface roughness and temperature or heat flux has been developed based on boundary layer approximations, the Businger-Dyer hypotheses for the non-dimensional wind shear and heat flux and a mixing length hypothesis.Results have been obtained for several situations, in particular, airflow with neutral upstream conditions encountering a step change in surface temperature or heat flux with no roughness change. In these cases large increases in shear stress at the outer edge of the internal boundary layer are predicted. The case of unstable upstream flow encountering a step change to zero heat flux is also considered.Two situations that may be encountered near the shores of the Great Lakes are considered.Notation B Businger-Dyer constant (= 16.0) in form for M, H - c p Specific heat at constant pressure - g Acceleration due to gravity - H Upward vertical heat flux - H 0 , H 1 Surface heat fluxes for x < 0, x 0 - k von Kármán's constant ( = 0.4) - l Mixing length - L Monin-Obukhov length - L 0 Upstream value of L - m Ratio of roughness lengths (= z 1/z 0) - RL * Non-dimensional parameter, see Equations (20, 22 and 24) - RL 1 * Same as RL * but with z 1 scaling (= mRL *) - T Scaled temperature - T 0 (z) Upstream temperature profile - u 0, u 1(x) Surface friction velocities for x < 0, x 0 - U, W Horizontal and vertical mean velocities - U 0 (z) Upstream velocity profile - x, z Horizontal and vertical coordinates - z i Local roughness length  相似文献   

3.
A common parametrization over snow-covered surfaces that are undergoing saltation is that the aerodynamic roughness length for wind speed (z 0) scales as au*2/g{\alpha u_\ast^2/g}, where u * is the friction velocity, g is the acceleration of gravity, and α is an empirical constant. Data analyses seem to support this scaling: many published plots of z 0 measured over snow demonstrate proportionality to u*2{u_\ast^2 }. In fact, I show similar plots here that are based on two large eddy-covariance datasets: one collected over snow-covered Arctic sea ice; another collected over snow-covered Antarctic sea ice. But in these and in most such plots from the literature, the independent variable, u *, was used to compute z 0 in the first place; the plots thus suffer from fictitious correlation that causes z 0 to unavoidably increase with u * without any intervening physics. For these two datasets, when I plot z 0 against u * derived from a bulk flux algorithm—and thus minimize the fictitious correlation—z 0 is independent of u * in the drifting snow region, u * ≥ 0.30 ms−1. I conclude that the relation z0 = au*2/g{z_0 = \alpha u_\ast^2/g} when snow is drifting is a fallacy fostered by analyses that suffer from fictitious correlation.  相似文献   

4.
It is frequently observed in field experiments that the eddy covariance heat fluxes are systematically underestimated as compared to the available energy. The flux imbalance problem is investigated using the NCAR’s large-eddy simulation (LES) model imbedded with an online scheme to calculate Reynolds-averaged fluxes. A top–down and a bottom–up tracer are implemented into the LES model to quantify the influence of entrainment and bottom–up diffusion processes on flux imbalance. The results show that the flux imbalance follows a set of universal functions that capture the exponential decreasing dependence on u */w *, where u * and w * are friction velocity and the convective velocity scale, respectively, and an elliptic relationship to z/z i , where z i is the mixing-layer height. The source location in the boundary layer is an important factor controlling the imbalance magnitude and its horizontal and vertical distributions. The flux imbalance of heat and the bottom–up tracer is tightly related to turbulent coherent structures, whereas for the top–down diffusion, such relations are weak to nonexistent. Our results are broadly consistent with previous studies on the flux imbalance problem, suggesting that the published results are robust and are not artefacts of numerical schemes.  相似文献   

5.
Summary In this paper the results of an urban measurement campaign are presented. The experiment took place from July 1995 to February 1996 in Basel, Switzerland. A total of more than 2000 undisturbed 30-minute runs of simultaneous measurements of the fluctuations of the wind vector u′, v′, w′ and the sonic temperature θ s ′ at three different heights (z=36, 50 and 76 m a.g.l.) are analysed with respect to the integral statistics and their spectral behaviour. Estimates of the zero plane displacement height d calculated by the temperature variance method yield a value of 22 m for the two lower levels, which corresponds to 0.92 h (the mean height of the roughness elements). At all three measurement heights the dimensionless standard deviation σ w /u * is systematically smaller than the Monin-Obukhov similarity function for the inertial sublayer, however, deviations are smaller compared to other urban turbulence studies. The σθ/θ* values follow the inertial sublayer prediction very close for the two lowest levels, while at the uppermost level significant deviations are observed. Profiles of normalized velocity and temperature variances show a clear dependence on stability. The profile of friction velocity u * is similar to the profiles reported in other urban studies with a maximum around z/h=2.1. Spectral characteristics of the wind components in general show a clear dependence on stability and dimensionless measurement height z/h with a shift of the spectral peak to lower frequencies as thermal stability changes from stable to unstable conditions and as z/h decreases. Velocity spectra follow the −2/3 slope in the inertial subrange region and the ratios of spectral energy densities S w (f)/S u (f) approach the value of 4/3 required for local isotropy in the inertial subrange. Velocity spectra and spectral peaks fit best to the well established surface layer spectra from Kaimal et al. (1972) at the uppermost level at z/h=3.2. Received September 26, 1997 Revised February 15, 1998  相似文献   

6.
For 390 ten-minute samples of turbulent flux, made with a trivane above a lake, the vertical alignment is determined within 0.1 ° through azimuth-dependent averaging. One degree of instrumental misalignment is found to produce an average tilt error of 9 ± 4% for momentum flux, and 4 ± 2% for heat flux. The tilt error in the vertical momentum flux depends mainly ons u/u*, and cannot be much diminished with impunity by high-pass pre-filtering of the turbulence signals. The effects of rain on trivane measurements of vertical velocity are shown to be negligible at high wind speeds, and adaptable to correction in any case.The normalized vertical velocity variance,s w/u*, appears to be proportional to the square root ofz/L for unstable stratification. For a wind speed range of 2 to 15 m s–1, the eddy correlation stresses measured at 4- and 8-m heights can be reasonably well estimated by using a constant drag coefficientC d=1.3 X 10-3, while cup anemometer profile measurements give an overestimate of eddy stress at high wind speeds. A good stress estimate is also obtained from the elevation variance; it is suggested that trivane measurement of this variance might be made from a mobile platform, e.g., a moderately stabilized spar buoy.  相似文献   

7.
In November 1993 an airborne field study was performed in order to investigate the microphysical and radiative properties of cooling tower water clouds initiated by water vapour emissions and polluted by the exhaust from coal-fired power plants. The number-median diameter of the droplet size distributions of these artificial clouds was in the range of 13 μm. The concentration of smaller droplets (diameters dD < 10 μm) increased with height and horizontal distance from the cooling towers. Close to the cooling towers, bimodal spectra were found with a second mode at 19 μm. The liquid water content (LWC) ranged between 2 and 5 g/m3 and effective droplet radii (Re) between 6 and 9 μm were measured. LWC and Re decreased with altitude, whereas the droplet concentration (ND) remained approximately constant (about 2000 cm−3 ). An enrichment of interstitial aerosol particles with particle diameters (dp) smaller 0.2 μm compared to the power plant plume in the vicinity of the clouds was observed. Particle activation for dm > 0.3 μm. was evident, especially in cooling tower clouds further apart and separated from their sources. Furthermore, radiation measurements were performed, which revealed differences in the vertical profiles of downwelling solar and UV radiation flux densities inside the clouds.The effective droplet radius Re was parameterized in terms of LWC and ND using equations known from literature. The close agreement between measured and parameterized Re indicates a similar coupling of Re, LWC and ND as in natural clouds.By means of Mie calculations, volume scattering coefficients and asymmetry factors are derived for both the cloud droplets and the aerosol particles. For the cloud droplets, the optical parameters were described by parameterizations from the literature. The results show, that the link between radiative and microphysical properties of natural clouds is not changed by the extreme pollution of the artificial clouds.  相似文献   

8.
Data collected during the SHEBA and CASES-99 field programs are employed to examine the flux–gradient relationship for wind speed and temperature in the stably stratified boundary layer. The gradient-based and flux-based similarity functions are assessed in terms of the Richardson number Ri and the stability parameter z/Λ*, z being height and Λ* the local Obukhov length. The resulting functions are expressed in an analytical form, which is essentially unaffected by self-correlation, when thermal stratification is strong. Turbulence within the stably stratified boundary layer is classified into four regimes: “nearly-neutral” (0 < z/Λ* < 0.02), “weakly-stable” (0.02 < z/Λ* < 0.6), “very-stable” (0.6 < z/Λ* < 50), and “extremely-stable” (z/Λ* > 50). The flux-based similarity functions for gradients are constant in “nearly-neutral” conditions. In the “very-stable” regime, the dimensionless gradients are exponential, and proportional to (z/Λ*)3/5. The existence of scaling laws in “extremely-stable” conditions is doubtful. The Prandtl number Pr decreases from 0.9 in nearly-neutral conditions and to about 0.7 in the very-stable regime. The necessary condition for the presence of steady-state turbulence is Ri < 0.7.  相似文献   

9.
The Ekman boundary-layer model is extended analytically for a gradually varying eddy diffusivity K(z) ≥ 0, z ≥ 0. A solution for the Ekman layer is provided having similar structure to the constant-K case; that is, exponentially decaying sine functions for the two horizontal wind components. The analytical asymptotic solution compares well with its numerical counterpart for various K(z). The result can be useful in theoretical studies such as Ekman pumping, for efficient estimation of the Ekman layer profiles in various analyses with near-neutral stratifications, or for a rapid initialization of mesoscale models.  相似文献   

10.
With observational data collected and interpreted by Crane et al. (1977), the adequacy of the O'Brien polynomial to represent the exchange profile of heat and pollution in a convective boundary layer is examined and a refinement suggested. Also, it is shown that the height of the surface layer, h=0.04 z i , developed by Blackadar and Tennekes (1968) for a neutrally stratified boundary layer (with z z =0.25u */f) appears to be equally valid for the convective boundary layer where z i , defined as the top of the mixed layer, is used.  相似文献   

11.
The roughness length, z 0u , and displacement height, d 0u , characterise the resistance exerted by the roughness elements on turbulent flows and provide a conventional boundary condition for a wide range of turbulent-flow problems. Classical laboratory experiments and theories treat z 0u and d 0u as geometric parameters independent of the characteristics of the flow. In this paper, we demonstrate essential stability dependences—stronger for the roughness length (especially in stable stratification) and weaker but still pronounced for the displacement height. We develop a scaling-analysis model for these dependences and verify it against experimental data.  相似文献   

12.
Turbulence statistics, including higher order moments, in the surface layer over plant canopies were compared with those observed over several different surfaces, using a nondimensional height (z – d)/z 0: The values of (z – d)/z 0extend over a very wide range from 10 over plant canopies to 107 over the ocean. Several properties such as intensities of turbulence and skewness factors show a remarkable height-dependency in the air layer below (z – d)/z 0 = 102, which is supposed to be much influenced by the underlying surface. In that layer, some peculiar phenomena, such as a downward energy transport and positive flux of shear stress, are frequently observed.  相似文献   

13.
A wind tunnel investigation of the wind erosion of uranium mine-tailings material typical of a northern Ontario site has been carried out. The aim of the study was to measure the effects of various parameters, including mean and turbulent wind characteristics of the boundary layer and surface moisture content, upon the erosion process. The analysis of experimental data has yielded a mathematical model for predicting the net vertical mass fluxes. The results show that the dry vertical flux is proportional to u * 2.3and the wet flux to u * 5.0 Partical size analysis was also carried out.  相似文献   

14.
We examine the unsteady response of a neutral atmospheric boundary layer (ABL) of depth h and friction velocity u * when a uniform surface heat flux is applied abruptly or decreased rapidly over a time scale t<inf>θ</inf> less than about h /(10u *). Standard Monin–Obukhov (MO) relationships are used for the perturbed eddy viscosity profile in terms of the changes to the heat flux and mean shear. Analytical solutions for changes in temperature, mean wind and shear stress profile are obtained for the surface layer, when there are small changes in h /|LMO| over the time scale tMO~|L MO|/(10u*) (where L MO and t MO are the length and time scales, respectively). They show that a maximum in the wind speed profile occurs at the top of the thermal boundary layer for weak surface cooling, i.e. a wind jet, whereas there is a flattening of the profile and no marked maximum for weak surface heating. The modelled profiles are approximately the same as those obtained from the U.K. Met Office Unified Model when operating as a mesoscale model at 12-km horizontal resolution. The theoretical model is modified when strong surface heating is suddenly applied, resulting in a large change in h /|L MO| (>>1), over the time scale t MO. The eddy structure is predicted to change significantly and the addition of convective turbulence increases the shear turbulence at the ground. A low-level wind jet can form, with convective turbulence adding to the mean momentum of the flow. This was verified by our laboratory experiment and direct numerical simulations. Additionally, it is shown that the effects of Coriolis acceleration diminish (rather than as suggested in the literature, amplify) the formation of the wind jets in the situations considered here. Hence, only when the surface heat flux changes over time scales greater than 1/f (where f is the Coriolis parameter) does the ABL adjust monotonically between its equilibrium states. These results are also applicable to the ABL passing over spatially varying surface heat fluxes.  相似文献   

15.
A laboratory study of scalar diffusion in the convective boundary layer has found results that are consistent with a 1999 large-eddy simulation (LES) study by Jonker, Duynkerke and Cuijpers. For bottom-up and top-down scalars (introduced as ‘infinite’ area sources of passive tracer at the surface and inversion, respectively) the dominant length scale was found to be much larger than the length scale for density fluctuations, the latter being equal to the boundary-layer depth h. The variance of the normalized passive scalar grew continuously with time and its magnitude was about 3–5 times larger for the top-down case than for the bottom-up case. The vertical profiles of the normalized passive scalar variance were found to be approximately constant through the convective boundary layer (CBL) with a value of about 3–8c*2 for bottom-up and 10–50c*2 for top-down diffusion. Finally, there was some evidence of a minimum in the variance and dominant length scale for scalar flux ratios (top-down to bottom-up flux) close to −0.5. All these convection tank results confirm the LES results and support the hypothesis that there is a distinct difference in behaviour between the dynamic and passive variables in the CBL.  相似文献   

16.
This paper discusses the importance of the aerodynamic characteristics of forest and other similar canopies to modelling of boundary-layer flow and to estimating the diffusivity coefficients of turbulence transfer mechanisms over such canopies.The hypothesis of Marunich (1971) reported by Tajchman (1981) that the zero-plane displacement, d, equals the upward displacement of the flow trajectory, is critically examined. It is concluded that Marunich's hypothesis is conceptually incorrect and that calculations of d based on Marunich's hypothesis are inherently in error.This paper presents a method based on the mass conservation principle and uses wind profiles in and above a forest canopy as the sole input for determining d, z 0 and u *.Sensitivities of calculated results to measurements errors of wind profile data are evaluated. It is found that an error of less than 1% in wind in the logarithmic regime above the canopy can introduce up to 100% errors in calculated values of d, z 0 and u *. It is also found that the high sensitivity to wind data accuracy, characteristic of the present method, can be used as a guide for the selection of high quality canopy wind data.  相似文献   

17.
We have devised a partial differential equation for the prediction of dust concentration in a thin layer near the ground. In this equation, erosion (detachment), transport, deposition and source are parameterised in terms of known quantities. The interaction between a wind prediction model in the boundary layer and this equation affects the evolution of the dust concentration at the top of the surface layer. Numerical integrations are carried out for various values of source strength, ambient wind and particle size. Comparison with available data shows that the results appear very reasonable and that the model should be subjected to further development and testing.Notation (x, y, z, t) space co-ordinates and time (cm,t) - u, v components of horizontal wind speed (cm s–1) - u g, vg components of the geostrophic wind (cm s–1) - V=(u2+v2)1/2 (cm s–1) - (û v)= 1/(h – k) k h(u, v)dz(cm s–1) - V * friction velocity (cm s–1) - z 0 roughness length (cm) - k 1 von Karman constant =0.4 - V d deposition velocity (cm s–1) - V g gravitational settling velocity (cm s–1) - h height of inversion (cm) - k height of surface layer (cm) - potential temperature (°K) - gr potential temperature at ground (°K) - K potential temperature at top of surface layer (°K) - P pressure (mb) - P 0 sfc pressure (mb) - C p/Cv - (t)= /z lapse rate of potential temperature (°K cm–1) - A(z) variation of wind with height in transition layer - B(z) variation of wind with height in transition layer - Cd drag coefficient - C HO transfer coefficient for sensible heat - C dust concentration (g m–3) - C K dust concentration at top of surface layer (g m–3) - D(z) variation with height of dust concentration - u, v, w turbulent fluctuations of the three velocity components (cm s–1) - A 1 constant coefficient of proportionality for heat flux =0.2 - Ri Richardson number - g gravitational acceleration =980 cm s–2 - Re Reynolds number = - D s thickness of laminar sub-layer (cm) - v molecular kinematic viscosity of air - coefficient of proportionality in source term - dummy variable - t time step (sec) - n time index in numerical equations On sabbatical leave at University of Aberdeen, Department of Engineering, September 1989–February 1990.  相似文献   

18.
An expression is derived relating the critical flux Richardson number with the critical (gradient) Richardson number. In contrast to an earlier analysis by Townsend (1958), which is restricted to the atmosphere well outside the earth's boundary layer, the present treatment is intended specifically for turbulent flow in the lower atmosphere and it takes account of the effect of evaporation on the stability. The effect of radiation on the rate of destruction of the mean square of the temperature fluctuations is obtained by considering the radiative flux divergence in a stratified atmosphere and by using a simple functional relationship to represent empirical emissivity data.It was found that evaporation and radiation increase the critical Richardson number by a sensible amount depending on the atmospheric conditions, mainly temperature, humidity and the gradients. There is no definite critical Richardson number but rather a range between 0.25, below which turbulence is very likely, and somewhat higher than 0.5, above which turbulence is improbable. The value of the critical Richardson number can be expressed in terms of evaporation, radiation and the ratio ( w /u *) which also appears not to have a definite critical value. Evaporation and radiation cause the ratio ( w /u *) to be larger than unity under neutral conditions. These results, based on the assumption of Reynolds' analogy,K H =K M , are consistent with the available experimental evidence.  相似文献   

19.
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.  相似文献   

20.
The scintillation method tested over a dry vineyard area   总被引：8，自引：1，他引：8  
Measurements of a scintillometer device mounted at 4 m above a dry vineyard area in La Mancha, Spain, are used to obtain an average sensible heat flux densityH. Averaging is over a rectangular area determined by the distance between the scintillometer light source and receptor (875 m) and some upwind distance governed by the horizontal wind speed perpendicular to that line. Using similarity relations obtained from La Crau, a good correspondence betweenH measured with the scintillometer and an eddy-correlation device in the centre of a vineyard is obtained. The friction velocityu * was either measured directly using a sonic anemometer or obtained indirectly from two wind speeds and known values of the roughness length zo and displacementd. The free convection formulation underestimates the sensible heat flux by about 30%. This is due to a significant contribution of mechanically generated turbulence to the total turbulent transport, which was caused by relatively strong winds and rough terrain.  相似文献   

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

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