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次律有系统偏离;这一偏离与烟流半径和高度的非线性关系相一致。上述事实与国外的一些观测结果均表明,环境湍流对烟气抬升的累积作用不容忽视。为此,提出了环境湍流与自生湍流的联合作用模式,由此导得的抬升路径方程可以满意地解释上述观测事实。模式还显示,卷挟速度除取决于切变速度及湍流强度以外,还是烟流半径的函数,因而在物理上更加合理。由联合作用模式导出的终极抬升公式比不考虑环境湍流累积作用的公式更符合观测结果。     

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.     

Scaling of the Vertical Spreading of a Plume of a Passive Tracer in a Rough-Wall Neutral Boundary Layer

The influence of surface roughness on the dispersion of a passive scalar in a rough wall turbulent boundary layer has been studied using wind-tunnel experiments. The surface roughness was varied using different sizes of roughness elements, and different spacings between the elements. Vertical profiles of average concentration were measured at different distances downwind of the source, and the vertical spread of the plume was computed by fitting a double Gaussian profile to the data. An estimate of the integral length scale is derived from the turbulence characteristics of the boundary layer and is then used to scale the measured values of plume spread. This scaling reduces the variability in the data, confirming the validity of the model for the Lagrangian integral time scale, but does not remove it entirely. The scaled plume spreading shows significant differences from predictions of theoretical models both in the near and in the far field. In the region immediately downwind of the source this is due to the influence of the wake of the injector for which we have developed a simple model. In the far field we explain that the differences are mainly due to the absence of large-scale motions. Finally, further downwind of the source the scaled values of plume spread fall into two distinct groups. It is suggested that the difference between the two groups may be related to the lack of dynamical similarity between the boundary-layer flows for varying surface roughness or to biased estimates of the plume spread.     

Macro- and micro-scale mixing in chemical reactive plumes

The average dispersion of a plume in the atmospheric boundary layer is strongly influenced by atmospheric turbulence. Atmospheric turbulence determines also concentration fluctuations due to turbulent meandering by large scale turbulent eddies and in-plume fluctuations, due to smaller scale eddies. Conversion of NO to NO₂ in a plume is influenced by micro-scale mixing, due to the concentration fluctuation correlation % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaa0aaaeaaca% qGobGaae4tamaaCaaaleqabaGaaeymaaaakiaab+eadaqhaaWcbaGa% ae4maaqaaiaabgdaaaaaaaaa!3AF4!\[\overline {{\rm{NO}}^{\rm{1}} {\rm{O}}_{\rm{3}}^{\rm{1}} } \] and macro-scale mixing, the mixing in of ambient air containing O₃ into the plume.The study of turbulent meandering, in-plume fluctuations, microscale and macro-scale mixing will contribute to a better understanding of concentration fluctuations in general.     

Buoyant plume rise in a Lagrangian framework

For the dispersion of buoyant material, the interaction with the environment by entrainment forms a serious obstacle for a formulation in a Lagrangian framework. Nevertheless an outline is given here on how buoyant plume rise in a Lagrangian sense could be described. Though the method contains a number of heuristic elements, it has all the advantages of a Lagrangian formulation. It is shown that it is possible to formulate a Lagrangian model which both is able to recover the classical formulations for plume rise in a calm environment and to accomodate more recent Eulerian formulations in a turbulent environment. Moreover, the method offers excellent possibilities to include the turbulent characteristics of the plume's environment and arbitrary stratifications of the boundary layer. These facts make it attractive for various practical applications. Some examples are given which illustrate this.     

Simulation of a Summer Urban Breeze Over Paris

Numerical simulations for an anticyclonic summer episode in the Paris area have been performed at the meso- scale for a 48-hour period, and compared to observations from a dense operational observational network. The meteorological stations have been classified, according to the extent of urbanization of their surroundings, into four classes (central Paris, urban, suburban, and rural). The atmospheric model, coupled with an urban surface scheme, correctly reproduces the temperature (within 1 K from the observations) and humidity. The intense urban heat island during the night is also well represented.Following the validation, the model is used to quantify atmospheric effects of Paris on the boundary layer, through a comparison with a purely rural simulation. At night, the model simulates a neutral or even slightly unstable boundary layer to a depth of 200 m over the city. In contrast, a very stable layer formed in the countryside. During the day, the boundary layer was more turbulent and 500 m deeper over Paris; vertical velocities of up to 1 m s^-1 were created over the city. This leads to an urban breeze with convergence at low levels (with winds around 5 to 7 m s^-1), and divergence at the boundary-layer top (with similar wind speeds). The horizontal extent of the breeze reaches for more than 50 km from the city centre, and could have an important impact on pollutant diffusion in the area for calm days.Finally, three other spring cases are presented briefly. These show that an urban breeze develops if the synoptic wind is weak enough or disorganized; an urban plume develops otherwise.     

