Comparing Two Implementations of a Micromixing Model. Part II: Canopy Flow

The sequential particle micromixing model (SPMMM) is used to estimate concentration fluctuations in plumes dispersing into a canopy flow. SPMMM uses the familiar single-particle Lagrangian stochastic (LS) trajectory framework to pre-calculate the required conditional mean concentrations, which are then used by an interaction by exchange with the conditional mean (IECM) micromixing model to predict the higher-order fluctuations of the scalar concentration field. The predictions are compared with experimental wind-tunnel dispersion data for a neutrally stratified canopy flow, and with a previously reported implementation using simultaneous particle trajectories. The two implementations of the LS–IECM model are shown to be largely consistent with one another and are able to simulate dispersion in a canopy flow with fair to good accuracy.     

Comparing Two Implementations of a Micromixing Model. Part I: Wall Shear-Layer Flow

A Lagrangian stochastic (LS) micromixing model is used for estimating concentration fluctuations in plumes of a passive, non-reactive tracer dispersing from elevated and ground-level compact sources into a neutral wall shear-layer flow. SPMMM (for sequential particle micromixing model) implements the familiar IECM (interaction by exchange with the conditional mean) micromixing scheme. The parametrization of the scalar micromixing time scale is identical to that proposed in a previously reported LS–IECM model (Cassiani et al., Atmos Environ 39:1457–1469, 2005a). However, while SPMMM is mathematically equivalent to the previously reported model, it differs in its numerical implementation: SPMMM releases N independent particles sequentially, whereas the previously reported model releases N independent particles simultaneously. In both implementations, the trajectories of the N particles are governed by single-point velocity statistics. The sequential particle implementation is computationally efficient, but cannot be applied to the case of reacting species. Results from both implementations are compared to experimental wind-tunnel dispersion data and to each other.     

First-Order Inconsistencies Caused by Rogue Trajectories

A theoretical requirement of the Interaction by Exchange with the Conditional Mean (IECM) micromixing model is that the mean concentration field produced by it must be consistent with the mean concentration field produced by a traditional Lagrangian stochastic (LS) marked particle model. We examine the violation of this requirement that occurs in a coupled LS–IECM model when unrealistically high particle velocities occur. No successful strategy was found to mitigate the effects of these rogue trajectories. It is our hope that this work will provide renewed impetus for investigation into rogue trajectories and methods to eliminate them from LS models.     

Micromixing modelling of concentration fluctuations in inhomogeneous turbulence in the convective boundary layer

The micromixing technique, widely used in engineering calculations of mixing and chemical reaction, is extended to atmospheric boundary-layer flows. In particular, a model based on the interaction-by-exchange-with-the-conditional-mean (IECM) micromixing approach is formulated to calculate concentration fluctuation statistics for a line source and a point source in inhomogeneous and non-Gaussian turbulence in the convective boundary layer. The mixing time scale is parameterised as a linear function of time with the intercept value determined by the source size at small times. Good agreement with laboratory data for the intensity of concentration fluctuations is obtained with a value of 0.9 for the coefficient of the linear term in the time-scale parameterisation for a line source, and a value of 0.6 for a point source. Calculation of higher-order moments of the concentration field for a line source shows that non-Gaussian effects persist into the vertically well-mixed region. The cumulative distribution function predicted by the model for a point source agrees reasonably well with laboratory data, especially in the far field. In the limit of zero mixing time scale, the model reduces to a meandering plume model, thus enabling the concentration variance to be partitioned into meandering and relative components. The meandering component is shown to be more persistent for a point source than for a line source.     

Concentration fluctuations according to fluctuating plume models in one and two dimensions

