首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 890 毫秒
1.
Lagrangian and Eulerian statistics are obtained from a water-channel experiment of an idealized two-dimensional urban canopy flow in neutral conditions. The objective is to quantify the Eulerian \((T^{\mathrm{E}})\) and Lagrangian \((T^{\mathrm{L}})\) time scales of the turbulence above the canopy layer as well as to investigate their dependence on the aspect ratio of the canopy, AR, as the latter is the ratio of the width (W) to the height (H) of the canyon. Experiments are also conducted for the case of flat terrain, which can be thought of as equivalent to a classical one-directional shear flow. The values found for the Eulerian time scales on flat terrain are in agreement with previous numerical results found in the literature. It is found that both the streamwise and vertical components of the Lagrangian time scale, \(T_\mathrm{u}^\mathrm{L} \) and \(T_\mathrm{w}^\mathrm{L} \), follow Raupach’s linear law within the constant-flux layer. The same holds true for \(T_\mathrm{w}^\mathrm{L} \) in both the canopies analyzed \((AR= 1\) and \(AR= 2\)) and also for \(T_\mathrm{u}^\mathrm{L} \) when \(AR = 1\). In contrast, for \(AR = 2\), \(T_\mathrm{u}^\mathrm{L} \) follows Raupach’s law only above \(z=2H\). Below that level, \(T_\mathrm{u}^\mathrm{L} \) is nearly constant with height, showing at \(z=H\) a value approximately one order of magnitude greater than that found for \(AR = 1\). It is shown that the assumption usually adopted for flat terrain, that \(\beta =T^{\mathrm{L}}/T^{\mathrm{E}}\) is proportional to the inverse of the turbulence intensity, also holds true even for the canopy flow in the constant-flux layer. In particular, \(\gamma /i_\mathrm{u} \) fits well \(\beta _\mathrm{u} =T_\mathrm{u}^\mathrm{L} /T_\mathrm{u}^\mathrm{E} \) in both the configurations by choosing \(\gamma \) to be 0.35 (here, \(i_\mathrm{u} =\sigma _\mathrm{u} / \bar{u} \), where \(\bar{u} \) and \(\sigma _\mathrm{u} \) are the mean and the root-mean-square of the streamwise velocity component, respectively). On the other hand, \(\beta _\mathrm{w} =T_\mathrm{w}^\mathrm{L} /T_\mathrm{w}^\mathrm{E} \) follows approximately \(\gamma /i_\mathrm{w} =0.65/\left( {\sigma _\mathrm{w} /\bar{u} } \right) \) for \(z > 2H\), irrespective of the AR value. The second main objective is to estimate other parameters of interest in dispersion studies, such as the eddy diffusivity of momentum \((K_\mathrm{{T}})\) and the Kolmogorov constant \((C_0)\). It is found that \(C_0\) depends appreciably on the velocity component both for the flat terrain and canopy flow, even though for the latter case it is insensitive to AR values. In all the three experimental configurations analyzed here, \(K_\mathrm{{T}}\) shows an overall linear growth with height in agreement with the linear trend predicted by Prandtl’s theory.  相似文献   