Concentration fluctuation measurements in a plume dispersing in a stable surface layer

A set of tracer experiments studying concentration fluctuations in a pollutant plume dispersing near the surface in a stably stratified nocturnal boundary layer is described, and the results are compared with those obtained in near-neutral stability conditions by Mylne and Mason (1991). The results highlight the importance of slow meandering of the plume which is characteristic of stable conditions. This meandering makes it impossible to conduct experiments under near-stationary conditions, resulting in considerable statistical variability in the results, but is important in reducing time-averaged concentrations. Spectral characteristics of the plume and general fluctuation statistics are qualitatively similar to those in near-neutral stability, but there are significant quantitative differences. Fluctuation time scales are shown to be substantially longer under stable conditions. This difference cannot be fully explained by the reduced windspeed alone, indicating that the length scale of plume elements is also longer. Some of the differences observed in stable conditions, particularly the longer time scales, are shown to substantially increase the potential hazard due to fluctuations in practical applications. A conceptual model of plume dispersion is described, which explains the observed plume structure under different conditions by relating it to the turbulent velocity spectra.     

基于大涡模拟评估GRAPES模式对对流边界层的模拟性能

为检验GRAPES半拉格朗日动力框架在大涡尺度上的模拟性能,为未来发展千米及其以下高分辨尺度的数值模式奠定基础,并构造GRAPES大涡模式以检验和发展边界层湍流参数化提供科学工具。通过在GRAPES模式中加入Smagorinsky-Lilly小尺度湍涡参数化,并将模式分辨率提高至50 m,构建GRAPES大涡模式(GRAPES_LES),以便分析GRAPES模式在大涡尺度上的适用性。同时利用广泛应用的已有大涡模式UCLA_LES作为参考,通过对干对流边界层湍流的模拟分析及与UCLA_LES模拟结果的对比,得出如下主要结论:GRAPES半拉格朗日动力框架能够模拟出与已有的大涡模式相似的边界层湍流特征;同时,通过分析也证明GRAPES存在由于采用半拉格朗日平流计算而带来过度耗散的问题:当使用相同的滤波尺度(Smagorinsky 常数)时,GRAPES_LES模拟出的速度场更为平滑,小尺度湍流结构过于光滑,通过对湍流能量的能谱分析更清楚地表明了这一点。进一步,对不同的Smagorinsky常数(对应不同的滤波尺度)进行了敏感性试验,表明可以通过改变滤波尺度,有效地缓解半拉格朗日框架隐含的耗散问题,得到更接近UCLA_LES所模拟的湍流特征。     

On the choice of diffusion schemes in numerical simulation of crystallizing aerosol propagation in the cloud medium

Considered is the parameterization of turbulent distribution of aerosol in the cloud medium, basing on three-stage diffusion model, in which the first, second, and third stages describe diffusion at small, intermediate (comparable with the Lagrangian scale), and large time intervals, respectively. Used are semiempirical diffusion equations which take into account the features of diffusion at different stages. The effects of temperature and humidity stratification on vertical turbulent regime are suggested to be described by means of the modified Richardson number. The distribution of crystallizing aerosol in a thick convective cloud is numerically studied from the point of view of hail-protection technology.     

The dispersion of chemically reactive species in the atmospheric boundary layer

Summary The role of turbulence in the dispersion of atmospheric pollutants that react with linear (decay) and nonlinear (second-order) chemical reactions is examined. The most relevant processes that drive the reactivity of species emitted in a surface area or released by a point source are studied by deriving the dimensionless scaling numbers from equations for the atmospheric turbulent reacting flow. The first number is the ratio of the time scale of turbulence to the time scale of the chemical reaction, namely the Damköhler number. The second number is the ratio of the concentrations of the species present in the chemical transformations. In this paper, model results and experimental studies of turbulent reacting flows in the atmospheric boundary layer are presented to show the modifications and control exerted by turbulence on the atmospheric chemistry as a function of these numbers and processes. We also discuss how the chemical transformation is affected when species are in a state of chemical equilibrium.By studying the plume dispersion of a reactant, that decays with a simple chemical reaction, one can analyse the dependence of concentration fluctuations on the Damköhler number. The study is extended to plumes that react nonlinearly. In such reacting systems, the large gradients and segregation of species result in a significant reduction in the reaction rates. Because of this modification, the chemistry of species related to NO_x and HO_x can be very different from the chemistry in conditions where the species are uniformly mixed. The lack of complete observational evidence is hampering our understanding of these processes and our evaluation of numerical modelling results. Finally, we discuss briefly how to represent, in the form of a parameterization, the effect that turbulence can have on the reactivity of species emitted by a point source or an area source.     

