首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
This study focuses on the behaviour of the turbulent Prandtl number, Pr t , in the stable atmospheric boundary layer (SBL) based on measurements made during the Surface Heat Budget of the Arctic Ocean experiment (SHEBA). It is found that Pr t increases with increasing stability if Pr t is plotted vs. gradient Richardson number, Ri; but at the same time, Pr t decreases with increasing stability if Pr t is plotted vs. flux Richardson number, Rf, or vs. ζ = z/L. This paradoxical behaviour of the turbulent Prandtl number in the SBL derives from the fact that plots of Pr t vs. Ri (as well as vs. Rf and ζ) for individual 1-h observations and conventional bin-averaged values of the individual quantities have built-in correlation (or self-correlation) because of the shared variables. For independent estimates of how Pr t behaves in very stable stratification, Pr t is plotted against the bulk Richardson number; such plots have no built-in correlation. These plots based on the SHEBA data show that, on the average, Pr t decreases with increasing stability and Pr t < 1 in the very stable case. For specific heights and stabilities, though, the turbulent Prandtl number has more complicated behaviour in the SBL.  相似文献   

2.
We present results of a technique for examining the scale-dependence of the gradient Richardson number, Ri, in the nighttime residual layer. The technique makes use of a series of high-resolution, in situ, vertical profiles of wind speed and potential temperature obtained during CASES-99 in south-eastern Kansas, U.S.A. in October 1999. These profiles extended from the surface, through the nighttime stable boundary layer, and well into the residual layer. Analyses of the vertical gradients of both wind speed, potential temperature and turbulence profiles over a wide range of vertical scale sizes are used to estimate profiles of the local Ri and turbulence structure as a function of scale size. The utility of the technique lies both with the extensive height range of the residual layer as well as with the fact that the sub-metre resolution of the raw profiles enables a metre-by-metre ‘sliding’ average of the scale-dependent Richardson number values over hundreds of metres vertically. The results presented here show that small-scale turbulence is a ubiquitous and omnipresent feature of the residual layer, and that the region is dynamic and highly variable, exhibiting persistent turbulent structure on vertical scales of a few tens of metres or less. Furthermore, these scales are comparable to the scales over which the Ri is less than or equal to the critical value of Ri c of 0.25, although turbulence is also shown to exist in regions with significantly larger Ri values, an observation at least consistent with the concept of hysteresis in turbulence generation and maintenance. Insofar as the important scale sizes are comparable to or smaller than the resolution of current models, it follows that, in order to resolve the observed details of small Ri values and the concomitant turbulence generation, future models need to be capable of significantly higher resolutions.  相似文献   

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

4.
Based on data obtained during the Hualhe River Basin Experiment （HUBEX） in 1999, this study intends to detect the quantitative discrepancies in the momentum （τ0）, sensible heat （H0） and latent heat （E0） fluxes among six sets of similarity functions with the aerodynamic method. It also aims to clarify the applicability of the functions under stable conditions. The relative discrepancy was studied with the normalized transfer coefficients for τ0, H0 and E0, namely CD, CH and CQ, respectively. Except for one set of functions that adopted a rather small von Kármán＇s constant （0.365）, the relative discrepancy in τ0 among the other functions was less than 10%, while that in H0（E0） sometimes reached 25% when the bulk Richardson number （R/B） was less than 0.07. The absolute discrepancy in the fluxes was studied with statistical computations. Among the six sets of functions, the discrepancy in τ0, H0 and E0 sometimes reached 0.03 kg m^-1 s^-2, 4 W m^-2 and 10 W m^-2, respectively, and the discrepancy in the energy balance ratio sometimes exceeded 0.1. Furthermore, when RiB exceeded the critical value （Ric） for a specific set of functions, no fluxes could be derived with the functions. It is therefore suggested that RiB be compared with Ric before computing the fluxes if RiB is less than Ric. Finally, two sets of nonlinear similarity functions are recommended, due to their unlimited applicability in terms of RiB.  相似文献   