2.
Evaporation from wet-canopy (\(E_\mathrm{C}\)) and stem (\(E_\mathrm{S}\)) surfaces during rainfall represents a significant portion of municipal-to-global scale hydrologic cycles. For urban ecosystems, \(E_\mathrm{C}\) and \(E_\mathrm{S}\) dynamics play valuable roles in stormwater management. Despite this, canopy-interception loss studies typically ignore crown-scale variability in \(E_\mathrm{C}\) and assume (with few indirect data) that \(E_\mathrm{S}\) is generally \({<}2\%\) of total wet-canopy evaporation. We test these common assumptions for the first time with a spatially-distributed network of in-canopy meteorological monitoring and 45 surface temperature sensors in an urban Pinus elliottii tree row to estimate \(E_\mathrm{C}\) and \(E_\mathrm{S}\) under the assumption that crown surfaces behave as “wet bulbs”. From December 2015 through July 2016, 33 saturated crown periods (195 h of 5-min observations) were isolated from storms for determination of 5-min evaporation rates ranging from negligible to 0.67 \(\hbox {mm h}^{-1}\). Mean \(E_\mathrm{S}\) (0.10 \(\hbox {mm h}^{-1}\)) was significantly lower (\(p < 0.01\)) than mean \(E_\mathrm{C}\) (0.16 \(\hbox {mm h}^{-1}\)). But, \(E_\mathrm{S}\) values often equalled \(E_\mathrm{C}\) and, when scaled to trunk area using terrestrial lidar, accounted for 8–13% (inter-quartile range) of total wet-crown evaporation (\(E_\mathrm{S}+E_\mathrm{C}\) scaled to surface area). \(E_\mathrm{S}\) contributions to total wet-crown evaporation maximized at 33%, showing a general underestimate (by 2–17 times) of this quantity in the literature. Moreover, results suggest wet-crown evaporation from urban tree rows can be adequately estimated by simply assuming saturated tree surfaces behave as wet bulbs, avoiding problematic assumptions associated with other physically-based methods.  相似文献   

3.
The Nieuwstadt closed-form solution for the stationary Ekman layer is generalized for katabatic flows within the conceptual framework of the Prandtl model. The proposed solution is valid for spatially-varying eddy viscosity and diffusivity (O’Brien type) and constant Prandtl number (Pr). Variations in the velocity and buoyancy profiles are discussed as a function of the dimensionless model parameters \(z_0 \equiv \hat{z}_0 \hat{N}^2 Pr \sin {(\alpha )} |\hat{b}_\mathrm{s} |^{-1}\) and \(\lambda \equiv \hat{u}_{\mathrm{ref}}\hat{N} \sqrt{Pr} |\hat{b}_\mathrm{s} |^{-1}\), where \(\hat{z}_0\) is the hydrodynamic roughness length, \(\hat{N}\) is the Brunt-Väisälä frequency, \(\alpha \) is the surface sloping angle, \(\hat{b}_\mathrm{s}\) is the imposed surface buoyancy, and \(\hat{u}_{\mathrm{ref}}\) is a reference velocity scale used to define eddy diffusivities. Velocity and buoyancy profiles show significant variations in both phase and amplitude of extrema with respect to the classic constant \(\textit{K}\) model and with respect to a recent approximate analytic solution based on the Wentzel-Kramers-Brillouin theory. Near-wall regions are characterized by relatively stronger surface momentum and buoyancy gradients, whose magnitude is proportional to \(z_0\) and to \(\lambda \). In addition, slope-parallel momentum and buoyancy fluxes are reduced, the low-level jet is further displaced toward the wall, and its peak velocity depends on both \(z_0\) and \(\lambda \).  相似文献   

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

5.
Mountain-top observations of greenhouse gas mixing ratios may be an alternative to tall-tower measurements for regional scale source and sink estimation. To investigate the equivalence or limitations of a mountain-top site as compared to a tall-tower site, we used the unique opportunity of comparing in situ measurements of methane (\(\hbox {CH}_{4}\)) and carbon dioxide (\(\hbox {CO}_{2}\)) mixing ratios at a mountain top (986 m above sea level, a.s.l.) with measurements from a nearby (distance 28.4 km) tall tower, sampled at almost the same elevation (1009 m a.s.l.). Special attention was given to, (i) how local wind statistics and greenhouse gas sources and sinks at the mountain top influence the observations, and (ii) whether mountain-top observations can be used as for those from a tall tower for constraining regional greenhouse gas emissions. Wind statistics at the mountain-top site are clearly more influenced by local flow systems than those at the tall-tower site. Differences in temporal patterns of the greenhouse gas mixing ratios observed at the two sites are mostly related to the influence of local sources and sinks at the mountain-top site. Major influences of local sources can be removed by applying a statistical filter (\(5{\mathrm{th}}\) percentile) or a filter that removes periods with unfavourable flow conditions. In the best case, the bias in mixing ratios between the mountain-top and the tall-tower sites after the application of the wind filter was \({-}0.0005\pm 0.0010\) ppm for methane (September, 0000–0400 UTC) and \(0.11\pm 0.18\) ppm for \(\hbox {CO}_{2}\) (February, 1200–1600 UTC). Temporal fluctuations of atmospheric \(\hbox {CH}_{4}\) and \(\hbox {CO}_{2}\) mixing ratios at both stations also showed good agreement (apart from \(\hbox {CO}_{2}\) during summertime) as determined by moving bi-weekly Pearson correlation coefficients (up to 0.96 for \(\hbox {CO}_{2}\) and 0.97 for \(\hbox {CH}_{4}\)). When only comparing mixing ratios minimally influenced by local sources (low bias and high correlation coefficients), our measurements indicate that mountain-top observations are comparable to tall-tower observations.  相似文献   

