Pollutant Dispersion in Boundary Layers Exposed to Rural-to-Urban Transitions: Varying the Spanwise Length Scale of the Roughness

Both large-eddy simulations (LES) and water-tunnel experiments, using simultaneous stereoscopic particle image velocimetry and laser-induced fluorescence, have been used to investigate pollutant dispersion mechanisms in regions where the surface changes from rural to urban roughness. The urban roughness was characterized by an array of rectangular obstacles in an in-line arrangement. The streamwise length scale of the roughness was kept constant, while the spanwise length scale was varied by varying the obstacle aspect ratio l / h between 1 and 8, where l is the spanwise dimension of the obstacles and h is the height of the obstacles. Additionally, the case of two-dimensional roughness (riblets) was considered in LES. A smooth-wall turbulent boundary layer of depth 10h was used as the approaching flow, and a line source of passive tracer was placed 2h upstream of the urban canopy. The experimental and numerical results show good agreement, while minor discrepancies are readily explained. It is found that for \(l/h=2\) the drag induced by the urban canopy is largest of all considered cases, and is caused by a large-scale secondary flow. In addition, due to the roughness transition the vertical advective pollutant flux is the main ventilation mechanism in the first three streets. Furthermore, by means of linear stochastic estimation the mean flow structure is identified that is responsible for street-canyon ventilation for the sixth street and onwards. Moreover, it is shown that the vertical length scale of this structure increases with increasing aspect ratio of the obstacles in the canopy, while the streamwise length scale does not show a similar trend.     

Effects of a Fence on Pollutant Dispersion in a Boundary Layer Exposed to a Rural-to-Urban Transition

Simultaneous particle-image velocimetry and laser-induced fluorescence combined with large-eddy simulations are used to investigate the flow and pollutant dispersion behaviour in a rural-to-urban roughness transition. The urban roughness is characterized by an array of cubical obstacles in an aligned arrangement. A plane fence is added one obstacle height h upstream of the urban roughness elements, with three different fence heights considered. A smooth-wall turbulent boundary layer with a depth of 10h is used as the approaching flow, and a passive tracer is released from a uniform line source 1h upstream of the fence. A shear layer is formed at the top of the fence, which increases in strength for the higher fence cases, resulting in a deeper internal boundary layer (IBL). It is found that the mean flow for the rural-to-urban transition can be described by means of a mixing-length model provided that the transitional effects are accounted for. The mixing-length formulation for sparse urban canopies, as found in the literature, is extended to take into account the blockage effect in dense canopies. Additionally, the average mean concentration field is found to scale with the IBL depth and the bulk velocity in the IBL.     

Effects of the Upstream-Flow Regime and Canyon Aspect Ratio on Non-linear Interactions Between a Street-Canyon Flow and the Overlying Boundary Layer

Large-scale structures within a rough-wall boundary layer generated over a cube array have recently been linked to small-scale fluctuations close to the roughness through a dynamical mechanism similar to amplitude modulation. Demonstrating the existence of this mechanism for different roughness types is a crucial step towards the development of a generic model for wind fluctuations in the urban canopy. Here the influence of the upstream roughness geometry (two-dimensional (2D) and three-dimensional (3D)) and planform packing density (\( \lambda_{p} \)) and street-canyon aspect ratio on the non-linear interactions between large-scale momentum regions and the small scales induced by the presence of the roughness is studied within a wind tunnel using combined particle-image velocimetry and hot-wire anemometry. A multi-time delay linear stochastic estimation is used to decompose the flow into large scales that participate in modulation and the remaining small scales. Using three different upstream roughness configurations composed of either 3D cubes or 2D rectangular blocks it is shown that the upstream roughness configuration has an influence on the non-linear interactions in the rough-wall boundary layer. Analysis of the turbulence skewness decomposition shows a change in the location of the maximum of the term \( \overline{{u_{L}^{\prime} u_{S}^{\prime 2}}} \), which represents the influence of the large-scale momentum regions on the small scales, whilst the temporal correlation shows a modification of the interaction located closer to the roughness with a change from 3D to 2D roughness. Furthermore, a two-point spatio–temporal correlation demonstrates that the non-linear relationship is significantly modified in the wake-interference-flow regime compared to the skimming-flow regime. Through skewness decomposition and temporal correlations the canyon aspect ratio is shown to have no influence on the non-linear interactions, indicating that the mechanism depends only on the flow developing upstream. Finally, although the upstream roughness configuration is shown to influence the non-linear interactions, the nature of the mechanism remains the same in all configurations.     