5.
Extensive eddy-correlation datasets are analyzed to examine the influence of nonstationarity of the mean flow on the flux–gradient relationship near the surface. This nonstationarity is due to wavelike motions, meandering of the wind vector, and numerous unidentified small-scale mesoscale motions. While the data do not reveal an obvious critical gradient Richardson number, the maximum downward heat flux increases approximately linearly with increasing friction velocity for significant stability. The largest of our datasets is chosen to more closely examine the influence of stability, nonstationarity, distortion of the mean wind profile and self-correlation on the flux-gradient relationship. Stability is expressed in terms of z/L, the gradient Richardson number or the bulk Richardson number over the tower layer. The efficiency of the momentum transport systematically increases with increasing nonstationarity and attendant distortion of the mean wind profile. Enhancement of the turbulent momentum flux associated with nonstationarity is examined in terms of the nondimensional shear, Prandtl number and the eddy diffusivity.  相似文献   

6.
Here we advance the physical background of the energy- and flux-budget turbulence closures based on the budget equations for the turbulent kinetic and potential energies and turbulent fluxes of momentum and buoyancy, and a new relaxation equation for the turbulent dissipation time scale. The closure is designed for stratified geophysical flows from neutral to very stable and accounts for the Earth’s rotation. In accordance with modern experimental evidence, the closure implies the maintaining of turbulence by the velocity shear at any gradient Richardson number Ri, and distinguishes between the two principally different regimes: “strong turbulence” at ${Ri \ll 1}$ typical of boundary-layer flows and characterized by the practically constant turbulent Prandtl number Pr T; and “weak turbulence” at Ri > 1 typical of the free atmosphere or deep ocean, where Pr T asymptotically linearly increases with increasing Ri (which implies very strong suppression of the heat transfer compared to the momentum transfer). For use in different applications, the closure is formulated at different levels of complexity, from the local algebraic model relevant to the steady-state regime of turbulence to a hierarchy of non-local closures including simpler down-gradient models, presented in terms of the eddy viscosity and eddy conductivity, and a general non-gradient model based on prognostic equations for all the basic parameters of turbulence including turbulent fluxes.  相似文献   

7.
A new method for obtaining instantaneous vertical profiles of two components of velocity and temperature in thermally stratified turbulent shear flows is presented. In this report, the design and construction of the traversing system will be discussed and results to date will be presented. The method is based on rapid vertical sampling whereby probe sensors are moved vertically at a high speed such that the measurement is approximately instantaneous. The system is designed to collect many measurements for the calculation of statistics such as vertical wave number spectra, mean square vertical gradients, and Thorpe scales. Results are presented for vertical profiles of temperature and compared to vertical profiles measured by single-point Eulerian time averages. The quality of the vertical profiles is found to be good over many profiles. Some comparisons are made between vertical measurements and standard single-point Eulerian measurements for three cases of stably stratified turbulent shear flow in which the initial microscale Reynolds number, Reλ≈30. In case 1, the mean conditions are characterized by a gradient Richardson number, Rig=0.015, for which the flow is “unstable”, meaning the spatially evolving turbulent kinetic energy (Ek) grows. In case 2, Rig=0.095, for which the evolving turbulent kinetic energy is almost constant. In case 3, the flow is highly stable, where Rig=0.25 and Ek decays with spatial evolution. The measurements indicate anisotropy in the small scales for all cases. In particular, it is found that the ratio grows initially to a maximum and then decays with further evolution. Maximum Thorpe displacements are measured and compared to single-point measures of the vertical scales. It is found that vertical length scales derived from single-point measurements, such as the Ozmidov scale, LO=(ε/N3)1/2 and the overturn scale, Lt=θ′/(dT/dz), do not represent well the wide range of overturning scales which are actually present in the turbulence.  相似文献   

8.
The parameterization of the dimensionless entrainment rate (w e /w *) versus the convective Richardson number (Ri δθ ) is discussed in the framework of a first-order jump model (FOM). A theoretical estimation for the proportionality coefficient in this parameterization, namely, the total entrainment flux ratio, is derived. This states that the total entrainment flux ratio in FOM can be estimated as the ratio of the entrainment zone thickness to the mixed-layer depth, a relationship that is supported by earlier tank experiments, and suggesting that the total entrainment flux ratio should be treated as a variable. Analyses show that the variability of the total entrainment flux ratio is actually the effect of stratification in the free atmosphere on the entrainment process, which should be taken into account in the parameterization. Further examination of data from tank experiments and large-eddy simulations demonstrate that the different power laws for w e /w * versus Ri δθ can be interpreted as the variability of the total entrainment flux ratio. These results indicate that the dimensionless entrainment rate depends not only on the convective Richardson number but also upon the total entrainment flux ratio.  相似文献   