Concentration probability density functions (pdfs) calculated according to fluctuating plume models in one- and two-dimensions, representing the limiting cases of one-dimensional dispersion from a line source or a point source in strongly anisotropic turbulence and of axisymmetric dispersion from a point source in isotropic turbulence, are discussed and analyzed in terms of the location of the sampling point within the mean plume and of the ratio, s/m, of the standard deviations for relative dispersion and meandering.In both cases, the pdfs cover the finite concentration range from zero to C ₀, the centreline concentration of the instantaneous plume. The main difference between them is that whereas the 2-D pdf is always unimodal, the 1-D pdf has a singularity at C ₀ which under some circumstances results in a bimodal form. However, the probability associated with this singularity is not always significant. Differences of practical importance in the shape of the pdfs occur mainly for centreline or near-centreline sampling locations when meandering is not too much larger than relative dispersion (1 < m ²/s² < 10) and for sampling locations a distance of order s from the centreline when relative dispersion is not too much larger than meandering (1 < s ²/m² < 5).Comparison against wind tunnel measurements not too far downstream of a line source in grid turbulence shows that the 1-D model reproduces the essential features and trends of the measurements. Under appropriate circumstances the measurements show the bimodal pdf predicted by the 1-D model (but not by the 2-D model) confirming that the effect of the anisotropy in the source distribution is observable.Present address: School of Mechanical Engineering, Aristotle University, Thessaloniki, 54006 Thessaloniki, Greece.     

A wind tunnel study of turbulent flow around single and multiple windbreaks, part I: Velocity fields

This paper describes wind-tunnel experiments on the flow around single and multiple porous windbreaks (height H), sheltering a model plant canopy (height H/3). The mean wind is normal to the windbreaks, which span the width of the wind tunnel. The incident turbulent flow simulates the adiabatic atmospheric surface layer. Five configurations are examined: single breaks of three solidities (low, medium, high; solidity = 1 - porosity), and medium-solidity multiple breaks of streamwise spacing 12H and 6H. The experimental emphases are on the interactions of the windbreak flow with the underlying plant canopy; the effects of solidity; the differences in shelter between single and multiple windbreaks; and the scaling properties of the flow. Principal results are: (1) the "quiet zones" behind each windbreak are smaller in multiple than single arrays, because of the higher turbulence level in the very rough-wall internal boundary layer which develops over the multiple arrays. Nevertheless, the overall shelter effectiveness is higher for multiple arrays than single windbreaks because of the "nonlocal shelter" induced by the array as a whole. (2) The flow approaching the windbreak decelerates above the canopy but accelerates within the canopy, particularly when the windbreak solidity is high. (3) A strong mixing layer forms just downwind of the top of each windbreak, showing some of the turbulence and scaling properties of the classical mixing layer formed between uniform, coflowing streams. (4) No dramatic increase in turbulence levels in the canopy is evident at the point where the deepening mixing layer contacts the canopy (around x/H = 3) but the characteristic inflection in the canopy wind profile is eliminated at this point.     

A Note on k - ε Modelling of Vegetation Canopy Air-Flows

The k - turbulence model is a standard of computational software packages for engineering, yet its application to canopy turbulence has not received comparable attention. This is probably due to the additional source (and/or sink) terms, whose parameterization remained uncertain. This model must include source terms for both turbulent kinetic energy (k) and the viscous dissipation rate (), to account for vegetation wake turbulence budget. In this note, we show how Kolmogorov's relation allows for an analytical solution to be calculated within the portion of a dense and homogeneous canopy where the mixing length does not vary. By substitution within model equations, this solution allows for a set of constraints on source term model coefficients to be derived.Those constraints should meet both Reynolds averaged Navier–Stokes equationsand large-eddy simulation sub-grid scale turbulence modelling requirements.Although originating from within a limited portion of the canopy, the predictedcoefficients values must be valid elsewhere in order to make the model capable of predicting the whole canopy-layer flow with a single set of constants.     

The Volumetric Particle Approach for Concentration Fluctuations and Chemical Reactions in Lagrangian Particle and Particle-grid Models