Estimation of atmospheric surface layer parameters and numerical simulation using MM5 at Mangalore, West Coast of India

Atmospheric surface layer meteorological observations obtained from 20-m-high meteorological tower at Mangalore, situated along the west coast of India are used to estimate the surface layer scaling parameters of roughness length (z _o) and drag coefficient (C _D), surface layer fluxes of sensible heat and momentum. These parameters are computed using the simple flux–profile relationships under the framework of Monin–Obukhov (M–O) similarity theory. The estimated values of z _o are higher (1.35–1.54 m) than the values reported in the literature (>0.4–0.9 m) probably due to the undulating topography surrounding the location. The magnitude of C _D is high for low wind speed (<1.5 m s^?1) and found to be in the range 0.005–0.03. The variations of sensible heat fluxes (SHF) and momentum fluxes are also discussed. Relatively high fluxes of heat and momentum are observed during typical days on 26–27 February 2004 and 10–11 April 2004 due to the daytime unstable atmospheric conditions. Stable or near neutral conditions prevail after 1700 h IST with negative SHF. A mesoscale model PSU/NCAR MM5 is run using a high-resolution (1 km) grid over the study region to examine the influence of complex topography on the surface layer parameters and the simulated fluxes are compared with estimated values. Spatial variations of the frictional velocity (u _*), C _D, surface fluxes, planetary boundary layer (PBL) height and surface winds are noticed according to the topographic variations in the simulation.     

Impact of surface variations on the momentum flux above the urban canopy

Based on the momentum flux–wind profile relationship of the Monin–Obukhov Similarity (MOS) theory, the observational data from the urban boundary layer field campaign in Nanjing are used to calculate the friction velocity ( $ {u_*} $ ) at the top of the urban canopy and the calculated results are evaluated. The urban surface roughness parameters (the roughness length z ₀ and zero-plane displacement height z _d) are estimated with the Ba method (Bottema’s morphological method). Two different regimes are employed for the calculations. In the homogeneous approach, z ₀ and z _d are averagely derived from the surface elements in the whole study area; while in the heterogeneous approach, z ₀ and z _d are locally derived from the surface elements in the corresponding upwind fetches (or source areas). The calculated friction velocities are compared to the measurement data. The results show that the calculated friction velocities from the heterogeneous approach are in better agreement with the observed values than those from the homogeneous approach are. This study implies that the local roughness parameters can properly represent the dynamical heterogeneity of urban surface, and its application can significantly improve the performance of parameterizations based on the MOS theory in the urban roughness sublayer.     

Velocity Adjustment and Passive Scalar Diffusion in and Above an Urban Canopy in Response to Various Approach Flows

We used wind-tunnel experiments to investigate velocity-field adjustment and scalar diffusion behaviour in and above urban canopies located downwind of various roughness elements. Staggered arrays of rectangular blocks of various heights H and plan area ratios λ_p were used to model the urban canopies. The velocity field in the roughness sublayer (height \({z \lesssim 2H}\)) reached equilibrium at distances proportional to \({\sqrt{L_{\rm c}H}}\) where L _c is the canopy-drag length scale determined as a function of λ_p and the block side length L. A distance of about \({20\sqrt{L_{\rm c}H}}\) was required for adjustment at z = H/2 (in the canopy), and a distance of about \({10\sqrt{L_{\rm c}H}}\) was required at z = 2H (near the top of the roughness sublayer). Diffusion experiments from a ground emission source revealed that differences in upwind roughness conditions had negligible effects on the plume growth near the source (up to a few multiples of L from the source) if the source was located at a fetch F larger than about \({10\sqrt{L_{\rm c}H}}\) from the upwind edge of the canopy. However, at locations farther downwind (more than several multiples of L from the source), upwind conditions had considerable effects on the plume growth. For a representative urban canopy, it was shown that a much larger fetch than required for velocity-field adjustment in the roughness sublayer was necessary to eliminate the effects of upwind conditions on plume widths at 24L downwind from the source.     

