首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 62 毫秒
1.
2.
Large-scale structures within a rough-wall boundary layer generated over a cube array have recently been linked to small-scale fluctuations close to the roughness through a dynamical mechanism similar to amplitude modulation. Demonstrating the existence of this mechanism for different roughness types is a crucial step towards the development of a generic model for wind fluctuations in the urban canopy. Here the influence of the upstream roughness geometry (two-dimensional (2D) and three-dimensional (3D)) and planform packing density (\( \lambda_{p} \)) and street-canyon aspect ratio on the non-linear interactions between large-scale momentum regions and the small scales induced by the presence of the roughness is studied within a wind tunnel using combined particle-image velocimetry and hot-wire anemometry. A multi-time delay linear stochastic estimation is used to decompose the flow into large scales that participate in modulation and the remaining small scales. Using three different upstream roughness configurations composed of either 3D cubes or 2D rectangular blocks it is shown that the upstream roughness configuration has an influence on the non-linear interactions in the rough-wall boundary layer. Analysis of the turbulence skewness decomposition shows a change in the location of the maximum of the term \( \overline{{u_{L}^{\prime} u_{S}^{\prime 2}}} \), which represents the influence of the large-scale momentum regions on the small scales, whilst the temporal correlation shows a modification of the interaction located closer to the roughness with a change from 3D to 2D roughness. Furthermore, a two-point spatio–temporal correlation demonstrates that the non-linear relationship is significantly modified in the wake-interference-flow regime compared to the skimming-flow regime. Through skewness decomposition and temporal correlations the canyon aspect ratio is shown to have no influence on the non-linear interactions, indicating that the mechanism depends only on the flow developing upstream. Finally, although the upstream roughness configuration is shown to influence the non-linear interactions, the nature of the mechanism remains the same in all configurations.  相似文献   

3.
This study presents spatiotemporally-resolved measurements of surface shear-stress τ s in live plant canopies and rigid wooden cube arrays to identify the sheltering capability against sediment erosion of these different roughness elements. Live plants have highly irregular structures that can be extremely flexible and porous resulting in considerable changes to the drag and flow regimes relative to rigid imitations mainly used in other wind-tunnel studies. Mean velocity and kinematic Reynolds stress profiles show that well-developed natural boundary layers were generated above the 8 m long wind-tunnel test section covered with the roughness elements at four different roughness densities (λ = 0, 0.017, 0.08, 0.18). Speed-up around the cubes caused higher peak surface shear stress than in experiments with plants at all roughness densities, demonstrating the more effective sheltering ability of the plants. The sheltered areas in the lee of the plants are significantly narrower with higher surface shear stress than those found in the lee of the cubes, and are dependent on the wind speed due to the plants ability to streamline with the flow. This streamlining behaviour results in a decreasing sheltering effect at increasing wind speeds and in lower net turbulence production than in experiments with cubes. Turbulence intensity distributions suggest a suppression of horseshoe vortices in the plant case. Comparison of the surface shear-stress measurements with sediment erosion patterns shows that the fraction of time a threshold skin friction velocity is exceeded can be used to assess erosion of, and deposition on, that surface.  相似文献   

4.
Prediction of windthrow risk to individual or groups of retained trees in harvested stands requires an improved understanding of canopy airflow dynamics. Large-eddy simulations were used to simulate wind-tunnel experiments in two and three dimensions to compare with observations for model validation and to address parameter space considerations for the design of subsequent retention pattern experiments. The three-dimensional simulations were similar to the observed wind-tunnel data for the statistical profiles for but there were greater differences in skewness and kurtosis. These results were obtained using a common leaf-area drag formulation without either skin friction or speed dependent drag that enables scaling with U 0 (ambient wind speed) and h (height of the canopy). This scaling results in a single non-dimensional parameter h/h c where h c (x, y, z) is the momentum range resulting from the canopy drag. The validity of the model scaling was tested using two-dimensional simulations. The irrotational component of the flow (potential flow) was found to be important when defining vertical domain limitations and has significant implications for time dependent flow (i.e. turbulent conditions) when considering retention pattern design. The sudden onset of drag associated with the isolated stand presents some unexpected challenges. The horizontal scales of the shearing instabilities were simulated in two dimensions and found to range between 2h for early times to 7h for later times. The early-time horizontal scales are in the range of logical retention pattern scales and as such need to be taken into account as part of the parameter space, i.e. a range of retention pattern lengths need consideration.  相似文献   

