DISSIPATION EFFECT ON THE FLOW OVER MOUNTAINS*

In the paper,the turbulent dissipation is considered in the model for studying the flow over mountains.The governing equation is a first order ordinary differential equation derived from an algebraic equation without dissipation case.The solution is sensitive to the upstream condition of dissipation.The dissipation not only reduces the strength of discontinuity but also changes the properties of the governing equation.In the paper,the qualitative characteristic features of the governing equation are discussed.The numerical results with super and sub-critical cases are also discussed in detail.The results show that the turbulent dissipation is an important factor and is not negligible.     

湍流的粘性和频散效应

本文应用Prandtl混合长理论导得了既包含湍流粘性又包含湍流频散的新的Reynolds平均运动方程组.分析指出:(1)湍流不仅存在粘性效应,而且存在频散效应;(2)正是湍流的频散效应可能导致能量逆转(即常说的负粘性现象,实质为一定条件下的负频散现象),并给出了能量逆转的必要条件和充分条件;(3)给出了湍流的KdV-Burgers方程模型.     

Turbulence structure in temperature inversions and in convection fields as observed by doppler sodar

The structure of turbulence in an inversion layer and in an homogeneous convective field of the planetary boundary layer is described. In the first part of the paper, we validate the sodar estimates of turbulent dissipation, by using measurements with an hot-wire anemometric system in situ. Limitations of an ε measurement technique using structure function calculations are given, taking account of atmospheric properties and acoustic Doppler instrumental effects. By comparison between isopleths of backscattering intensity and of turbulent dissipation rates, we observe that in the early morning, turbulence is advected by mechanical turbulence generated by wind shear. The same mechanism seems to be operating in the case of an inversion layer capping thermal instability, when the convective activity is not too greatly developed. A turbulent kinetic energy budget is examined using aircraft, sodar, and tower measurements. This indicates a constant turbulent dissipation profile through a deep convective layer.     

Note on the viscous heating term in the temperature variance equation

Experimental results in the atmospheric surface layer indicate that the viscous heating term in the equation for the mean-square value of the turbulent temperature fluctuation is small, but the correlation between temperature and turbulent energy dissipation fluctuations is significant.     

大气运动的动能在湍流边界层中的耗散

在大气能量学中,大气平均运动的动能通过湍流摩擦力转化为脉动运动的能量并最终转化为内能的“变性能量”是大气动能的“汇”,正确估计其大小对研究动能的收支,以至天气系统的生衰,发展均有一定意义。此项由于与湍流通量有关,因而难于从大尺度气象参数计算。以往的研究中,或是将动能平衡方程中其余各项算出,而将此项作为余项估计出,或用地转风由粗浅的边界层模式算出。     

Second-Order Modelling of Turbulence in Katabatic Flows

A complete one-dimensional second-order closure model is used to simulate katabatic flows observed on glaciers and ice caps. The model is tested with two different closure assumptions for the viscous dissipation, one based on a prognostic equation for and the other on a diagnostic buoyant length scale. Both formulations give quite similar results. Model simulations are compared to observations made over sloping ice surfaces during periods dominated by katabatic flow. In general, good agreement is found for both mean wind and temperature profiles as well as eddy correlation measurements. It is also found that the turbulent transport terms play an important role in katabatic flows as opposed to the classical stable boundary layer where these terms are usually ignored. Even the turbulent transport of temperature variance, which leads to the well-known countergradient term in unstable boundary layers, is relatively important for modelling the observed temperature profiles. The effect of these terms on the flux-profile relationships, using observed and simulated profiles, is also discussed.     

NUMERICAL SIMULATIONS OF TURBULENT FLOW IN AND ABOVE PLANT CANOPIES

The turbulent flow in and above plant canopies is of fundamental importance to the understanding oftransport processes of momentum,heat and mass between plant canopies and atmosphere,and to microme-teorology.The Reynolds stress equation model(RSM)has been applied to calculate the turbulence in cano-pies in this paper.The calculated mean wind velocity profiles,Reynolds stress,turbulent kinetic energy andviscous dissipation rate in a corn canopy and a spruce forest are compared with field observed data and withWilson's and Shaw's model.The velocity profiles and Rynolds stress calculated by both models are in goodagreement,and the length scale of turbulence appears to be similar.     

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

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