Scalar fluxes from urban street canyons. Part I: Laboratory simulation

Flow over urban surfaces depends on surface morphology and interaction with the boundary layer above. However, the effect of the flow on scalar fluxes is hard to quantify. The naphthalene sublimation technique was used to quantify scalar vertical fluxes out of a street canyon under neutral conditions. For an array of eight canyons with aspect ratio H/W=0.75 (here, H is building height and W is the street width), increased flux was observed in the first two or three canyons for moderate and low roughness upstream. This is consistent with predictions of the length scale for initial adjustment of flow to an urban canopy. The flux was constant after the initial adjustment region and thus dependent only on local geometry. For a street canyon in the equilibrium part of the array, each facet of the street canyon was coated with naphthalene to simulate scalar release from street, walls and roof, to evaluate the effect of street canyon geometry on fluxes for H/W=0.25, 0.6, 1 and 2. Fluxes from the roof and downstream wall were considerably larger than fluxes from the street and upstream wall, and only the flux from the downstream wall exhibited a simple decrease with H/W. For each H/W there was a monotonic decrease between downstream wall, street and upstream wall transfer. This suggests that flow decelerates around the recirculation region in the lee of the upstream building, i.e. a recirculating jet rather than a symmetrical vortex. The addition of a second source within the street canyon resulted in reduced fluxes from each facet for H/W>0.25, due to increased concentration of naphthalene in the canyon air.     

应用城市冠层模式研究建筑物形态对城市边界层的影响

文中将城市冠层模式耦合到南京大学城市尺度边界层模式中,通过模拟对比发现,耦合模式对城市地区气温模拟结果更接近于观测值,尤其是对城市地区夜间气温模拟的改进.运用改进耦合模式通过多个敏感性试验的模拟,从城市面积扩张、建筑物高度增加、建筑物分布密度变化等角度研究城市建筑物三维几何形态变化对城市边界层及城市气象环境的影响,试验结果表明:(1)城市面积扩张使得城市下垫面的热通量增大,热力湍流活动增强,动量通量输送增强,城市湍能增大,湍流扩散系数变大,城市气温升高,且对不同时刻城市区域大气层结稳定度均有不同程度的影响.(2)建筑物高度增加增大了城市下垫面的粗糙度和零平面位移.同时也增大了城市街渠高宽比.城市建筑物越高,白天城市地区地表热通量越小,城市上空大气温度越低,平均风速减小,湍能减小;夜间由于高大建筑物释放储热比低矮建筑物要多,其热力湍流相对活跃,地表热通量增大,使得城市区域气温较高.(3)建筑物密度增大,会减小城市下垫面的粗糙度同时增强街渠对辐射的影响.建筑物密度增大在白天会减小地表热通量和动量通量,使城市气温降低,平均风速增大,城市湍流活动能力减弱;夜间城市释放较多储热使得气温较高.     

Momentum- and Heat-Flux Parametrization at Dome C,Antarctica: A Sensitivity Study

