SOME MEASUREMENTS OF TURBULENCE CHARACTERISTICS OVER COMPLEX TERRAIN

Results are presented from a windturbulence measurements campaign at the summit of ahill in complex terrain in Andros island (Greece)where a wind turbine park has been constructed. Meanturbulence parameters and power spectra from three 30 mmasts located at the summit, upwind and downwind ofthe hill are analysed to ascertain the differencesfrom respective parameters in flat terrain and toassess the influence of the irregular topography. Thevariances of the horizontal (vertical) wind components are found tobe smaller (larger) than the onesfrom flat terrain. Of the available correlations forthe power spectra, the Von Karman one was found togive better results, although the low frequency partsshow pronounced lags. The asymmetry of the upwindterrain for different wind directions and especiallyits slope is clearly felt both by variances and byspectra. This influence, as expected, is found toincrease with increasing stability. These results havean important significance for wind turbine design andwind energy utilisation.     

Turbulence in Complex Terrain

Theoretical results and data are presented in order to describethe spectral modifications induced by the effects of surfaceinhomogeneity on the features of turbulent boundary-layer flows.The hill case is used by way of example, and atmospheric andlaboratory data in neutral conditions are discussed. Shear andblocking effects with different characteristics are shown to bepresent both upwind and over the hill top.     

Landscape Roughness Parameters for Sherwood Forest – Experimental Results

The aim of this work is to present experimentally evaluated effective roughnesses (zoe) of a partly forested landscape. Although the ratio of boundary-layer height to obstacle size was only of the order of 50, there still seemed to exist a height range of 75–200 m where surface-layer similarity was approximately valid. Attempts were made to use conventional wind profile analysis to evaluate zoe, but the small height range and the large number of variables initially led to unacceptable uncertainties. Fixing the displacement height zd, rather than fitting it, reduced the data scatter to an acceptable level. The profile-derived roughness lengths zop obtained in this way were in good agreement with previous work, and with an alternative roughness length estimate zof for which flux-derived profile parameters u* and * were used. This implies that the profile-derived roughnesses were consistent with the measured surface-layer momentum flux. Comparison of both roughness estimates also yielded an improved estimate of the displacement height. Besides this, the authors tested a landscape roughness evaluation method which makes use of the gustiness parameter Tu = u/U in the surface layer. The results obtained by this method were in fair agreement with the profile-derived data. In previous work, the gustiness method was advocated because it could be used at relatively low levels, perhaps even within the roughness sub-layer. At the present measuring site, this was not the case as the gustiness method was only valid in an approximate way, and for a limited height range.     

Ship-board measurements and estimations of air-sea fluxes in the western tropical pacific during TOGA COARE

Direct air-sea flux measurements were made on R/V Kexue #1 at 4 ° S, 156 ° E during the Tropical Ocean Global Atmosphere (TOGA) Coupled Ocean-Atmospheric Response Experiment (COARE) Intensive Observation Period (IOP). An array of six accelerometers was used to measure the motion of the anchored ship, and a sonic anemometer and Lyman-α hygrometer were used to measure the turbulent wind vector and specific humidity. The contamination of the turbulent wind components by ship motion was largely removed by an improvement of a procedure due to Shao based on the acceleration signals. The scheme of the wind correction for ship motion is briefly outlined. Results are presented from data for the best wind direction relative to the ship to minimize flow distortion effects. Both the time series and the power spectra of the sonic-measured wind components show swell-induced ship motion contamination, which is largely removed by the accelerometer correction scheme. There was less contamination in the longitudinal wind component than in the vertical and transverse components. The spectral characteristics of the surface-layer turbulence properties are compared with those from previous land and ocean results. Momentum and latent heat fluxes were calculated by eddy correlation and compared to those estimated by the inertial dissipation method and the TOGA COARE bulk formula. The estimations of wind stress determined by eddy correlation are smaller than those from the TOGA COARE bulk formula, especially for higher wind speeds, while those from the bulk formula and inertial dissipation technique are generally in agreement. The estimations of latent heat flux from the three different methods are in reasonable agreement. The effect of the correction for ship motion on latent heat fluxes is not as large as on momentum fluxes.     

On momentum flux and velocity spectra over waves

