Turbulent Intensities and Velocity Spectra for Bare and Forested Gentle Hills: Flume Experiments

To investigate how velocity variances and spectra are modified by the simultaneous action of topography and canopy, two flume experiments were carried out on a train of gentle cosine hills differing in surface cover. The first experiment was conducted above a bare surface while the second experiment was conducted within and above a densely arrayed rod canopy. The velocity variances and spectra from these two experiments were compared in the middle, inner, and near-surface layers. In the middle layer, and for the canopy surface, longitudinal and vertical velocity variances () were in phase with the hill-induced spatial mean velocity perturbation (Δu) around the so-called background state (taken here as the longitudinal mean at a given height) as predicted by rapid distortion theory (RDT). However, for the bare surface case, and remained out of phase with Δu by about L/2, where L is the hill half-length. In the canopy layer, wake production was a significant source of turbulent energy for , and its action was to re-align velocity variances with Δu in those layers, a mechanism completely absent for the bare surface case. Such a lower ‘boundary condition’ resulted in longitudinal variations of to be nearly in phase with Δu above the canopy surface. In the inner and middle layers, the spectral distortions by the hill remained significant for the background state of the bare surface case but not for the canopy surface case. In particular, in the inner and middle layers of the bare surface case, the effective exponents derived from the locally measured power spectra diverged from their expected  − 5/3 value for inertial subrange scales. These departures spatially correlated with the hill surface. However, for the canopy surface case, the spectral exponents were near  − 5/3 above the canopy though the minor differences from  − 5/3 were also correlated with the hill surface. Inside the canopy, wake production and energy short-circuiting resulted in significant departures from  − 5/3. These departures from  − 5/3 also appeared correlated with the hill surface through the wake production contribution and its alignment with Δu. Moreover, scales commensurate with Von Karman street vorticies well described wake production scales inside the canopy, confirming the important role of the mean flow in producing wakes. The spectra inside the canopy on the lee side of the hill, where a negative mean flow delineated a recirculation zone, suggested that the wake production scales there were ‘broader’ when compared to their counterpart outside the recirculation zone. Inside the recirculation zone, there was significantly more energy at higher frequencies when compared to regions outside the recirculation zone.     

A Flume Experiment on the Adjustment of the Mean and Turbulent Statistics to a Transition from Short to Tall Sparse Canopies

Water-flume experiments are conducted to study the structure of turbulent flow within and above a sparse model canopy consisting of two rigid canopies of different heights. This difference in height specifies a two-dimensional step change from a rough to a rougher surface, as opposed to a smooth-to-rough transition. Despite the fact that the flow is in transition from a rough to a rougher surface, the thickness of the internal boundary layer scales as x ^4/5, consistent with smooth-to-rough boundary layer adjustment studies, where x is the downstream distance from the step change. However, the analogy with smooth-to-rough transitions no longer holds when the flow inside the canopy and near the canopy top is considered. Results show that the step change in surface roughness significantly increases turbulence intensities and shear stress. In particular, there is an adjustment of the mean horizontal velocity and shear stress as the flow passes over the rougher canopy, so that their vertical profiles adjust to give maximum values at the top of this canopy. We also observe that the magnitude and shape of the inflection in the mean horizontal velocity profile is significantly affected by the transition. The horizontal and vertical turbulence spectra compare well with Kolmogorov’s theory, although a small deviation at high frequencies is observed in the horizontal spectrum within the canopy. Here, for relatively low leaf area index, shear is found to be a more effective mechanism for momentum transfer through the canopy structure than vortex shedding.     

A Numerically Efficient Parametrization of Turbulent Wind-Turbine Flows for Different Thermal Stratifications

The wake characteristics of a wind turbine in a turbulent atmospheric boundary layer under different thermal stratifications are investigated by means of large-eddy simulation with the geophysical flow solver EULAG. The turbulent inflow is based on a method that imposes the spectral energy distribution of a neutral boundary-layer precursor simulation, the turbulence-preserving method. This method is extended herein to make it applicable for different thermal stratification regimes (convective, stable, neutral) by including suitable turbulence assumptions, which are deduced from velocity fields of a diurnal-cycle precursor simulation. The wind-turbine-wake characteristics derived from simulations that include the parametrization result in good agreement with diurnal-cycle-driven wind-turbine simulations. Furthermore, different levels of accuracy are tested in the parametrization assumptions, representing the thermal stratification. These range from three-dimensional matrices of the precursor-simulation wind field to individual values. The resulting wake characteristics are similar, even for the simplest parametrization set-up, making the diurnal-cycle precursor simulation non-essential for the wind-turbine simulations. Therefore, the proposed parametrization results in a computationally fast, simple, and efficient tool for analyzing the effects of different thermal stratifications on wind-turbine wakes by means of large-eddy simulation.     