An extensive meteorological observational dataset at Dome C, East Antarctic Plateau, enabled estimation of the sensitivity of surface momentum and sensible heat fluxes to aerodynamic roughness length and atmospheric stability in this region. Our study reveals that (1) because of the preferential orientation of snow micro-reliefs (sastrugi), the aerodynamic roughness length \(z_{0}\) varies by more than two orders of magnitude depending on the wind direction; consequently, estimating the turbulent fluxes with a realistic but constant \(z_{0}\) of 1 mm leads to a mean friction velocity bias of \(24\,\%\) in near-neutral conditions; (2) the dependence of the ratio of the roughness length for heat \(z_{0t}\) to \(z_{0}\) on the roughness Reynolds number is shown to be in reasonable agreement with previous models; (3) the wide range of atmospheric stability at Dome C makes the flux very sensitive to the choice of the stability functions; stability function models presumed to be suitable for stable conditions were evaluated and shown to generally underestimate the dimensionless vertical temperature gradient; as these models differ increasingly with increases in the stability parameter z / L, heat flux and friction velocity relative differences reached \(100\,\%\) when \(z/L > 1\); (4) the shallowness of the stable boundary layer is responsible for significant sensitivity to the height of the observed temperature and wind data used to estimate the fluxes. Consistent flux results were obtained with atmospheric measurements at heights up to 2 m. Our sensitivity study revealed the need to include a dynamical parametrization of roughness length over Antarctica in climate models and to develop new parametrizations of the surface fluxes in very stable conditions, accounting, for instance, for the divergence in both radiative and turbulent fluxes in the first few metres of the boundary layer.     

城市建筑物动力冠层方案的引入及应用研究

文中将建筑物动力冠层方案引入到城市边界层精细模式中,该方案描叙了城市地Ⅸ建筑物对气流的拖曳作用以及建筑物形态对城市地区湍流活动动量输送的影响.采用建筑物拖曳法和建筑物动力冠层方案两种不同建筑物动力学处理方法对北京地区东南部进行模拟,通过城市地区水平风速的模拟结果与自动气象站实测资料对比发现模拟结果与实测相当吻合,建筑物动力冠层方案的引入能够更好地实施对城市地区的水平风速分布的数值模拟.分别采用建筑物拖曳法与建筑物动力冠层方案对一个实际小规模城市进行模拟试验,分析表明引入建筑物动力冠层方案可以模拟出小规模城市地区水平风速偏低的现象;还可以模拟出城市地区建筑物动力作用对湍流活动的影响;对湍流动能模拟结果表明比较符合实际分布情形.采用建筑物动力冠层方案对小规模城市地区建筑物高度变化对城市地区的流场及湍流活动的影响进行模拟分析,结果表明:城市地区建筑物高度增加,风速变小,但是高大建筑物底部风速略高于低矮建筑物底部;城市地区大气湍能增加,高大建筑物底部湍能较小,冠层高度内湍能随高度增加而增加;同时城市地Ⅸ的垂直扩散系数也随建筑物高度增加而增加;城市地区污染物排放高度处污染物浓度较低,下游乡村地区地面污染物浓度较低,但高空污染物浓度却较高.     

The Effects of Inner- and Outer-Layer Turbulence in a Convective Boundary Layer on the Near-Neutral Inertial Sublayer Over an Urban-Like Surface

A large-eddy simulation of the atmospheric boundary layer, large enough to contain both an urban surface layer and a convective mixed layer, was performed to investigate inner-layer and outer-layer scale motions. The objective was to determine the applicability of Monin–Obukhov similarity theory to inner-layer motions, to investigate the influence of outer-layer motions on surface-layer structure, as well as to assess the interaction of the two scales of motion. The urban surface roughness consisted of square-patterned cubic buildings of dimension H (40 m). A spatial filter was used to decompose the two scales in the inertial sublayer. The horizontal square filter of size 10H was effective in separating the inner-layer (surface-layer height ≈ 2 H) and outer-layer scales (boundary-layer height δ ≈ 30H), where the Reynolds stress contribution of the inner layer dominates in the logarithmic layer (depth 2H). Similarity coefficients for velocity fluctuations were successfully determined for inner-layer motions in the surface layer, proving the robustness of Monin–Obukhov similarity for surface-layer turbulence. The inner-layer structures exhibit streaky structures that have similar streamwise length but narrower spanwise width relative to the streamwise velocity fluctuation field, consistent with observations from an outdoor scale model. The outer-layer motions to some extent influence the location of ejections and sweeps through updraft and downdraft motions, respectively, thus, disturbing the homogeneity and similarity of inner-layer motions. Although the horizontal averages of the variances and covariance of motions reveal that the Reynolds stresses are dominated by inner-layer structures, the localized influence of the interaction of outer-layer horizontal and inner-layer vertical motions on the Reynolds stress is not insignificant.     