5.
The air–sea transfer velocity of $\mathrm{CO}_{2}\, (k_{\mathrm{CO}_{2}})$ was investigated in a shallow estuary in March to July 2012, using eddy-covariance measurements of $\mathrm{CO}_{2}$ fluxes and measured air–sea $\mathrm{CO}_{2}$ partial-pressure differences. A data evaluation method that eliminates data by nine rejection criteria in order to heighten parametrization certainty is proposed. We tested the data evaluation method by comparing two datasets: one derived using quality criteria related solely to the eddy-covariance method, and the other derived using quality criteria based on both eddy-covariance and cospectral peak methods. The best parametrization of transfer velocity normalized to a Schmidt number of 600 $(k_{600})$ was determined to be: $k_{600} = 0.3\,{U_{10}}^{2.5}$ where $U_{10}$ is the wind speed in m $\mathrm{s}^{-1}$ at 10 m; $k_{600}$ is based on $\mathrm{CO}_{2}$ fluxes calculated by the eddy-covariance method and including the cospectral peak method criteria. At low wind speeds, the transfer velocity in the shallow water estuary was lower than in other coastal waters, possibly a symptom of low tidal amplitude leading to low intensity water turbulence. High transfer velocities were recorded above wind speeds of 5 m $\mathrm{s}^{-1}$ , believed to be caused by early-breaking waves and the large fetch (6.5 km) of the estuary. These findings indicate that turbulence in both air and water influences the transfer velocity.  相似文献   

6.
We present surface-layer measurements of temperature fluctuation variance from a site characterized by small-scale inhomogeneities. Periods of marked radiative forcing are selected. The data characterized by diabatic conditions and vertical heat flux larger than some threshold (here, chosen to be 0.01 K ms−1) agree quite well with the convective scaling in unstable cases, and with the z-less parameterisation (with a large scatter) in stable cases. For near-neutral cases, the similarity function diverges because of the loss of significance of the temperature scale. Departures from similarity are highlighted in cases with smaller thermal fluxes, because horizontal heterogeneity and unsteadiness become important as production terms.  相似文献   

7.
The efficacy of vegetation dynamics simulations in offline land surface models (LSMs) largely depends on the quality and spatial resolution of meteorological forcing data. In this study, the Princeton Global Meteorological Forcing Data (PMFD) and the high spatial resolution and upscaled China Meteorological Forcing Data (CMFD) were used to drive the Simplified Simple Biosphere model version 4/Top-down Representation of Interactive Foliage and Flora Including Dynamics (SSiB4/TRIFFID) and investigate how meteorological forcing datasets with different spatial resolutions affect simulations over the Tibetan Plateau (TP), a region with complex topography and sparse observations. By comparing the monthly Leaf Area Index (LAI) and Gross Primary Production (GPP) against observations, we found that SSiB4/TRIFFID driven by upscaled CMFD improved the performance in simulating the spatial distributions of LAI and GPP over the TP, reducing RMSEs by 24.3% and 20.5%, respectively. The multi-year averaged GPP decreased from 364.68 gC \begin{document}${\mathrm{m}}^{-2}\;{\mathrm{y}\mathrm{r}}^{-1}$\end{document} to 241.21 gC \begin{document}${\mathrm{m}}^{-2}\;{\mathrm{y}\mathrm{r}}^{-1}$\end{document} with the percentage bias dropping from 50.2% to –1.7%. When using the high spatial resolution CMFD, the RMSEs of the spatial distributions of LAI and GPP simulations were further reduced by 7.5% and 9.5%, respectively. This study highlights the importance of more realistic and high-resolution forcing data in simulating vegetation growth and carbon exchange between the atmosphere and biosphere over the TP.  相似文献   

