首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 83 毫秒
1.
2.
A Similarity Theory for Saltation and Application to Aeolian Mass Flux   总被引:7,自引:0,他引:7  
A similarity theory for saltation is presented. To derive the theory, a saltation model is developed and applied to simulating particle motion in turbulent flows. The numerical results are then fitted to simple universal expressions, or similarity functions. This approach allows the determination of saltation mass flux and other quantities. While the theory can be general, we focus on studying the saltation mass flux profile q(z). It is shown that q(z) is determined by friction velocity, aerodynamic roughness length and soil particle size distribution. There are two limiting situations, i.e., the saltation of large particles in weak turbulence and that of small particles in strong turbulence, for these two cases, q(z) is respectively exponential and Gaussian. Modified saltation has an intermediate saltation mass flux profile. For multi-sized particles, q(z) is a weighted superposition of many different profiles. The theory is compared with wind-tunnel observations, and uncertainties of the theory are discussed.  相似文献   

3.
Based upon comparisons between published experimental data and simulated results on the vertical sand flux distribution in the saltation layer, Shao’s similarity saltation model has been greatly improved by correcting the average vertical particle lift-off velocity and using a more suitable universal roughness length. By the improved model, the vertical sand flux profile over the bare, dry and loose uniform sandy surface, which is quite representative of real desert surfaces, can be reproduced very well. Meanwhile, the surface transport rate and the characteristic and average saltation heights have been simulated and analyzed in detail, disclosing their relationships with friction velocity, particle size and roughness length, and the possible underlying mechanisms. Besides, the average particle lift-off velocity and the average mean vertical aerodynamic action upon the ascending particle, which determine the saltation process, are explicitly expressed by parameters involved in the similarity model, and their relationships with friction velocity, particle size and roughness length are also described concisely. The corrected average particle lift-off velocity makes it possible to investigate the characteristic particle trajectory, whose initial velocity equals the average lift-off velocity, so as to estimate the average particle against surface impacting velocity and the average aerodynamic action upon the saltation process.  相似文献   

4.
Observations of drifting snow on small scales have shown that, in spite of nearly steady winds, the snow mass flux can strongly fluctuate in time and space. Most drifting snow models, however, are not able to describe drifting snow accurately over short time periods or on small spatial scales as they rely on mean flow fields and assume equilibrium saltation. In an attempt to gain understanding of the temporal and spatial variability of drifting snow on small scales, we propose to use a model combination of flow fields from large-eddy simulations (LES) and a Lagrangian stochastic model to calculate snow particle trajectories and so infer snow mass fluxes. Model results show that, if particle aerodynamic entrainment is driven by the shear stress retrieved from the LES, we can obtain a snow mass flux varying in space and time. The obtained fluctuating snow mass flux is qualitatively compared to field and wind-tunnel measurements. The comparison shows that the model results capture the intermittent behaviour of observed drifting snow mass flux yet differences between modelled turbulent structures and those likely to be found in the field complicate quantitative comparisons. Results of a model experiment show that the surface shear-stress distribution and its influence on aerodynamic entrainment appear to be key factors in explaining the intermittency of drifting snow.  相似文献   

5.
Heterogeneous Saltation: Theory, Observation and Comparison   总被引:7,自引:0,他引:7  
In the theory of saltation, under development since the 1940s, it is often assumed that saltation is homogeneous, i.e., saltating particles are uniform in size and follow identical trajectories. This assumption is a limitation to the development of saltation theory. In this paper, as in some other studies made since the 1980s, we are concerned with the saltation of multi-sized particles in turbulent flows, a process that we refer to as heterogeneous saltation. The theory deals with several questions, including the variation of particle size distribution with height, the entrainment rates of particles in different size ranges, and the associated profiles of saltation flux, particle momentum flux and particle concentration. It is hypothesised that saltation is dynamically similar and ‘universal’ similarity functions can be established. The similarity function for saltation flux is presented. Field observations are carried out at the southern fringe of the Takla Makan Desert in April 2002. Measurements of streamwise saltation flux of 32 particle size groups are made using a sand particle counter, together with measurements of wind speed and other atmospheric variables. These data are used to validate the saltation theory, by examining whether the observed saltation flux and particle size distribution can be reproduced. It is shown that the theory is promising in predicting these quantities.  相似文献   