The Influence Of Vertical Wind Direction Shear On Dispersion In The Convective Boundary Layer, And Its Incorporation In Coastal Fumigation Models

The mean concentration distributionwithin a plume released from a point source in the atmosphericboundary layer can be greatly influenced by the systematic turningof wind with height (i.e. vertical wind direction shear). Such aninfluence includes a deflection of the plume centroid, with anassociated shearing of the vertical plume cross-section, and anenhancement of dispersion, in the horizontal plane. Wind directionshear is normally not accounted for in coastal fumigation models,although dispersion observations with shear acting as acontrolling parameter are not uncommon. A three-dimensionalLagrangian stochastic model is used to investigate the influenceof uniform wind direction shear on the diffusion of a point-sourceplume within the horizontally homogeneous convective boundarylayer, with the source located at the top of the boundary layer.Parameterisations are developed for the plume deflection andenhanced dispersion due to shear within the framework of aprobability density function (PDF) approach, and compared with theLagrangian model results. These parameterisations are thenincorporated into two applied coastal fumigation models: a PDFmodel, and a commonly used model that assumes uniform andinstantaneous mixing in the vertical direction. The PDF modelrepresents the vertical mixing process more realistically. A moreefficient version of the PDF model, which assumes a well-mixedconcentration distribution in the vertical at large times, isapplied to simulate sulfur dioxide data from the Kwinana CoastalFumigation Study. A comparison between the model results and thedata show that the model performs much better when the wind-sheareffects are included.     

Dispersion of a Point-Source Release of a Passive Scalar Through an Urban-Like Array for Different Wind Directions

The dispersion of a point-source release of a passive scalar in a regular array of cubical, urban-like, obstacles is investigated by means of direct numerical simulations. The simulations are conducted under conditions of neutral stability and fully rough turbulent flow, at a roughness Reynolds number of Re _τ  = 500. The Navier–Stokes and scalar equations are integrated assuming a constant rate release from a point source close to the ground within the array. We focus on short-range dispersion, when most of the material is still within the building canopy. Mean and fluctuating concentrations are computed for three different pressure gradient directions (0°, 30°, 45°). The results agree well with available experimental data measured in a water channel for a flow angle of 0°. Profiles of mean concentration and the three-dimensional structure of the dispersion pattern are compared for the different forcing angles. A number of processes affecting the plume structure are identified and discussed, including: (i) advection or channelling of scalar down ‘streets’, (ii) lateral dispersion by turbulent fluctuations and topological dispersion induced by dividing streamlines around buildings, (iii) skewing of the plume due to flow turning with height, (iv) detrainment by turbulent dispersion or mean recirculation, (v) entrainment and release of scalar in building wakes, giving rise to ‘secondary sources’, (vi) plume meandering due to unsteady turbulent fluctuations. Finally, results on relative concentration fluctuations are presented and compared with the literature for point source dispersion over flat terrain and urban arrays.     

A parameterization of the stable atmospheric boundary layer

Two formulations of the stable atmospheric boundary layer are proposed for use in weather forecasting or climate models. They feature the log-linear profile near the surface, but are free from the associated critical Richardson number. The diffusion coefficients in the Ekman layer are a natural extension of the surface layer. They are locally determined using wind shear in one case and turbulent kinetic energy in the other. The parameterizations are tested in a one-dimensional model simulating the evolution of the nocturnal boundary layer with and without radiative cooling. Both formulations give very similar results, except near the top of the boundary layer where the transition to the free atmosphere is smoother with the wind shear formulation. A distinctive feature of these schemes is that they retain their simulating skill when resolution is reduced. This is verified for a wide range of situations. In practice, this means that there is no need for a large-scale model to have a level below 50 m or so.     

On the dispersion parameters of plumes from tall stacks in a shoreline environment

Experimental data for buoyant plumes released from high sources into layers having little ambient turbulence show that plume dispersion parameters vary in a manner similar to that during initial plume rise. This is consistent with general plume rise theory. Dispersion of plumes from tall stacks in a shoreline environment where a thermal internal boundary layer is formed often demonstrates this behaviour.     

