首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 625 毫秒
1.
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.  相似文献   

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

3.
The empirical dependence of turbulence Prandtl number (Pr) on gradient Richardson number (Ri) is presented, derived so as to avoid the effects of self-correlation from common variables. Linear power relationships between the underlying variables that constitute both Pr and Ri are derived empirically from flux and profile observations. Pr and Ri are then reconstructed from these power laws, to indicate their interdependence whilst avoiding self-correlation. Data are selected according to the stability range prior to regression, and the process is iterated from neutral to higher stability until error analysis indicates the method is no longer valid. A Butterworth function is fitted to the resulting Pr −1(Ri) regression to give an empirical summary of the analysis. The form suggests that asymptotically Pr −1 decreases as Ri 3/2. Scatter in the data increases above Ri ~ 1, however, indicating additional constraints to Pr are not captured by Ri alone in this high stability regime. The Butterworth function is analytic for all Ri > 0, and may be included in suitable boundary-layer parameterisation schemes where the turbulent diffusivity for heat is derived from the turbulent diffusivity for momentum.  相似文献   

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

5.
利用兰州大学半干旱气候与环境观测站(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,湍流平稳时间长度反映了中尺度运动的发生频率.  相似文献   

6.
We advance our prior energy- and flux-budget (EFB) turbulence closure model for stably stratified atmospheric flow and extend it to account for an additional vertical flux of momentum and additional productions of turbulent kinetic energy (TKE), turbulent potential energy (TPE) and turbulent flux of potential temperature due to large-scale internal gravity waves (IGW). For the stationary, homogeneous regime, the first version of the EFB model disregarding large-scale IGW yielded universal dependencies of the flux Richardson number, turbulent Prandtl number, energy ratios, and normalised vertical fluxes of momentum and heat on the gradient Richardson number, Ri. Due to the large-scale IGW, these dependencies lose their universality. The maximal value of the flux Richardson number (universal constant ≈0.2–0.25 in the no-IGW regime) becomes strongly variable. In the vertically homogeneous stratification, it increases with increasing wave energy and can even exceed 1. For heterogeneous stratification, when internal gravity waves propagate towards stronger stratification, the maximal flux Richardson number decreases with increasing wave energy, reaches zero and then becomes negative. In other words, the vertical flux of potential temperature becomes counter-gradient. Internal gravity waves also reduce the anisotropy of turbulence: in contrast to the mean wind shear, which generates only horizontal TKE, internal gravity waves generate both horizontal and vertical TKE. Internal gravity waves also increase the share of TPE in the turbulent total energy (TTE = TKE + TPE). A well-known effect of internal gravity waves is their direct contribution to the vertical transport of momentum. Depending on the direction (downward or upward), internal gravity waves either strengthen or weaken the total vertical flux of momentum. Predictions from the proposed model are consistent with available data from atmospheric and laboratory experiments, direct numerical simulations and large-eddy simulations.  相似文献   

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

9.
  相似文献   

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

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

12.
An analytical model is developed for the initial stage of surface wave generation at an air–water interface by a turbulent shear flow in either the air or in the water. The model treats the problem of wave growth departing from a flat interface and is relevant for small waves whose forcing is dominated by turbulent pressure fluctuations. The wave growth is predicted using the linearised and inviscid equations of motion, essentially following Phillips [Phillips, O.M., 1957. On the generation of waves by turbulent wind. J. Fluid Mech. 2, 417–445], but the pressure fluctuations that generate the waves are treated as unsteady and related to the turbulent velocity field using the rapid-distortion treatment of Durbin [Durbin, P.A., 1978. Rapid distortion theory of turbulent flows. PhD thesis, University of Cambridge]. This model, which assumes a constant mean shear rate Γ, can be viewed as the simplest representation of an oceanic or atmospheric boundary layer.For turbulent flows in the air and in the water producing pressure fluctuations of similar magnitude, the waves generated by turbulence in the water are found to be considerably steeper than those generated by turbulence in the air. For resonant waves, this is shown to be due to the shorter decorrelation time of turbulent pressure in the air (estimated as ∝ 1/Γ), because of the higher shear rate existing in the air flow, and due to the smaller length scale of the turbulence in the water. Non-resonant waves generated by turbulence in the water, although being somewhat gentler, are still steeper than resonant waves generated by turbulence in the air. Hence, it is suggested that turbulence in the water may have a more important role than previously thought in the initiation of the surface waves that are subsequently amplified by feedback instability mechanisms.  相似文献   

13.
Our study employed large-eddy simulation (LES) based on a one-equation subgrid-scale model to investigate the flow field and pollutant dispersion characteristics inside urban street canyons. Unstable thermal stratification was produced by heating the ground of the street canyon. Using the Boussinesq approximation, thermal buoyancy forces were taken into account in both the Navier–Stokes equations and the transport equation for subgrid-scale turbulent kinetic energy (TKE). The LESs were validated against experimental data obtained in wind-tunnel studies before the model was applied to study the detailed turbulence, temperature, and pollutant dispersion characteristics in the street canyon of aspect ratio 1. The effects of different Richardson numbers (Ri) were investigated. The ground heating significantly enhanced mean flow, turbulence, and pollutant flux inside the street canyon, but weakened the shear at the roof level. The mean flow was observed to be no longer isolated from the free stream and fresh air could be entrained into the street canyon at the roof-level leeward corner. Weighed against higher temperature, the ground heating facilitated pollutant removal from the street canyon.  相似文献   

14.
Turbulence in the Stable Boundary Layer at Higher Richardson Numbers   总被引：1，自引：1，他引：0  
We present some algebraic and numerical simulations of the stable boundary layer. We also discuss the problem of the existence of a critical Richardson number (Ri), beyond which the turbulence is suppressed. We compare the results of a second-order algebraic model with those of a third-order numerical model and, to this purpose, numerical simulations of a wind-tunnel flow, which is characterized by various Richardson numbers, were performed. As far as the second-order model is concerned, solutions, for the Richardson number greater than any critical value, can be obtained by modifying the time scales of the second-order equation pressure correlation terms in order to account for a buoyancy damping factor. We show that using a third-order model allows the same results (no critical Richardson number) to be obtained without modifications to the time scales. It is suggested that the non-locality, accounted for by the third-order moments, could allow the turbulence to persist also for Ri > 1.  相似文献   

15.
Turbulent mixing across heat-stratified density interfaces was studied in the laboratory using oscillating-grid generated turbulence. The aim was to study the transition between the entrainment regimes dominated by interfacial wave-breaking and molecular diffusion, and to study the characteristics of the latter. It was observed that, above a critical Richardson number Ric, which depends on the Peclet number Pe, the mixing due to wave breaking disappears and that Ric Pe−n, where the mean value of the exponent n is approximately . Above Ric, the entrainment is molecular-diffusion dominated and takes place through a sequence of events: the buoyancy gradient of the initially sharp density interface is weakened by molecular diffusion until the mixed-layer eddies can engulf a portion of the interfacial layer wherefore the interface sharpens again. Thus, the entrainment events are recurrent with a rate-controlling diffusion stage between them. An entrainment law of the form E Ri−2Pe−2, where E is the entrainment coefficient and Ri is the Richardson number, is suggested for the diffusion-dominated entrainment regime.  相似文献   

16.
Observations obtained over a glacier surface in a predominantlykatabatic flow and with a distinctwind maximum below 13-m height are presented. The data werecollected using a 13-m high profilemast and two sonic anemometers (at about 2.5-m and 10-m heights).The spectra at frequencies belowthat of the turbulence range appear to deviate considerably fromthe curves obtained by Kaimal andco-workers during the 1968 Kansas experiment. The characteristicsof these deviations are compared tothe observations of others in surface-layers disturbed by anykind of large-scale outer-layer (orinactive) turbulence. In our case the disturbances arelikely to be induced by the highmountain ridges that surround the glacier. Moreover, the deviationsobserved in the cospectra seemto result from an, as yet, unspecified interaction between theinactive outer-layer turbulenceand the local surface-layer turbulence. Near the distinctwind maximum turbulence production ceasedwhile turbulence itself did not, probably the result ofturbulence transport from other levels. Consequently, we studied thelocal similarity relations using w instead of u* as an alternative velocity scale. Wellbelow the wind maximum, and for relatively low stability(0< Rig <0.2), the flow behaves accordingto well established local-scaling similarity relationshipsin the stable boundary layer. For higherstability (Rig > 0.2), and near or above the wind maximum, the boundary-layer structure conforms tothat of z-less stratification suggesting that the eddy sizeis restricted by the local stability ofthe flow. In line with this we observed that the sensibleheat fluxes relate remarkably well to thelocal flow parameters.  相似文献   

17.
We have analyzed eddy covariance data collected within open canopy to investigate the influence of non-flat terrain and wind direction shear on the canopy turbulence. The study site is located on non-flat terrain with slopes in both south-north and east-west directions. The surface elevation change is smaller than the height of roughness element such as building and tree at this site. A variety of turbulent statistics were examined as a function of wind direction in near-neutral conditions. Heterogeneous surface characteristics results in significant differences in measured turbulent statistics. Upwind trees on the flat and up-sloping terrains yield typical features of canopy turbulence while upwind elevated surface with trees yields significant wind direction shear, reduced u and w skewness, and negligible correlation between u and w. The directional dependence of turbulence statistics is due that strong wind blows more horizontally rather than following terrain, and hence combination of slope related momentum flux and canopy eddy motion decreases the magnitude of Sk w and r uw for the downslope flow while it enhances them for the upslope flow. Significant v skewness to the west indicates intermittent downdraft of northerly wind, possibly due to lateral shear of wind in the presence of significant wind direction shear. The effects of wind direction shear on turbulent statistics were also examined. The results showed that correlation coefficient between lateral velocities and vertical velocity show significant dependence on wind direction shear through change of lateral wind shear. Quadrant analysis shows increased outward interaction and reduced role of sweep motion for longitudinal momentum flux for the downslope flow. Multi-resolution analysis indicates that uw correlation shows peak at larger averaging time for the upslope flow than for the downslope flow, indicating that large eddy plays an active role in momentum transfer for the upslope flow. On the other hand, downslope flow shows larger velocity variances than other flows despite similar wind speed. These results suggest that non-flatness of terrain significantly influences on canopy-atmosphere exchange.  相似文献   

18.
Application of the E-ε turbulence model to the atmospheric boundary layer   总被引：1，自引：0，他引：1  
  相似文献   

19.
We analyzed the structure and evolution of turbulent transfer and the wind profile in the atmospheric boundary layer in relation to aerosol concentrations during an episode of heavy haze pollution from 6 December 2016 to 9 January 2017. The turbulence data were recorded at Peking University’s atmospheric science and environment observation station. The results showed a negative correlation between the wind speed and the PM2.5 concentration. The turbulence kinetic energy was large and showed obvious diurnal variations during unpolluted (clean) weather, but was small during episodes of heavy haze pollution. Under both clean and heavy haze conditions, the relation between the non-dimensional wind components and the stability parameter z/L followed a 1/3 power law, but the normalized standard deviations of the wind speed were smaller during heavy pollution events than during clean periods under near-neutral conditions. Under unstable conditions, the normalized standard deviation of the potential temperature σ θ /|θ*| was related to z/L, roughly following a –1/3 power law, and the ratio during pollution days was greater than that during clean days. The three-dimensional turbulence energy spectra satisfied a –2/3 power exponent rate in the high-frequency band. In the low-frequency band, the wind velocity spectrum curve was related to the stability parameters under clear conditions, but was not related to atmospheric stratification under polluted conditions. In the dissipation stage of the heavy pollution episode, the horizontal wind speed first started to increase at high altitudes and then gradually decreased at lower altitudes. The strong upward motion during this stage was an important dynamic factor in the dissipation of the heavy haze.  相似文献   

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

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

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