A Method of Analysis for Turbulent Flows Using the Streamline Coordinate System

A procedure is described for analysing the transport equations for Reynolds stresses written in a streamline coordinate system, starting from the fields of first- and second-order moments of wind velocity measured in a terrain-following system over topography. In the analysis, the equations are split into two parts: the first contains the terms that can be calculated directly from measurements; the second contains third-order moments that are parameterized using suitable models. To evaluate the error associated with both parts, a Monte-Carlo technique that takes into account the experimental errors is proposed. An example of the application of this method for the Reynolds shear stress equation, using wind-tunnel data for non-separating flow over a two-dimensional valley, is reported. The comparison between the measured and modelled parts is fair near the surface, while at higher levels, the modelled part can be shown to miss a correct treatment of the third-order moments. In the frame of this analysis, the use of the correct derivative transformation has been found to be significant even for moderately sloping topography.     

Direct Numerical Simulation of Turbulent Katabatic Slope Flows with an Immersed-Boundary Method

We investigate a Cartesian-mesh immersed-boundary formulation within an incompressible flow solver to simulate laminar and turbulent katabatic slope flows. As a proof-of-concept study, we consider four different immersed-boundary reconstruction schemes for imposing a Neumann-type boundary condition on the buoyancy field. Prandtl’s laminar solution is used to demonstrate the second-order accuracy of the numerical solutions globally. Direct numerical simulation of a turbulent katabatic flow is then performed to investigate the applicability of the proposed schemes in the turbulent regime by analyzing both first- and second-order statistics of turbulence. First-order statistics show that turbulent katabatic flow simulations are noticeably sensitive to the specifics of the immersed-boundary formulation. We find that reconstruction schemes that work well in the laminar regime may not perform as well when applied to a turbulent regime. Our proposed immersed-boundary reconstruction scheme agrees closely with the terrain-fitted reference solutions in both flow regimes.     

Large-Eddy-Simulation Study of the Effects of Building-Height Variability on Turbulent Flows over an Actual Urban Area

Large-eddy simulation (LES) is used to investigate the effects of building-height variability on turbulent flows over an actual urban area, the city of Kyoto, which is reproduced using a 2-m resolution digital surface dataset. Comparison of the morphological characteristics of Kyoto with those of European, North American, and other Japanese cities indicates a similarity to European cities but with more variable building heights. The performance of the LES model is validated and found to be consistent with turbulence observations obtained from a meteorological tower and from Doppler lidar. We conducted the following two numerical experiments: a control experiment using Kyoto buildings, and a sensitivity experiment in which all the building heights are set to the average height over the computational region \(h_{all}\). The difference of Reynolds stress at height \(z=2.5h_{all}\) between the control and sensitivity experiments is found to increase with the increase in the plan-area index (\(\lambda _p\)) for \(\lambda _p > 0.32\). Thus, values of \(\lambda _p\approx 0.3\) can be regarded as a threshold for distinguishing the effects of building-height variability. The quadrant analysis reveals that sweeps contribute to the increase in the Reynolds stress in the control experiment at a height \(z= 2.5h_{all}\). The exuberance in the control experiment at height \(z=0.5h_{all}\) is found to decrease with increase in the building-height variability. Although the extreme momentum flux at height \(z=2.5h_{all}\) in the control experiment appears around buildings, it contributes little to the total Reynolds stress and is not associated with coherent motions.     

Large-eddy Simulation of Turbulent Flow Across a Forest Edge. Part I: Flow Statistics

The statistics of turbulent flow across a forest edge have been examined using large-eddy simulation, and results compared with field and wind-tunnel observations. The moorland-to-forest transition is characterized by flow deceleration in the streamwise direction, upward distortion of the mean flow, formation of a high pressure zone immediately in front of the edge, suppression of the standard deviations and covariance of velocity components, and enhancement of velocity skewnesses. For the selected forest density, it is observed that the maximum distortion angle is about 8 degrees from the horizontal. Instead of approaching a downwind equilibrium state in a monotonic manner, turbulence (standard deviations and covariances of velocity components) and mean streamwise velocity undershoot in the transition zone behind the edge. Evolution of flow statistics clearly reveals the growth of an internal boundary layer, and the establishment of an equilibrium layer downwind of the edge. It is evident that lower-order moments generally adjust more quickly over the new rough surface than do higher-order moments. We also show that the streamwise velocity standard deviation at canopy height starts its recovery over the rough surface sooner than does the vertical velocity standard deviation, but completes full adjustment later than the latter. Despite the limited domain size upstream of the edge, large-eddy simulation has successfully reproduced turbulent statistics in good agreement with field and wind-tunnel measurements.     

The Analysis of Turbulence Intensity and Reynolds Shear Stress in Wall-Bounded Turbulent Flows at High Reynolds Numbers