9.
In this study, we evaluate four different parameterizations of the turbulent Prandtl (Schmidt) number Prt = νt/Γt where νt is the eddy viscosity and Γt is the scalar eddy diffusivity, for stably stratified flows. All four formulations of Prt are strictly functions of the gradient Richardson number Ri, which provides a measure of the strength of the stratification. A zero-equation (i.e. no extra transport equations are required) turbulence model for νt in a one-dimensional, turbulent channel flow is considered to evaluate the behavior of the different formulations of Prt. Both uni-directional and oscillatory flows are considered to simulate conditions representative of practical flow problems such as atmospheric boundary layer flows and tidally driven estuarine flows, to quantify the behavior of each of the four formulations of Prt. We perform model-to-model comparisons to highlight which of the models of Prt allow for a higher rate of turbulent mixing and which models significantly inhibit turbulent mixing in the presence of buoyancy forces resulting from linear (continuous) stratification as well as two-layer stratification. The basis underlying the formulation of each model in conjunction with the simulation results are used to emphasize the considerable variability in the different formulations and the importance of choosing an appropriate parameterization of Prt given a model for νt in stably stratified flows.  相似文献   

10.
Summary Sublayer-Stanton numbers, Bi, of heat and matter for the interfacial sublayer over aerodynamically smooth surfaces determined for forced convective conditions by elementary and numerical integration are reviewed and evaluated. The results are based on Roth’s modified Heisenberg model for the spectral energy transfer in the equilibrium range under locally isotropic conditions and the approaches of Reichardt, Elser, Deissler, van Driest, Rannie, Sheppard, and Spalding for the normalized eddy diffusivity Km/v. The results substantiate that with the exception of Sheppard’s Km/v-approach all formulations are appropriate to provide sublayer-Stanton numbers with a sufficient degree of accuracy. From a theoretical point of view the Km/v-relationships of Roth, Reichardt, van Driest, and Spalding are to be preferred, when a turbulent Prandtl number Prt = 1 is presumed. Since within the framework of mesoscale meteorological modelling numerical integration techniques would consume too much CPU-time because of the large number of near-wall grid points, a parameterization formula for the sublayer-Stanton number is proposed and evaluated. Compared to the Bi −1-results obtained by numerical integration, this kind of parameterization leads to a relative error of less than 5 percent for roughness Reynolds numbers, ηr, ranging from 30 to 600. Received January 2001 Revised November 5, 2001  相似文献   

11.
The gradient-based similarity approach removes turbulent fluxes as governing parameters and replaces them with vertical gradients of mean wind speed and potential temperature. As a result, the gradient Richardson number, Ri, appears as a stability parameter instead of the Monin–Obukhov stability parameter z/L (L is the Obukhov length). The gradient-based scaling is more appropriate for moderate and very stable conditions when the gradients are large and their errors are relatively small whereas z/L becomes ambiguous in these conditions because turbulent fluxes are small. However, the gradient-based formulation is faced with a problem related to the influence of Ri outliers: outliers with high values of Ri can exist in conditions that are really near-neutral. These outliers are mapped into the very stable range in plots in which Ri is the independent variable and may lead to spurious dependencies for bin-averaged data (spurious bin-averaging). This effect is quite large for functions that are steep for the gradient-based scaling. The present study uses the Surface Heat Budget of the Arctic Ocean (SHEBA) data to examine the problem and proposes two methods, conditional analysis and independent binning, to limit the influence of outliers on bin-averaging. A disadvantage of the conditional analysis is associated with eliminating outliers based on criteria that could be considered as subjective. The independent bin-averaging method does not have this disadvantage, but the scatter of the bin-averaged points is higher than for the conditional analysis, rendering data analysis and interpretation difficult.  相似文献   

12.
Summary Field observations were carried out in order to determine the transfer coefficeient of sensible heat flux above a melting snow surface at the Moshiri experimental site. The coefficient is calculated as the ratio of the sensible heat flux determined by the eddy correlation method using a sonic anemometer to the product of the wind speed and the temperature difference between the air and snow surface. The sensible heat fluxes are also compared with the result of the precise heat balance observations. The nondimensional transfer coefficienth shows a good correlation with the atmospheric stability,Ri. The value ofh=2.3×10–3 is obtained in the range of 0<Ri<0.1, however, it is smaller and scattered under stronger stable condition (Ri>0.1). The dimensional transfer coefficent for sensible heat flux is calculated, and a linear relationship is obtained as a function of the logarithm of atmospheric stability.With 11 Figures  相似文献   