6.
A large-eddy simulation (LES) with the dynamic Smagorinsky-Germano subgrid-scale (SGS) model is used to study the dispersion of solid particles in a turbulent boundary layer. Solid particles are tracked in a Lagrangian way. The instantaneous velocity of the surrounding fluid is considered to have a large-scale part (directly computed by the LES) and a small-scale part. The SGS velocity of the surrounding fluid is given by a three-dimensional Langevin model written in terms of SGS statistics at a mesh level. An appropriate Lagrangian correlation time scale is considered in order to include the influences of gravity and inertia of the solid particle. Inter-particle collisions and the influence of particles on the mean flow are also taken into account. The results of the LES are compared with the wind-tunnel experiments of Nalpanis et al. (1993 J Fluid Mech 251: 661–685) and of Tanière et al. (1997 Exp in Fluids 23:463–471) on sand particles in saltation and in modified saltation, respectively.  相似文献   

7.
There are many geometrical factors than can influence the aerodynamic parameters of urban surfaces and hence the vertical wind profiles found above. The knowledge of these parameters has applications in numerous fields, such as dispersion modelling, wind loading calculations, and estimating the wind energy resource at urban locations. Using quasi-empirical modelling, we estimate the dependence of the aerodynamic roughness length and zero-plane displacement for idealized urban surfaces, on the two most significant geometrical characteristics; surface area density and building height variability. A validation of the spatially-averaged, logarithmic wind profiles predicted by the model is carried out, via comparisons with available wind-tunnel and numerical data for arrays of square based blocks of uniform and heterogeneous heights. The model predicts two important properties of the aerodynamic parameters of surfaces of heterogeneous heights that have been suggested by experiments. Firstly, the zero-plane displacement of a heterogeneous array can exceed the surface mean building height significantly. Secondly, the characteristic peak in roughness length with respect to surface area density becomes much softer for heterogeneous arrays compared to uniform arrays, since a variation in building height can prevent a skimming flow regime from occurring. Overall the simple model performs well against available experimental data and may offer more accurate estimates of surface aerodynamic parameters for complex urban surfaces compared to models that do not include height variability.  相似文献   

8.
In many land-surface models using bulk transfer (one-source) approaches, the application of radiometric surface temperature observations in energy flux computations has given mixed results. This is due in part to the non-unique relationship between the so-called aerodynamic temperature, which relates to the efficiency of heat exchange between the land surface and overlying atmosphere, and a surface temperature measurement from a thermal-infrared radiometer, which largely corresponds to a weighted soil and canopy temperature as a function of radiometer viewing angle. A number of studies over the past several years using multi-source canopy models and/or experimental data have developed simplified methods to accommodate radiometric–aerodynamic temperature differences in one-source approaches. A recent investigation related the variability in the radiometric–aerodynamic relation to solar radiation using experimental data from a variety of landscapes, while another used a multi-source canopy model combined with measurements over a wide range in vegetation density to derive a relationship based on leaf area index. In this study, simulations by a detailed multi-source soil–plant–environment model, Cupid, which considers both radiative and turbulent exchanges across the soil–canopy–air interface, are used to explore the radiometric–aerodynamic temperature relations for a semi-arid shrubland ecosystem under a range of leaf area/canopy cover, soil moisture and meteorological conditions. The simulated radiometric-aerodynamic temperatures indicate that, while solar radiation and leaf area both strongly affect the magnitude of this temperature difference, the relationships are non-unique, having significant variability depending on local conditions. These simulations also show that soil–canopy temperature differences are highly correlated with variations in the radiometric–aerodynamic temperature differences, with the slope being primarily a function of leaf area. This result suggests that two-source schemes with reliable estimates of component soil and canopy temperatures and associated resistances may be better able to accommodate variability in the radiometric–aerodynamic relation for a wider range in vegetated canopy cover conditions than is possible with one-source schemes. However, comparisons of sensible heat flux estimates with Cupid using a simplified two-source model and a one-source model accommodating variability in the radiometric-aerodynamic relation based on vegetation density gave similar scatter. On the other hand, with experimental data from the shrubland site, the two-source model generally outperformed the one-source scheme. Clearly, vegetation density/leaf area has a major effect on the radiometric–aerodynamic temperature relation and must be considered in either one-source or two-source formulations. Hence these adjusted one-source models require similar inputs as in two-source approaches, but provide as output only bulk heat fluxes; this is not as useful for monitoring vegetation conditions.  相似文献   