A new approach is proposed to predict concentration fluctuations in the framework of one-particle Lagrangian stochastic models. The approach is innovative since it allows the computation of concentration fluctuations in dispersing plumes using a Lagrangian one-particle model with micromixing but with no need for the simulating of background particles. The extension of the model for the treatment of chemically reactive plumes is also accomplished and allows the computation of plume-related chemical reactions in a Lagrangian one-particle framework separately from the background chemical reactions, accounting for the effect of concentration fluctuations on chemical reactions in a general, albeit approximate, manner. These characteristics should make the proposed approach an ideal tool for plume-in-grid calculations in chemistry transport models. The results are compared to the wind-tunnel experiments of Fackrell and Robins (J Fluid Mech, 117:1–26, 1982) for plume dispersion in a neutral boundary layer and to the measurements of Legg et al. (Boundary-Layer Meteorol, 35:277–302, 1986) for line source dispersion in and above a model canopy. Preliminary reacting plume simulations are also shown comparing the model with the experimental results of Brown and Bilger (J Fluid Mech, 312:373–407, 1996; Atmos Environ, 32:611–628, 1998) to demonstrate the feasibility of computing chemical reactions in the proposed framework.     

A numerical study of the air flow within and around a single tree

A three-dimensional non-hydrostatic numerical model has been used to investigate the air flow and turbulence around a single tree.The results for velocity and turbulence distributions have been compared with available data from windtunnel experiments; the agreement is satisfactory.Simulations have been carried out for different meteorological conditions (wind speed, thermal stratification) as well as for different canopy characteristics (stem height, crown diameter, crown height, porosity).Dedicated to Prof. F. Wippermann on the occasion of his 65th birthday.     

A two-dimensional trajectory-simulation model for non-Gaussian,inhomogeneous turbulence within plant canopies

We report a two-dimensional (alongwind u, vertical w) trajectory-simulation model, consistent with Thomson's (1987) well-mixed criteria, that allows for the non-Gaussian turbulence typical of flow within a plant canopy. The effect of non-Gaussian turbulence was examined by formulating a non-Gaussian u, w joint probability density function (PDF) as the sum of two Gaussian joint-PDFs. The resultant PDF reproduced the desired means, variances, skewnesses, and kurtoses, and the correct covariance. In prediction of the location of maximum concentration downwind of a line source in homogeneous, slightly non-Gaussian turbulence, it proved advantageous to incorporate skewness and kurtosis. However, in the case of inhomogeneous, highly non-Gaussian turbulence, the addition of skewness and kurtosis in the model resulted in substantially worse agreement with measurements than the results of the model using Gaussian PDFs. This may be due to inaccuracy in our PDF formulation. Dispersion predictions from the model with Gaussian PDFs were generally not statistically different from measurements. These results indicate that a two-dimensional Gaussian trajectory-simulation approach is adequate to predict mean concentrations and fluxes resulting from sources within plant canopies.     

Large-Eddy Simulation of Flow Over a Vegetation-Like Canopy Modelled as Arrays of Bluff-Body Elements

Turbulent flow over a vegetation canopy under neutral atmospheric conditions is investigated using large-eddy simulation. Each model tree, which consists of a sphere-shaped tree crown and a cylindrical trunk, is fully resolved. The resulting turbulence statistics and the drag force on the vegetation agree well with measurements from the corresponding wind-tunnel experiment described by Böhm et al. (Boundary-Layer Meteorol, 146:393–419, 2013). Statistically, this kind of model canopy exhibits both vegetation and bluff-body-flow characteristics. The time-averaged flow skims over the top of the underlying canopy, forming a low-momentum recirculation zone on the lee-side of the bluff elements, which causes significant dispersive stress within the canopy layer. Two other numerical representations of vegetation canopies, referred to as the drag-element and drag-crown approaches, have also been developed to assess the performance of simulations. Turbulence statistics suggest that the canopy shear layer interferes with wakes behind stems and crowns. The drag-crown approach yields better agreement between numerical results and experimental measurements than does the traditional drag-element approach, thus providing a promising numerical model for simulating canopy turbulence.     