植被内部及其上方湍流场的数值模拟

植被内部及其上方的湍流流场对于了解植被与大气之间的动量、热量和质量交换过程极其重要。本文把高阶湍流封闭模型的Reynolds应力方程模型(RSM)应用于植被湍流的计算,得到了风速、湍流动能、Reynolds应力及能量耗散率的垂直分布,与现场观测数据比较,甚为满意。     

Analysis of atmospheric flow over a surface protrusion using the turbulence kinetic energy equation

Flow over surface obstructions can produce significantly large wind shears such that adverse flying conditions can occur for aeronautical systems (helicopters, STOL vehicles, etc.) Atmospheric flow fields resulting from a semi-elliptical surface obstruction in an otherwise horizontally homogeneous statistically stationary flow are modelled with the boundary-layer / Boussinesq-approximation of the governing equation of fluid mechanics. The turbulence kinetic energy equation is used to determine the dissipative effects of turbulent shear on the mean flow. Mean-flow results are compared with those given in a previous paper where the same problem was attacked using a Prandtl mixing-length hypothesis. The diffusion and convection of turbulence kinetic energy not accounted for in the Prandtl mixing-length concept cause departures of the mean wind profiles from those previously computed, primarily in the regions of strong pressure gradients. Iso-lines of turbulence kinetic energy and turbulence intensity are plotted in the plane of the flow. They highlight regions of high turbulence intensity in the stagnation zone and sharp gradients in intensity along the transition from adverse to favourable pressure gradient.     

THEORETICAL STUDY OF COLD-AIR DRAINAGE

The formation and characteristics of cold-air drainage are studied in this paper.The terrain is considered as an ideal two-dimensional slope which is infinitely extended in y-direction.In general situations,the regularity of cold-air drainage in temporal and spatial variations is obtained by solving numerically the atmospheric equations of motion,the continuity equation,the thermodynamic equation,the surface thermal balance equation and the soil heat conduction equation.We also test the changes of cold-air drainage with cloudiness,the coefficient of turbulent exchange,slope and its length,and the heat flux and thermal characteristics of soil.The calculated results are in agreement with observations.     

A modified turbulent energy model for diffusion from elevated and ground point sources in neutral boundary layers

A simple model for turbulent diffusion is proposed based on a turbulent energy model modified by Galperin and Hassid (1986) wherein local equilibrium in the concentration flux equation is assumed resulting in a turbulent diffusivity tensor. The early stage of plume development is simulated through adjustment of the integral length scale of diffusion. A simple semi-empirical formulation is also proposed for the dissipation length scale of concentration variance. Although it is substantially simpler than that of Sykes et al. (1984), the model is shown to compare well with the measurements of Fackrell and Robins (1982a, b).     

考虑动量和浮力通量耗散的烟气抬升解析模式

环流湍流对烟流抬升有重要的作用。本文考虑了环境湍流引起的烟流的动量和浮力通量耗散，导出了烟流抬升的控制方程，给出了烟流的轨迹方程和中性条件下热浮力烟流的终极抬升高度。在烟囱附近，轨迹方程接近2/3次律；当下风距离增大时，与2/3次律有明显偏离并逐渐变平。与外场试验资料的比较表明，本模式能较好地模拟浮力烟流的轨迹和终极高度。     

采用隐式立方样条计算平流过程的数值模式及理想实验

本文介绍了一种用隐式立方样条求解平流方程的数值方法,并从理论上对其无条件稳定性进行了证明,在此基础上建立了一个在地形坐标系下的两维原始方程模式,模式在行星边界层参数化中引入了湍流动能方程,在模式顶部引入了吸收层.数值实验表明:模式有较好的计算稳定性,对较高的模式水平分辨率和复杂地形均有较强的适应能力;对复杂地形和下垫面非均匀热源条件下中尺度系统的模拟能获得合理的结果,并具有较高的精确度.     