8.
The spatial peak surface shear stress tS¢¢{\tau _S^{\prime\prime}} on the ground beneath vegetation canopies is responsible for the onset of particle entrainment and its precise and accurate prediction is essential when modelling soil, snow or sand erosion. This study investigates shear-stress partitioning, i.e. the fraction of the total fluid stress on the entire canopy that acts directly on the surface, for live vegetation canopies (plant species: Lolium perenne) using measurements in a controlled wind-tunnel environment. Rigid, non-porous wooden blocks instead of the plants were additionally tested for the purpose of comparison since previous wind-tunnel studies used exclusively artificial plant imitations for their experiments on shear-stress partitioning. The drag partitioning model presented by Raupach (Boundary-Layer Meteorol 60:375–395, 1992) and Raupach et al. (J Geophys Res 98:3023–3029, 1993), which allows the prediction of the total shear stress τ on the entire canopy as well as the peak (tS ¢¢/t)1/2{(\tau _S ^{\prime\prime}/\tau )^{1/2}} and the average (tS/t)1/2{(\tau _S^{\prime}/\tau )^{1/2}} shear-stress ratios, is tested against measurements to determine the model parameters and the model’s ability to account for shape differences of various roughness elements. It was found that the constant c, needed to determine the total stress τ and which was unspecified to date, can be assumed a value of about c = 0.27. Values for the model parameter m, which accounts for the difference between the spatial surface average tS{\tau _S^{\prime}} and the peak tS ¢¢{\tau _S ^{\prime\prime}} shear stress, are difficult to determine because m is a function of the roughness density, the wind velocity and the roughness element shape. A new definition for a parameter a is suggested as a substitute for m. This a parameter is found to be more closely universal and solely a function of the roughness element shape. It is able to predict the peak surface shear stress accurately. Finally, a method is presented to determine the new a parameter for different kinds of roughness elements.  相似文献   

9.
We examine the effect of varying roughness-element aspect ratio on the mean velocity distributions of turbulent flow over arrays of rectangular-prism-shaped elements. Large-eddy simulations (LES) in conjunction with a sharp-interface immersed boundary method are used to simulate spatially-growing turbulent boundary layers over these rough surfaces. Arrays of aligned and staggered rectangular roughness elements with aspect ratio >1 are considered. First the temporally- and spatially-averaged velocity profiles are used to illustrate the aspect-ratio effects. For aligned prisms, the roughness length (\(z_\mathrm{o}\)) and the friction velocity (\(u_*\)) increase initially with an increase in the roughness-element aspect ratio, until the values reach a plateau at a particular aspect ratio. The exact value of this aspect ratio depends on the coverage density. Further increase in the aspect ratio changes neither \(z_\mathrm{o}\), \(u_*\) nor the bulk flow above the roughness elements. For the staggered cases, \(z_\mathrm{o}\) and \(u_*\) continue to increase for the surface coverage density and the aspect ratios investigated. To model the flow response to variations in roughness aspect ratio, we turn to a previously developed phenomenological volumetric sheltering model (Yang et al., in J Fluid Mech 789:127–165, 2016), which was intended for low to moderate aspect-ratio roughness elements. Here, we extend this model to account for high aspect-ratio roughness elements. We find that for aligned cases, the model predicts strong mutual sheltering among the roughness elements, while the effect is much weaker for staggered cases. The model-predicted \(z_\mathrm{o}\) and \(u_*\) agree well with the LES results. Results show that the model, which takes explicit account of the mutual sheltering effects, provides a rapid and reliable prediction method of roughness effects in turbulent boundary-layer flows over arrays of rectangular-prism roughness elements.  相似文献   

