首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 410 毫秒
1.
Increased heat fluxes near a forest edge   总被引：1，自引：0，他引：1  
Summary ?Observations of sensible and latent heat flux above forest downwind of a forest edge show these fluxes to be larger than the available energy over the forest. The enhancement averages to 56 W m−2, or 16% of the net radiation, at fetches less than 400 m, equivalent to fetch to height ratios less than 15. The enhancement of turbulent energy fluxes is explained by advection and increases with the difference in temperature and humidity of the air over the upwind area as compared to the forest. The relatively high temperature and humidity of the upwind air are not caused by high surface heat fluxes, but are explained by the relatively low aerodynamic roughness of the upwind surface. Although the heat fluxes over forest are enhanced, the momentum fluxes are almost adjusted to the underlying forest. The different behaviour of heat and momentum fluxes is explained by absorption of momentum by pressure gradients near the forest edge. It is concluded that fetch requirements to obtain accurate surface fluxes from atmospheric observations need to be more stringent for scalar fluxes as compared to momentum fluxes. Received November 23, 2001; accepted May 13, 2002  相似文献   

2.
We propose a new turbulence closure model based on the budget equations for the key second moments: turbulent kinetic and potential energies: TKE and TPE (comprising the turbulent total energy: TTE = TKE + TPE) and vertical turbulent fluxes of momentum and buoyancy (proportional to potential temperature). Besides the concept of TTE, we take into account the non-gradient correction to the traditional buoyancy flux formulation. The proposed model permits the existence of turbulence at any gradient Richardson number, Ri. Instead of the critical value of Richardson number separating—as is usually assumed—the turbulent and the laminar regimes, the suggested model reveals a transitional interval, , which separates two regimes of essentially different nature but both turbulent: strong turbulence at ; and weak turbulence, capable of transporting momentum but much less efficient in transporting heat, at . Predictions from this model are consistent with available data from atmospheric and laboratory experiments, direct numerical simulation and large-eddy simulation.  相似文献   

3.
A modification of the most popular two-equation (E–φ) models, taking into account the plant drag, is proposed. Here E is the turbulent kinetic energy (TKE) and φ is any of the following variables: El (product of E and the mixing length l), (dissipation rate of TKE), and ω (specific dissipation of TKE, ). The proposed modification is due to the fact that the model constants estimated experimentally for ‘free-air’ flow do not allow for adequate reconstruction of the ratio between the production and dissipation rates of TKE in the vegetation canopy and have to be adjusted. The modification is universal, i.e. of the same type for all E–φ models considered. The numerical experiments carried out for both homogeneous and heterogeneous plant canopies with E–φ models (and with the E–l model taken as a kind of reference) show that the modification performs well. They also suggest that E– and E–ω schemes are more promising than the E–El scheme for canopy flow simulation since they are not limited by the need to use a wall function.In addition, a new parameterization for enhanced dissipation within the plant canopy is derived. It minimizes the model sensitivity to C μ, the key parameter for two-equation schemes, and whose estimates unfortunately vary considerably from experiment to experiment. The comparison of results of new modified E– and E –ω models with observations from both field and wind-tunnel experiments shows that the proposed parameterization is quite robust. However, because of uncertainties with the turbulence Prandtl and Schmidt numbers for the E– model within the canopy, the E–ω model is recommended for future implementation, with the suggested modifications.  相似文献   

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

5.
The flux footprint, that is the contribution per unit emission from each element of the upwind surface area to measurement of the vertical flux of a passive scalar, is calculated for fluxes estimated by micrometeorological profile techniques. It is found that the upwind extent of the footprint for concentration-profile flux estimates is similar to that of the footprint for eddy-covariance flux measurements, when the eddy-covariance measurement is made at a height equal to the arithmetic mean of the highest and lowest profile measurement heights for stable stratification or the geometric mean for unstable stratification. The concentration-profile flux footprint depends on the ratio of the highest to the lowest measurement height, but is insensitive to the number of measurement levels. The concentration-profile flux footprint extends closer to the measurement location than does the 'equivalent eddy-covariance flux footprint, and the difference becomes more pronounced as the ratio of the profile measurement heights increases. The flux footprint for the Bowen-ratio technique is identical to that for a two-level profile measurement only for very limited circumstances. In the more general case, a flux footprint cannot be defined for the Bowen-ratio technique and the uniform upwind fetch required for representative flux measurements depends on the specific spatial distribution of surface fluxes.  相似文献   