13.
Air temperature T a , specific humidity q, CO2 mole fraction χ c , and three-dimensional winds were measured in mountainous terrain from five tall towers within a 1 km region encompassing a wide range of canopy densities. The measurements were sorted by a bulk Richardson number Ri b . For stable conditions, we found vertical scalar differences developed over a “transition” region between 0.05 < Ri b < 0.5. For strongly stable conditions (Ri b > 1), the vertical scalar differences reached a maximum and remained fairly constant with increasing stability. The relationships q and χ c have with Ri b are explained by considering their sources and sinks. For winds, the strong momentum absorption in the upper canopy allows the canopy sublayer to be influenced by pressure gradient forces and terrain effects that lead to complex subcanopy flow patterns. At the dense-canopy sites, soil respiration coupled with wind-sheltering resulted in CO2 near the ground being 5–7 μmol mol−1 larger than aloft, even with strong above-canopy winds (near-neutral conditions). We found Ri b -binning to be a useful tool for evaluating vertical scalar mixing; however, additional information (e.g., pressure gradients, detailed vegetation/topography, etc.) is needed to fully explain the subcanopy wind patterns. Implications of our results for CO2 advection over heterogenous, complex terrain are discussed.  相似文献   

14.
Measurements of atmospheric turbulence made over the Arctic pack ice during the Surface Heat Budget of the Arctic Ocean experiment (SHEBA) are used to determine the limits of applicability of Monin–Obukhov similarity theory (in the local scaling formulation) in the stable atmospheric boundary layer. Based on the spectral analysis of wind velocity and air temperature fluctuations, it is shown that, when both the gradient Richardson number, Ri, and the flux Richardson number, Rf, exceed a ‘critical value’ of about 0.20–0.25, the inertial subrange associated with the Richardson–Kolmogorov cascade dies out and vertical turbulent fluxes become small. Some small-scale turbulence survives even in this supercritical regime, but this is non-Kolmogorov turbulence, and it decays rapidly with further increasing stability. Similarity theory is based on the turbulent fluxes in the high-frequency part of the spectra that are associated with energy-containing/flux-carrying eddies. Spectral densities in this high-frequency band diminish as the Richardson–Kolmogorov energy cascade weakens; therefore, the applicability of local Monin–Obukhov similarity theory in stable conditions is limited by the inequalities Ri < Ri cr and Rf < Rf cr. However, it is found that Rf cr  =  0.20–0.25 is a primary threshold for applicability. Applying this prerequisite shows that the data follow classical Monin–Obukhov local z-less predictions after the irrelevant cases (turbulence without the Richardson–Kolmogorov cascade) have been filtered out.  相似文献   

15.
利用兰州大学半干旱气候与环境观测站(Semi-Arid Climate and Environment Observatory of Lanzhou University,简称SACOL)2008年12月观测资料,研究了稳定边界层湍流特征.使用涡动相关资料研究湍流通量时,定义湍流的平均时间τ内的中尺度运动是造成湍流统计量变化范围大的主要原因,稳定情形? τ取几十秒至几分钟.对梯度理查森数大于0.3的强稳定情形的湍流尺度分解(MRD)谱分析表明,感热通量在112.4～449.9 s存在谱隙,尺度大于谱隙的中尺度运动造成了通量观测资料离散性大,甚至有支配性影响.动量通量的谱隙在112.4～224.9 s之间.弱风时,中尺度运动的影响更大,垂直风速标准差以0.1的比率随中尺度风速变化;垂直风速标准差同广义风速表现出很好的相关性,并随着广义风速消失而消失.三维风速标准差与摩擦速度呈很好的线性关系,垂直、水平、横风风速的无量纲标准差分别为1.35、2.54、2.21.对湍流动能的研究发现,在梯度理查森数大于0.3的条件下,仍然存在连续的湍流.以湍动能为依据,分析了湍流的平稳时间长度,其长度随稳定度变化而变化,2008年12月7～11日从133.5 s变化到856.2 s,湍流平稳时间长度反映了中尺度运动的发生频率.  相似文献   