Data from a research tower in Lake Ontario are used to study the validity of Monin--Obukhov scaling in the marine atmospheric boundary layer under various wave conditions. It is found that over pure wind seas, the velocity spectra and cospectra follow established universal scaling laws. However, in the presence of swells outrunning weak winds, velocity spectra and cospectra no longer satisfy universal spectral shapes. Here, Monin–Obukhov similarity theory, and the classical logarithmic boundary layers, are no longer valid. It is further shown that, in the presence of such swells, the momentum flux can be significantly modified in comparison to pure wind sea values. The implications of these findings for bulk flux estimations and on the inertial dissipation method for calculating fluxes are discussed.     

Eddy Covariance Measurements On Mountain Slopes: The Advantage Of Surface-Normal Sensor Orientation Over A Vertical Set-Up

The measurement of scalar fluxes employing the eddy covariance method is a widely used experimental approach,for which the flow distortion due to obstacles (e.g., sensor mounts and mast)is a well-known but not fully solved problem. In order to reduce flow distortion we installed a sonic anemometer in a surface-normal orientationrelative to the terrain slope, and a second instrument in a verticalposition at a horizontal distance of 1.54 m from the first instrumentWe found a significant reduction in the rotation angle necessary for the coordinaterotation procedure in the x-z plane whencomputing 30-minute flux averages with the surface-normal orientation. In 91% of all cases this rotation angleremained within the angle of incidence of ±10° recommended bythe manufacturer. In contrast, only 24% of the measurements taken with the vertically mounted anemometer were obtained at an angle of incidencewithin ±10°, and 3% were outside the ±30° range specified for an acceptable operation.A data quality test based on the variance of vertical windspeed normalized with friction velocity (_w/u_*) revealed problems for application under stable conditions due to large uncertainties in the determination of the Monin–Obukhov stability parameter z/L. An alternative test using the bulk drag coefficient C_D revealed other problems related to the dependence of C_D on z/z₀, the measuring height normalized by the roughness length, which do not appear to be constantin complex terrain. With both tests, a tendency for a slightly improved dataquality was found for the surface normal set-up, which, however, proved statistically insignificant.It is concluded that the surface-normal set-up of a sonic anemometer significantly reduces flow distortion by thesensor head. Although the surface-normal mounting position therefore appears to be the preferred one, with decreased flow distortion and a slightly improved data quality, no significant differences in turbulent quantities were found between the two set-uppositions. Hence, the consequences for short-term measurements of massand energy fluxes with a surface-normal set-up in complex terrain appearto be relevant only if single flux events are to be inspected, while for long-term measurements of integrated fluxes both the surface-normaland vertical installation of the sonic anemometer are adequate,indicating that eddy covariance measurements in complex terrain are lessdelicate than expected.     

Evaluation of Haney-Type Surface ThermalBoundary Conditions Using a CoupledAtmosphere and Ocean Model

A coupled atmosphere-ocean model developed at the Institute for Space Studies at NASA Goddard Space Flight Center (Russell et al., 1995) was used to verify the validity of Haney-type surface thermal boundary condition, which linearly connects net downward surface heat flux Q to air / sea temperature difference △T by a relaxation coefficient k. The model was initiated from the National Centers for Environmental Prediction (NCEP) atmospheric observations for 1 December 1977, and from the National Ocean Data Center (NODC) global climatological mean December temperature and salinity fields at 1° ×1° resolution. The time step is 7.5 minutes. We integrated the model for 450 days and obtained a complete model-generated global data set of daily mean downward net surface flux Q, surface air temperature TA,and sea surface temperature To. Then, we calculated the cross-correlation coefficients (CCC) between Q and △T. The ensemble mean CCC fields show (a) no correlation between Q and △T in the equatorial regions, and (b) evident correlation (CCC≥ 0.7) between Q and △T in the middle and high latitudes.Additionally, we did the variance analysis and found that when k= 120 W m-2K-1, the two standard deviations, σQ and σk△T, are quite close in the middle and high latitudes. These results agree quite well with a previous research (Chu et al., 1998) on analyzing the NCEP re-analyzed surface data, except that a smaller value of k (80 W m-2K-1) was found in the previous study.     

Surface Characteristics Observed over the Central Tropical Indian Ocean During Indoex IFP99