10.
Data from the Antarctic winter at Halley Base have been used in order to evaluate qualitatively and quantitatively how the stratification in the low atmosphere (evaluated with the gradient Richardson number, Ri) influences the eddy transfers of heat and momentum. Vertical profiles of wind and temperature up to 32 m, and turbulent fluxes ( , and ) measured from three ultrasonic thermo-anemometers installed at 5, 17 and 32 m are employed to calculate Ri, the friction velocity (u *) and the eddy diffusivities for heat (K h ) and momentum (K m ). The results show a big dependence of stability onK m ,K h andu *, with a sharp decrease of these turbulent parameters with increasing stability. The ratio of eddy diffusivities (K h /K m ) is also analyzed and presents a decreasing tendency as Ri increases, reaching values even less than 1, i.e., there were situations where the turbulent transfer of momentum was greater than that of heat. Possible mechanisms of turbulent mixing are discussed.  相似文献   

11.
In 2005 the Study of Stable Boundary Layer Environment at Dome C (STABLEDC) experimental campaign was conducted at the plateau station of Concordia at Dome C, Antarctica. Temperature profiles measured with a microwave radiometer were used to study the characteristics of surface-based temperature inversions over the course of a year. Statistics of temperature profiles for every month are discussed; the difference between daytime and nocturnal cases observed during the summer months disappears during winter. Surface-based temperature inversions occurred in 70 % of the time during summer, and almost all of the time during winter. During winter the occurrence of warming events leads to a decrease in the temperature difference between the top and the base of the inversion (i.e. the inversion strength). The inversion strength maxima ranged between $3\,^{\circ }\mathrm{C}$ (December) and $35\,^{\circ }\mathrm{C}$ (August) corresponding to gradients of 0.1 and $0.3\,^{\circ }\mathrm{C}\, \mathrm{m}^{-1}$ , respectively. The average surface-based inversion height presents a daily cycle during the summer months with values up to 200 m in the morning hours, while it affects a layer always deeper than 100 m during the winter months. The relationships between inversion strength and the downward longwave radiative flux, absolute temperature, and wind speed are examined. The inversion strength decreases as the longwave radiation increases. A clear anti-correlation between inversion strength and near-surface temperature is evident throughout the year. During the winter, the largest inversion strength values were observed under low wind-speed conditions; in contrast, a clear dependence was not found during the summer.  相似文献   

12.
Spatial Variability of Flow Statistics within Regular Building Arrays   总被引:2,自引:2,他引:0  
Turbulence statistics obtained by direct numerical simulations are analysed to investigate spatial heterogeneity within regular arrays of building-like cubical obstacles. Two different array layouts are studied, staggered and square, both at a packing density of . The flow statistics analysed are mean streamwise velocity (), shear stress (), turbulent kinetic energy (k) and dispersive stress fraction (). The spatial flow patterns and spatial distribution of these statistics in the two arrays are found to be very different. Local regions of high spatial variability are identified. The overall spatial variances of the statistics are shown to be generally very significant in comparison with their spatial averages within the arrays. Above the arrays the spatial variances as well as dispersive stresses decay rapidly to zero. The heterogeneity is explored further by separately considering six different flow regimes identified within the arrays, described here as: channelling region, constricted region, intersection region, building wake region, canyon region and front-recirculation region. It is found that the flow in the first three regions is relatively homogeneous, but that spatial variances in the latter three regions are large, especially in the building wake and canyon regions. The implication is that, in general, the flow immediately behind (and, to a lesser extent, in front of) a building is much more heterogeneous than elsewhere, even in the relatively dense arrays considered here. Most of the dispersive stress is concentrated in these regions. Considering the experimental difficulties of obtaining enough point measurements to form a representative spatial average, the error incurred by degrading the sampling resolution is investigated. It is found that a good estimate for both area and line averages can be obtained using a relatively small number of strategically located sampling points.  相似文献   

