Ekman动量近似下中间边界层模式中的风场结构

发展了一个准三维的、中等复杂的边界层动力学模式，该模式包含了EKman动量近似下的惯性加速度和Blackadar的非线性湍流粘性系数，它进一步改进了Tan和Wu(1993）提出的边界层理论模型。该模型在数值计算复杂性上与经典Ekman模式相类似，但由于包含了Ekman动量近似下的惯性项，使得该模式比传统Ekman模式更近于实际过程。中详细地比较了该模式与其他简化边界层模式在动力学上的差异，结果表明：在经典的Ekman模式中，由于忽略了流动的惯性项作用，导致在气旋性切变气流（反气旋性切变气流）中风速和边界层顶部的垂直速度的高估（低估），而在半地转边界层模式中，由于高估了流动惯性项的作用，结果与经典Ekman模式相反。同样，该模式可以应用于斜压边界层，对于Ekman动量下的斜压边界层风场同时具有经典斜压边界层和Ekman动量近似边界层的特征。     

Near Wall Flow over Urban-like Roughness

In this study, comprehensive measurements over a number of urban-type surfaces with the same area density of 25% have been performed in a wind tunnel. The experiments were conducted at a free stream velocity of 10 m s^-1 and the main instrumentation was 120 ° x-wire anemometry, but measurement accuracy was checked using laser Doppler anemometry.The results haveconfirmed the strong three-dimensionalityof the turbulent flow inthe roughness sublayer and the depths of the inertial sublayer (log-law region) and roughness sublayer for each surface have been determined. Spatial averaging has been used to remove the variability of the flow in the roughness sublayer due to individual obstacles and it is shown that the spatially averaged mean velocity in the inertial sublayer and roughness sublayer can,together, be described by a single log-law with a mean zero-plane displacement and roughness length for the surface, provided that the proper surface stress is known. The spatially averaged shear stresses in the inertial sublayer and roughness sublayer are compared with the surface stress deduced from form drag measurements on the roughness elements themselves.The dispersive stress arising from the spatial inhomogeneity in the mean flow profiles was deduced from the data and is shown to be negligible compared with the usual Reynolds stresses in the roughness sublayer. Comparisons have been made between a homogeneous (regular element array) surface and one consisting of random height elements of the same total volume. Although the upper limits of the inertial sublayer for both surfaces were almost identical at equivalent fetch, the roughness sublayer was much thicker for the random surface than for the uniform surface, the friction velocity and the roughness length were significantly larger and the `roughness efficiency' was greater. It is argued that the inertial sublayer may not exist at all in some of the more extreme rough urban areas. These results will provide fundamental information for modelling urban air quality and forecasting urban wind climates.     

Near-Wall Flow Development After A Step Change In Surface Roughness

An experimental study of the initial flow field downstream of a step change in surface roughness is presented. The roughness length of the downstream surface was approximately tenfold that of the upstream roughness and, unlike all previous studies, attention was concentrated on the roughness sublayer region beneath the inertial (log-law) region. The experiments were conducted at a boundary layer Reynolds number of about 6 × 10⁴ (based on layer thickness andfree-stream velocity) and around a longitudinal location where the (downstream) roughness length, z_o2, was about 1% of the boundary-layer thickness atthe roughness change point.The thickness of the roughness sublayer was found for the two roughness. It was observed that the vertical profiles of mean velocity and turbulence characteristics started to show similarity after about 160z₀₂ downstream of the roughness change. The presence of a shear stress overshoot is shown to depend strongly on the precise location (with respect to the roughness elements) at which the measurements are made and the thickness of the equilibrium layer is shown to be very sensitive to the way it is defined. It is demonstrated that the growing equilibrium layer has first to encompass the roughness sublayer before any thickness of inertial sublayer can be developed. It follows that, in somepractical cases, like flows across some urban environments, the latter(log-law) region may never exist at all.     

Laboratory Studies Of Wind Stress Over Surface Waves

Simultaneous laboratory observations of wind speed, wind stress, and surfacewind-wave spectra are made under a variety of wind forcing patterns using cleanwater as well as water containing an artificial surfactant. Under typical experimentalconditions, more than half of the total stress is supported by the wave-induced stressrather than by the surface viscous stress. When the surfactant reduces the shortwind-wave spectra, the wind stress also decreases by as much as 20–30% at agiven wind speed. When the wind forcing is modulated in time, the wind stresstends to be higher under decreasing wind than under increasing wind at a givenwind speed, mainly because the response of short wind-wave spectra to varyingwind forcing is delayed in time. These examples clearly demonstrate that therelationship between the wind speed and the wind stress can be significantlymodified if the surface wave field is not in equilibrium with the wind forcing.Next, we examine whether the wind stress is estimated accurately if the wave-inducedstress by all surface wave components is explicitly evaluated by linear superpositionand is added to the surface viscous stress. It is assumed that the surface viscous stressis uniquely related to the wind speed, and that the wind input rate is determined by thelocal, reduced turbulent stress rather than the total stress. Our wind stress estimatesincluding the wave contributions agree well with observed wind stress values, evenif the surface wave field is away from its equilibrium with the wind in the presenceof surface films and/or under time-transient wind forcing. These observations stronglysuggest that the wind stress is accurately evaluated as a sum of the wave-induced stressand the surface viscous stress. At very high winds, our stress estimates tend to be lowerthan the observations. We suspect that this is because of the enhancement of wind stressover very steep (or breaking) short wind-waves.     

Wind velocity and sand transport on a barchan dune

We present measurements of wind velocity and sand flux performed on the windward side of a large barchan dune in Jericoacoara, northeastern Brazil. From the measured profile, we calculate the air shear stress using an analytical approximation and treat the problem of flow separation by an heuristic model. We find that the results from this approach agree well with our field data. Moreover, using the calculated shear velocity, we predict the sand flux according to well-known equilibrium relations and with a phenomenological continuum saltation model that includes saturation transients and thus allows for nonequilibrium conditions. Based on the field data and theoretical predicted results, we indicate the principal differences between saturated and nonsaturated sand flux models. Finally, we show that the measured dune moves with invariant shape and predict its velocity from our data and calculations.     

泰山顶部近地层风的方差与谱分析

风在速度与方向上的不稳定性可引发结构振动并可导致损坏。山顶风较少受局地作用影响, 流场较稳定，作为典型风况对山顶风湍流脉动情形进行的观测和分析，对于研究建筑结构风载振动情形有着重要意义。在不同风况下对泰山气象站近地层风的风速和风向进行了同时、逐秒测记，所得样本经检验符合正态分布。计算表明，风向样本方差显著大于风速样本方差。绘制了各风速、风向样本自、互谱密度曲线图。谱图显示诸样本具有红噪声序列特征，不同样本自、互谱密度曲线具有相近形状，而以风速自谱曲线吻合程度最好，谱密度曲线在周期为4秒及2秒处有峰值。此后的研究应结合记录仪器的改进增加采样频率，对周期为2秒以下谱曲线作进一步分析。     

Numerical Simulation of Drag Partition over Rough Surfaces

We present a numerical simulation of drag partition over rough surfaces. A computational fluid dynamics model is applied with high resolution to simulatingturbulent flows over arrays of roughness elements positioned on asmooth surface. The skin drag on the surface and the pressure drag on the roughnesselements are computed. The simulated drag partition compares well with wind-tunnelmeasurements and theoretical estimates for similar rough surfaces. This confirms that the computational approach offers an alternative to wind-tunnel and field experiments in studying drag and drag partition. The model is then applied to studying drag partition over rough surfaces with various roughness configurations. It is shown that drag partition depends not only on the magnitude of the roughness frontal area but also on the sizes and arrangement of roughness elements, because (1) the pressure drag coefficient is sensitive to roughness-element dimensions and (2) the arragement of roughness elements lead to different interferences of turbulent wakes. The impact ofthe latter factor is not insignificant.     

Splash-saltation trajectories of soil particles under wind-driven rain

It is usually recognized that relatively large amounts of soil particles cannot be transported by raindrop splashes under windless rain. However, the splash-saltation process can cause net transportation in the prevailing wind direction since variations in splash-saltation trajectory due to the wind are expected in wind-driven rain. Therefore, determining the combined effect of rain and wind on the process should enable improvement of the estimation of erosion for any given prediction technique. This paper presents experimental data on the effects of slope aspect, slope gradient, and horizontal wind velocity on the splash-saltation trajectories of soil particles under wind-driven rain. In a wind tunnel facility equipped with a rainfall simulator, the rains driven by horizontal wind velocities of 6, 10, and 14 m s⁻¹ were allowed to impact three agricultural soils packed into 20×55 cm soil pans placed at both windward and leeward slopes of 7%, 15%, and 20%. Splash-saltation trajectories were measured by trapping the splashed particles at distances downwind on a 7-m uniform slope segment in the upslope and downslope directions, respectively, for windward and leeward slopes. Exponential decay curves were fitted for the mass distribution of splash-saltation sediment as a function of travel distance, and the average splash-saltation trajectory was derived from the average value of the fitted functions. The results demonstrated that the average trajectory of a raindrop-induced and wind-driven soil particle was substantially affected by the wind shear velocity, and it had the greatest correlation (r=0.96 for all data) with the shear velocity; however, neither slope aspect nor slope gradient significantly predicted the splash-saltation trajectory. More significantly, a statistical analysis conducted with nonlinear regression model of C₁(u_*²/g) showed that average trajectory of splash saltation was approximately three times greater than that of typical saltating sand grain.     

A parameterisation for the threshold shear velocity to initiate deflation of dry and wet sediment

A new parameterisation for the threshold shear velocity to initiate deflation of dry and wet particles is presented. It is based on the balance of moments acting on particles at the instant of particle motion. The model hence includes a term for the aerodynamic forces, including the drag force, the lift force and the aerodynamic-moment force, and a term for the interparticle forces. The effect of gravitation is incorporated in both terms. Rather than using an implicit function for the effect of the aerodynamic forces as reported earlier in literature, a constant aerodynamic coefficient was introduced. From consideration of the van der Waals force between two particles, it was further shown that the effect of the interparticle cohesion force between two dry particles on the deflation threshold should be inversely proportional to the particle diameter squared. The interparticle force was further extended to include wet bonding forces. The latter were considered as the sum of capillary forces and adhesive forces. A model that expresses the capillary force as a function of particle diameter squared and the inverse of capillary potential was deduced from consideration of the well-known model of Fisher and the Young–Laplace equation. The adhesive force was assumed to be equal to tensile strength, and a function which is proportional to particle diameter squared and the inverse of the potential due to adhesive forces was derived. By combining the capillary-force model and the adhesive force model, the interparticle force due to wet bonding was simplified and written as a function of particle diameter squared and the inverse of matric potential. The latter was loglinearly related to the gravimetric moisture content, a relationship that is valid in the low-moisture content range that is important in the light of deflation of sediment by wind. By introducing a correction to force the relationship to converge to zero moisture content at oven dryness, the matric potential–moisture content relationship contained only one unknown model parameter, viz. moisture content at −1.5 MPa. Working out the model led to a rather simple parameterisation containing only three coefficients. Two parameters were incorporated in the term that applies to dry sediment and were determined by using experimental data as reported by Iversen and White [Sedimentology 29 (1982) 111]. The third parameter for the wet-sediment part of the model was determined from wind-tunnel experiments on prewetted sand and sandy loam aggregates. The model was validated using data from wind-tunnel experiments on the same but dry sediment, and on data obtained from simulations with the model of Chepil [Soil Sci. Soc. Am. Proc. 20 (1956) 288]. The experiments showed that soil aggregates should be treated as individual particles with a density equal to their bulk density. Furthermore, it was shown that the surface had to dry to a moisture content of about 75% of the moisture content at −1.5 MPa before deflation became sustained. The threshold shear velocities simulated with our model were found to be in good agreement with own observations and with simulations using Chepil's model.