Large-eddy simulation of turbulent flow and dispersion over a complex urban street canyon

Turbulent flow and dispersion characteristics over a complex urban street canyon are investigated by large-eddy simulation using a modified version of the Fire Dynamics Simulator. Two kinds of subgrid scale (SGS) models, the constant coefficient Smagorinsky model and the Vreman model, are assessed. Turbulent statistics, particularly turbulent stresses and wake patterns, are compared between the two SGS models for three different wind directions. We found that while the role of the SGS model is small on average, the local or instantaneous contribution to total stress near the surface or edge of the buildings is not negligible. By yielding a smaller eddy viscosity near solid surfaces, the Vreman model appears to be more appropriate for the simulation of a flow in a complex urban street canyon. Depending on wind direction, wind fields, turbulence statistics, and dispersion patterns show very different characteristics. Particularly, tall buildings near the street canyon predominantly generate turbulence, leading to homogenization of the mean flow inside the street canyon. Furthermore, the release position of pollutants sensitively determines subsequent dispersion characteristics.     

On the heating environment in street canyon

This study investigates the impact of building aspect ratio (building-height-to-street-canyon-width-ratio), wind speed and surface and air-temperature difference (Δθ_s−a) on the heating environment within street canyon. The Reynolds-averaged Navier-Stokes (RANS) and energy transport equations were solved with Renormalization group (RNG) theory version of k-e{\varepsilon} turbulence model. The validation process demonstrated that the model could be trusted for simulating air-temperature and velocity trends. The temperature and velocity patterns were discussed in idealized street canyons of different aspect ratios (0.5–2.0) with varying ambient wind speeds (0.5–1.5 m/s) and Δθ_s−a (2–8 K). Results show that air-temperatures are directly proportional to bulk Richardson number (R _b ) for all but ground heating situation. Conversely, air-temperatures increase significantly across the street canyon with a decrease in ambient wind speed; however, the impact of Δθ_s−a was negligible. Clearly, ambient wind speed decreases significantly as it passes over higher AR street canyons. Notably, air-temperatures were the highest when the windward wall was heated and the least during ground heating. Conversely, air-temperatures were lower along the windward side but higher within the street canyon when the windward wall was heated.     

Numerical simulations of the effect of building configurations and wind direction on fine particulate matters dispersion in a street canyon

To explore the effect of traffic emissions on air quality within street canyon, the wind flow and pollutant dispersion distribution in urban street canyons of different H/W, building gap and wind direction are studied and discussed by 3D computational fluid dynamics simulations. The largest PM_2.5 concentrations are 46.4, 37.5, 28.4 µg/m³ when x = ? 88, ? 19.3, ? 19.3 m in 1.5 m above the ground level and the ratio of H/W is 1:1, 1:2 and 2:1, respectively. The flow around the top of the building and clearance flow between the buildings in street canyon influence by different H/W, which affected the diffusion of fine particulate matters. The largest PM_2.5 concentrations are 88.1, 31.6 and 33.7 µg/m³ when x = 148.0, ? 92.3 and ? 186.7 m above the ground level of 1.5 m height and the building gap of 0, 20 and 40%, respectively. The air flows are cut by the clearance in the street canyons, and present the segmental characteristics. The largest PM_2.5 concentrations are 10.6, 11.2 and 16.0 µg/m³ when x = 165.3 m, x = 58.0 and 1.5 m above the ground level of 1.5 m height and wind direction of the parallel to the street, perpendicular to the street and southwest, respectively. Modelled PM_2.5 concentrations are basic agreement with measured PM_2.5 concentrations for southwest wind direction. These results can help analyze the difussion of PM_2.5 concentration in street canyons and urban planning.

     

A semi-coupled model for determining the distribution of two-dimensional coherent structures in open channel flow

Environmental Fluid Mechanics - Mixing layers and associated large-scale turbulent structures are some of the most common features of turbulent shallow water flow. Due to their...     

Evaluation of the QUIC-URB fast response urban wind model for a cubical building array and wide building street canyon

This paper describes the QUIC-URB fast response urban wind modeling tool and evaluates it against wind tunnel data for a 7 × 11 cubical building array and wide building street canyon. QUIC-URB is based on the Röckle diagnostic wind modeling strategy that rapidly produces spatially resolved wind fields in urban areas and can be used to drive urban dispersion models. Röckle-type models do not solve transport equations for momentum or energy; rather, they rely heavily on empirical parameterizations and mass conservation. In the model-experiment comparisons, we test two empirical building flow parameterizations within the QUIC-URB model: our implementation of the standard Röckle (SR) algorithms and a set of modified Röckle (MR) algorithms. The MR model attempts to build on the strengths of the SR model and introduces additional physically based, but simple parameterizations that significantly improve the results in most regions of the flow for both test cases. The MR model produces vortices in front of buildings, on rooftops and within street canyons that have velocities that compare much more favorably to the experimental results. We expect that these improvements in the wind field will result in improved dispersion calculations in built environments.     

