首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 234 毫秒
1.
A Second-Order Closure for Neutrally Stratified Vegetative Canopy Flows   总被引：1，自引：1，他引：0  
An existing second-order closure model is modified to include the effects on mean and turbulent motions of form and viscous drag in vegetative canopies. The additional physical mechanisms represented by the closure are viscous and pressure drag on canopy elements, their role in momentum absorption, in the creation of fine scale turbulent eddies and in enhancing the total viscous dissipation in the canopy airspace. Viscous dissipation is split into a standard 'isotropic contribution associated with the spectral eddy cascade and a foliage contribution associated with work against pressure and viscous drag on the foliage. Changes in the turbulent time scale that result from these mechanisms are included in the standard parameterisations of third moments and of the eddy cascade contribution to dissipation. The model is tested against a wind- tunnel 'wheat canopy, a corn canopy and a eucalypt canopy, a height range from 50 mm to 12.6 m. Model results show that the parameterisations of foliage interaction used in the closure are sufficiently robust to reproduce second-moment profiles within and above vegetative canopies to a high degree of accuracy without resorting to 'tuning of the model constants. The model also shows the natural emergence of two length scales, one associated with the familiar eddy cascade isotropic contribution to total dissipation and the other associated with the length scales of the canopy elements.  相似文献   