Footprints in Homogeneously and Heterogeneously Driven Boundary Layers Derived from a Lagrangian Stochastic Particle Model Embedded into Large-Eddy Simulation

A Lagrangian stochastic (LS) model, which is embedded into a parallelised large-eddy simulation (LES) model, is used for dispersion and footprint evaluations. For the first time an online coupling between LES and LS models is applied. The new model reproduces concentration patterns, which were obtained in prior studies, provided that subgrid-scale turbulence is included in the LS model. Comparisons with prior studies show that the model evaluates footprints successfully. Streamwise dispersion leads to footprint maxima that are situated less far upstream than previously reported. Negative flux footprints are detected in the convective boundary layer (CBL). The wide range of applicability of the model is shown by applying it under neutral and stable stratification. It is pointed out that the turning of the wind direction with height leads to a considerable dependency of source areas on height. First results of an application to a heterogeneously heated CBL are presented, which emphasize that footprints are severely affected by the inhomogeneity.     

Improving the Representation of the Nocturnal Near-Neutral Surface Layer in the Urban Environment with a Mesoscale Atmospheric Model

A new approach to improve the representation of surface processes in the Global Environmental Multiscale (GEM) atmospheric model associated with the exchanges between the urban canopy and the atmosphere is presented. Effects of the urban canopy on the evolution of surface-layer wind, temperature, moisture, and turbulence are directly parametrized in order to allow realistic interactions between the canopy elements (i.e., roofs, roads, and walls) and the atmosphere at GEM’s multiple vertical levels that are positioned inside the canopy. Surface energy budgets as implemented in the Town Energy Balance (TEB) scheme have been used to determine temperatures of the urban canopy elements for the proposed multilayer scheme. Performance of the multilayer scheme is compared against standard implementations of the TEB scheme for one nighttime intensive observation period of the Joint Urban 2003 experiment held in Oklahoma City, USA. Although the new approach is found to have a negligible impact on urban surface-layer wind profiles, it improves the prediction of near-neutral nocturnal profiles of potential temperature close to the surface. The urban heat island effect is simulated with a better accuracy by the multilayer approach. The horizontal temperature gradient across the central business district of the city along the direction of flow is also reasonably well captured by the proposed scheme.     

Large-eddy simulation of plant canopy flows using plant-scale representation

Turbulent flow in a corn canopy is simulated using large-eddy simulation (LES) with a Lagrangian dynamic Smagorinsky model. A new numerical representation of plant canopies is presented that resolves approximately the local structure of plants and takes into account their spatial arrangement. As a validation, computational results are compared with experimental data from recent field particle image velocimetry (PIV) measurements and two previous experimental campaigns. Numerical simulation using the traditional modelling method to represent the canopy (field-scale approach) is also conducted as a comparison to the plant-scale approach. The combination of temporal PIV data, LES and spatial PIV data allows us to couple a wide range of relevant turbulence scales. There is good agreement between experimental data and numerical predictions using the plant-scale approach in terms of various turbulence statistics. Within the canopy, the plant-scale approach also allows the capture of more details than the field-scale approach, including instantaneous gusts that penetrate deep inside the canopy.     

A Lagrangian Stochastic Model for Heavy Particle Dispersion in the Atmospheric Marine Boundary Layer

The dispersion of heavy particles and pollutants is often simulated with Lagrangian stochastic (LS) models. Although these models have been employed successfully over land, the free surface at the air-sea interface complicates the implementation of traditional LS models. We present an adaptation of traditional LS models to the atmospheric marine boundary layer (MBL), where the bottom boundary is represented by a realistic wavy surface that moves and deforms. In addition, the correlation function for the turbulent flow following a particle is extended to the anisotropic, unsteady case. Our new model reproduces behaviour for Lagrangian turbulence in a stratified air flow that departs only slightly from the expected behaviour in isotropic turbulence. When solving for the trajectory of a heavy particle in the air flow, the modelled turbulent forcing on the particle also behaves remarkably well. For example, the spectrum of the turbulence at the particle location follows that of a massless particle for time scales approximately larger than the Stokes’ particle response time. We anticipate that this model will prove especially useful in the context of sea-spray dispersion and its associated momentum, sensible and latent heat, and gas fluxes between spray droplets and the atmosphere.     