Sensitivity and robustness of Lagrangian coherent structures in coastal water systems

Environmental Fluid Mechanics - In coastal water systems, horizontal chaotic dispersion plays a significant role in the distribution and fate of pollutants. Lagrangian Coherent Structures (LCSs)...     

Spatial evolution of coherent motions in finite-length vegetation patch flow

A number of experimental studies on submerged canopy flows have focused on fully-developed flow and turbulent characteristics. In many natural rivers, however, aquatic vegetation occurs in patches of finite length. In such vegetated flows, the shear layer is not formed at the upstream edge of the vegetation patch and coherent motions develop downstream. Therefore, more work is neededz to reveal the development process for large-scale coherent structures within vegetation patches. For this work, we considered the effect of a limited length vegetation patch. Turbulence measurements were intensively conducted in open-channel flows with submerged vegetation using Particle Image Velocimetry (PIV). To examine the transition from boundary-layer flow upstream of the vegetation patch to a mixing-layer-type flow within the patch, velocity profiles were measured at 33 positions in a longitudinal direction. A phenomenological model for the development process in the vegetation flow was developed. The model decomposed the entire flow region into four zones. The four zones are the following: (i) the smooth bed zone, (ii) the diverging flow zone, (iii) the developing zone and (iv) the fully-developed zone. The PIV data also confirmed the efficiency of the mixing-layer analogy and provided insight into the spatial evolution of coherent motions.     

Three-dimensional analysis of coherent turbulent flow structure around a single circular bridge pier

The coherent turbulent flow around a single circular bridge pier and its effects on the bed scouring pattern is investigated in this study. The coherent turbulent flow and associated shear stresses play a major role in sediment entrainment from the bed particularly around a bridge pier where complex vortex structures exist. The conventional two-dimensional quadrant analysis of the bursting process is unable to define sediment entrainment, particularly where fully three-dimensional flow structures exist. In this paper, three-dimensional octant analysis was used to improve understanding of the role of bursting events in the process of particle entrainment. In this study, the three-dimensional velocity of flow was measured at 102 points near the bed of an open channel using an Acoustic Doppler Velocity meter (Micro-ADV). The pattern of bed scouring was measured during the experiment. The velocity data were analysed using the Markov process to investigate the sequential occurrence of bursting events and to determine the transition probability of the bursting events. The results showed that external sweep and internal ejection events were an effective mechanism for sediment entrainment around a single circular bridge pier. The results are useful in understanding scour patterns around bridge piers.     

Fourier analysis of the roll-up and merging of coherent structures in shallow mixing layers

The current study investigates the role of nonlinearity in the development of two-dimensional coherent structures (2DCS) in shallow mixing layers. A nonlinear numerical model based on the depth-averaged shallow water equations is used to investigate temporal shallow mixing layers, where the mapping from temporal to spatial results is made using the velocity at the center of the mixing layer. The flow is periodic in the stream-wise direction and the transmissive boundary conditions are used in the cross-stream boundaries to prevent reflections. The numerical results are examined with the aid of Fourier decomposition. Results show that the previous success in applying local linear theory to shallow mixing layers does not imply that the flow is truly linear. Linear stability theory is confirmed to be only valid within a short distance from the inflow boundary. Downstream of this linear region, nonlinearity becomes important for the roll-up and merging of 2DCS. While the energy required for the merging of 2DCS is still largely provided by the velocity shear, the merging mechanism is one where nonlinear mode interaction changes the velocity field of the subharmonic mode and the gradient of the along-stream velocity profile which, in turn, changes the magnitude of the energy production of the subharmonic mode by the velocity shear implicitly. The nonlinear mode interaction is associated with energy up-scaling and is consistent with the inverse energy cascade which is expected to occur in shallow shear flows. Current results also show that such implicit nonlinear interaction is sensitive to the phase angle difference between the most unstable mode and its subharmonic. The bed friction effect on the 2DCS is relatively small initially and grows in tandem with the size of the 2DCS. The bed friction also causes a decrease in the velocity gradient as the flow develops downstream. The transition from unstable to stable flow occurs when the bed friction balances the energy production. Beyond this point, the bed friction is more dominant and the 2DCS are progressively damped and eventually get annihilated. The energy production by the velocity shear plays an important role from the upstream end all the way to the point of transition to stable flow. The fact that linear stability theory is valid only for a short distance from the inflow boundary suggests that some elements of nonlinearity is incorporated in the mean velocity profile in experiments by the averaging process. The implicit nature of nonlinear interaction in shallow mixing layers and the sensitivity of the nonlinear interaction to phase angle difference between the most unstable mode and its subharmonic allows local linear theory to be successful in reproducing features of the instability such as the dominant mode of the 2DCS and its amplitude.     