6.
Wind and tracer data from the Oklahoma City Joint Urban 2003 (JU2003) and the Manhattan Madison Square Garden 2005 (MSG05) urban field experiments are being analyzed to aid in understanding air flow and dispersion near street-level in built-up downtown areas. The mean winds are separately calculated for groups of anemometers having similar exposures such as “near street level” and “on building top”. Several general results are found, such as the scalar wind speed at street level is about 1/3 of that at building top. Turbulent standard deviations of wind speed components and temperature, and vertical fluxes of momentum and sensible heat, are calculated from sonic anemometers near street level at 20 locations in JU2003 and five locations in MSG05, and from two rooftop locations in MSG05. The turbulence observations are consistent with observations in the literature at other cities, although the JU2003 and MSG05 data are unique in that many data are available near street level. For example, it is found that the local (i.e., at the measuring height) averages about 1.5 and the local averages about 0.25 in the two cities, where is the standard deviation of vertical velocity fluctuations, is the friction velocity, and u is the wind speed. The ratio of temperature fluctuations to temperature scale, , averages about −3 in both cities, consistent with similarity theory for slightly unstable conditions, where is the standard deviation of temperature fluctuations, and is the temperature scale. The calculated Obukhov length, L, is also consistent with slightly unstable conditions near street level, even at night during JU2003. The SF6 tracer concentration observations from JU2003 are analyzed. Values of for the continuous releases are calculated for each release and arc distance, where is the 30-min average arc maximum concentration, Q is the continuous source emission rate, and u is the spatial-averaged wind speed in the downtown area. The basic characteristics of the JU2003 plot of averaged agree reasonably well with similar plots for other urban experiments in Salt Lake City and London (i.e., at . A is found to be about 3 during the day and about 10 during the night.  相似文献   

7.
Numerical analysis of flux footprints for different landscapes   总被引：5，自引：0，他引：5  
Summary A model for the canopy – planetary boundary layer flow and scalar transport based on E- closure was applied to estimate footprint for CO2 fluxes over different inhomogeneous landscapes. Hypothetical heterogeneous vegetation patterns – forest with clear-cuts as well as hypothetical heterogeneous relief – a bell-shaped valley and a ridge covered by forest were considered. The distortions of airflow caused by these heterogeneities are shown – the upwind deceleration of the flow at the ridge foot and above valley, acceleration at the crest and the flow separation with the reversed flow pattern at lee slopes of ridge and valley. The disturbances induce changes in scalar flux fields within the atmospheric surface layer comparing to fluxes for homogeneous conditions: at a fixed height the fluxes vary as a function of distance to disturbance. Correspondingly, the flux footprint estimated from model data depends on the location of the point of interest (flux measurement point) and may significantly deviate from that for a flat terrain. It is shown that proposed method could be used for the choice of optimal sensor position for flux measurements over complex terrain as well as for the interpretation of data for existing measurement sites. To illustrate the latter the method was applied for experimental site in Solling, Germany, taking into account the complex topography and vegetation heterogeneities. Results show that in certain situations (summer, neutral stratification, south or north wind) and for a certain sensor location the assumptions of idealized air flow structure could be used for measurement interpretation at this site, though in general, extreme caution should be applied when analytical footprint models are used in the interpretation of flux measurements over complex sites.  相似文献   

8.
We describe one-dimensional (1D) simulations of the countergradient zone of mean potential temperature observed in the convective boundary layer (CBL). The method takes into account the third-order moments (TOMs) in a turbulent scheme of relatively low order, using the turbulent kinetic energy equation but without prognostic equations for other second-order moments. The countergradient term is formally linked to the third-order moments and , and a simple parameterization of these TOMs is proposed. It is validated for several cases of a dry CBL, using large-eddy simulations that have been realized from the MESO-NH model. The analysis of the simulations shows that TOMs are responsible for the inversion of the sign of in the higher part of the CBL, and budget analysis shows that the main terms responsible for turbulent fluxes and variances are now well reproduced.  相似文献   

9.
Most of our knowledge on forest-edge flows comes from numerical and wind-tunnel experiments where canopies are horizontally homogeneous. To investigate the impact of tree-scale heterogeneities (\({>}1\) m) on the edge-flow dynamics, the flow in an inhomogeneous forest edge on Falster island in Denmark is investigated using large-eddy simulation. The three-dimensional forest structure is prescribed in the model using high resolution helicopter-based lidar scans. After evaluating the simulation against wind measurements upwind and downwind of the forest leading edge, the flow dynamics are compared between the scanned forest and an equivalent homogeneous forest. The simulations reveal that forest inhomogeneities facilitate flow penetration into the canopy from the edge, inducing important dispersive fluxes in the edge region as a consequence of the flow spatial variability. Further downstream from the edge, the forest inhomogeneities accentuate the canopy-top turbulence and the skewness of the wind-velocity components while the momentum flux remains unchanged. This leads to a lower efficiency in the turbulent transport of momentum within the canopy. Dispersive fluxes are only significant in the upper canopy. Above the canopy, the mean flow is less affected by the forest inhomogeneities. The inhomogeneities induce an increase in the mean wind speed that was found to be equivalent to a decrease in the aerodynamic height of the canopy. Overall, these results highlight the importance of forest inhomogeneities when looking at canopy–atmosphere exchanges in forest-edge regions.  相似文献   

10.
Large-eddy simulations of the neutrally stratified flow over the Askervein Hill were performed, to improve the knowledge of the flow obtained from field measurements and numerical simulations with Reynolds averaged Navier-Stokes (RANS) methods. A Lagrangian dynamic subgrid model was used but, to avoid the underdissipative character near the ground, it was merged with a damped Smagorinsky model. Simulations of a flat boundary-layer flow with this subgrid model showed that the turbulent vertical motions and shear stress were better resolved using grids with a stream to spanwise aspect ratio Δx / Δy = 2 than with an aspect ratio Δx / Δy = 1. Regarding the flow over the Askervein Hill, it was found that large-eddy simulations provide an acceptable solution for the mean-velocity field and better predictions of the turbulent kinetic energy in the upstream side of the hill than the model. However, as with the model, grid convergence was not achieved in the lee side and the size of the zone with reversed flow increased with the grid refinement. Nevertheless, the existence of the intermittent separation predicted with unsteady RANS in part one of this work seems unquestionable, due to the deceleration of the flow. In our opinion, a better modelling of the decelerating boundary layer in the lee side is required to improve the results obtained using equilibrium assumptions and achieve grid convergence.  相似文献   

11.
This paper reports on measurements of sensible and latent heat and CO2 fluxes made over an irrigated potato field, growing next to a patch of desert. The study was conducted using two eddy correlation systems. One measurement system was located within the equilibrium boundary layer 800 m downwind from the edge of the potato field. The other measurement system was mobile and was placed at various downwind positions to probe the horizontal transition of vertical scalar fluxes. Latent (LE) and sensible (H) heat fluxes, measured at 4 m above the surface, exhibited marked variations with downwind distance over the field. Only after the fetch to height ratio exceeded 75 to 1 didLE andH become invariant with downwind distance. When latent and sensible heat fluxes were measured upwind of this threshold, significant advection of humidity-deficit occurred, causing a vertical flux divergence ofH andLE.The measured fluxes of momentum, heat, and moisture were compared with predictions from a second-order closure two-dimensional atmospheric boundary layer model. There is good agreement between measurements and model predictions. A soil-plant-atmosphere model was used to examine nonlinear feedbacks between humidity-deficits, stomatal conductance and evaporation. Data interpretation with this model revealed that the advection of hot dry air did not enhance surface evaporation rates near the upwind edge of the potato field, because of negative feedbacks among stomatal conductance, humidity-deficits, andLE. This finding is consistent with results from several recent studies.  相似文献   

12.
The turbulent flow inside dense canopies is characterized by wake production and short-circuiting of the energy cascade. How these processes affect passive scalar concentration variability in general and their spectral properties in particular remains a vexing problem. Progress on this problem is frustrated by the shortage of high resolution spatial concentration measurements, and by the lack of simplified analytical models that connect spectral modulations in the turbulent kinetic energy (TKE) cascade to scalar spectra. Here, we report the first planar two-dimensional scalar concentration spectra (ϕ cc ) inside tall canopies derived from flow visualization experiments. These experiments were conducted within the deeper layers of a model canopy composed of densely arrayed cylinders welded to the bottom of a large recirculating water channel. We found that in the spectral region experiencing wake production, the ϕ cc exhibits directional scaling power laws. In the longitudinal direction (x), or the direction experiencing the largest drag force, the ϕ cc (k x ) was steeper than and followed an approximate at wavenumbers larger than the injection scale of wake energy, where k x is the longitudinal wavenumber. In the lateral direction (y), the spectra scaled as up to the injection scale, and then decayed at an approximate power law. This departure from the classical inertial subrange scaling (i.e., k −5/3) was reproduced using a newly proposed analytical solution to a simplified scalar spectral budget equation. Near the velocity viscous dissipation range, the scalar spectra appear to approach an approximate k −3, a tantalizing result consistent with dimensional analysis used in the inertial-diffusive range. Implications to subgrid modelling for large-eddy simulations (LES) inside canopies are briefly discussed.  相似文献   

13.
Time Scales in the Unstable Atmospheric Surface Layer   总被引：2，自引：2，他引：0  
Calculation of eddy covariances in the atmospheric surface layer (ASL) requires separating the instantaneous signal into mean and fluctuating components. Since the ASL is not statistically stationary, an inherent ambiguity exists in defining the mean quantities. The present study compares four methods of calculating physically relevant time scales in the unstable ASL that may be used to remove the unsteady mean components of instantaneous time signals, in order to yield local turbulent fluxes that appear to be statistically stationary. The four mean-removal time scales are: (t c ) based on the location of the maximum in the ogive of the heat flux cospectra, () the location of the zero crossing in the multiresolution decomposition of the heat flux, (t *) the ratio of the mixed-layer depth over the convective velocity, and () the convergence time of the vertical velocity and temperature variances. The four time scales are evaluated using high quality, three-dimensional sonic anemometry data acquired at the Surface Layer Turbulence and Environmental Science Test (SLTEST) facility located on the salt flats of Utah’s western desert. Results indicate that and , with t c achieving values about 2–3 times greater than t *. The sensitivity of the eddy covariances to the mean-removal time scale (given a fixed 4-h averaging period during midday) is also demonstrated.  相似文献   

14.
A comprehensive model for the prediction of concentration fluctuations in plumes dispersing in the complex and highly disturbed wind flows in an urban environment is formulated. The mean flow and turbulence fields in the urban area are obtained using a Reynolds-averaged Navier-Stokes (RANS) flow model, while the standard k-ϵ turbulence model (k is the turbulence kinetic energy and ϵ is the viscous dissipation rate) is used to close the model. The RANS model provides a specification of the velocity statistics of the highly disturbed wind flow in the urban area, required for the solution of the transport equations for the mean concentration and concentration variance (both of which are formulated in the Eulerian framework). A physically-based formulation for the scalar dissipation time scale t d , required for the closure of the transport equation for , is presented. This formulation relates t d to an inner time scale corresponding to “internal” concentration fluctuation associated with relative dispersion, rather than an outer time scale associated with the entire portion of the fluctuation spectrum. The two lowest-order moments of concentration ( and ) are used to determine the parameters of a pre-chosen functional form for the concentration probability density function (clipped-gamma distribution). Results of detailed comparisons between a water-channel experiment of flow and dispersion in an idealized obstacle array and the model predictions for mean flow, turbulence kinetic energy, mean concentration, concentration variance, and concentration probability density function are presented.  相似文献   

15.
We address some of the methodological challenges associated with the measurement of turbulence and use of scintillometers in the urban roughness sublayer (RSL). Two small-aperture scintillometers were located near the roof interface in a densely urbanized part of Basel, Switzerland, as part of the Basel Urban Boundary-Layer Experiment (BUBBLE) in the summer of 2002. Eddy correlation instruments were co-located near the mid-point of each scintillometer path for data verification purposes. The study presents the first values of the inner length scale of turbulence (l 0) and the refractive index structure parameter of air for a city and demonstrates the influence of mechanical driven turbulence on dissipation. Comparison of dissipation values determined from the two approaches show large scatter that is possibly due to the spatial inhomogeneity of the turbulence statistics within the RSL. Velocity and temperature spectra display a −2/3 slope in the inertial subrange, although the spectral ratio is less than the theoretical prediction of 4/3 expected for isotropy. Conventional Monin–Obukhov equations used to calculate fluxes from the scintillometer were replaced with urban forms of the equations. The results suggest that the scintillometer may be an appropriate tool for the measurement of sensible heat flux (Q H ) above the rooftops given a suitable determination of the effective measurement height.  相似文献   

16.
To investigate how velocity variances and spectra are modified by the simultaneous action of topography and canopy, two flume experiments were carried out on a train of gentle cosine hills differing in surface cover. The first experiment was conducted above a bare surface while the second experiment was conducted within and above a densely arrayed rod canopy. The velocity variances and spectra from these two experiments were compared in the middle, inner, and near-surface layers. In the middle layer, and for the canopy surface, longitudinal and vertical velocity variances () were in phase with the hill-induced spatial mean velocity perturbation (Δu) around the so-called background state (taken here as the longitudinal mean at a given height) as predicted by rapid distortion theory (RDT). However, for the bare surface case, and remained out of phase with Δu by about L/2, where L is the hill half-length. In the canopy layer, wake production was a significant source of turbulent energy for , and its action was to re-align velocity variances with Δu in those layers, a mechanism completely absent for the bare surface case. Such a lower ‘boundary condition’ resulted in longitudinal variations of to be nearly in phase with Δu above the canopy surface. In the inner and middle layers, the spectral distortions by the hill remained significant for the background state of the bare surface case but not for the canopy surface case. In particular, in the inner and middle layers of the bare surface case, the effective exponents derived from the locally measured power spectra diverged from their expected  − 5/3 value for inertial subrange scales. These departures spatially correlated with the hill surface. However, for the canopy surface case, the spectral exponents were near  − 5/3 above the canopy though the minor differences from  − 5/3 were also correlated with the hill surface. Inside the canopy, wake production and energy short-circuiting resulted in significant departures from  − 5/3. These departures from  − 5/3 also appeared correlated with the hill surface through the wake production contribution and its alignment with Δu. Moreover, scales commensurate with Von Karman street vorticies well described wake production scales inside the canopy, confirming the important role of the mean flow in producing wakes. The spectra inside the canopy on the lee side of the hill, where a negative mean flow delineated a recirculation zone, suggested that the wake production scales there were ‘broader’ when compared to their counterpart outside the recirculation zone. Inside the recirculation zone, there was significantly more energy at higher frequencies when compared to regions outside the recirculation zone.  相似文献   

17.
A linear model for neutral surface-layer flow over orography is presented. The Reynolds-Averaged Navier-Stokes and E– turbulence closure equations are expressed in a terrain-following coordinate system created from a simple analytical expression in the Fourier domain. The perturbation equations are solved spectrally horizontally and by numerical integration vertically. Non-dimensional solutions are stored in look-up tables for quick re-use. Model results are compared to measurements, as well as other authors’ flow models in three test cases. The model is implemented and tested in two-dimensional space; the equations for a full three-dimensional version are presented.  相似文献   

18.
The note presents a rational approach to modelling the source/sink due to vegetation or buoyancy effects that appear in the turbulent kinetic energy, E, equation and a supplementary equation for a length-scale determining variable, φ, when two-equation closure is applied to canopy and atmospheric boundary-layer flows. The approach implements only standard model coefficients C φ1 and C φ2 in the production and destruction terms of the φ equation, respectively. Numerical tests illustrate the practical applicability of the method, where, for example, simulations with the E–ω model (where is the specific dissipation and is the dissipation rate of E) properly reproduce both the surface-layer wind profile estimated from the Monin-Obukhov similarity theory and the mixing-height evolution observed above forested terrain in Southern Finland.  相似文献   

19.
The importance of soil moisture inputs and improved model physics in the prediction of the daytime boundary-layer structure during the Southern Great Plains Hydrology Experiment 1997 (SGP97) is investigated using the non-hydrostatic fifth-generation Pennsylvania State University/National Center for Atmospheric Research (PSU/NCAR) Mesoscale Model MM5. This is Part II of a two-part study examining the relationship of surface heterogeneity to observed boundary-layer structure. Part I focuses on observations and utilizes a simple model while Part II uses observations and MM5 modelling. Soil moisture inputs tested include a lookup table based on soil type and season, output from an offline land-surface model (LSM) forced by atmospheric observations, and high-resolution ( 800 m) airborne microwave remotely sensed data. Model physics improvements are investigated by comparing an LSM directly coupled with the MM5 to a simpler force-restore method at the surface. The scale of land surface heterogeneities is compared to the scale of their effects on boundary-layer structure.The use of more detailed soil moisture fields allowed the MM5 to better represent the large-scale (hundreds of km) and small-scale (tens of km) horizontal gradients in surface-layer weather and, to a lesser degree, the atmospheric boundary-layer (ABL) height, which was evaluated against observations measured by differential absorption lidar (DIAL). The benefits of coupling an LSM to the MM5 were not readily evident in this summertime case, with the model having particular difficulty simulating the timing of maximum surface fluxes while underestimating the depth of the mixed layer.  相似文献   

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

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

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