The present study is based on the observations carried out over the IndianOcean from the Indian research vessel ORV Sagar Kanya during the intensive field phase of the Indian Ocean Experiment in January–March 1999. The study area spanned from 15°N to 20°S in the central Indian Ocean. Near surface variations and surface fluxes along the cruise track are presented. A comparison of near surface characteristics over the Indian Ocean and tropical west Pacific has been made. It is observed that the average difference between the sea surface temperature and air temperature at 10 m height was 0.7 °C over the study area, nearly half of that observed over the tropical west Pacific. A comparison between observed and NCEP reanalysissurface data has been made. We find good agreement between ship measured andNCEP reanalysis surface pressure, specific humidity and wind fields.On the other hand, surface air temperature in the reanalysis tends to be lowcompared to observations. The components of the net surface heat flux comparebetter in the north Indian Ocean than in the southern Indian Ocean.     

Optimization of Upstream Profiles in Modelled Flow Over Complex Terrain

This paper describes an adjoint method for data assimilation intoupstream boundary conditions of numerical modelsusing optimal control theory. Mathematical formalisms are given along with the numerical implementation of the schemein a column model of the atmospheric boundary layer. The optimized mean and turbulence profiles are used as an upstream solutionin a model of turbulent flow in complex terrain. To contrast thiswith other methods, two solutions for flow over an isolatedhill are calculated, one with an optimized upstream solution andone with a simple surface-layer formulation for the upstream solution.These two solutions are compared to observations and analytical theory. The adjoint optimization method is shown to producesolutions of flow in complex terrain that are substantively differentat the two solutions, with the optimized solution giving more accurate results.     

Numerical Simulation of Flow over Two-Dimensional Hills Using a Second-Order Turbulence Closure Model

We study turbulent flow over two-dimensional hills. The Reynolds stresses are represented by a second-order closure model, where advection, diffusion, production and dissipation processes are all accounted for. We solve a full set of primitive non-hydrostatic dynamic equations for mean flow quantities using a finite-difference numerical method. The model predictions for the mean velocity and Reynolds stresses are compared with the measured data from a wind-tunnel experiment that simulates the atmospheric boundary layer. The agreement is good. The performance of the second-order closure model is also compared withthat of lower level turbulence models, including the eddy-viscositymodel and algebraic Reynolds stress models. It is concluded that thepresent closure is a considerable improvement over the other modelsin representing various physical effects in flow over hills. Thefeasibility of running a finite-difference numerical simulationincorporating a full second-order closure model on an IBM workstationis also demonstrated.     

Boundary-Layer Turbulence Over The Nebraska Sandhills

Data from National Centre for Atmospheric Research (NCAR) Queen Air boundary-layer flights over the Nebraska Sandhills are analyzed to investigate the effects of these low hills on boundary-layer turbulence. The Sandhills are an area of anisotropic rolling terrain with characteristic wavelengths of order 2km and rms height variations of order 25m. The biggest impact is found in early morning flight data where horizontal velocity perturbations appear at the same wavelengths as the terrain and variances (normalised by u ² , where u is the local friction velocity) are significantly enhanced relative to standard flat terrain values. By contrast the vertical velocity variance seems less affected and terrain effects are much less evident in data from the afternoon convective boundary layer.     

Turbulence Closures In Neutral Boundary Layer Over Complex Terrain

In this work, three turbulence closure models, Mellor andYamada level 2.5, E - l and E - implemented in a circulation model, are compared in neutral condition over complex terrain. They are firstly applied to a one-dimensional case on flat terrain and then to a schematic two-dimensional valley. The simulation results, in terms of wind field and turbulent kinetic energy, are tested against measurements from a wind-tunnel experiment. The empirical constants defining the characteristic length scales of the closures are modified based on turbulence parameters estimated in the experiment. The formulation of the diffusion coefficients is analysed to explain the differences among the various closures in the simulation results. Regarding the mean flow, both on flat and complex terrain, all the closures yield satisfactory results. Concerning the turbulent kinetic energy, the best results are obtained by E - l and E - closures.     

Daytime Turbulence Statistics above a Steep Forested Slope