6.
Methane (\(\mathrm {CH}_{4}\)) is known to be emitted from lakes to the atmosphere via processes such as diffusion and ebullition (i.e., bubble emission). We developed a practical method for partitioning eddy-covariance \(\mathrm {CH}_{4}\) fluxes from a shallow lake into diffusive and ebullitive fluxes using a wavelet analysis based on local scalar similarity between the \(\mathrm {CH}_{4}\) concentration and other reference scalars, such as the air temperature or water vapour concentration, in the wavelet time-scale domain, with the hypothesis that similar and dissimilar fluctuation components are related to diffusive and ebullitive \(\mathrm {CH}_{4}\) fluxes, respectively. Our method is applied to approximately two weeks of data obtained at a shallow mid-latitude lake. The partitioned diffusive flux has a physically sound relationship with wind speed, supporting the validity of the method. The ratio of the diffusive flux to the total flux is typically 0.11 with flow from an area of steady bubble emission, but otherwise 0.36. Further validation is required using a larger dataset and data from other lakes. The proposed method can be easily applied to historical data because it requires only 10-Hz data of \(\mathrm {CH}_{4}\) concentration and other reference scalars, along with an empirical parameter.  相似文献   

7.
Wind-tunnel experiments were carried out on fully-rough boundary layers with large roughness (\(\delta /h \approx 10\), where h is the height of the roughness elements and \(\delta \) is the boundary-layer thickness). Twelve different surface conditions were created by using LEGO? bricks of uniform height. Six cases are tested for a fixed plan solidity (\(\lambda _\mathrm{P}\)) with variations in frontal density (\(\lambda _\mathrm{F}\)), while the other six cases have varying \(\lambda _\mathrm{P}\) for fixed \(\lambda _\mathrm{F}\). Particle image velocimetry and floating-element drag-balance measurements were performed. The current results complement those contained in Placidi and Ganapathisubramani (J Fluid Mech 782:541–566, 2015), extending the previous analysis to the turbulence statistics and spatial structure. Results indicate that mean velocity profiles in defect form agree with Townsend’s similarity hypothesis with varying \(\lambda _\mathrm{F}\), however, the agreement is worse for cases with varying \(\lambda _\mathrm{P}\). The streamwise and wall-normal turbulent stresses, as well as the Reynolds shear stresses, show a lack of similarity across most examined cases. This suggests that the critical height of the roughness for which outer-layer similarity holds depends not only on the height of the roughness, but also on the local wall morphology. A new criterion based on shelter solidity, defined as the sheltered plan area per unit wall-parallel area, which is similar to the ‘effective shelter area’ in Raupach and Shaw (Boundary-Layer Meteorol 22:79–90, 1982), is found to capture the departure of the turbulence statistics from outer-layer similarity. Despite this lack of similarity reported in the turbulence statistics, proper orthogonal decomposition analysis, as well as two-point spatial correlations, show that some form of universal flow structure is present, as all cases exhibit virtually identical proper orthogonal decomposition mode shapes and correlation fields. Finally, reduced models based on proper orthogonal decomposition reveal that the small scales of the turbulence play a significant role in assessing outer-layer similarity.  相似文献   