13.
On the location and orientation of the South Pacific Convergence Zone   总被引:2,自引:0,他引:2  
Three semi-permanent cloud bands exist in the Southern Hemisphere extending southeastward from the equator, through the tropics, and into the subtropics. The most prominent of these features occurs in the South Pacific and is referred to as the South Pacific Convergence Zone (SPCZ). Similar bands, with less intensity, exist in the South Indian and Atlantic oceans. We attempt to explain the physical mechanisms that promote the diagonal orientation of the SPCZ and the processes that determine the timescales of its variability. It is argued that the slowly varying sea surface temperature patterns produce upper tropospheric wind fields that vary substantially in longitude. Regions where 200?hPa zonal winds decrease with longitude (i.e., negative zonal stretching deformation, or $ {{\partial \overline{U} } \mathord{\left/ {\vphantom {{\partial \overline{U} } {\partial x}}} \right. \kern-0em} {\partial x}} < 0 $ ) reduce the group speed of the eastward propagating synoptic (3?C6?day period) Rossby waves and locally increase the wave energy density. Such a region of wave accumulation occurs in the vicinity of the SPCZ, thus providing a physical basis for the diagonal orientation and earlier observations that the zone acts as a ??graveyard?? of propagating synoptic disturbances. In essence, $ {{\partial \overline{U} } \mathord{\left/ {\vphantom {{\partial \overline{U} } {\partial x}}} \right. \kern-0em} {\partial x}} = 0 $ demarks the boundary of the graveyard while regions where $ {{\partial \overline{U} } \mathord{\left/ {\vphantom {{\partial \overline{U} } {\partial x}}} \right. \kern-0em} {\partial x}} < 0 $ denote the graveyard itself. Composites of the life cycles of synoptic waves confirm this hypothesis. From the graveyard hypothesis comes a more general theory accounting for the SPCZ??s spatial orientation and its longer term variability influenced by the El Ni?o-Southern Oscillation (ENSO), or alternatively, the changing background SST associated with different phases of ENSO.  相似文献   

14.
In order to quantitatively investigate the role of leads and sea-ice in air-mass modification, aircraft observations were conducted over the partially ice-covered Sea of Okhotsk. We investigated two cold-air outbreak events with different sea-ice concentrations. In both cases, the difference between the temperatures of surface air and the sea surface (ΔT) dropped rapidly with the accumulated fetch-width of leads up to about 35-40 km, and then decreased very slowly. The surface sensible heat flux originating from open water was about 300 W m−2 within a few kilometres from the coast and decreased with increasing accumulated fetch-width. The sensible heat flux was about 100 W m−2 on average. These results indicate that the downwind air-mass modification depends mainly on the total (accumulated) extent of open water. The total buoyancy flux calculated by the joint frequency distribution method correlated very well with ice concentration. Such a relationship was not clear in the case of the moisture flux . The ratio between rising thermals and cold downdrafts differed significantly between upwind and downwind regions; that is, the buoyancy flux was dominated by in the developing stage of the boundary layer, while also became important after the development of the boundary layer.  相似文献   

15.
The present work provides a new methodology to determine onset dates of the rainy season (ONR) in central Amazon (CAM) using the antisymmetric in relation to the equator outgoing longwave radiation (OLR) (AOLR) data, for the 1979–2006 period. Spatial averages of the AOLR ( $\overline {AOLR} $ ) over the CAM for the ONR periods are obtained. These periods correspond to 25 pentads centered on the mean pentad for the ONR. The sign changes from positive to negative of the $\overline {AOLR} $ for the ONR periods indicate the transition from dry to wet season. Composites of several variables are done for pentads before and after the ONR dates. These composites show physically consistent features. The potential of the $\overline {AOLR} $ time series as an index for monitoring tasks is analyzed. The results here show that the $\overline {AOLR} $ for the ONR period captures the transition from dry to wet conditions in the CAM area during 2006. The advantages of this method are discussed. The new simple method proposed here seems to be efficient in determining the ONR in the CAM.  相似文献   

