细网格非静力二维高阶闭合模式对山体流场特征的模拟分析

本文采用二阶矩湍流闭合方案，分别就非静力和准静力两种条件建立ＰＢＬ数值模式，并模拟了二维山脊流场的结构，计算了二维山体对气流平均场和湍流场的影响。计算结果表明，无论是平均场还是湍流场，采用准静力假设都会造成模拟结果的较大偏差，且对平均场的影响比对湍流场的影响更为显著。两种模式计算结果的差别还随地形坡度和背景风增大而增大。     

大气近地层湍流能谱特重的再分析

根据Kolmogorov局地均匀性湍流假设，大气近地层湍流能谱在惯性副区符合“－5／3律”。但是作者通过翻阅大量文献以及实测资料分析发现，大气近地层湍流能谱符合“－5／3律”的频率范围比较狭窄，而用幂函数＋指数函数形式“fS(f)∝f^-ae^-bf”拟合更恰当一些，这种能谱函数形式是连续时间混沌系统所特有的。利用HEIHE实验湍流观测资料分析结果表明，对于径向风分量，a=0.447,b=0.228；切向风分量：a=0.489,b=0.190；垂直风分量：a=0.551,b=0.124；温度：a=0.588,b=0.123。     

大气近地层湍流能谱特征的再分析

根据Kolmogorov局地均匀性湍流假设,大气近地层湍流能谱在惯性副区符合"-5/3律".但是作者通过翻阅大量文献以及实测资料分析发现,大气近地层湍流能谱符合"-5/3律"的频率范围比较狭窄,而用幂函数+指数函数形式"fS(f)∝F-ae-bf"拟合更恰当一些,这种能谱函数形式是连续时间混沌系统所特有的.利用HEIHE实验湍流观测资料分析结果表明,对于径向风分量,a＝0.447,b＝0.228;切向风分量:a＝0.489,b＝0.190;垂直风分量:a＝0.551,b＝0.124;温度:a＝0.588,b＝0.123.     

谐波分析方法在沿海风速数值预报订正的应用

利用福建漳浦沿海一座测风塔2009年12月—2010年11月的观测资料和中央气象台发布的MM5风速数值预报资料,统计分析了沿海近地层风速和风速数值预报误差的变化特征,统计结果表明:全年和各季节的风速变化具有明显的日变化特征,日最小风速出现在上午至午后,日最大风速则出现在旁晚至凌晨;同时,风速数值预报与测风塔实测风速的误差也存在较明显的日变化特征,这种变化特征与近地层湍流的发生发展有密切联系。在统计分析基础之上,提出一种谐波分析与人工神经网络(ANN)相结合的24 h短期风速数值预报的订正方法,订正试验结果表明:该方法可以每24 h更新一次预报结果,相对于风速数值预报的精度,两次独立样本检验的平均绝对误差分别减小25.6%、28.8%,且订正后使得风速数值预报的系统性偏差有明显下降。     

森林近地层大气湍流特性观测分析

利用长白山森林生态系统定位研究站观测资料,及2003年8月和9月涡旋相关资料,分析和比较了该地区近地层包括风速、风向、大气稳定度在内的平均场特征,以及湍流强度、无量纲化风脉动方差相似性和地表通量变化特征。结果表明:(1)8月和9月稳定度都基本集中在0附近;(2)风速2 m·s-1的环境中,湍流发展最为旺盛,随着风速的增大湍流强度先迅速减小,当风速增大到3 m·s-1后,湍流强度偏离0值变大了一些,再继续增大到一定风速大小以后,湍流强度基本不随风速变化;(3)无量纲化三维风脉动方差符合Monin-Obukhov相似理论的"1/3"定律,其最佳通用相似函数在稳定和不稳定条件下都可以拟合得到;(4)地表通量均表现出明显的日变化特征,8月以潜热为主,感热较小;9月仍以潜热为主,但是相比8月明显偏小,感热变化不大。     

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

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

六种近地层湍流动量输送系数计算方案对比分析