8.
The assumption that the roughness Reynolds number \(( Re_{*})\) can be used as a basis for quantifying the boundary-layer property \({ kB}^{-1} (= \ln (z_{0}/z_{0T}))\) as in some modern numerical models is questioned. While \({ Re}_{*}\) is a useful property in studies of pipe flow, it appears to have only marginal applicability in the case of treeless terrain, as studied in the two experimental situations presented here. For both the daytime and night-time cases there appears to be little correlation between \({ kB}^{-1}\) and \({ Re}_{*}\). For daytime, the present studies indicate that the assumption \({ kB}^{-1} \approx 2\) is acceptable, while for night-time, the scatter involved in relating \({ kB}^{-1}\) to \({ Re}_{*}\) suggests there is little reason to assume a direct relationship. However, while the scatter affecting all of the night-time results is large, there remains a significant correlation between the heat and momentum fluxes upon which an alternative methodology for describing bulk air–surface exchange at night could be constructed. The friction coefficient (\(C_{f}\)) and the turbulent Stanton number \(({ St}_{*})\) are discussed as possible alternatives for describing bulk properties of the air layer adjacent to the surface. While describing the surface roughness in terms of the friction coefficient provides an attractive simplification relative to the conventional methodologies based on roughness length and stability considerations, use of the Stanton number shares many of uncertainties that affect \({ kB}^{-1}\). The transitions at dawn and dusk remain demanding situations to address.  相似文献   

9.
The influence of random errors in turbulence measurements on scalar similarity for temperature, water vapour, \(\hbox {CO}_{2}\), and \(\hbox {NH}_{3}\) is investigated using two eddy-covariance datasets collected over a lake and a cattle feedlot. Three measures of scalar similarity, namely, the similarity constant in the flux–variance relationship, the correlation coefficient between two scalars and the relative transport efficiency, are examined. The uncertainty in the similarity constant \(C_{s}\) in the flux–variance relationship resulting from random errors in turbulence measurements is quantified based on error propagation analyses and a Monte-Carlo sampling method, which yields a distribution instead of a single value for \(C_{s}\). For different scalars, the distributions of \(C_{s}\) are found to significantly overlap, implying that scalars are transported similarly under strongly unstable conditions. The random errors in the correlation coefficients between scalars and the relative transport efficiencies are also quantified through error propagation analyses, and they increase as the atmosphere departs from neutral conditions. Furthermore, the correlation coefficients between three scalars (water vapour, \(\hbox {CO}_{2}\), and \(\hbox {NH}_{3}\)) are statistically different from unity while the relative transport efficiencies are not, which highlights the difference between these two measures of scalar similarity. The results suggest that uncertainties in these measures of scalar similarity need to be quantified when using them to diagnose the existence of dissimilarity among different scalars.  相似文献   

10.
The influence of wave-associated parameters controlling turbulent \(\hbox {CO}_2\) fluxes through the air–sea interface is investigated in a coastal region. A full year of high-quality data of direct estimates of air–sea \(\hbox {CO}_2\) fluxes based on eddy-covariance measurements is presented. The study area located in Todos Santos Bay, Baja California, Mexico, is a net sink of \(\hbox {CO}_2\) with a mean flux of \(-1.3\, \upmu \hbox {mol m}^{-2}\hbox {s}^{-1}\) (\(-41.6\hbox { mol m}^{-2}\hbox {yr}^{-1}\)). The results of a quantile-regression analysis computed between the \(\hbox {CO}_2\) flux and, (1) wind speed, (2) significant wave height, (3) wave steepness, and (4) water temperature, suggest that the significant wave height is the most correlated parameter with the magnitude of the flux but the behaviour of the relation varies along the probability distribution function, with the slopes of the regression lines presenting both positive and negative values. These results imply that the presence of surface waves in coastal areas is the key factor that promotes the increase of the flux from and into the ocean. Further analysis suggests that the local characteristics of the aqueous and atmospheric layers might determine the direction of the flux.  相似文献   

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

