用拉氏动力学理论研究PBL非均匀湍流扩散

将Gifford和Legg使用的拉氏动力学技巧从地面层延伸到行星边界层(PBL),讨论了平均风和湍流的垂直非均匀性以及大气稳定度对中尺度扩散特征的影响,得到了它们之间的解析表达式,理论结果与野外实测资料之间有较好的一致性.本文结果表明,经典的湍流统计理论不能直接用来研究中尺度范围的扩散,要合理地考虑动力学因子的影响.     

关于不同尺度大气运动中的雷诺交换

本文从多尺度分解概念出发,建立了多尺度Reynolds方程组。证明了平均运动的Reynolds交换项应为各级子平均运动非线性项的平均和。文中将子尺度运动处理为描述质点个别运动的Langevin形式以避免Euler方式导致的高阶闭合困难,简单地得到了K闭合表示,给出了动量、热量、质点交换系数的表达式并阐明了它们之间差异的形成原因。此外本文还给出了表现无序脉动量在导向力(如浮力、科氏力)作用下建立有序平均量梯度的负K值及其存在条件,解释了大气中大尺度运动的一些负粘性现象。     

REYNOLDS EXCHANGE IN THE ATMOSPHERIC MULTI-SCALE MOTIONS

In this paper the system of Reynolds equations of the multi-scaled atmospheric motions is set up based on the con-cept of decomposing the meteorological elements into multi-scale disturbances.It is proved to be true that the Reynoldsexchange term in the averaged motion is equal to the sum of averaged nonlinear terms in all sub-averaged motions.Inorder to avoid the higher order closure in Eulerian approaches,a new K-theory based on the multi-scaled Reynoldsequations is given in which the subscale motions are described by Langevin equation as the air particles are moving inthe Eulerian average background.From the new K-theory are derived the momentum,heat and mass exchangecoefficients as the functions of statistical variables such as variances and Lagrangian time scales of velocity,temperatureand other meteorological elements in disturbances.The new K-theory also expounds the causes for the differences be-tween the exchange coefficients of one element and another and gives the ambient conditions in which the buoyancyand/or Coriolis force Will build the chaotic disturbances into the orderly gradient of mean values of the correspondingelements.In consequence the K-theory can be used to explain some of negative viscosity phenomena in atmospheric mo-tions.     

合力随时间线性变化的热带气旋路径分析

热带气旋在移动过程中,环境气压场和其自身内力等都在发生变化。假设这些影响热带气旋的力之和随时间的变化在某一时段内是张性的,从而导出热带气旋的路径曲线方程。分析方程,得逞以热带气旋移动的三个基本规律。通过实例计算,可以看到方程曲线与实际路径几乎一致的个例验证结果,并有意义地证实一些人们以前的研究成果。     

LAGRANGIAN METHOD USED IN THE RESEARCH OF ATMOSPHERIC CIRCULATION

Based on the conservation of entropy and potential vorticity in adiabatic atmospheric motion without theconsideration of friction,calculation is made of the trajectory of a particle on an isentropic surface by use ofthe data of FGGE III-b.Results of several calculation schemes of the trajectory discussed show that thelocal data interpolation and Runge-Kutta time-integral scheme is the best.The calculated trajectory reflectsthe large-scale atmospheric motion only and the small-scale motion emerges as a deviation term of thecalculated trajectory.And then the outbreak and propagation of planetary wave are studied by means ofthe deformation of a material line,with the result showing that the material line can be tracked in the trop-osphere only in a few days,beyond which the interaction between the small-scale waves and large-scale motionleads to its dramatical twisting and deformation.Therefore,the Lagrangian method is assumed to be an effectivemeans of diagnostic research in the nonlinear interaction in atmospheric circulation,in addition to the generalstudy of the atmospheric circulation.     

A Lagrangian Decorrelation Time Scale in the Convective Boundary Layer