16.
The turbulent characteristics of the neutral boundary layer developing over rough surfaces are not well predicted with operational weather-forecasting models. The problem is attributed to inadequate mixing-length models, to the anisotropy of the flow and to a lack of controlled experimental data against which to validate numerical studies. Therefore, in order to address directly the modelling difficulties for the development of a neutral boundary layer over rough surfaces, and to investigate the turbulent momentum transfer of such a layer, a set of hydraulic flume experiments were carried out. In the experiments, the mean and turbulent quantities were measured by a particle image velocimetry (PIV) technique. The measured velocity variances and fluxes \({(\overline{{u_{i}^{\prime}}{u_{j}^{\prime}}})}\) in longitudinal vertical planes allowed the vertical and longitudinal gradients (?/?z and ?/?x) of the mean and turbulent quantities (fluxes, variances and third-order moments) to be evaluated and the terms of the evolution equations for ?e/?t, \({\partial \overline{u^{\prime 2}}/\partial t}\), \({\partial \overline{w^{\prime 2}}/\partial t}\) and \({\partial \overline{{u^{\prime}}{w^{\prime}}}/\partial t}\) to be quantified, where e is the turbulent kinetic energy. The results show that the pressure-correlation terms allow the turbulent energy to be transferred equitably from \({\overline{{u^{\prime}}^{2}}}\) to \({\overline{{w^{\prime}}^{2}}}\). It appears that the repartition between the constitutive terms of the budget of e, \({\overline{{u^{\prime}}^{2}}}\), \({\overline{{w^{\prime}}^{2}}}\) and \({\overline{{u^{\prime}}{w^{\prime}}}}\) is not significantly affected by the development of the rough neutral boundary layer. For the whole evolution, the transfers of energy are governed by the same terms that are also very similar to the smooth-wall case. The PIV measurements also allowed the spatial integral scales to be computed directly and to be compared with the dissipative and mixing length scales, which were also computed from the data.  相似文献   

17.
The performance of an atmospheric single-column model (SCM) is studied systematically for stably-stratified conditions. To this end, 11 years (2005–2015) of daily SCM simulations were compared to observations from the Cabauw observatory, The Netherlands. Each individual clear-sky night was classified in terms of the ambient geostrophic wind speed with a \(1\hbox { m} \hbox { s}^{-1}\) bin-width. Nights with overcast conditions were filtered out by selecting only those nights with an average net radiation of less than \(-\,30\hbox { W }\hbox {m}^{-2}\). A similar procedure was applied to the observational dataset. A comparison of observed and modelled ensemble-averaged profiles of wind speed and potential temperature and time series of turbulent fluxes showed that the model represents the dynamics of the nocturnal boundary layer (NBL) at Cabauw very well for a broad range of mechanical forcing conditions. No obvious difference in model performance was found between near-neutral and strongly-stratified conditions. Furthermore, observed NBL regime transitions are represented in a natural way. The reference model version performs much better than a model version that applies excessive vertical mixing as is done in several (global) operational models. Model sensitivity runs showed that for weak-wind conditions the inversion strength depends much more on details of the land-atmosphere coupling than on the turbulent mixing. The presented results indicate that in principle the physical parametrizations of large-scale atmospheric models are sufficiently equipped for modelling stably-stratified conditions for a wide range of forcing conditions.  相似文献   

18.
19.
A Reynolds-averaged Navier–Stokes model is used to investigate the evolution of the sectional drag coefficient and turbulent length scales with the layouts of aligned arrays of cubes. Results show that the sectional drag coefficient is determined by the non-dimensional streamwise distance (sheltering parameter), and the non-dimensional spanwise distance (channelling parameter) between obstacles. This is different than previous approaches that consider only plan area density $(\lambda _\mathrm{p})$ . On the other hand, turbulent length scales behave similarly to the staggered case (e. g. they are function of $\lambda _\mathrm{p}$ only). Analytical formulae are proposed for the length scales and for the sectional drag coefficient as a function of sheltering and channelling parameters, and implemented in a column model. This approach demonstrates good skill in the prediction of vertical profiles of the spatially-averaged horizontal wind speed.  相似文献   

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

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

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