2.
Second-order closure models for the canopy sublayer (CSL) employ aset of closure schemes developed for `free-air' flow equations andthen add extra terms to account for canopy related processes. Muchof the current research thrust in CSL closure has focused on thesecanopy modifications. Instead of offering new closure formulationshere, we propose a new mixing length model that accounts for basicenergetic modes within the CSL. Detailed flume experiments withcylindrical rods in dense arrays to represent a rigid canopy areconducted to test the closure model. We show that when this lengthscale model is combined with standard second-order closureschemes, first and second moments, triple velocity correlations,the mean turbulent kinetic energy dissipation rate, and the wakeproduction are all well reproduced within the CSL provided thedrag coefficient (CD) is well parameterized. The maintheoretical novelty here is the analytical linkage betweengradient-diffusion closure schemes for the triple velocitycorrelation and non-local momentum transfer via cumulant expansionmethods. We showed that second-order closure models reproducereasonably well the relative importance of ejections and sweeps onmomentum transfer despite their local closure approximations.Hence, it is demonstrated that for simple canopy morphology (e.g.,cylindrical rods) with well-defined length scales, standard closureschemes can reproduce key flow statistics without much revision.When all these results are taken together, it appears that thepredictive skills of second-order closure models are not limitedby closure formulations; rather, they are limited by our abilityto independently connect the drag coefficient and the effectivemixing length to the canopy roughness density. With rapidadvancements in laser altimetry, the canopy roughness densitydistribution will become available for many terrestrialecosystems. Quantifying the sheltering effect, the homogeneity andisotropy of the drag coefficient, and more importantly, thecanonical mixing length, for such variable roughness density isstill lacking.  相似文献   

3.
The aim of this study is to quantify the impact of turbulence closure on the simulation of surface air temperature at screen height (1.5 m) over Belgium. The mesoscale model MAR (Modèle Atmosphérique Régional), developed at the Université catholique de Louvain, is used to examine one-dimensional situations. A new second-order closure (level 2.5) is implemented containing prognostic equations for all three velocity variances, and diagnostic or prognostic formulations for the dissipation. This closure is compared with first and one-and-a-half order closures. Idealized nearly-neutral and convective cases underline the differences between first and second-order closures, and between diagnostic and prognostic equations for the dissipation. The one-and-a-half and second-order closures give satisfying results, but the first-order closure produces generally less appropriate vertical diffusion. Observed clear sky and weak horizontal advection situations have shown the sensitivity of 24 h temperature evolution to the choice of the turbulent closure.  相似文献   

4.
An Investigation of Higher-Order Closure Models for a Forested Canopy   总被引：11，自引：10，他引：1  
Simultaneous triaxial sonic anemometer velocity measurements vertically arrayed at six levels within and above a uniform pine forest were used to examine two parameterization schemes for the triple-velocity correlation tensor employed in higher-order closure models. These parameterizations are the gradient-diffusion approximation typically used in second-order closure models, and the full budget for the triple-velocity correlation tensor typically employed in third-order closure models. Both second- and third-order closure models failed to reproduce the measured profiles of the triple-velocity correlation within and above the canopy. However, the Reynolds stress tensor profiles (including velocity variances) deviated greatly from the measurements only within the lower levels of the canopy. It is shown that the Reynolds stresses are most sensitive to the parameterization of the triple-velocity correlation in these lower canopy regions where local turbulent production is negligible and turbulence is mainly sustained by the flux transport term. The failure of the third-order closure model to reproduce the measured third moments in the upper layers of the canopy-top contradicts conclusions from a previous study over shorter vegetation but agrees with another study for a deciduous forest. Whether the third-order closure model failure is due to the zero-fourth-cumulant closure approximation is therefore considered. Comparisons between measured and predicted quadruple velocity correlations suggest that the zero-fourth-cumulant approximation is valid close to the canopy-atmosphere in agreement with recent experiments.  相似文献   

5.
Turbulence structure in a deciduous forest   总被引：5，自引：2，他引：5  
Three-dimensional wind velocity components were measured at two levels above and at six levels within a fully-leafed deciduous forest. Greatest shear occurs in the upper 20% of the canopy, where over 70% of the foliage is concentrated. The turbulence structure inside the canopy is characterized as non-Gaussian, intermittant and highly turbulent. This feature is supported by large turbulence intensities, skewness and kurtosis values and by the large infrequent sweeps and ejections that dominate tangential momentum transfer. Considerable day/night differences were observed in the vertical profiles of the mean streamwise wind velocity and turbulence intensities since the stability of the nocturnal boundary layer dampens turbulence above and within the canopy.  相似文献   

6.
An extended Lagrangian stochastic dispersion model that includes time variations of the turbulent kinetic energy dissipation rate is proposed. The instantaneous dissipation rate is described by a log-normal distribution to account for rare and intense bursts of dissipation occurring over short durations. This behaviour of the instantaneous dissipation rate is consistent with field measurements inside a pine forest and with published dissipation rate measurements in the atmospheric surface layer. The extended model is also shown to satisfy the well-mixed condition even for the highly inhomogeneous case of canopy flow. Application of this model to atmospheric boundary-layer and canopy flows reveals two types of motion that cannot be predicted by conventional dispersion models: a strong sweeping motion of particles towards the ground, and strong intermittent ejections of particles from the surface or canopy layer, which allows these particles to escape low-velocity regions to a high-velocity zone in the free air above. This ejective phenomenon increases the probability of marked fluid particles to reach far regions, creating a heavy tail in the mean concentration far from the scalar source.  相似文献   

7.
The two-scalar covariance budget is significant within the canopy sublayer (CSL) given its role in modelling scalar flux budgets using higher-order closure principles and in estimating the segregation ratio for chemically reactive species. Despite its importance, an explicit expression describing how the two-scalar covariance is modified by inhomogeneity in the flow statistics and in the vertical variation in scalar emission or uptake rates within the canopy volume remains elusive even for passive scalars. To progress on a narrower version of this problem, an analytical solution to the two-scalar covariance budget in the CSL is proposed for the most idealized flow conditions: a stationary and planar homogeneous flow inside a uniform and dense canopy with a constant leaf area density distribution. The foliage emission (or uptake) source strengths are assumed to vary exponentially with depth while the forest floor emission is represented as a scalar flux. The analytical solution is a superposition of a homogeneous part that describes how the two-scalar covariance at the canopy top is transported and dissipated within the canopy volume, and an inhomogeneous part governed by local production mechanisms of the two-scalar covariance. The homogeneous part is primarily described by the canopy adjustment length scale, and the attenuation coefficients of the turbulent kinetic energy and the mean velocity. Conditions for which the vertical variation of the two-scalar covariance is controlled by the rapid attenuation in the mean velocity and turbulent kinetic energy profiles, vis-à-vis the vertical variation of the scalar source strength, are explicitly established. This model also demonstrates how dissimilarity in the emissions from the ground, even for the extreme binary case with one scalar turned ‘on’ and the other scalar turned ‘off’, modifies the vertical variation of the two-scalar covariance within the CSL. To assess its applicability to field conditions, the analytical model predictions were compared with observations made at two different forest types—a sparse pine forest at the Hyytiälä SMEAR II-station (in Finland) and a dense alpine hardwood forest at Lavarone (in Italy). While the model assumptions do not represent the precise canopy morphology, attenuation properties of the turbulent kinetic energy and the mean velocity, observed mixing length, and scalar source attenuation properties for these two forest types, good agreement was found between measured and modelled two scalar covariances for multiple scalars and for the triple moments at the Hyytiälä site.  相似文献   

8.
A semi-analytical method for describing the mean wind profile and shear stress within plant canopies and for estimating the roughness length and the displacement height is presented. This method incorporates density and vertical structure of the canopy and includes simple parameterizations of the roughness sublayer and shelter factor. Some of the wind profiles examined are consistent with first-order closure techniques while others are consistent with second-order closure techniques. Some profiles show a shearless region near the base of the canopy; however, none displays a secondary maximum there. Comparing several different analytical expressions for the canopy wind profile against observations suggests that one particular type of profile (an Airy function which is associated with the triangular foliage surface area density distribution) is superior to the others. Because of the numerical simplicity of the methods outlined, it is suggested that they may be profitably used in large-scale models of plant-atmosphere exchanges.  相似文献   

9.
Observations of wind statistics within and above a Scots pine forest are comparedwith those predicted from an analytical second-order closure model. The roughnesssublayer (RSL) effects, and the influence of stability on similarity functions, arestudied using observations. The commonly accepted forms of similarity functionsdescribe the influence of diabatic effects above the RSL well. According to earlierstudies they are expected also to apply within the RSL. As an exception, the averagewind speed normalised with friction velocity was found to be invariant with stabilityclose to the canopy top under unstable conditions. Lagrangian stochastic trajectorysimulations were used to evaluate the influence of canopy turbulence profiles onfootprint prediction. The main uncertainty was found to arise from parameterisationof the random forcing term in the Lagrangian velocity equation. The influence ofdiabatic conditions was studied, and it was found that thermal stability affectssignificantly the footprint function above the forest canopy, but significantuncertainty exists because of uncertainties in the formulation of stability functions.  相似文献   

10.
Second- and third-order turbulence closure models have met with mixed success when applied to the prediction of turbulent flows within and above plant canopies and model predictions are typically no better than those obtained using simpler two-equation (k-e){(k-\varepsilon)} models. This is because the local gradient diffusion approximation (a crucial model requirement) cannot represent accurately turbulent transport that is dominated by the presence of ejections and sweeps whose length scales are comparable with the canopy height. To make progress, turbulent transport must be treated without approximation, as in Lagrangian probability density function (PDF) models. This study is the first to develop and to validate a PDF model of horizontally-homogeneous canopy flow. The model relies upon a prescribed length-scale that has been used elsewhere in the modelling of turbulent flows. Model predictions compare favourably with measurements of neutrally stratified turbulent flows within and above canopies of mature corn and forested eucalypt trees.  相似文献   

11.
Turbulence Spectra And Dissipation Rates Above And Within A Forest Canopy   总被引：4，自引：0，他引：4  
Three velocity componentsand temperature were measured usingthree-dimensional sonic anemometers/thermometers attwo levels, above and within a forest canopy, in theChangbai Mountains of northeast China. Turbulencespectral structure, local isotropy anddissipation rates above and within the forest canopywere calculated using the eddy correlation method.Results show that the normalized turbulent spectralcurves have -2/3 slopes in the inertial subrange.While the shapes of the spectra are in good agreementwith the Kansas flat terrain results, the atmosphericturbulence is anisotropic above the forest canopy. Dueto breaking down of large eddies by the foliage,branches and trunks, the spectral peak frequencies forvelocity and temperature are higher withinthan above the forest canopy. Compared withmeasurements from previous studies over flat terrain,the velocity and temperature spectra above andinside the forest canopy appear to shift toward higherfrequencies. The turbulence is approximately isotropicin the inertial subrange within the forest canopy, and isanisotropic above the forest canopy. The turbulentkinetic energy and heat energy dissipation rates aboveand inside the forest canopy are much larger thanthose obtained by Kaimal and Hogstrom over grasslandand grazing land. The distinct features in the resultsof the present experiment may be attributed to thedynamic forcing caused by the rough surface of the forestcanopy.  相似文献   

12.
The knowledge of the concentration probability density function (pdf) is of importance in a number of practical applications, and a Lagrangian stochastic (LS) pdf model has been developed to predict statistics and concentration pdf generated by continuous releases of non-reactive and reactive substances in canopy generated turbulence. Turbulent dispersion is modelled using a LS model including the effects of wind shear and along-wind turbulence. The dissipation of concentration fluctuations associated with turbulence and molecular diffusivity is simulated by an Interaction by Exchange with the Conditional Mean (IECM) micromixing model. A general procedure to obtain the micromixing time scale needed in the IECM model useful in non-homogeneous conditions and for single and multiple scalar sources has been developed. An efficient algorithm based on a nested grid approach with particle splitting, merging techniques and time averaging has been used, thus allowing the calculation for cases of practical interest. The model has been tested against wind-tunnel experiments of single line and multiple line releases in a canopy layer. The approach accounted for chemical reactions in a straightforward manner with no closure assumptions, but here the validation is limited to non-reacting scalars.  相似文献   

13.
Particle image velocimetry (PIV) data obtained in a wind-tunnel model of a canopy boundary layer is used to examine the characteristics of mean flow and turbulence. The vector spacing varies between 1.7 and 2.5 times the Kolmogorov scales. Conditional sampling based on quadrants, i.e. based on the signs of velocity fluctuations, reveals fundamental differences in flow structure, especially between sweep and ejection events, which dominate the flow. During sweeps, the downward flow generates a narrow, highly turbulent, shear layer containing multiple small-scale vortices just below canopy height. During ejections, the upward flow expands this shear layer and the associated small-scale flow structures to a broad region located above the canopy. Consequently, during sweeps the turbulent kinetic energy (TKE), Reynolds stresses, as well as production and dissipation rates, have distinct narrow peaks just below canopy height, whereas during ejections these variables have broad maxima well above the canopy. Three methods to estimate the dissipation rate are compared, including spectral fits, measured subgrid-scale (SGS) energy fluxes at different scales, and direct measurements of slightly underresolved instantaneous velocity gradients. The SGS energy flux is 40–60% of the gradient-based (direct) estimates for filter sizes inside the inertial range, while decreasing with scale, as expected, within the dissipation range. The spectral fits are within 5–30% of the direct estimates. The spectral fits exceed the direct estimates near canopy height, but are lower well above and below canopy height. The dissipation rate below canopy height increases with velocity magnitude, i.e. it has the highest values during sweep and quadrant 1 events, and is significantly lower during ejection and quadrant 3 events. Well above the canopy, ejections are the most dissipative. Turbulent transport during sweep events acts as a source below the narrow shear layer within the canopy and as a sink above it. Transport during ejection events is a source only well above the canopy. The residual term in the TKE transport equation, representing mostly the effect of pressure–velocity correlations, is substantial only within the canopy, and is dominated by sweeps.  相似文献   

14.
Extensive turbulence measurements from the Limagne and Beauce experiments were used to compute a characteristic time scale of the turbulence field (Τ = second moment/dissipation rate) for turbulent kinetic energy, temperature and humidity variances, and temperature-humidity covariance. The height variations of these time scales were analysed. The characteristic half-time scale Τ/2 of the turbulent velocity field was found, as expected, to be of the same order of magnitude as the large-eddy time scale Τ L = Zi/w*, showing that the turbulence structure is controlled by large eddies in the bulk of the mixed layer. The increase of Τ/2 above z/Z i ～- 0.7 implies, however, that this time scale is no longer relevant to destruction of turbulent kinetic energy in the statically stable region with negative heat fluxes. An effective time scale Τeff, introduced by Zeman (1975), has been computed and its behaviour discussed. The scales for θ′ 2, q′2, and θ′q′ were found to be much shorter than Τ. Furthermore, a significant difference in behaviour was also revealed between the characteristic time scales of temperature and humidity fields in the stable layer. By using these experimental estimates, we tested some of the models for molecular dissipations, which are currently in use in higher order closure atmospheric boundary-layer models. The parameterized dissipation rates for θ′ 2, and q′ 2 agree well qualitatively with experimental estimates in the bulk of the mixed layer. In the stable layer, however, the parameterized dissipation rate ε θ tends to become larger than the experimental ones although the parameterized dissipation rate ε q still agrees with the experimental ones. For the molecular dissipation of θ′q′, this current model becomes physically inconsistent in the middle part of the mixed layer, because this term may become a production term for temperature-humidity covariance.  相似文献   

15.
The multifractality of energy and thermal dissipation of fully developed intermittent turbulence is investigated in the urban canopy layer under unstable conditions by the singularity spectrum for the fractal dimensions of sets of singularities characterizing multifractals. In order to obtain high-order moment properties of smallscale turbulent dissipation in the inertial range, an ultrasonic anemometer with a high sampling frequency of 100 Hz was used. The authors found that the turbulent signal could be singular everywhere. Moreover, the singular exponents of energy and thermal dissipation rates are most frequently encountered at around 0.2, which is significantly smaller than the singular exponents for a wind tunnel at a moderate Reynolds number. The evidence indicates a higher intermittency of turbulence in the urban canopy layer at a high Reynolds number, which is demonstrated by the data with high temporal resolution. Furthermore, the temperature field is more intermittent than the velocity field. In addition, a large amount of samples could be used for verification of the results.  相似文献   

16.
A one-dimensional analytical model that predicts foliage CO2 uptake rates, turbulent fluxes, and mean concentration throughout the roughness sub-layer (RSL), a layer that extends from the ground surface up to 5h, where h is canopy height, is proposed. The model combines the mean continuity equation for CO2 with first-order closure principles for turbulent fluxes and simplified physiological and radiative transfer schemes for foliage uptake. This combination results in a second-order ordinary differential equation in which soil respiration (R) and CO2 concentration well above the RSL are imposed as lower and upper boundary conditions, respectively. An inverse version of the model was tested against datasets from two contrasting ecosystems: a tropical forest (h = 40m) and a managed irrigated rice canopy (h = 0.7m), with good agreement noted between modelled and measured mean CO2 concentration profiles within the entire RSL. Sensitivity analysis on the model parameters revealed a plausible scaling regime between them and a dimensionless parameter defined by the ratio between external (R) and internal (stomatal conductance) characteristics controlling the CO2 exchange process. The model can be used to infer the thickness of the RSL for CO2 exchange, the inequality in zero-plane displacement between CO2 and momentum, and its consequences on modelled CO2 fluxes. A simplified version of the solution is well suited for being incorporated into large-scale climate models. Furthermore, the model framework here can be used to a priori estimate relative contributions from the soil surface and the atmosphere to canopy-air CO2 concentration, thereby making it synergetic to stable isotopes studies.  相似文献   

17.
  相似文献   

18.
The processes influencing turbulence in a deciduous forest and the relevant length and time scales are investigated with spectral and cross-correlation analysis. Wind velocity power spectra were computed from three-dimensional wind velocity measurements made at six levels inside the plant canopy and at one level above the canopy. Velocity spectra measured within the plant canopy differ from those measured in the surface boundary layer. Noted features associated with the within-canopy turbulence spectra are: (a) power spectra measured in the canopy crown peak at higher wavenumbers than do those measured in the subcanopy trunkspace and above the canopy; (b) peak spectral values collapse to a relatively universal value when scaled according to a non-dimensional frequency comprised of the product of the natural frequency and the Eulerian time scale for vertical velocity; (c) at wavenumbers exceeding the spectral peak, the slopes of the power spectra are more negative than those observed in the surface boundary layer; (d) Eulerian length scales decrease with depth into the canopy crown, then increase with further depth into the canopy; (e) turbulent events below crown closure are more correlated with turbulent events above the canopy than are those occurring in the canopy crown; and (f) Taylor's frozen eddy hypothesis is not valid in a plant canopy. Interactions between plant elements and the mean wind and turbulence alter the processes that produce, transport and remove turbulent kinetic energy and account for the noted observations.  相似文献   

19.
An integrated canopy micrometeorological model is described for calculating CO2, water vapor and sensible heat exchange rates and scalar concentration profiles over and within a crop canopy. The integrated model employs a Lagrangian random walk algorithm to calculate turbulent diffusion. The integrated model extends previous Lagrangian modelling efforts by employing biochemical, physiological and micrometeorological principles to evaluate vegetative sources and sinks. Model simulations of water vapor, CO2 and sensible heat flux densities are tested against measurements made over a soybean canopy, while calculations of scalar profiles are tested against measurements made above and within the canopy. The model simulates energy and mass fluxes and scalar profiles above the canopy successfully. On the other hand, model calculations of scalar profiles inside the canopy do not match measurements.The tested Lagrangian model is also used to evaluate simpler modelling schemes, as needed for regional and global applications. Simple, half-order closure modelling schemes (which assume a constant scalar profile in the canopy) do not yield large errors in the computation of latent heat (LE) and CO2 (F c ) flux densities. Small errors occur because the source-sink formulation of LE andF c are relatively insensitive to changes in scalar concentrations and the scalar gradients are small. On the other hand, complicated modelling frames may be needed to calculate sensible heat flux densities; the source-sink formulation of sensible heat is closely coupled to the within-canopy air temperature profile.  相似文献   

20.
The turbulent flow in and above plant canopies is of fundamental importance to the understanding oftransport processes of momentum,heat and mass between plant canopies and atmosphere,and to microme-teorology.The Reynolds stress equation model(RSM)has been applied to calculate the turbulence in cano-pies in this paper.The calculated mean wind velocity profiles,Reynolds stress,turbulent kinetic energy andviscous dissipation rate in a corn canopy and a spruce forest are compared with field observed data and withWilson's and Shaw's model.The velocity profiles and Rynolds stress calculated by both models are in goodagreement,and the length scale of turbulence appears to be similar.  相似文献   

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

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