9.
赵建华  张强  隆霄 《气象学报》2013,71(4):731-742
局地地形起沙对沙尘模式的预报精度具有重要影响。基于狭管与下滑效应以及气流过山的已有研究成果,研究了中性层结与局地平缓地形条件下风速在这两种效应下次网格参数化问题,然后基于平坦地形的起沙公式、引入局地地形的摩擦速度,得到了相应的次网格起沙公式,最后对狭管效应、下滑效应以及两效应的叠加进行了理想试验。试验表明:(1)对于狭管效应的次网格参数化问题,风从狭管的阔口进(出)、窄口出(进)时,风速、跃移通量与起沙量均增大(减小);(2)下滑效应致使风速、跃移通量与起沙量均增大,且山脚增幅大于山坡;(3)对于两种效应的叠加,当气流从狭管的阔口进、窄口出时,风速、跃移通量和起沙量同时增大,山脚风速大于山坡,大粒径沙粒跃移通量增幅大于小粒径沙粒。反之,窄口进、阔口出时,出入口宽度之比存在一临界值,大于此临界值,风速减小;反之增大。而且,除了在山脚的小粒径沙粒在一定风速条件下跃移通量会增大外,其余情况均减少。最后,提出了本参数化方法与数值模式耦合的分算法与整算法。  相似文献   

10.
In wind tunnel studies of aeolian transport, the number and position of pitot tubes are decided by the researcher, so that there are important variations in the computation ofU * between studies. Velocity measurements seldom are made very close to mobile surfaces because the tubes become blocked by drifting sand grains. This practice is fortuitous as demonstrated by recent selfregulatory models of saltation which indicate that fluid and grain-borne shear stress vary substantially within the lowest 0.01 m and application of the logarithmic law is therefore unsound. This study reports detailed velocity measurements which further suggest that no single logarithmic expression, based on fixed values of κ and τ, adequately represents the full wind profile which includes the inner saltation cloud above 0.01 m and the outer grain-free region of the boundary layer. A much improved fit over the logarithmic wind profile model is achieved with a square root relation, although there is no known physical basis for this specific form of power model. Relatively shallow boundary-layer development in wind tunnels forces the velocity gradient above the region of momentum extraction to attain exceptionally large values, uncommon in natural settings.  相似文献   

11.
The traditional method for computing the steady oceanic circulation has been by stepping an oceanic model forward in time until transients are damped by friction. An alternative method, which has the potential for being more economical is to minimize the sum of the squares of the residuals of the steady model equations. A variety of algorithms might be considered for computing the minimum; attention here is focused on preconditioned conjugate-gradient descent with the gradient computed using an adjoint model. The choice of varibles, i.e. the preconditioning transformation used in the optimization process, is found to be critical to the efficiency of the method. An appropriate preconditioning transformation can be suggested by a heuristic analysis similar to that commonly used to test the stability of numerical models. The method is demonstrated within the context of the barotropic vorticity equation.  相似文献   

12.
A three-dimensional Lagrangian stochastic (LS) model to evaluate pollutant dispersion in the atmospheric boundary layer has been developed. The model satisfies the well-mixed criterion of Thomson and allows for inhomogeneous, skew turbulence. Making use of the spherical reference frame, one of the possible solutions has been obtained. A skewed joint probability density function (PDF), which reproduces the given velocity moments (means, variances, skewness and covariances), has been built-up by a linear combination of eight Gaussian PDFs. In order to verify consistency with the well-mixed criterion, the long term results have been compared with the theoretical behaviour. A comparison between our model and Thomson's published algorithms was also carried out. By comparing wind-tunnel data and numerical predictions, a further validation of our LS model has been obtained. From an analysis of the numerical results, we can state that our model is able to evaluate dispersion in the case of complex flows where the application of previous models is unsuccessful.  相似文献   

13.
利用区域边界层模式对杭州市热岛的模拟研究   总被引:31,自引:11,他引:20  
利用一个三维非静力区域边界层数值模式,对杭州地区城市热岛现象进行了数值模拟。数值模拟结果表明,杭州地区存在明显的城市热岛现象。夏季城市热岛现象较强,春秋季次之,冬季较弱。利用资源卫星遥感资料反演所得的地表温度,以及与地面气象站观测资料和数值模拟结果相比较,数值模拟结果和实测结果吻合得较好。在此基础上,模拟并探讨了人为热源、风速、地面粗糙度等因素对城市热岛强度的影响。  相似文献   

14.
本文提出了一种用于双向套网格模式的变格距差分计算方案。该方案在不同格距的网格区采用不同精度的差分格式,它自然地连接粗细网格,避免了一般套网格方案在粗细网格相重合点上进行的重复计算。用解析法和数值试验证明了:它与其它一些变格距差分格式相比,对短波的穿透能力有明显改进,虚假的反射也较小。应用该方案建立了正压原始方程双向套网格模式,并采用空间分解和时间分解计算方法。这不仅使二维问题转化为二个一维问题,而且二维套网格也可简化为一维均匀网格和一维套网格两部分,从而使计算和程序简化。用理想场为初值所做的一系列数值试验表明,该模式中的波可以自由进出粗细网格区,计算稳定。最后,还用该模式做了台风路径预报试验,给出了一些试验结果。  相似文献   