A new method for deriving the Lagrangian decorrelation time scales for inhomogeneous turbulence is described. The expression for the time scales here derived for the convective boundary layer is compared to those estimated by Hanna during the Phoenix experiment. Then the values of C₀, the Lagrangian velocity structure function constant, and of B_i, the Lagrangian velocity spectrum constant, were evaluated from the Eulerian velocity spectra and from the Lagrangian time scales derived, under unstable conditions, from Taylor's statistical diffusion theory. The numerical coefficient of the lateral and vertical Lagrangian spectra in the inertial subrange was found equal to 0.21, in good agreement with previous experimental estimates.     

ON THE PDF MODELS FOR ATMOSPHERIC DIFFUSION IN BOUNDARY LAYER

In this paper, taking its turbulent exchange coefficient as a function of the Lagrangian timescale and standard variance of the turbulence in atmosphere, the atmospheric dispersion PDFmodels are obtained on the basis of atmospheric diffusion K-theory. In the model the statistics ofwind speed are directly used as its parameters instead of classic dispersion parameters. The bi-Gaussian PDF is derived in convective boundary layer (CBL), from the statistics of verticalvelocity in both of the downdraft and updraft regions that are investigated theoretically in the otherpart of this paper. Giving the driven parameters of the CBL (including the convective velocity scalew* and the mixing depth h_i) and the time-averaged wind speed at release level, the PDF model isable to simulate the distribution of concentration released at any levels in the CBL. The PDF'ssimulations are fairly consistent with the measurements in CONDORS experiment or the resultsbrought out by some numerical simulations.     

AN ATMOSPHERIC MOTION EQUATION BUILT ON THE CONSERVATIVE RELATIONSHIP OF THE ANGULAR MOMENTUM EXCHANGE BETWEEN THE SOLID EARTH AND THE ATMOSPHERE ON SEASONAL-ANNUAL TIMESCALE*

Within the seasonal-annual timeseale,there exists an angular momentum conservative exchange relationship between the solid earth and the atmosphere,and their angular momentum exchange not only can cause variations in length-of-day(LOD) but also can express anomalies in atmospheric general circulation.Therefore,their angular momentum exchange mechanism should be introduced into the general circulation model.Considering the angular momentum anomalous exchange caused by the air-earth interface friction effect,a whole-layer atmospheric motion equation is derived in this paper including the earth spin anomalous friction force parameterized by using the change in the earth rotation rate.Through analysing the equation,it shows that the magnitude of the earth spin anomalous friction force is the same as that of Coriolis force on seasonal-annual timescale.     

Large-Eddy Simulation of air Pollution Dispersion in the Nocturnal Cloud-Topped Atmospheric Boundary Layer

Effects of stratocumulus clouds on the dispersion of contaminants are studied in the nocturnal atmospheric boundary layer. The study is based on a large-eddy simulation (LES) model with a bulk parametrization of clouds. Computations include Lagrangian calculations of atmospheric dispersion of a passive tracer released from point sources at various heights above the ground. The results obtained show that the vertical diffusion is non-Gaussian and depends on the location of a source in the boundary layer.     

UNIVERSALITY OF THE LAGRANGIAN VELOCITY STRUCTURE FUNCTION CONSTANT (C0) ACROSS DIFFERENT KINDS OF TURBULENCE

In this paper, we evaluate the Lagrangian velocity structure function constant, C0, in the inertial subrange by comparing experimental diffusion data and simulation results obtained with applicable Lagrangian stochastic models. We find in several different flows (grid turbulence, laboratory boundary-layer flow and the atmospheric surface layer under neutral stratification) the value for C0 is 3.0 ± 0.5. We also identify the reasons responsible for earlier studies having not reached the present result.     

AN ATMOSPHERIC MOTION EQUATION BUILT ON THE CONSERVATIVE RELATIONSHIP OF THE ANGULAR MOMENTUM EXCHANGE BETWEEN THE SOLID EARTH AND THE ATMOSPHERE ON SEASONAL-ANNUAL TIMESCALE