选取30多年来近地层湍流通量研究中具有代表性的六种参数化方案, 应用GAME/Tibet试验中那曲通量观测站的实测资料, 对比分析了各方案计算所得的湍流动量输送系数 (C_M) 之间的差异。结果表明:六种参数化方案计算得到的湍流动量输送系数之间存在较大差异。对于那曲观测站稀疏短草下垫面而言, 稳定条件下当理查孙数小于0.1时, 除Businger71方案存在显著低估以外, 其他各方案均能较好估算湍流动量输送系数; 不稳定条件下, Dyer74方案对湍流动量输送系数的估算效果最好, 其次为Wang02, Launiainen95和Louis82方案, Businger71方案误差较大。     

珠穆朗玛峰绒布河谷近地层大气湍流及能量输送特征分析

利用中国科学院2005年珠穆朗玛峰地区科学考察期间 (4月2日至6月7日) 收集的大气观测资料，分析了珠峰绒布河谷近地层水平风速、温度、湍流强度、湍流通量日变化及地表能量平衡特征。通过分析得出近地层三维风速方差与稳定度的关系基本满足1/3次方规律；珠峰绒布河谷近地层大气水平风速、温度、动量通量、感热通量和潜热通量均存在明显的日变化；地表获得的能量很大一部分以感热形式散失掉了，潜热所占比重很小。另外，还发现绒布河谷地区地表能量通量各分量并不满足能量平衡方程Rn=Hs+Le+G。通过对地面加热场的分析发现珠峰地表白天是强热源，晚上转变为弱冷源。     

雷达定量测量降水中的异方差性及加权判别函数处理

利用Z-I关系法,通过2006年两次积层混合云降水过程的回波和雨量资料分析A,b值的变化规律.针对Z-I关系的非线性和不稳定性,发现Z-I统计关系拟合中存在着异方差性,提出了对传统的最小二乘法和最优判别函数法进行加权处理的方法,并与其他方法的计算结果作了对比分析.研究结果表明:使用加权判别函数法,可以在符合Z-I关系实际变化规律的前提下抑制异方差性在计算过程中的不稳定,同时保持回波反演降水与实测降水之间的偏差最小化.     

白天谷地逆温层生消过程的数值模拟

本文用二维数值模式模拟了谷地白天边界层的发展过程.通过求解大气动力.热力学方程组.模拟出了谷地上空的温度场、流场以及在两种不同层结条件下的湍流动能q~2的分布特征.从而对白天谷地上稳定层结(或逆温层)的生消机制进行了初步探讨.模拟出的温度分布廓线等与有关实测资料基本一致.     

Improvement of the Mellor–Yamada–Nakanishi–Niino Planetary Boundary-Layer Scheme Based on Observational Data in China

The Mellor–Yamada–Nakanishi–Niino (MYNN) planetary boundary-layer (PBL) scheme is a second-order turbulence closure model that is an improved version of the Mellor–Yamada scheme based on large-eddy simulation data. It simulates PBL structure and evolution well, particularly over the ocean surface. However, when used with various underlying surfaces in China, the scheme overestimates the turbulent momentum flux and the sensible heat flux. Based on observations of surface fluxes in China, we attempt to improve the MYNN model by modifying the parameters and representation of the turbulence scale. Closure constants and empirical expressions in the diagnostic equation are chosen first, and an additional component of the turbulent heat flux is considered in the potential temperature prognostic equation to improve the surface heat-flux modelling. The modified MYNN scheme is incorporated into a three-dimensional mesoscale model and is evaluated using various underlying surface observations. Amelioration of the surface turbulent fluxes is confirmed at five observational sites in China over different land-use types.     

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.     

Impact of Turbulence Closures on Diurnal Temperature evolution for Clear Sky Situations Over Belgium

The aim of this study is to quantify the impact of turbulence closure on the simulation of surface air temperature at screen height (1.5 m) over Belgium. The mesoscale model MAR (Modèle Atmosphérique Régional), developed at the Université catholique de Louvain, is used to examine one-dimensional situations. A new second-order closure (level 2.5) is implemented containing prognostic equations for all three velocity variances, and diagnostic or prognostic formulations for the dissipation. This closure is compared with first and one-and-a-half order closures. Idealized nearly-neutral and convective cases underline the differences between first and second-order closures, and between diagnostic and prognostic equations for the dissipation. The one-and-a-half and second-order closures give satisfying results, but the first-order closure produces generally less appropriate vertical diffusion. Observed clear sky and weak horizontal advection situations have shown the sensitivity of 24 h temperature evolution to the choice of the turbulent closure.     

An assessment of mixing-length closure schemes for models of turbulent boundary layers over complex terrain