Observations and Modelling of the Atmospheric Boundary Layer Over Sea-Ice in a Svalbard Fjord

Sonic anemometer and profile mast measurements made in Wahlenbergfjorden, Svalbard Arctic archipelago, in May 2006 and April 2007 were employed to study the atmospheric boundary layer over sea-ice. The turbulent surface fluxes of momentum and sensible heat were calculated using eddy correlation and gradient methods. The results showed that the literature-based universal functions underestimated turbulent mixing in strongly stable conditions. The validity of the Monin-Obukhov similarity theory was questionable for cross-fjord flow directions and in the presence of mesoscale variability or topographic effects. The aerodynamic roughness length showed a dependence on the wind direction. The mean roughness length for along-fjord wind directions was (2.4 ± 2.6) × 10⁻⁴ m, whereas that for cross-fjord directions was (5.4 ± 2.8) × 10⁻³ m. The thermal stratification and turbulent fluxes were affected by the synoptic situation with large differences between the 2 years. Channelling effects and drainage flows occurred especially during a weak large-scale flow. The study periods were simulated applying the Weather Research and Forecasting (WRF) model with 1-km horizontal resolution in the finest domain. The results for the 2-m air temperature and friction velocity were good, but the model failed to reproduce the spatial variability in wind direction between measurement sites 3 km apart. The model suggested that wind shear above the stable boundary layer provided a non-local source for the turbulence observed.     

Variability of energy fluxes in relation to the net-radiation of urban and suburban areas: a case study

This paper addresses the relation between the net-radiation (Q ^*) and the ground heat flux (Q _G), the energy stored in the soil ( $\Updelta Q_{\rm S}$ ), and the residual of the energy partition (R = Q ^* ? Q _H  ? Q _E ) of urban and suburban areas of Oklahoma City, USA. These three forms of energy were observed or estimated from observations taken during Joint Urban 2003 Campaign. The database includes net-radiation, soil temperature, ground heat flux, and turbulent fluxes. In most cases the estimates of the energy stored in the soil were obtained by assuming roughly a certain type of soil and an effective soil depth. From the residuals it seems to be possible to distinguish the urban boundary layer from the suburban boundary layer when plotted as a function of net-radiation. Hysteresis coefficients were computed for fits of net-radiation against R, $\Updelta Q_{\rm S}$ and Q _G. In particular, the hysteresis patterns show that Q ^* vs. R represents clearer urban areas or suburban areas under the influence of an urban “plume”. On the other hand, hysteresis curves obtained from $\Updelta Q_{\rm S}$ or Q _G account for better the ground composition. A possible consequence is that the land use of urban areas could be roughly inferred from curve shapes such as Q ^* vs. R, or Q ^* versus another input variable representing the storage term. The objective is to show the variability of the subsurface-related energy fluxes across an urban area using these three different quantities and also to show that $\Updelta Q_{\rm S}, \,Q_{\rm G}$ , or R (and their corresponding hysteresis curves) are likely to be quantitatively different, which have not been clearly stated in the literature.     

Wind profiles,momentum fluxes and roughness lengths at Cabauw revisited