Probability Density Function Integrals And Small-Scale Mixing In Turbulent Plumes

We analyse cross-wind-integrated statistics of theconcentration field of a conserved scalar for pointand line sources in grid turbulence. In particular,using wind-tunnel measurements we calculate thecross-wind integrated probability density function(pdf) for the scalar concentration. We then use thatquantity in the exact evolution equation for the pdfto calculate the cross-wind integrated mean of therate of dissipation of scalar variance, conditional onthe scalar concentration. Much of the variation ofthese statistics with distance downstream is accountedfor by scaling with concentration, length and timescales based on the development of the mean plume.This scaling thus suggests some simple practicalparameterisations of these statistics in terms ofmean-field quantities. One of the motivations for thiswork is to find a simple parameterisation for thescalar dissipation that can be used for modellingchemical reactions in plumes.We also consider the cross-wind integral of the firstfew moments of the concentration field and show thatthe integration greatly simplifies the budgets forthese moments. Thus the first moment is just thedownstream flux of the scalar, which is constant. Thesecond moment budget provides a check on the meandissipation estimated directly from the pdf evolutionequation.     

一方程模型在植被冠层中的应用

利用建立的一方程模式对植被气象场和湍流场进行了模拟计算。结果表明：在植被叶面积密度最大值处，由植被叶面积产生的阻力使风速急剧减小，而后平缓接近地面风速值。Reynolds应力从植被顶部向下剧烈减小，在Z／Hc值为0．4～0．6之间衰减迅速；冠层下部风速小且切变弱，湍流未能充分发展，因而动量输送甚微。由于在植被内部环境中存在着湍流通量的辐散或辐合的现象，湍流强度从植被底部开始由下至上逐渐增强，在冠层上方则基本保持不变。     

对流边界层中地面源的铅直扩散模拟

对流边界层(CBL)中高架源浓度轴线下倾而地面源浓度轴线上升。本文分析了形成这种相反的几何图象的机理，从而认为:只要能够恰当模拟CBL中铅直湍流结构的特征，一个粒子随机扩散模式应既能模拟高架源的扩散，也能模拟地面源的扩散。因此，用作者早先模拟高架源的模式模拟了地面源的铅直扩散，同样获得较好的结果，进一步证实了模式的有效性。本文还应用Thomson准则检验了模式，讨论了它在理论上的合理性和适用范围。     

Characteristics of canopy turbulence over a deciduous forest on complex terrain

Canopy turbulence plays an important role in mass and energy exchanges at the canopy-atmosphere interface. Despite extensive studies on canopy turbulence over a flat terrain, less attention has been given to canopy turbulence in a complex terrain. The purpose of this study is to scrutinize characteristics of canopy turbulence in roughness sublayer over a hilly forest terrain. We investigated basic turbulence statistics, conditionally sampled statistics, and turbulence spectrum in terms of different atmospheric stabilities, wind direction and vertical structures of momentum fluxes. Similarly to canopy turbulence over a homogeneous terrain, turbulence statistics showed coherent structure. Both quadrant and spectrum analysis corroborated the role of intermittent and energetic eddies with length scale of the order of canopy height, regardless of wind direction except for shift of peak in vertical wind spectrum to relatively high frequency in the down-valley wind. However, the magnitude of the momentum correlation coefficient in a neutral condition was smaller than typical value over a flat terrain. Further scrutiny manifested that, in the up-valley flow, temperature skewness was larger and the contribution of ejection to both momentum and heat fluxes was larger compared to the downvalley flow, indicating that thermal instability and weaker wind shear in up-valley flow asymmetrically affect turbulent transport within the canopy.