We investigated the turbulent intensities and Reynolds shear stress at high Reynolds number $({Re_\tau = 5 \times 10^{6}})$ in the atmosphere surface layer (ASL) through analyzing observations in near-neutral stratified conditions. The results show that with increasing Reynolds number the streamwise turbulent intensity increases linearly, and the peak of the Reynolds shear stress extends to a higher non-dimensional height, which means that the thickness of the logarithmic region increases. Furthermore, our results provide evidence for treating the ASL as a canonical turbulent boundary layer, the results of which can be extended and applied to higher Reynolds number wall turbulence in the ASL.     

On Heavy Particle Dispersion in Turbulent Shear Flows: 3-D Analysis of the Effects of Crossing Trajectories

In the approaches used to predict the dispersion of discrete particles moving in a turbulent flow, the effects of crossing trajectories due to gravity (or any other external force field) are generally accounted for by modifying the integral time scales according to the well-known analysis of Csanady (J Atmos Sci 20:201–208, 1963). Here, an alternative theoretical analysis of the time correlation of the fluid velocity fluctuations along a particle trajectory is presented and applied in a turbulent shear flow. The study is carried out in the frame of three-dimensional Langevin-type stochastic models, where the main unknowns are the drift tensor components rather than the conventional integral time scales of the fluid seen by the particles. Starting from a model for the space-time velocity covariance tensor of the turbulence under the assumption of homogeneous shear flow, the various components of the time correlation tensor of the fluid seen are expressed in the asymptotic case of large mean relative velocity (between the particles and the flow) compared to the particle velocity fluctuations. In order to provide comparison with the generally used expressions arising from isotropic turbulence assumption, we examine also the conventional integral time scales of the fluid seen in the directions parallel and perpendicular to the mean relative velocity. The most prominent deviations from isotropic turbulence are observed when the external force field is in the direction of the mean velocity gradient: in this case the loss of correlation in the mean flow direction is significantly lower than expected in a uniform flow, an observation that is in qualitative agreement with the few available data.     

人工增雨作业试验分析

该文利用多普勒雷达回波和天气形势等资料,对2006-06-28人工增雨作业进行统计分析,结果表明:46次增雨作业效果明显,在中、低纬度地区,低空有弱切变,高空温度≤-3℃时进行人工催化能有效增加降雨量。     

Evening and Morning Transition of Katabatic Flows

An experimental investigation of the evening and morning transition phases of katabatic slope flows has been conducted to identify the mechanisms for their development and destruction over an isolated slope. The momentum and energy equations of the flow have been used to describe these mechanisms for the particular topographic features of the studied slope, and to outline the differences from the dynamics of well-developed simple slope flows. In the lowest portion of the slope, frontal characteristics have been identified in early evening periods when the local pre-existing near-surface thermal structure does not impose a katabatic acceleration. The frontal shape is determined by the near-surface thermal stability and ambient wind. The flow initiation is distinctly different when it is linked to the local surface cooling, in which case it develops gradually and produces a slight local warming.The erosion of the katabatic layer at mid-slope precedes that at the foot and is closely linked to dilution of the local surface inversion. The flow erosion at the foot is often delayed, as the warming of air proceeds uniformly at all heights near the ground, so maintaining the inversion due to warming produced by mixing and advective processes linked to the upslope flow development. The latter initiates first at mid-slope and then at the foot, where for a non-negligible time period it flows over the persistent katabatic flow. The prerequisite for the development of this structure is the maintenance of a shallow inversion in the first 2–3 m above the ground surface.The morning dilution of the katabatic flow is apparently different from common experience over simple slopes and may be attributed to the steep upper portion of the slope in association with its easterly orientation, which results in strong non-uniformity of the solar heating along the slope.     

总辐射表夜间零点偏移试验与分析

总辐射表夜间零点偏移试验在夜间晴天条件下,采用室外平行对比法对不同型号的总辐射表的零点进行测试。通过数据分析,分别比较了黑白型和全黑型总辐射表、国产与进口总辐射表夜间零点偏移的大小,并分析了形成零点偏移的原因。试验结果表明,在夜间晴天条件下国产总辐射表的零点偏移小于10 W·m~(-2)通过总辐射表感应面向下水平安装和加盖试验也验证了零点偏移确实是由"冷天空"所致。     

对流层大气的负熵流和熵产生

将原始熵平衡方程应用于大气系统，详细计算了对流层大气平均状态下的负熵流和熵产生及其各项分量，得出了它们之间的一些数量关系，并应用熵理论初步解释了大气运动的某些规律，对于天气预报和气候变化的研究具有一定参考意义。     

Large-eddy Simulation of Turbulent Flow across a Forest Edge. Part II: Momentum and Turbulent Kinetic Energy Budgets