Six levels of simultaneously sampled ultrasonic data are used to analyse the turbulence structure within a mixed forest of 13 m height on a steep slope (35°) in an alpine valley. The data set is compared to other studies carried out over forests in more ideal, flat terrain. The analysis is carried out for 30-min mean data, joint probability distributions, length scales and spectral characteristics.Thermally induced upslope winds and cold air drainage lead to a wind speed maximum within the trunk space. Slope winds are superimposed on valley winds and the valley-wind component becomes stronger with increasing height. Slope and valley winds are thus interacting on different spatial and time scales leading to a quite complex pattern in momentum transport that differs significantly from surface-layer characteristics. Directional shear causes lateral momentum transports that are in the same order or even larger than the longitudinal ones. In the canopy, however, a sharp attenuation of turbulence is observed. Skewed distributions of velocity components indicate that intermittent turbulent transport plays an important role in the energy distribution.Even though large-scale pressure fields lead to characteristic features in the turbulent structure that are superimposed on the canopy flow, it is found that many statistical properties typical of both mixing layers and canopy flow are observed in the data set.     

复杂山地近地层强风特性分析

利用在贵州省西南部复杂山地上获取的近地层梯度风观测资料和三维超声测风仪观测资料,从中筛选出具有该地气候特征的强风样本,利用数量统计和谱分析等方法,计算分析了由于复杂地形影响而导致的局地低层强风的平均和脉动特征.平均风场主要表现在受当地主要特征地形(西北-东南走向的深切峡谷)影响,其全年的主导风向和最大风速出现的方向几乎完全转为沿峡谷走向,即使在符合中性大气层结条件下,风的垂直廓线也完全不满足幂指数分布,风攻角远远大于规范推荐的值,并且不同风向的强风攻角因地形影响其差异可达20°;强风条件下脉动风场的主要特征是:不同风向强风的湍流强度有所不同,在纵、横和垂直方向的湍流强度比值与现行设计规范给出的三维湍流强度比值有明显差异,其中垂直方向的湍流强度显著偏大是突出特征;湍流积分尺度偏大,其中纵向值偏大20%-60%,横向值在某些风向上可偏大3倍以上,垂直方向则普遍较平坦地形偏大一个量级左右;在桥梁结构较为敏感的频域范围内,各风向的湍流谱密度值有显著差异,其中不同风向在纵向上最大差值可达8倍,横向和垂直向可相差6倍,但无论哪个方向的湍流谱密度值均比台风中心要小1-2个量级.     

Vertical Structure Of Turbulence In Offshore Flow During Rasex

The adjustment of the boundary layer immediately downstream froma coastline is examined based on two levels of eddy correlation data collected on a mast at the shore and six levels of eddy correlation data and profiles of mean variables collected from a mast 2 km offshore during the Risø Air-Sea Experiment. The characteristics of offshore flow are studied in terms of case studies and inter-variable relationships for the entire one-month data set. A turbulent kinetic energy budget is constructed for each case study.The buoyancy generation of turbulence is small compared to shear generation and dissipation. However, weakly stable and weakly unstable cases exhibit completely different vertical structure. With flow of warm air from land over cooler water, modest buoyancy destruction of turbulence and reduced shear generation of turbulence over the less rough sea surface cause the turbulence to rapidly weaken downstream from the coast. The reduction of downward mixing of momentum by the stratification leads to smaller roughness lengths compared to the unstable case. Shear generation at higher levels and advection of stronger turbulence from land often lead to an increase of stress and turbulence energy with height and downward transport of turbulence energy toward the surface.With flow of cool air over a warmer sea surface, a convective internal boundary layer develops downstream from the coast. An overlying relatively thick layer of downward buoyancy flux (virtual temperature flux) is sometimes maintained by shear generation in the accelerating offshore flow.     

新源县复杂地形下的风电场风场结构特征分析

利用新源县风电场周边气象站及测风塔观测资料、ERA5-Land再分析资料和数值模拟结果,分析新源县风电场的风场结构特征。结果表明：（1）风电场周边气象站的逐小时10m风速均呈现出早、晚偏低,中午偏高的变化规律。喀拉布拉镇站和公安农场站的风向主要以东风和南风为主,肖尔布拉克沟站的风向则以南风为主。测风塔4302的风速随着高度变化不明显。测风塔4301-4306均存在南风和北风,但各风向占比有一定差异。随高度升高,测风塔南风风向呈现出东南转西南的趋势。（2）ERA5-Land资料不能很好地再现研究区域风场变化。（3）数值模拟的风场变化具有一定山谷风特征。夜间的风向以东南风为主,白天则低海拔河谷地带风吹向山顶,北部谷风逐渐主导东北风向。20时,风电场区域主要以西北风为主。     

AN EVALUATION OF SEVERAL TURBULENCE SCHEMES FOR THE PREDICTION OF MEAN AND TURBULENT FIELDS IN COMPLEX TERRAIN