15.
A new non-iterative bulk parametrization for surface-layer transfer coefficients for momentum and heat is presented. It is applicable for a wide range of aerodynamic and thermal roughness lengths, and includes the effect of the roughness sublayer. As a consequence, the non-iterative method is suitable for every surface type, especially for urban surfaces for which existing non-iterative parametrizations fail. The analytical approximation compares very well with an iterative approach. Our method can be easily implemented in atmospheric numerical models that already employ a non-iterative approach.  相似文献   

16.
A simple analytical model of the effect of a hill on potential evaporation (i.e., evaporation from well-watered vegetation) is presented. It is shown to reproduce the results from a two-dimensional non-linear numerical model. The analytical model is used to demonstrate the sensitivity of the changes in evaporation to surface resistance, aerodynamic resistance, temperature, slope and sun angle. Methods used by MORECS (Meteorological Office Regional Evaporation Calculation System) to estimate the potential evaporation at high elevations, and by Numerical Weather Prediction models to represent sub-grid scale hills, take account of the change in altitude but neglect the slope effects. A slope of 20 degrees, which is typical of upland terrain in the U.K., increases the potential evaporation by about 5%.  相似文献   

17.
 A simple moments model used in studying the large-scale thermally driven ocean circulation, in one hemisphere, is extended with a set of evolution equations for the basin-averaged salinity gradients. Natural formulations of the boundary conditions for the heat flux and the (virtual) salt flux are given, the latter based on the SST-evaporation feedback. Stommel’s box model result, a coexisting thermal and saline solution, is retrieved in the limit of no rotation. Including rotation in a salt-dominated setting, a steady circulation is found which bifurcates for higher Rayleigh numbers in a periodic solution which becomes chaotic through a cascade of subharmonic bifurcations. Periodic motion results from two different mechanisms. First, the stable stationary state bifurcates into a periodic solution where anomalously saline water is advected by the overturning circulation. Second, this periodic solution bifurcates into a state which is dominated, during the larger part of the cycle, by diffusion and inertia, characterized by a decreasing overturning rate, and, during the subsequent shorter part of the cycle, by rapid advection and restratification of the entire basin. The basin-averaged vertical density field is stably stratified in the steady and the periodic regimes and remains statically stable in the chaotic regime. Received: 16 July 1996 / Accepted: 29 December 1997  相似文献   

18.
风沙流中风速廓线的数值模拟与实验验证   总被引:2,自引:0,他引:2  
如何描述风沙流中被风沙运动改变了的风速廓线是风沙相互作用研究中的关键问题之一.该文中将跃移风沙流视为一种颗粒拟流体,将跃移颗粒对气流产生的阻力用颗粒流的阻力系数来表达,建立了描写两场相互作用的数学模型.颗粒流的阻力系数采用了前人在液态流化床研究中得出的阻力系数表达形式,通过引入一个修正系数,使其适用于风沙流(气-固两相流).将风沙边界层划分为跃移颗粒所产生的阻力不可忽略的内边界层和跃移颗粒阻力可以忽略但受内边界层影响的外边界层,分别建立了内边界层和外边界层的风速廓线表达式.应用所建立的数学模型,根据由风洞实验测定的跃移风沙流的浓度分布和速度分布资料,计算了跃移风沙流中的风速廓线,并与风洞实测结果进行了对比.结果表明,计算风速廓线与实测风速廓线吻合得比较好,在半对数图上均为上凸的曲线,有别于无风沙运动时的直线.跃移边界层外风速分布可较好地用对数函数来描述.对风沙流中风速廓线的进一步分析证实了风沙物理学奠基人Bagnold在其早期观测风沙流中的风速廓线时提出的"结点现象"(Bagnold结),该结点的高度随风速的增大而升高,随颗粒粒径的增大而降低.根据数值模拟和模拟实验,可以认为有风沙运动的动床剪切风速是综合反映风场与跃移层以及地表之间相互作用的物理量.  相似文献   

19.
20.
The stability properties of a source-sink driven flow in a two-layer, rotating fluid is considered. A steady flow in an annulus has been determined previously (Maeland, 1982) by Ekman-layer dynamics at the interface of the two layers; in this paper, the stability will be studied by linear, inviscid theory. For a given configuration, the two parameters which govern the stability of the basic flow are the internal Froude numbers of each layer. The growth rates for the instabilities are calculated by a simple numerical method (matrix method).  相似文献   

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

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