We describe the results of an experiment focusing on wind speed and momentum fluxes in the atmospheric boundary layer up to 200 m. The measurements were conducted in 1996 at the Cabauw site in the Netherlands. Momentum fluxes are measured using the K-Gill Propeller Vane. Estimates of the roughness length are derived using various techniques from the wind speed and flux measurements, and the observed differences are explained by considering the source area of the meteorological parameters. A clear rough-to-smooth transition is found in the wind speed profiles at Cabauw. The internal boundary layer reaches the lowest k-vane (20 m) only in the south-west direction where the obstacle-free fetch is about 2 km. The internal boundary layer is also reflected in the roughness lengths derived from the wind speed profiles. The lower part of the profile (< 40 m) is not in equilibrium and no reliable roughness analysis can be given. The upper part of the profile can be linked to a large-scale roughness length. Roughness lengths derived from the horizontal wind speed variance and gustiness have large footprints and therefore represent a large-scale average roughness. The drag coefficient is more locally determined but still represents a large-scale roughness length when it is measured above the local internal boundary layer. The roughness length at inhomogeneous sites can therefore be determined best from drag coefficient measurements just above the local internal boundary layers directly, or indirectly from horizontal wind speed variance or gustiness. In addition, the momentum and heat fluxes along the tower are analysed and these show significant variation with height related to stability and possibly surface heterogeneity. It appears that the dimensionless wind speed gradients scale well with local fluxes for the variety of conditions considered, including the unstable cases.     

Pollutant Plume Dispersion in the Atmospheric Boundary Layer over Idealized Urban Roughness

The Gaussian model of plume dispersion is commonly used for pollutant concentration estimates. However, its major parameters, dispersion coefficients, barely account for terrain configuration and surface roughness. Large-scale roughness elements (e.g. buildings in urban areas) can substantially modify the ground features together with the pollutant transport in the atmospheric boundary layer over urban roughness (also known as the urban boundary layer, UBL). This study is thus conceived to investigate how urban roughness affects the flow structure and vertical dispersion coefficient in the UBL. Large-eddy simulation (LES) is carried out to examine the plume dispersion from a ground-level pollutant (area) source over idealized street canyons for cross flows in neutral stratification. A range of building-height-to-street-width (aspect) ratios, covering the regimes of skimming flow, wake interference, and isolated roughness, is employed to control the surface roughness. Apart from the widely used aerodynamic resistance or roughness function, the friction factor is another suitable parameter that measures the drag imposed by urban roughness quantitatively. Previous results from laboratory experiments and mathematical modelling also support the aforementioned approach for both two- and three-dimensional roughness elements. Comparing the UBL plume behaviour, the LES results show that the pollutant dispersion strongly depends on the friction factor. Empirical studies reveal that the vertical dispersion coefficient increases with increasing friction factor in the skimming flow regime (lower resistance) but is more uniform in the regimes of wake interference and isolated roughness (higher resistance). Hence, it is proposed that the friction factor and flow regimes could be adopted concurrently for pollutant concentration estimate in the UBL over urban street canyons of different roughness.     

Optimal nudging strategies in regional climate modelling: investigation in a Big-Brother experiment over the European and Mediterranean regions

The objective of this work is to gain a general insight into the key mechanisms involved in the impact of nudging on the large scales and the small scales of a regional climate simulation. A “Big Brother experiment” (BBE) approach is used where a “reference atmosphere” is known, unlike when regional climate models are used in practice. The main focus is on the sensitivity to nudging time, but the BBE approach allows to go beyond a pure sensitivity study by providing a reference which model outputs try to approach, defining an optimal nudging time. Elaborating upon previous idealized studies, this work introduces key novel points. The BBE approach to optimal nudging is used with a realistic model, here the weather research and forecasting model over the European and Mediterranean regions. A winter simulation (1 December 1989–28 February 1990) and a summer simulation (1 June 1999–31 August 1999) with a 50 km horizontal mesh grid have been performed with initial and boundary conditions provided by the ERA-interim reanalysis of the European Center for Medium-range Weather Forecast to produce the “reference atmosphere”. The impacts of spectral and indiscriminate nudging are compared all others things being equal and as a function of nudging time. The impact of other numerical parameters, specifically the domain size and update frequency of the large-scale driving fields, on the sensitivity of the optimal nudging time is investigated. The nudged simulations are also compared to non-nudged simulations. Similarity between the reference and the simulations is evaluated for the surface temperature, surface wind and for rainfall, key variables for climate variability analysis and impact studies. These variables are located in the planetary boundary layer, which is not subject to nudging. Regarding the determination of a possible optimal nudging time, the conclusion is not the same for indiscriminate nudging (IN) and spectral nudging and depends on the update frequency of the driving large-scale fields τ _a . For IN, the optimal nudging time is around τ = 3 h for almost all cases. For spectral nudging, the best results are for the smallest value of τ used for the simulations (τ = 1 h) for frequent update of the driving large-scale fields (3 and 6 h). The optimal nudging time is 3 for 12 h interval between two consecutive driving large-scale fields due to time sampling errors. In terms of resemblance to the reference fields, the differences between the simulations performed with IN and spectral nudging are small. A possible reason for this very similar performance is that nudging is active only above the planetary boundary layer where small-scale features are less energetic. As expected from previous studies, the impact of nudging is weaker for a smaller domain size. However the optimal nudging time itself is not sensitive to domain size. The proposed strategy ensures a dynamical consistency between the driving field and the simulated small-scale field but it does not ensure the best “observed” fine scale field because of the possible impact of incorrect driving large-scale field. This type of downscaling provides an upper bound on the skill possible for recent historical past and twenty-first century projections. The optimal nudging strategy with respect to dynamic downscaling could add skill whenever the parent global model has some level of skill.     