A prognostic three-dimensional mesoscale model has been developed andused in one- and two-dimensional modes to evaluate ten local turbulenceclosure schemes. The schemes ranged from first-order to the two-equationprognostic schemes. Predictions by the models were compared for aone-dimensional convective boundary layer using mixed layer scaling andmeasurements to interpret the results. Two-dimensional simulations were alsoperformed for a sea-breeze flow and for flow over a hill. The results showedthat for all of the models considered, minor differences were produced in themean meteorological fields and in the vertical scalar fluxes, but majordifferences were apparent in the velocity variances and dissipation rate.Predicted tracer concentrations were very sensitive to the turbulence modelformulation for dispersion from a point source in the convective boundarylayer, particularly for the prediction of maximum concentrations. Predictedtracer concentrations from a surface volume source for the two-dimensionalsimulations were similar for all models, although the degree of mixing in themorning growth period produced some differences. Generally, good results forthe mean meteorological fields can be obtained with first-order schemes, evenif they underpredict the magnitude of turbulence in the convective boundarylayer, and reasonable tracer concentrations can also be obtained with thesemodels provided near-source effects are not important. The two-equationprognostic models performed best for the prediction of turbulence in theconvective boundary layer.     

秋季珠峰复杂地形下地表能量通量卫星遥感研究

依据目前较为成熟的地表能量通量遥感参数化方案,利用2006年9～10月珠峰地区秋季近地面大气梯度实测资料,结合中分辨率成像光谱仪MODIS卫星资料,遥感估算了包括珠峰地区在内的青藏高原南部地区地表能量通量的空间分布情况.结果表明:地表净辐射峰值约为420.0 W·m-2,介于200.0～620.0 W·m-2之间;地热通量峰值约为110.0 W·m-2,介于50.0～180.0 W·m-2之间;感热通量峰值约为270.0 W·m-2,介于120.0～280.0 W·m-2之间;潜热通量峰值约为90.0 W·m-2,介于0.0～250.0 W·m-2之间.由此可见,在高原季风期后的秋季,该地区白天地面加热场主要通过感热传输加热大气,地气之间水汽传输耗热是地表净辐射平衡中的小量.     

复杂地形风场的精细数值模拟

风能是一种重要气候资源,随着我国风电规模的迅速增大,发展风能资源评估系统和风功率预测系统已成为一项重要的研究内容。国内外对复杂地形风场结构的数值模拟有大量研究,随着计算机能力增强,以往用于空气动力学精细流场计算的计算流体力学(Computational Fluid Dynamics,CFD)模式越来越多地在气象领域得到应用,人们开始研究用中尺度预报模式和CFD模式结合进行复杂地形风场的数值模拟。本文的耦合模式系统采用中尺度气象模式(WRF),通过嵌套网格到内层尺度(一般是几公里),然后通过耦合CFD模式Fluent软件获得高分辨率(水平30~100 m,垂直150 m高度以下10 m)的风速分布资料,得到精细化的风场信息。通过对鄱阳湖北部区域和云南杨梅山复杂地形的风场模拟,提供了风能评估和预报的一种可行的方法。     

Measured And Simulated Latent And Sensible Heat Fluxes At Two Marine Sites In The Baltic Sea

In this study, turbulent heat flux data from two sites within the Baltic Sea are compared with estimates from two models. The main focus is on the latent heat flux. The measuring sites are located on small islands close to the islands of Bornholm and Gotland. Both sites have a wide wind direction sector with undisturbed over-water fetch. Mean parameters and direct fluxes were measured on masts during May to December 1998.The two models used in this study are the regional-scale atmospheric model HIRLAM and the ocean model PROBE-Baltic. It is shown that both models overestimate the sensible and latent heat fluxes. The overestimation can, to a large extent, be explained by errors in the air-water temperature and humidity differences. From comparing observed and modelled data, the estimated 8-month mean errors in temperature and humidity are up to 1 °C and 1 g kg^-1, respectively. The mean errors in the sensible and latent heat fluxes for the same period are approximately 15 and 30 W m^-2, respectively.Bulk transfer coefficients used for calculating heat and humidity fluxes at the surface were shown to agree rather well with the measurements, at least for the unstable data. For stable stratification, the scatter in data is generally large, and it appears that the bulk formulation chosen overestimates turbulent heat fluxes.