Significance of small-scale spatio-temporal heterogeneity in phytoplankton abundance for energy flow in Mytilus edulis

The effect of small-scale heterogeneity of phytoplankton abundance on mussel (Mytilus edulis L.) energy flow was studied by comparing filtration rate, consumption and assimilation efficiency in two groups of mussels fed under constant (control) and fluctuating (24 min) seston regimes of the same average. The experiment was carried out twice. Mussels appeared to regulate their filtration rate in response to varying seston concentration in one experiment, but not in the other. This behavioral difference may have been caused by differences in mean food level. Mean consumption was not different in the two feeding regimes, but assimilation efficiency was probably lower under fluctuating seston concentration. This suggests that small-scale heterogeneity in phytoplankton abundance would not affect the total transfer of energy between plankton and benthos. However, since a higher proportion of consumed energy was probably lost as faeces and/or pseudofaeces under the fluctuating seston regime, it may be expected that growth rate of mussels is depressed by small-scale heterogeneity in phytoplankton availability.Contribution to the programme GIROQ (Groupe interuniversitaire de recherches océanographiques du Québec)     

The role of coherent turbulent structures in explaining scalar dissimilarity within the canopy sublayer

Scalar similarity is widely assumed in models and interpretation of micro-meteorological measurements. However, in the air space within and just above the canopy (the so-called canopy sublayer, CSL) scalar similarity is generally violated. The scalar dissimilarity has been mainly attributed to differences in the distribution of sources and sinks throughout the canopy. Since large-scale coherent structures in the CSL (e.g. double roller and sweep/ejection) arise from the instabilities generated by the interaction between the mean flow and the canopy, they may encode key dynamical features about the production term responsible for the source–sink dissimilarity of scalars. Therefore, it is reasonable to assume that the geometric attributes of coherent structures are tightly coupled to the onset and the vertical extent of scalar dissimilarity within the CSL. Large-eddy simulation (LES) runs were used to investigate the role of coherent structures in explaining scalar dissimilarity among three scalars (potential air temperature, water vapour and $\text{ CO }_2$ concentration) within the CSL under near-neutral conditions for horizontally uniform but vertically varying vegetation leaf area density. It was shown that coherent structures, when identified from the first mode of a novel proper orthogonal decomposition (POD) approach, were able to capture some features of the scalar dissimilarity in the original LES field. This skill was quantified by calculating scalar–scalar correlation coefficients and turbulent Schmidt numbers of the original field and the coherent structures, respectively. However, coherent structures tend to magnify the magnitude of scalar–scalar correlation, particularly in cases where this correlation is already strong. The ability of coherent structures to describe more complex features such as the scalar sweep-ejection cycle was also explored. It was shown that the first mode of the POD does not capture the relative importance of sweeps to ejections in the original LES field. However, the superposition of few secondary coherent structures, derived from higher order POD modes, largely diminish the discrepancies between the original field and the POD expansion.     

LES validation of urban flow,part II: eddy statistics and flow structures

Time-dependent three-dimensional numerical simulations such as large-eddy simulation (LES) play an important role in fundamental research and practical applications in meteorology and wind engineering. Whether these simulations provide a sufficiently accurate picture of the time-dependent structure of the flow, however, is often not determined in enough detail. We propose an application-specific validation procedure for LES that focuses on the time dependent nature of mechanically induced shear-layer turbulence to derive information about strengths and limitations of the model. The validation procedure is tested for LES of turbulent flow in a complex city, for which reference data from wind-tunnel experiments are available. An initial comparison of mean flow statistics and frequency distributions was presented in part I. Part II focuses on comparing eddy statistics and flow structures. Analyses of integral time scales and auto-spectral energy densities show that the tested LES reproduces the temporal characteristics of energy-dominant and flux-carrying eddies accurately. Quadrant analysis of the vertical turbulent momentum flux reveals strong similarities between instantaneous ejection-sweep patterns in the LES and the laboratory flow, also showing comparable occurrence statistics of rare but strong flux events. A further comparison of wavelet-coefficient frequency distributions and associated high-order statistics reveals a strong agreement of location-dependent intermittency patterns induced by resolved eddies in the energy-production range. The validation concept enables wide-ranging conclusions to be drawn about the skill of turbulence-resolving simulations than the traditional approach of comparing only mean flow and turbulence statistics. Based on the accuracy levels determined, it can be stated that the tested LES is sufficiently accurate for its purpose of generating realistic urban wind fields that can be used to drive simpler dispersion models.     