The effectiveness of closure assumptions implemented in turbulent boundary-layer models is rather uncertain over complex terrain. Different closure schemes for Reynolds shear stress based on the mixing-length concept are compared with data from wind tunnel experiments over complex terrain and the results are analysed on the basis of second-order moment equations. A good estimation of the vertical momentum flux velocity scale turns out to be given by the standard deviation of the vertical velocity while the turbulent kinetic energy scaling gives less satisfactory results in regions where turbulence anisotropy is large. Fairly good results are given by closure models implementing a shear-limited mixing-length already proposed for non-logarithmic wind profiles, while large errors characterize traditional mixing-length formulations.     

A Hierarchy of Energy- and Flux-Budget (EFB) Turbulence Closure Models for Stably-Stratified Geophysical Flows

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.     

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.     

E − ε − 〈θ 2〉 turbulence closure model for an atmospheric boundary layer including the urban canopy

A modified three-parameter model of turbulence for a thermally stratified atmospheric boundary layer (ABL) is presented. The model is based on tensor-invariant parameterizations for the pressure–strain and pressure–temperature correlations that are more complete than the parameterizations used in the Mellor–Yamada model of level 3.0. The turbulent momentum and heat fluxes are calculated with explicit algebraic models obtained with the aid of symbol algebra from the transport equations for momentum and heat fluxes in the approximation of weakly equilibrium turbulence. The turbulent transport of heat and momentum fluxes is assumed to be negligibly small in this approximation. The three-parameter $E - \varepsilon - {\left\langle {\theta ^{2} } \right\rangle }$ model of thermally stratified turbulence is employed to obtain closed-form algebraic expressions for the fluxes. A computational test of a 24-h ABL evolution is implemented for an idealized two-dimensional region. Comparison of the computed results with the available observational data and other numerical models shows that the proposed model is able to reproduce both the most important structural features of the turbulence in an urban canopy layer near the urbanized ABL surface and the effect of urban roughness on a global structure of the fields of wind and temperature over a city. The results of the computational test for the new model indicate that the motion of air in the urban canopy layer is strongly influenced by mechanical factors (buildings) and thermal stratification.     

An Analytical one-Dimensional Second-Order Closure Model of Turbulence Statistics and the Lagrangian Time Scale Within and Above Plant Canopies of Arbitrary Structure

An analytical one-dimensional second-order closure model is developed to describe the within canopy velocity variances, turbulent intensities, dissipation rates, Lagrangian time scale and Lagrangian far field diffusivities for vegetation canopies of arbitrary structure and density. The model incorporates and extends the model of momentum transfer developed by Massman (1997) and the model of within canopy velocity variances developed by Weil (unpublished) from the second-order closure model of Wilson and Shaw (1977). Model predictions of within and above canopy velocity variances, turbulent intensities, dissipation rates and the Lagrangian time scale are in reasonable agreement with previously measured or estimated values for these parameters. The present model suggests that the Lagrangian time scale and the far field diffusivity could be strongly dependent upon foliage structure and density through the foliage effects on the velocity variances. A simple formulation for the Lagrangian time scale at canopy height is derived from model results. Taken as a whole, the present model may provide a relatively simple way to incorporate turbulence parameters into models of soil/canopy/atmosphere mass transfer.     

On the determination of the closure parameters in higher-order closure models

In higher-order closure models at least the pressure redistribution and the dissipation of the turbulent kinetic energy and the temperature variance have to be parameterized. Due to this, the introduction of proportionality coefficients — the so-called closure parameters — is forced, which have to be determined before the model is used. We compare a group of models which use the return-to-isotropy hypothesis (Rotta, 1951) to describe the pressure redistribution and assume local isotropy for the smallest eddies in order to parameterize the dissipation. Special concern is given to the method of Mellor and Yamada (1982). Some of the closure parameters are re-derived on the basis of sensitivity studies requiring that both shear production and buoyancy behave in a realistic way if pressure redistribution or dissipation is changed by varying the closure parameters. This set of parameters is compared with those obtained by fitting to experimental data, by use of the Monin-Obukhov similarity theory and by considering ratios of variances, covariances and mean flow gradients, respectively. It is shown that the various sets of closure parameters are at least of the same order. The differences give some insight into the advantages and disadvantages of the various determination procedures. However, the general accordance of the different parameter sets supports the assumption of universal constants.