Extension of a fluctuating plume model of tracer dispersion to a sheared boundary layer and to a large array of obstacles

Fluctuating plume models provide a useful conceptual paradigm in the understanding of plume dispersion in a turbulent flow. In particular, these models have enabled analytical predictions of higher-order concentration moments, and the form of the one-point concentration probability density function (PDF). In this paper, we extend the traditional formalism of these models, grounded in the theory of homogeneous and isotropic turbulent flow, to two cases: namely, a simple sheared boundary layer and a large array of regular obstacles. Some very high-resolution measurements of plume dispersion in a water channel, obtained using laser-induced fluorescence (LIF) line-scan techniques are utilised. These data enable us to extract time series of plume centroid position (plume meander) and dispersion in the relative frame of reference in unprecedented detail. Consequently, experimentally extracted PDFs are able to be directly compared with various theoretical forms proposed in the literature. This includes the PDF of plume centroid motion, the PDF of concentration in the relative frame, and a variety of concentration moments in the absolute and relative frames of reference. The analysis confirms the accuracy of some previously proposed functional forms of model components used in fluctuating plume models, as well as suggesting some new forms necessary to deal with the complex boundary conditions in the spatial domain.     

An Advanced Puff Model Based On A Mixed Eulerian/Lagrangian Approach For Turbulent Dispersion In The Convective Boundary Layer

An advanced model aimed at describing the problem of dispersion in theconvective boundary layer is proposed. The pollutant particles are groupedin clusters and modelled as Gaussian puffs. The expansion of each puff ismodelled according to the concept of relative dispersion and expressed interms of the spectral properties of the energy containing eddies of the turbulent field. The centre of mass of each puff is moved along a stochastic trajectory, obtained using a Lagrangian stochastic model and filtering the velocity with a recursive Kalman filter. At any instant, a filtering procedure, depending both on travel time and on puff size, acts to select spectral components involved in the expansion and in the meandering of the puff. Such an approach requires only a moderate number of puff releases, so that the proposed model is faster to run than a standard Lagrangian model. On the other hand, unlike the traditional puff model, it allows us to simulate both expansion and meandering of the puff. Therefore, it is well suited to simulate dispersion when the turbulent structures are larger thanthe plume dimensions, as for example in convective conditions. Being based onspectral formulations in both Eulerian and Lagrangian parts, the model is consistent in all the turbulent parameterizations utilised. Comparisons with a standard Lagrangian particle model as well as with a classical convective experimental dataset show good performance of the proposed model.     

Three-dimensional air pollutant modeling in the lower atmosphere

A steady-state, three-dimensional turbulent diffusion equation describing the concentration distribution of an air pollutant from an elevated point source in the lower atmosphere is solved analytically. The same formulation can be used to obtain solutions from line, area or other kinds of sources. The solution is developed for the cases in which the velocity, vertical and lateral diffusivities are given by the power law. The model preserves the beauty of analytical solution without sacrificing much on the accuracy of approximating the velocity and eddy diffusivities. Methods of evaluating the parameters, which are required for the model applications, are discussed. Results indicate that the ratio of the plume width to the plume length increases with decreasing stability and with increasing source height. These consequences are in response to the variations of the size of eddies in the vertical direction.     

A modeling study of an Antarctica polynya event

An event of polynya at Terra Nova Bay (TNB), occurring from 15 July to 17 July 2006, is simulated by a recent version of the mesoscale Eta model. Simulation results and observational data describe the surface conditions during the period. The spatial and temporal structure of the atmospheric boundary layer in response to the warm area of the polynya is also investigated. Numerical experiments show that the latter influences significantly the wind intensity, the temperature and the specific humidity of the air over Terra Nova Bay. The significant heating of the low atmosphere results in a three-dimensional anomaly in the baric field, which extends far in the Ross Sea, embedded in the complex pressure field obtained by the Eta model also without taking into account the polynya. A turbulent kinetic energy plume, indicating turbulent mixing, and an increased vertical diffusion of horizontal momentum are also simulated over the polynya. The downward flux of high momentum air and the modified pressure gradient force change the wind speed at low level over TNB.     

不同下垫面的大气边界层湍流扩散特征

利用海洋、戈壁、城郊和城市4种不同下垫面上的湍流观测资料,运用多尺度湍流概念和方法计算分析了风速谱和扩散参数等湍流特征量。结果表明:下垫面的粗糙度对大气边界层湍流特征影响很大,随着粗糙度的增大,速度谱的峰值频率向高频方向移动,扩散参数增大。