12.
Nine methods to determine local-scale aerodynamic roughness length \((z_{0})\) and zero-plane displacement \((z_{d})\) are compared at three sites (within 60 m of each other) in London, UK. Methods include three anemometric (single-level high frequency observations), six morphometric (surface geometry) and one reference-based approach (look-up tables). A footprint model is used with the morphometric methods in an iterative procedure. The results are insensitive to the initial \(z_{d}\) and \(z_{0}\) estimates. Across the three sites, \(z_{d}\) varies between 5 and 45 m depending upon the method used. Morphometric methods that incorporate roughness-element height variability agree better with anemometric methods, indicating \(z_{d}\) is consistently greater than the local mean building height. Depending upon method and wind direction, \(z_{0}\) varies between 0.1 and 5 m with morphometric \(z_{0}\) consistently being 2–3 m larger than the anemometric \(z_{0}\). No morphometric method consistently resembles the anemometric methods. Wind-speed profiles observed with Doppler lidar provide additional data with which to assess the methods. Locally determined roughness parameters are used to extrapolate wind-speed profiles to a height roughly 200 m above the canopy. Wind-speed profiles extrapolated based on morphometric methods that account for roughness-element height variability are most similar to observations. The extent of the modelled source area for measurements varies by up to a factor of three, depending upon the morphometric method used to determine \(z_{d}\) and \(z_{0}\).  相似文献   

13.
Urban morphology and inter-building shadowing result in a non-uniform distribution of surface heating in urban areas, which can significantly modify the urban flow and thermal field. In Part I, we found that in an idealized three-dimensional urban array, the spatial distribution of the thermal field is correlated with the orientation of surface heating with respect to the wind direction (i.e. leeward or windward heating), while the dispersion field changes more strongly with the vertical temperature gradient in the street canyon. Here, we evaluate these results more closely and translate them into metrics of “city breathability,” with large-eddy simulations coupled with an urban energy-balance model employed for this purpose. First, we quantify breathability by, (i) calculating the pollutant concentration at the pedestrian level (horizontal plane at \(z\approx 1.5\)–2 m) and averaged over the canopy, and (ii) examining the air exchange rate at the horizontal and vertical ventilating faces of the canyon, such that the in-canopy pollutant advection is distinguished from the vertical removal of pollution. Next, we quantify the change in breathability metrics as a function of previously defined buoyancy parameters, horizontal and vertical Richardson numbers (\(Ri_\text {h}\) and \(Ri_\text {v}\), respectively), which characterize realistic surface heating. We find that, unlike the analysis of airflow and thermal fields, consideration of the realistic heating distribution is not crucial in the analysis of city breathability, as the pollutant concentration is mainly correlated with the vertical temperature gradient (\(Ri_\text {v}\)) as opposed to the horizontal (\(Ri_\text {h}\)) or bulk (\(Ri_\text {b}\)) thermal forcing. Additionally, we observe that, due to the formation of the primary vortex, the air exchange rate at the roof level (the horizontal ventilating faces of the building canyon) is dominated by the mean flow. Lastly, since \(Ri_\text {h}\) and \(Ri_\text {v}\) depend on the meteorological factors (ambient air temperature, wind speed, and wind direction) as well as urban design parameters (such as surface albedo), we propose a methodology for mapping overall outdoor ventilation and city breathability using this characterization method. This methodology helps identify the effects of design on urban microclimate, and ultimately informs urban designers and architects of the impact of their design on air quality, human health, and comfort.  相似文献   