Near-Wall Flow Development After A Step Change In Surface Roughness

An experimental study of the initial flow field downstream of a step change in surface roughness is presented. The roughness length of the downstream surface was approximately tenfold that of the upstream roughness and, unlike all previous studies, attention was concentrated on the roughness sublayer region beneath the inertial (log-law) region. The experiments were conducted at a boundary layer Reynolds number of about 6 × 10⁴ (based on layer thickness andfree-stream velocity) and around a longitudinal location where the (downstream) roughness length, z_o2, was about 1% of the boundary-layer thickness atthe roughness change point.The thickness of the roughness sublayer was found for the two roughness. It was observed that the vertical profiles of mean velocity and turbulence characteristics started to show similarity after about 160z₀₂ downstream of the roughness change. The presence of a shear stress overshoot is shown to depend strongly on the precise location (with respect to the roughness elements) at which the measurements are made and the thickness of the equilibrium layer is shown to be very sensitive to the way it is defined. It is demonstrated that the growing equilibrium layer has first to encompass the roughness sublayer before any thickness of inertial sublayer can be developed. It follows that, in somepractical cases, like flows across some urban environments, the latter(log-law) region may never exist at all.     

Wind-tunnel Modelling of Dispersion from a Scalar Area Source in Urban-Like Roughness

A wind-tunnel study was conducted to investigate ventilation of scalars from urban-like geometries at neighbourhood scale by exploring two different geometries a uniform height roughness and a non-uniform height roughness, both with an equal plan and frontal density of λ _p = λ _f = 25%. In both configurations a sub-unit of the idealized urban surface was coated with a thin layer of naphthalene to represent area sources. The naphthalene sublimation method was used to measure directly total area-averaged transport of scalars out of the complex geometries. At the same time, naphthalene vapour concentrations controlled by the turbulent fluxes were detected using a fast Flame Ionisation Detection (FID) technique. This paper describes the novel use of a naphthalene coated surface as an area source in dispersion studies. Particular emphasis was also given to testing whether the concentration measurements were independent of Reynolds number. For low wind speeds, transfer from the naphthalene surface is determined by a combination of forced and natural convection. Compared with a propane point source release, a 25% higher free stream velocity was needed for the naphthalene area source to yield Reynolds-number-independent concentration fields. Ventilation transfer coefficients w _T /U derived from the naphthalene sublimation method showed that, whilst there was enhanced vertical momentum exchange due to obstacle height variability, advection was reduced and dispersion from the source area was not enhanced. Thus, the height variability of a canopy is an important parameter when generalising urban dispersion. Fine resolution concentration measurements in the canopy showed the effect of height variability on dispersion at street scale. Rapid vertical transport in the wake of individual high-rise obstacles was found to generate elevated point-like sources. A Gaussian plume model was used to analyse differences in the downstream plumes. Intensified lateral and vertical plume spread and plume dilution with height was found for the non-uniform height roughness.