Within the seasonal-annual timeseale,there exists an angular momentum conservative exchange relationship be-tween the solid earth and the atmosphere,and their angular momentum exchange not only can cause variations inlength-of-day(LOD)but also can express anomalies in atmospheric general circulation.Therefore,their angular mo-mentum exchange mechanism should be introduced into the general circulation model.Considering the angular momentum anomalous exchange caused by the air-earth interface friction effect,awhole-layer atmospheric motion equation is derived in this paper including the earth spin anomalous friction forceparameterized by using the change in the earth rotation rate.Through analysing the equation,it shows that the magni-tude of the earth spin anomalous friction force is the same as that of Coriolis force on seasonal-annual timescale.     

珠江三角洲东部低层大气流动与中尺度扩散特性

采用拉格朗日随机粒子扩散模式和三维风场诊断模式及实际气象观测资料,模拟分析以大亚湾核电厂厂址为中心的珠江三角洲东部地区低层大气流动和气载污染物的中尺度扩散情况。结果表明:(1)该区域低层大气流动受盛行风背景、天气系统、海陆风局地环流以及当地地形扰动的共同作用,季节变化明显。(2)冬-春-秋季,核电厂扩散的可能影响区域总体分布在偏北内陆方向和西南方的珠江出海口及外海地区,夏季总体分布于北方内陆地区。(3)污染物扩散输送出200km×200 km模拟区域的平均时间约12~20 h;污染物在模拟区域内滞留24 h以上的情况,冬-春-秋季三季的出现频率为1/5左右,夏季出现较少,约占1/20的比例。     

利用测温声雷达回波图估算大气扩散的一种参数化方法

本文提出了一种基于湍流扩散统计理论的利用声雷达回波图估算大气扩散参数的新方法。这种新方法,仅需要从声雷达回波图上读出边界层结构的信息并取得常规地面测风资料,即可利用参数化公式进行计算,不需进行Pasquill稳定度分类。 所应用的参数化公式用外场观测资料进行了检验,表明公式计算结果与观测值比较符合。这种方法计算得到的扩散参数也与P-G扩散曲线进行了比较。     

A similarity model for maximum ground-level concentration in a height-invariant,stably stratified atmospheric boundary layer

A model is presented for determining the location and magnitude of the maximum ground-level concentration arising from an elevated buoyant source in a very stable atmospheric boundary layer. The development combines the turbulent structure of such a boundary layer, Lagrangian similarity of the diffusion process, and similarity solutions of the conservation equations of the buoyant plume with mass conservation to produce a simple, experimentally verifiable formulation. Functional analogy with previous results for the constant flux layer and a deep convectively unstable layer suggest a heuristic model by which to visualize the process.     

A Velocity–Dissipation Lagrangian Stochastic Model for Turbulent Dispersion in Atmospheric Boundary-Layer and Canopy Flows

An extended Lagrangian stochastic dispersion model that includes time variations of the turbulent kinetic energy dissipation rate is proposed. The instantaneous dissipation rate is described by a log-normal distribution to account for rare and intense bursts of dissipation occurring over short durations. This behaviour of the instantaneous dissipation rate is consistent with field measurements inside a pine forest and with published dissipation rate measurements in the atmospheric surface layer. The extended model is also shown to satisfy the well-mixed condition even for the highly inhomogeneous case of canopy flow. Application of this model to atmospheric boundary-layer and canopy flows reveals two types of motion that cannot be predicted by conventional dispersion models: a strong sweeping motion of particles towards the ground, and strong intermittent ejections of particles from the surface or canopy layer, which allows these particles to escape low-velocity regions to a high-velocity zone in the free air above. This ejective phenomenon increases the probability of marked fluid particles to reach far regions, creating a heavy tail in the mean concentration far from the scalar source.     

An advection scheme based on the combination of particle mesh method and pure Lagrangian approach