14.
Comprehensive, ground-based observations from the US Department of Energy Atmospheric Radiation Measurements program Southern Great Plains site are used to study the variability of turbulence forcings and cloud-scale turbulence structures in a continental stratocumulus cloud. The turbulence observations are made from an upward facing cloud (35 GHz) Doppler radar. Cloud base and liquid water path are characterized using a lidar at the surface and a microwave radiometer. The turbulence characterizations are compared and contrasted with those observed in marine stratocumulus clouds. During the 16-h observation period used in this study the cloud-base and cloud-top heights evolve with time and changes in liquid water path observed by the radiometer are consistent with variations in cloud depth. Unlike marine stratocumulus clouds, a diurnal cycle of cloud thickness and liquid water path is not observed. The observed surface latent, sensible, and virtual sensible heat fluxes and the radiative fluxes exhibit a diurnal cycle with values increasing from sunrise to afternoon and decreasing afterwards. During the night, the sensible heat, virtual sensible heat and the net radiative fluxes at the surface are slightly negative. Solar radiative heating prevails in the cloud layer during the day and strong radiative cooling exists at cloud top even during the day. Unlike marine stratocumulus, surface heating described by the convective velocity scale \(W_\mathrm{s}^{*}\) and cloud-top cooling described by \(W_\mathrm{r}^{*}\) are both important in driving the in-cloud turbulence during the day, whereas cloud-top cooling is the exclusive contributor during the night. The combined \(W_\mathrm{s}^{*}\) and \(W_\mathrm{r}^{*}\) (the total velocity scale \(W_\mathrm{t}^{*})\) provides a useful way to track the evolution of the turbulence structure in the cloud. The variance of the radar-measured radial velocity, which is related to resolved turbulence, follows the diurnal cycle and is consistent with the total velocity scale \(W_\mathrm{t}^{*}\) variations. It is higher during the day and lower during the night, which is contrary to that in marine stratocumulus. The \(W_\mathrm{t}^{*}\) values are lowest around sunset when the radiative cooling is also small due to upper-level clouds observed above the low-level stratus. The vertical distribution of the variance results from the surface heating during the day and cloud-top cooling during the night. The squared spectrum width, which is related to turbulence structures within the radar sampling volume (unresolved turbulence) also follows the diurnal cycle. Its vertical distribution indicates that the unresolved turbulence more closely relates to the processes near cloud top. Turbulence in the cloud requires about an hour to respond to the external forcings of surface heating and cloud-top radiative cooling. Positive skewness prevails during the day and negative skewness prevails at night with a sharp transition around sunset. Resolved turbulence dominates near cloud base whereas unresolved turbulence dominates near cloud top. The turbulence characteristics and variability defined in this study can be used to evaluate the time evolution of turbulence structures in large eddy simulation forced by surface and cloud-top radiative forcings.  相似文献   

15.
The scale properties of anisotropic and isotropic turbulence in the urban surface layer are investigated. A dimensionless anisotropic tensor is introduced and the turbulent tensor anisotropic coefficient, defined as C, where \(C = 3d_{3}\,+\,1 (d_{3}\) is the minimum eigenvalue of the tensor) is used to characterize the turbulence anisotropy or isotropy. Turbulence is isotropic when \(C \approx 1\), and anisotropic when \(C \ll 1\). Three-dimensional velocity data collected using a sonic anemometer are analyzed to obtain the anisotropic characteristics of atmospheric turbulence in the urban surface layer, and the tensor anisotropic coefficient of turbulent eddies at different spatial scales calculated. The analysis shows that C is strongly dependent on atmospheric stability \(\xi = (z-z_{\mathrm{d}})/L_{{\textit{MO}}}\), where z is the measurement height, \(z_{\mathrm{d}}\) is the displacement height, and \(L_{{\textit{MO}}}\) is the Obukhov length. The turbulence at a specific scale in unstable conditions (i.e., \(\xi < 0\)) is closer to isotropic than that at the same scale under stable conditions. The maximum isotropic scale of turbulence is determined based on the characteristics of the power spectrum in three directions. Turbulence does not behave isotropically when the eddy scale is greater than the maximum isotropic scale, whereas it is horizontally isotropic at relatively large scales. The maximum isotropic scale of turbulence is compared to the outer scale of temperature, which is obtained by fitting the temperature fluctuation spectrum using the von Karman turbulent model. The results show that the outer scale of temperature is greater than the maximum isotropic scale of turbulence.  相似文献   