Momentum and turbulent kinetic energy (TKE) budgets across a forest edge have been investigated using large-eddy simulation (LES). Edge effects are observed in the rapid variation of a number of budget terms across this vegetation transition. The enhanced drag force at the forest edge is largely balanced by the pressure gradient force and by streamwise advection of upstream momentum, while vertical turbulent diffusion is relatively insignificant. For variance and TKE budgets, the most important processes at the forest edge are production due to the convergence (or divergence) of the mean flow, streamwise advection, pressure diffusion and enhanced dissipation by canopy drag. Turbulent diffusion, pressure redistribution and vertical shear production, which are characteristic processes in homogeneous canopy flow, are less important at the forest transition. We demonstrate that, in the equilibrated canopy flow, a substantial amount of TKE produced in the streamwise direction by the vertical shear of the mean flow is redistributed in the vertical direction by pressure fluctuations. This redistribution process occurs in the upper canopy layers. Part of the TKE in the vertical velocity component is transferred by turbulent and pressure diffusion to the lower canopy levels, where pressure redistribution takes place again and feeds TKE back to the streamwise direction. In this TKE cycle, the primary source terms are vertical shear production for streamwise velocity variance and pressure redistribution for vertical velocity variance. The evolution of these primary source terms downwind of the forest edge largely controls the adjustment rates of velocity variances.     

我国沙尘天气微气象学和湍流输送特征研究进展

基于风蚀起沙的物理机制,从微气象学角度回顾了国内沙尘天气的研究进展,介绍了沙尘天气过程近地面层微气象学特征,以及湍流通量演变规律和湍流作用对沙尘暴过程参数化的影响。同时,结合国外的相关研究,探讨了国内在沙尘天气研究方面存在的不足和今后的研究方向,如:开展沙尘暴过程数值模式湍流参数化方案的分析及与卫星和实验观测结果的比对研究;获取沙尘天气过程近地面层气象要素垂直廓线的精细分布,加深理解不同沙尘天气演变规律及特征;加强不同沙尘源区土壤湿度对起沙阈值影响的研究;开展不同粒径土壤的起沙率以及沙尘通量廓线演变规律的研究等。     

Effects of Fetch on Turbulent Flow and Pollutant Dispersion Within a Cubical Canopy

The effects of fetch on turbulent flow and pollutant dispersion within a canopy formed by regularly-spaced cubical objects is investigated using large-eddy simulation. Six tracer gases are simultaneously released from a ground-level continuous pollutant line source placed parallel to the spanwise axis at the first, second, third, fifth, seventh and tenth rows. Beyond the seventh row, the standard deviations of the fluctuations in the velocity components and the Reynolds shear stresses reach nearly equivalent states. Low-frequency turbulent flow is generated near the bottom surface around the first row and develops as the fetch increases. The turbulent flow eventually passes through the canopy at a near-constant interval. The mean concentration within the canopy reaches a near-constant value beyond the seventh row. In the first and second rows, narrow coherent structures frequently affect the pollutant escape from the top of the canopy. These structures increase in width as the fetch increases, and they mainly affect the removal of pollutants from the canopy.     

AREM模式对“05.7”四川大暴雨的敏感性试验

利用AREM中尺度数值模式,对2005年7月8日四川盆地大暴雨天气过程进行了数值模拟,比较分析了AREM模式提供的不同降水方案、地表通量方案、地表辐射参数化方案对此次降水过程的模拟,以及对该过程中西南低涡活动特征的模拟。试验结果     

Experiments on numerical modeling of intense convection

Results of numerical simulations using the WRF-ARW nonhydrostatic model are presented for eight episodes of intense convection over European Russia in the summer of 2007. The calculations were performed on four nested grids with horizontal grid meshes of 27, 9, 3, and 1 km. Convection was parametrized on the first two grids and explicitly resolved on the other two. It has been found that simulations on finer grids with explicit calculation of convective flows make it possible to reproduce heavy rainfalls and strong-wind zones in the areas of intense convection. A preliminary verification of the short-range predictions of convective systems shows that the maximum 12-h precipitation totals and the maximum winds at 10 m are close, in the order of magnitude, to the observed values. Prediction of convection centers is the weakest point. Difficulties in the model verification associated with the absence of data with high space-time resolution are discussed.     

CINRAD雷达产品“缺角”现象分析及处理

基于CINRAD/SA雷达图像产品出现"缺角"的现象,分析信号流程并通过BDAVC5软件观察天线运动轨迹,发现激磁电压过低导致旋转变压器有时辨识不出来正弦波信号,引起方位跳码,使得天线在未完成本仰角扫描的情况下就提前跳转到下一仰角进行扫描,从而造成了仰角剩余方位角度上的数据缺测。详细阐述了激磁电压过低而导致数据缺测的原因,指出此类故障多数发生在伺服系统,提出利用BDAVC5软件分析天线运行轨迹在伺服故障定位方面能起到很好的辅助作用。