Turbulent dissipation of cold air lake in a basin

Summary Problems of turbulent dissipation of a cold air lake (CAL) and the inversion layer bordering CAL on the upper boundary are presented and studied with the compound model. In wintertime such cold air lakes can persist for days even if rather strong winds are blowing above them. The required conditions for CAL dissipation are removed processes of its formation or maintenance, as well as a sufficiently strong invasion of turbulence in the inversion layer from above down-wards. By this, the inversion layer at first becomes stronger and dissipation is stopped, until the increase of turbulent kinetic energy of the upper flow enables further dissipation. Such turbulent dissipation process is shown by the model for typical conditions and for different initial values of the relevent variables.With 6 Figures     

Multifractal Characteristics of Intermittent Turbulence in the Urban Canopy Layer

The multifractality of energy and thermal dissipation of fully developed intermittent turbulence is investigated in the urban canopy layer under unstable conditions by the singularity spectrum for the fractal dimensions of sets of singularities characterizing multifractals. In order to obtain high-order moment properties of smallscale turbulent dissipation in the inertial range, an ultrasonic anemometer with a high sampling frequency of 100 Hz was used. The authors found that the turbulent signal could be singular everywhere. Moreover, the singular exponents of energy and thermal dissipation rates are most frequently encountered at around 0.2, which is significantly smaller than the singular exponents for a wind tunnel at a moderate Reynolds number. The evidence indicates a higher intermittency of turbulence in the urban canopy layer at a high Reynolds number, which is demonstrated by the data with high temporal resolution. Furthermore, the temperature field is more intermittent than the velocity field. In addition, a large amount of samples could be used for verification of the results.     

Conditional Analysis of an Internal Boundary Layer

Data from a convective internal boundary layer (IBL) are analyzed by focusing on the instantaneousstructure of the top of the IBL instead of the time-average structure.A conditional averaging technique is developed todiscriminate between air from above the IBL and air from below the IBL , which alternately invade some instrument levels due tosubstantial variation of the top of the convective IBL.Sensitivity to the conditional sampling criteria is examined.Inside the IBL , buoyant and mechanicalproduction and dissipation dominate the turbulent kinetic energy budget.The horizontal advection and turbulent transport terms are smaller, but not negligible. The inferred pressure correlation term is negligible.Above the IBL , buoyant production and dissipation, although weak,dominate the turbulent kinetic energy budget. Shear generation andturbulent transport are smaller but significant.     

Dissipation length in stable layers

An attempt is made, for an anticyclonic stable situation in the atmospheric boundary layer during November 1978, to estimate the dissipation length. The turbulent dissipation and the turbulent kinetic energy, as measured on a 100 m high mast, exhibit large variability, due to the sporadic turbulent structure of the stable layer, and associated variability of the Richardson number. A comparison with different parameterizations of the turbulent dissipation rate allows a validation of these parameterizations.     

Numerical simulation of January 28 cold air outbreak during GALE Part I: The model and sensitivity tests of turbulence closures

A mesoscale planetary boundary layer (PBL) numerical model has been developed to study airflow over complex topography. Turbulence closures using the turbulent kinetic energy (TKE) and dissipation () equations are investigated in combination with the level 2.5 scheme of Mellor and Yamada (1982) to determine eddy diffusivities for momentum and heat. This modified E- closure is simpler than the level 3 one which requires more prognostic equations for moist turbulent transport.One-dimensional (1-D) model results show that the PBL mean flows under various stability conditions are not significantly sensitive to the modified Blackadar and Kolmogorov eddy mixing-length formulations used in this E- model, although the latter yields excessively large mixing lengths in the entrainment region of the upper PBL. Eddy mixing lengths in the Kolmogorov-type formulation can be better defined by introducing background dissipation. Using the same prognostic TKE equation, the 1-D model results are not significantly affected by different diagnostic formulations in the closures. The simulated results compare well with large-eddy simulations and those obtained using higher-order closure schemes including the level 3 one. The results are found to be insensitive to eddy Prandtl number, in contrast to the 2-D model results (see Part II).