16.
Local Scales of Turbulence in the Stable Boundary Layer   总被引：1，自引：1，他引：0  
Local, gradient-based scales, which contain the vertical velocity and temperature variances, as well as the potential temperature gradient, but do not include fluxes, are tested using data collected during the CASES-99 experiment. The observations show that the scaling based on the temperature variance produces relatively smaller scatter of empirical points. The resulting dimensionless statistical moments approach constant values for sufficiently large values of the Richardson number Ri. This allows one to derive predictions for the Monin–Obukhov similarity functions φ m and φ h , the Prandtl number Pr and the flux Richardson number Rf in weak turbulence regime.  相似文献   

17.
Wind speed profiles above a forest canopy relate to scalar exchange between the forest canopy and the atmosphere. Many studies have reported that vertical wind speed profiles above a relatively flat forest can be classified by a stability index developed assuming wind flow above a flat plane. However, can such a stability index be used to classify vertical wind speed profiles observed above a sloping forest at nighttime, where drainage flow often occurs? This paper examines the use of the bulk Richardson number to classify wind speed profiles observed above a sloping forest at nighttime. Wind speed profiles above a sloping forest were classified by the bulk Richardson number Ri B . Use of Ri B distinguished between drainage flow, shear flow, and transitional flow from drainage flow to shear flow. These results suggest that Ri B is useful to interpret nighttime CO2 and energy fluxes above a sloping forest. Through clear observational evidence, we also show that the reference height should be high enough to avoid drainage-flow effects when calculating Ri B .  相似文献   

18.
Based on classic iterative computation results, new equations to calculate the surface turbulent transfer coefficients are proposed, which allow for large ratios of the momentum and heat roughness lengths. Compared to the Launiainen scheme, our proposed scheme generates results closer to classical iterative computations. Under unstable stratification, the relative error in the Launiainen scheme increases linearly with increasing instability, even exceeding 15%, while the relative error of the present scheme is always less than 8.5%. Under stable stratification, the Launiainen scheme uses two equations, one for 0 < Ri B ≤ 0.08 and another for 0.08 < Ri B ≤ 0.2, and does not consider the condition that Ri B > 0.2, while its relative errors in the region 0 < Ri B ≤ 0.2 exceed 31 and 24% for momentum and heat transfer coefficients, respectively. In contrast, the present scheme uses only one equation for 0 < Ri B ≤ 0.2 and another equation for Ri B > 0.2, and the relative error of the present scheme is always less than 14%.  相似文献   

19.
We give a new derivation of the familiar linear relation for the dimensionless velocity gradient in the stably stratified surface layer and provide physical and empirical grounds for its universal applicability in stationary homogeneous turbulence over the whole range of static stabilities from Ri =  0 to very large Ri. Combining this relation with the budget equation for the turbulent kinetic energy we obtain the “equilibrium formulation” of the turbulent dissipation length scale, and recommend it for use in turbulence closure models.  相似文献   

20.
A continuous measurement of number size distributions and chemical composition of aerosol particles was conducted in Beijing in a dust storm event during 21-26 March 2001. The number concentration of coarse particles （ 〉2μm） increased more significantly than fine particles （ 〈2μm） during the dust storm due to dust weather, while the anthropogenic aerosols collected during the non-dust-storm period tended to be associated with fine particles. Elemental compositions were analyzed by using proton-induced X-ray emission （PIXE）. The results show that 20 elements in the dust storm were much higher than in the non-dust-storm period. The calculated soil dust concentration during the dust storm was, on average, 251.8μg m^-3, while it was only 52.1μg m^-3 on non-dust-storm days. The enrichment factors for Mg, A1, P, K, Ca, Ti, Mn, Fe, C1, Cu, Pb, and Zn show small variations between the dust storm and the non-dust-storm period, while those for Ca, Ni and Cr in the dust storm were much lower than those in the non-dust-storm period due to significant local emission sources. A high concentration and enrichment factor for S were observed during the dust storm, which implies that the dust particles were contaminated by aerosol particles from anthropogenic emissions during the long-range transport. A statistical analysis shows that the elemental composition of particles collected during the dust storm in Beijing were better correlated with those of desert soil colleted from desert regions in Inner Mongolia. Air mass back-trajectory analysis further confirmed that this dust storm event could be identified as streaks of dust plumes originating from Inner Mongolia.  相似文献   

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

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