16.
The spatial structure of the marine atmospheric boundary layer (MABL) over the Aegean Sea is investigated using the Weather Research and Forecasting (WRF) mesoscale model. Two ‘first-order’ non-local and five ‘1.5-order’ local planetary boundary-layer (PBL) parametrization schemes are used. The predictions from the WRF model are evaluated against airborne observations obtained by the UK Facility for Airborne Atmospheric Measurements BAe-14 research aircraft during the Aegean-GAME field campaign. Statistical analysis shows good agreement between measurements and simulations especially at low altitude. Despite the differences between the predicted and measured wind speeds, they reach an agreement index of 0.76. The simulated wind-speed fields close to the surface differ substantially among the schemes (maximum values range from 13 to \(18\hbox { m s}^{-1}\) at 150-m height), but the differences become marginal at higher levels. In contrast, all schemes show similar spatial variation patterns in potential temperature fields. A warmer (1–2 K) and drier (2–3\(\hbox { g kg}^{-1})\) layer than is observed, is predicted by almost all schemes under stable conditions (eastern Aegean Sea), whereas a cooler (up to 2 K) and moister (1–2\(\hbox { g kg}^{-1})\) layer is simulated under near-neutral to nearly unstable conditions (western Aegean Sea). Almost all schemes reproduce the vertical structure of the PBL and the shallow MABL (up to 300 m) well, including the low-level jet in the eastern Aegean Sea, with non-local schemes being closer to observations. The simulated PBL depths diverge (up to 500 m) due to the different criteria applied by the schemes for their calculation. Under stable conditions, the observed MABL depth corresponds to the height above the sea surface where the simulated eddy viscosity reaches a minimum; under neutral to slightly unstable conditions this is close to the top of the simulated entrainment layer. The observed sensible heat fluxes vary from ?40 to \(25\hbox { W m}^{-2}\), while the simulated fluxes range from ?40 to \(40\hbox { W m}^{-2}\); however, all of the schemes’ predictions are close to the observations under unstable conditions. Finally, all schemes overestimate the friction velocity, although the simulated range (from 0.2 to \(0.5\hbox { m s}^{-1})\) is narrower than that observed (from 0.1 to \(0.7\hbox { m s}^{-1})\).  相似文献   

17.
We present a simple model based on already existing and widely used equations for estimating particle mass fluxes on surfaces sheltered by live vegetation. Wind-tunnel measurements of vertical profiles of mass flux in three different dense live plant canopies, and as a function of the spatially averaged skin friction velocity \({u_{\tau }}'\), provide the baseline set of data. For the bare-sand surface, the total mass flux Q shows the typical \(b({u_\tau }' - {u_{\tau t}}')^{3 }\) increase with increasing skin friction velocity \({u_{\tau }}'\), where b is a constant and \({u_{\tau t}}'\) is the threshold at the onset of particle erosion. Similar relations, however, with different values for b and \({u_{\tau t}}'\) compared to the bare-sand surface were found for experiments with 5.25 and 24.5 plants \(\hbox {m}^{-2}\) and can be explained by the spatial variations of \(u_{\tau }\) for the canopy cases. Based on the resulting parameters b and \({u_{\tau t}}'\), which are found to be functions of the roughness density \(\lambda \), we present a final simple relation \(Q(\lambda ,\, {u_{\tau }}')\) used for estimating the total mass flux for surfaces sheltered by live vegetation.  相似文献   

18.
Enhanced meteorological observations were made during the 2015 Pan and Parapan American Games in Toronto in order to measure the vertical and horizontal structure of lake-breeze events. Two scanning Doppler lidars (one fixed and one mobile), a C-band radar, and a network including 53 surface meteorological stations (mesonet) provided pressure, temperature, humidity, and wind speed and direction measurements over Lake Ontario and urban areas. These observations captured the full evolution (prior, during, and after) of 27 lake-breeze events (73% of observation days) in order to characterize the convective and dynamic processes driving lake breezes at the local scale and mesoscale. The dominant signal of a passing lake-breeze front (LBF) was an increase in dew-point temperature of \(2.3 \pm 0.3 \,^{\circ }\hbox {C}\), coinciding with a \(180^{\circ }\) shift in wind direction and a decrease in air temperature of \(2.1 \pm 0.2 \,^{\circ }\hbox {C}\). Doppler lidar observations over the lake detected lake breezes 1 hour (on average) before detection by radar and mesonet. On days with the synoptic flow in the offshore direction, the lidars observed wedge-shaped LBFs with shallow depths, which inhibited the radar’s ability to detect the lake breeze. The LBF’s ground speed and inland penetration distance were found to be well-correlated (\(r = 0.78\)), with larger inland penetration distances occurring on days with non-opposing (non-offshore) synoptic flow. The observed enhanced vertical motion \(({>} 1\hbox { m s}^{-1})\) at the LBF, observed by the lidar on 54% of lake-breeze days, was greater (at times \({>} 2.5\hbox { m s}^{-1}\)) than that observed in previous studies and longer-lasting over the lake than over land. The weaker and less pronounced lake-breeze structure over land is illustrated in two case studies highlighting the lifetime of the lake-breeze circulation and the impact of propagation distance on lake-breeze intensity.  相似文献   