Evaluation of fast atmospheric dispersion models in a regular street network

Environmental Fluid Mechanics - The need to balance computational speed and simulation accuracy is a key challenge in designing atmospheric dispersion models that can be used in scenarios where...     

Resident mysids: community structure,abundance and small-scale distributions in a coral reef lagoon

Seasonal and diel variations in community structure and abundance of coral-reef lagoon mysids were examined at Davies Reef in the central region of the Great Barrier Reef (GBR) between June 1980 and May 1981. Twenty-five mysid species belonging to three subfamilies of the family Mysidae were captured during the study, with six new records for the GBR. The epibenthic mysid community differed from that in the overlying water, was faunistically uniform, but formed characteristic seasonal and diel groupings. The dominant epibenthic species were Erythrops sp., Anisomysis pelewensis, Doxomysis littoralis, A. laticauda, Prionomysis stenolepis, A. lamellicauda, and A. australis, five of which formed schools. Total mysid abundances ranged between 110 and 790 m^-3 with peak abundance in October. Schooling species occurred at local densities of up to 500 000 m^-3. Mysids were absent from shallow and midwater depths during the day, but were distributed throughout all depths at night with peak abundances in mid-water and deep layers. The dominant species in the water column at night were Pseudanchialina inermis, A. laticauda and Gastrosaccus indicus, in descending order of abundance. Lagoonal mysids contribute little to the food of sessile reef planktivores, as all but three species remain concentrated near or on the lagoon floor both day and night. The contribution of resident lagoon mysids to reef trophodynamics is probably through remineralization of lagoon detritus. Given the vast reef areas comprised of sandy lagoons, the large populations and relatively large size of lagoon mysids, this trophodynamic role may be of considerable importance.A.I.M.S. Contribution No. 477     

Rare earth elements in street dust and associated health risk in a municipal industrial base of central China

The content levels, distribution characteristics, and health risks associated with 15 rare earth elements (REEs) in urban street dust from an industrial city, Zhuzhou, in central China were investigated. The total REE content (∑REE) ranged from 66.1 to 237.4 mg kg^?1, with an average of 115.9 mg kg^?1, which is lower than that of Chinese background soil and Yangtze river sediment. Average content of the individual REE in street dust decreased in the order Ce > La > Nd > Y > Pr > Sm > Gd > Dy > Er > Yb > Eu > Ho > Tb > Tm > Lu. The chondrite-normalized REE pattern indicated light REE (LREE) enrichment, a relatively steep LREE trend, heavy REE (HREE) depletion, a flat HREE trend, a Eu-negative anomaly and a Ce-positive anomaly. Foremost heavy local soil and to less degree anthropogenic pollution are the main sources of REE present in street dust. Health risk associated with the exposure of REE in street dust was assessed based on the carcinogenic and non-carcinogenic effect and lifetime average daily dose. The obtained cancer and non-cancer risk values prompt for no augmented health hazard. However, children had greater health risks than that of adults.     

小城镇街尘污染特征及控制对策

以重庆市永川区小城镇为例进行调查研究,结果表明,工业园区道路街尘和城乡结合部道路街尘以200μm以上粒径段为主(65%—72%),而商业区和双竹镇街道以200μm以下粒径段为主(60%—72%);有机质和总磷、总氮、锌、铅、镉的最高浓度出现在商业区街尘(分别为103.04 mg.g-1和785.48、1672.04、439.72、133.30、8.44 mg·kg-1),工业园区道路街尘含有最高的铜、镍浓度(分别为241.52和92.30 mg·kg-1).小城镇街道街尘污染评价结果表明,商业区街尘和工业园区道路街尘污染风险较高,其中商业区街尘中的镍处于重度污染等级,工业园区道路街尘中的镉达到极度污染最高等级,同时,小城镇街尘控制对策分析指出,在"源"-"传输过程"-"汇"三个环节分别采取措施控制小城镇街尘污染是可行的.     

Health risks from PAHs and potentially toxic elements in street dust of a coal mining area in India

Environmental Geochemistry and Health - Polycyclic aromatic hydrocarbons (PAHs) and potentially toxic elements (PTEs) (Ba, Zn, Pb, Cu, Cr, Ni, As, Co) were determined in the road dusts of a coal...