Possibility of using pure Lagrangian approach in modeling transport phenomena is described in this paper. The application of pure Lagrangian approach in real atmospheric field induces highly irregular spatial distribution of grid points, after only a few time steps. In order to avoid problems caused by that irregularity, a quasi interpolation procedure is proposed. Proposed interpolation procedure is similar to the radial basis functions interpolation and does not impose any demands about spatial distribution of the grid points or about continuity and differentiability of the field that needs to be interpolated. Besides that, proposed procedure is explicitly mass conserving. Combination of particle mesh method and pure Lagrangian approach creates efficient transport scheme that does not produce any new local maxima and minima in advected field. In proposed advection scheme motion of points are performed in Lagrangian manner while spatial derivatives are evaluated on the basis of values interpolated onto regular grid. Applicability of proposed advection scheme in an unambiguous way is proved by performing “standard” numerical tests with (i) the slotted cylinder under solid body rotation, (ii) the test with Doswell’s idealized cyclogenesis as well as (iii) integration of shallow water equations.     

The impact of different vertical diffusion schemes in a three-dimensional oil spill model in the Bohai Sea

Vertical transport is critical to the movement of oil spills in seawater. Breaking waves play an important role by developing a well-defined mixing layer in the upper part of the water column. A three-dimensional （3-D） Lagrangian random walk oil spill model was used here to study the influence of sea surface waves on the vertical turbulence movement of oil particles. Three vertical diffusion schemes were utilized in the model to compare their impact on oil dispersion and transportation. The first scheme calculated the vertical eddy viscosity semi-empirically. In the second scheme, the vertical diffusion coefficient was obtained directly from an Eulerian hydrodynamic model （Princeton Ocean Model, POM2k） while considering wave- caused turbulence. The third scheme was formulated by solving the Langevin equation. The trajectories, percentages of oil particles intruding into water, and the vertical distribution structures of oil particles were analyzed for a series of numerical experiments with different wind magnitudes. The results showed that the different vertical diffusion schemes could generate different horizontal trajectories and spatial distributions of oil spills on the sea surface. The vertical diffusion schemes caused different water-intruding and resurfacing oil particle behaviors, leading to different horizontal transport of oil particles at the surface and subsurface of the ocean. The vertical diffusion schemes were also applied to a realistic oil spill simulation, and these results were compared to satellite observations. All three schemes yielded acceptable results, and those of the third scheme most closely simulated the observed data.     

A Wind Tunnel and Numerical Investigation of Turbulent Dispersion in Coastal Atmospheric Boundary Layers

Thermal internal boundary layers in onshore air flows have a significant influence on pollutant diffusion in coastal areas. Although several models for this diffusion problem exist, measurements for model verification are scarce. In this paper, we present a set of wind tunnel observations and examine the performance of a Lagrangian stochastic model. The good agreement between the model simulation and the tunnel measurements confirms the usefulness of the Lagrangian stochastic model for practical purposes. Sensitivity tests of the model to turbulence statistics show that uncertainty in velocity skewness to the extent of observational scatter does not seem to have a significant influence on pollutant dispersion, while uncertainties in turbulence intensity (variance) significantly influence the dispersion pattern.     

Comments on the 'Universality of the Lagrangian Velocity Structure Function Constant (C0) ACROSS DIFFERENT KINDS OF TURBULENCE’

Recently Du ( Boundary-Layer Meteorology 83, 207–219, 1997) estimated the value of the Lagrangian velocity structure constant, C0, in the inertial subrange by comparing experimental diffusion data and simulation results obtained with the one-dimensional form of Thomson's model ( J. Fluid Mech. 180, 529–556, 1987). Du reported that for several different flows (grid turbulence, a wind-tunnel boundary layer and the atmospheric surface layer under neutral stratification) the value of C0 is 3.0±0.5. Here, it is shown that optimal model agreement with experimental diffusion data for the wind-tunnel boundary layer is, in fact, obtained when C0=5.0 ± 0.5. It is also shown that accounting for the skewness of velocity statistics and finite Reynolds number effects does not significantly change this estimate for the value of C0. It is suggested that one-dimensional Lagrangian stochastic models are inconsistent with the supposed universality of C0.