19.
Stereoscopic particle image velocimetry was used to provide a three-dimensional characterization of the flow around a simplified urban model defined by a 5 by 7 array of blocks, forming four parallel streets, perpendicular to the incoming wind direction corresponding to a zero angle of incidence. Channeling of the flow through the array under consideration was observed, and its effect increased as the incoming wind direction, or angle of incidence (AOI), was changed from \(0^{\circ }\) to \(15^{\circ }\), \(30^{\circ }\), and \(45^{\circ }\). The flow between blocks can be divided into two regions: a region of low turbulence kinetic energy (TKE) levels close to the leeward side of the upstream block, and a high TKE area close to the downstream block. The centre of the arch vortex is located in the low TKE area, and two regions of large streamwise velocity fluctuation bound the vortex in the spanwise direction. Moreover, a region of large spanwise velocity fluctuation on the downstream block is found between the vortex legs. Our results indicate that the reorientation of the arch vortex at increasing AOI is produced by the displacement of the different TKE regions and their interaction with the shear layers on the sides and top of the upstream and downstream blocks, respectively. There is also a close connection between the turbulent structure between the blocks and the wind gusts. The correlations among gust components were also studied, and it was found that in the near-wall region of the street the correlations between the streamwise and spanwise gusts \(R_{uv}\) were dominant for all four AOI cases. At higher wall-normal positions in the array, the \(R_{uw}\) correlation decreased with increasing AOI, whereas the \(R_{uv}\) coefficient increased as AOI increased, and at \({\textit{AOI}}=45^{\circ }\) all three correlations exhibited relatively high values of around 0.4.  相似文献   

20.
The three turbulent velocity components, water vapour (\(\text {H}_2\text {O}\)), carbon dioxide (\(\text {CO}_{2}\)), and methane (\(\text {CH}_{4}\)) concentration fluctuations are measured above a boreal peatland and analyzed using conditional sampling and quadrant analysis. The overarching question to be addressed is to what degree lower-order cumulant expansion methods describe transport efficiency and the relative importance of ejections and sweeps to momentum, \(\text {CH}_{4}\), \(\text {CO}_{2}\) and \(\text {H}_2\text {O}\) fluxes across a range of atmospheric flow regimes. The patchy peatland surface creates distinctly different source and sink distributions for the three scalars in space and time thereby adding to the uniqueness of the set-up. The measured and modelled fractional contributions to the momentum flux show that sweep events dominate over ejections in agreement with prior studies conducted in the roughness sublayer. For scalar fluxes, ejections dominate the turbulent fluxes over sweeps. While ejective motions persist longer for momentum transport, sweeping events persist longer for all three scalars. Third-order cumulant expansions describe many of the results detailed above, and the results are surprising given the highly non-Gaussian distribution of \(\text {CH}_{4}\) turbulent fluctuations. Connections between the asymmetric contributions of sweeps and ejections and the flux-transport term arising in scalar turbulent-flux-budget closure are derived and shown to agree reasonably well with measurements. The proposed model derived here is much simpler than prior structural models used to describe laboratory experiments. Implications of such asymmetric contributions on, (i) the usage of the now proliferating relaxed-eddy-accumulation method in turbulent flux measurements, (ii) the constant-flux assumption, and (iii) gradient-diffusion closure models are presented.  相似文献   

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

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