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.     

Computational fluid dynamics assessment of damaging wind loads on the One Indiana Square tower

In this paper, we simulated damaging wind loads on the One Indiana Square tower in Indianapolis due to the storm of April 2nd 2006. We followed recommended practice guidelines for this urban wind modeling. First, a test case, Aerodynamics of Commonwealth Advisory Aeronautical Council (CAARC) building were modeled and simulated to compare with a publicly available experiment and other computational studies. Based on the modeling parameters in the CAARC study, then, as a clean building configuration, we modeled the One Indiana tower alone without surrounding buildings. Finally, the flow field around the tower including nearby downtown buildings were simulated. We used the Fluent flow analysis software tools. The domain was meshed using unstructured grids, the first boundary layer grid element being 10 cm (4 in.) and 15 cm (6 in.) in height from the tower and the ground for the CAARC building and the One Indiana tower, respectively. Two different wind directions of 260 \(^\circ \) and 280 \(^\circ \) at 137 km/h (85 mph) speed were considered to estimate wind loads on the One Indiana tower façades. Simulated pressure distributions on the tower and flow patterns over the downtown buildings were discussed to draw conclusions about the mechanism of extreme wind load that caused the damage. The simulations revealed that suction forces are almost twice higher hence more damaging at the corners of the west façades than straight wind. It was also seen in the simulation results that upstream building topology, specifically Chase, One America, and some low-rise towers, augmented the actual wind load unfavorably on the One Indiana Square tower. Although this study presents a specific case, the applicability of its findings are of more general interest. Similar wind events are common especially during storm seasons both in urban and suburban areas. In similar incidents, one can follow the same procedure to analyze their problems as certain modeling guidelines were followed in this study.     

Application of CFD modelling to air quality in Kuwait City

High-resolution computational fluid dynamics (CFD) simulations have been performed to assess the dispersion of air pollutants (CO₂) emanating from traffic in a busy street and in the vicinity of a complex configuration of buildings located in Salmiya, Kuwait City. New buildings are planned for this area, and the work here includes predictions for the dispersion of pollutants after the buildings’ completion. The CFD simulations are based on calculated CO₂ concentration levels for traffic counts taken on location in Salmiya with the existing configuration of buildings. As the computer code used in this work has been evaluated previously, it will be applied here to predict with confidence any potential air pollution problem areas on the addition of the new buildings. It was found for very light wind, that the proposed new buildings help reduce pollution in the vicinity of residential buildings within the configuration of buildings, but as the wind becomes moderate to strong, there was a tendency for the pollutant to get trapped in the residential area. Results are given for both exceptionally high ambient temperatures and very light wind, which are not often reported in studies found in the literature.     

An approach to simulate wind fields around an urban environment for wind energy application

This work proposes an approach to simulate wind flow fields around an urban environment with the aim of evaluating the potential impact of buildings on the general wind patterns and power production using the current generation of commercial wind turbines. The simulation process was performed with the aid of accessible computational tools that can potentially render the proposed procedure applicable in other cases of interest. The roughness of the urban environment was defined as the association of roughness map, topography, and an alternative process for obtaining the volumetry of buildings. A case study was conducted in a region located at the district of Boa Viagem (Recife-PE) for assessing the applicability of the approach. Scenarios were designed in order to simulate wind flow patterns and pre-identify sites that have suitable wind energy potential for electric power production by investigating the combination of wind speed magnitude and turbulence intensity. From the results obtained, it was possible to identify zones of potential wind sources that are not detected in classical wind atlas probably due to the influence of the built environment on local wind flow patterns.     

Evaluation of a fast response pressure solver for flow around an isolated cube

This work describes and evaluates a pressure solver that has been incorporated into a fast response three-dimensional building-resolving diagnostic wind modeling system. The solver computes the three-dimensional pressure field around buildings and on exterior walls in terms of a coefficient of pressure by solving a simplified pressure Poisson equation (that neglects turbulence stresses in the Navier-Stokes) for incompressible flow. The input to the solver is the three-dimensional mean wind field obtained from a fast response empirical-diagnostic urban wind model. The present study is an evaluation of the pressure solver using wind-tunnel data for flow normal to and at a 45° angle to an isolated cubical building. Results for the normal incident wind angle case indicate that the model satisfactorily reproduces the general spatial patterns and the magnitude of the pressure difference around much of the cube. Details of the flow field that are not satisfactorily predicted include the spatial distribution of pressure on the roof and the lower half of the front side of the building and the magnitude along the sidewalls where pressures are over predicted. The results for the 45° case show reasonable agreement between the model and experiments on the front and the back walls, but over predict pressures on the leading edge of the rooftop. Regions with poor pressure predictions appear to be a result of unsatisfactory mean wind modeling.     

Wave hydrodynamics around a multi-functional artificial reef at Leirosa

This paper describes an application of the Boussinesq-type COULWAVE model to study the wave hydrodynamics in the vicinity of a multi-functional artificial reef (MFAR). This reef is under investigation and consists of a supplementary protection solution for the Leirosa sand dune system located at South of Figueira da Foz, on the Portuguese West coast. Such installation near the coastline is expected to contribute to enhance the surfing conditions in the area, protect the sand dune system in the surroundings of Leirosa beach, and increase its environmental value. Numerical calculations with the COULWAVE model were performed for four test cases, considering two reef geometries (differing in the reef angle) and two incident wave conditions (storm condition and a common wave condition). Comparisons between the results obtained, in terms of wave heights and breaking line positions allow us to assess the influence of the reef on the hydrodynamics near the beach and around the reef. Moreover, the reef performance was analysed in terms of surfability and coastal protection. The surfability parameters (breaker height, Iribarren number and peel angle) were calculated for each test case using the numerical wave heights, wave directions and wave breaking positions. Comparisons of parameters allow characterizing the most appropriate configuration of the reef to improve the surfing conditions in the study area. A methodology based on numerical free surface elevations and horizontal velocity components was developed to calculate wave directions, since this is not a direct output of the COULWAVE model. Concerning coastal protection, analyses of the mean currents around the reef were used together with observations of the velocity cells near the shoreline as an indication of the sediment transport.     

居住区绿色空间的生态规划与设计——以珠海华发新城为例

居住区绿色空间是城市绿色空间的重要组成部分，它对提高居民的生活环境质量，增进居民的身心健康至关重要，具有重要的生态服务功能。从生态学的角度出发，提出了居住区绿色空间生态规划的8项原则，并以珠海华发新城为例，从滨河公园、住区外围防护绿色空间、建筑组团外围和庭院绿色空间、建筑架空底部绿色空间、屋顶和垂直绿化等方面进行了规划设计。     

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.     

Exposure assessment of particulate matter and blood chromium levels in people living near a cement plant

This study evaluates the effect of air pollution caused by cement plants on nearby residential areas and performs an exposure assessment of particulate matter (PM) and total Cr, Cr⁶⁺, Pb, and Al. Further, the blood Cr levels of residents exposed to PM released by cement plants are also assessed. Nine buildings (eight residential and one elementary school building) close to cement plants were selected for this study, which were located in Pyeongtaek port, in west of Gyeonggi Province, South Korea. A total of 51 suspended particulate samples were collected at a flow rate of 2.0 L/min. Total Cr was more widely detected in residents’ houses and elementary schools. PM levels were higher at distances of 4.1 and 4.8 km than those at closer distances of 2.7 km. This was due to the influence of wind direction. The estimated mean blood level of Cr for the study participants was 3.80 μg/L, which is higher than levels estimated by other studies on Cr blood levels. Therefore, cement plants could cause an increase in total Cr and blood Cr levels in residential areas, and more continuous monitoring is necessary to better understand their impacts.     

Experimental surface flow visualization and numerical investigation of flow structure around an oblong stockpile

Emission factors are largely used to quantify particle emissions from industrial open storage piles. These factors are based on the knowledge of velocity distribution and flow patterns over the stockpile surface which still requires further research. The aim of the present work is to investigate the airflow characteristics over a single typical oblong pile and in its near-ground surroundings for various wind flow directions. Wind tunnel experiments using an oil-film surface coating technique were carried out for near-wall flow visualization. Numerical simulation results, favorably compared to PIV measurements, were used to allow comparison analysis of flow features. For the stockpile oriented 90° to the wind main direction, typical topology of flow around wall-mounted obstacles were observed, notably a wake zone downstream the pile including two main counter-rotating vortices. Further analysis of numerical wall shear stress distribution and streamlines indicates that two complex three-dimensional vortical flow structures develop downstream the pile. For other incoming wind flow directions (30 and 60°), the flow characteristics over the storage pile greatly differ as a single helical main vortex develops from the pile’s crest. Corresponding high values of wall shear stress are noticed downstream the storage pile. For each configuration studied, downwash and upwash zones are induced by the vortical structures developed. This near-wall flow topology combined with areas of high friction levels may be linked to potential dust emission from the ground surface surrounding industrial stockpiles.     

广州城市近地风场特征研究

利用实测广州市近地风速、风向数据资料,分析广州市近地风场特征。结果表明：（1）城市近地风速微弱,风速范围在0.1~0.71 m.s-1,平均风速约0.25 m.s-1,处于静风状况。城市对盛行风风速的削减达90%左右;（2）在盛行风被极度削弱的同时,珠江江风促成了城区一些区域风速的增强。同时,白云山等北部山林的山风以及大型绿地对周围测点的风场均有影响,这些局地的热力效应促进了大气的流动,增强近地风速,改善近地风速微弱甚至无风的状况。广州大道北联白云山,南通珠江,道路开敞,受到白云山山风和珠江江风的影响,广州大道监测点风速相对较大。有些地方则成为＂风影区＂,监测点基本处于无风状态。随着城市的扩展,城市近地风速微弱,通风不良,风向变化复杂,这对于汽车日益增多,尾气污染激增的城市是极为不利的。因此,城市的河流、山林以及大型绿地对城市近地风场的影响值得关注,值得进一步研究如何利用其布局改善城市的无风状况。     

A CFD study of wind patterns over a desert dune and the effect on seed dispersion

In the Namib Desert seed distribution is greatly influenced by wind patterns. Existing literature regarding wind patterns over dunes focuses on two-dimensional simulations of flow over simplified dune structures. The three-dimensional geometries of the sand dunes suggests far more complex flow features exist, which are not captured by two-dimensional simulations. Computational fluid dynamics (CFD) was used to reproduce the three-dimensional near surface wind patterns around a dune with the aim to learn more about seed distribution. Field work included terrain mapping, wind speed, direction and temperature metering. The CFD results show the expected two-dimensional flow features of high pressure at the dune toe, low pressure at the crest and flow acceleration up windward slope. Also observed are some three-dimensional flow features such as a spiral vortex near the crest and transverse flow due to crest-line curvature of the dune. It was also observed how the wall shear stress differs due to the three-dimensional shape of the dune. The wall shear stress suggests that seed accumulation is more likely to occur behind trailing (down-wind) crest edges. Particle tracking showed how seeds tend to move over the dune crest and recirculate towards the crest on the lee-side. The study showed that adding the third dimension makes the simulations more complex, adds to computational requirements and increases simulation time but also provides vital flow information which is not possible with two-dimensional simulations.     

Stratification and Circulation in a Shallow Turbid Waterbody

Shallow waterbodies are often assumed to be well mixed in the vertical. However, when they are characterised by high turbidity levels, absoption of solar heating within a relatively thin surface layer can produce thermal stratification. Results from an intensive monitoring program have been combined with three-dimensional circulation modelling to examine the diurnal stratification cycles in a small turbid waterbody. The waterbody, known as Rushy Billabong, is located in southeastern Australia and its high turbidity coupled with forcing by wind and solar radiation resulted in regular diurnal cycles of stratification and overturning. Under conditions of light wind and high solar radiation, the model results were generally consistent with the observed temperature field. However, under stronger winds, preferential cooling and sinking of shallow water around the edge of the lake began to contribute significantly to the interior stratification. Model estimates then became more sensitive to the detailed bathymetry and the choice of turbulence parameterisation. The level of stratification is also shown to influence the circulation in the billabong by trapping the wind-driven flow near the surface. Insights provided by the observations and modelling may have broader implications for the management of small turbid systems such as settling ponds, aquaculture ponds, and some natural wetlands.     

Aeolian erosion of storage piles yards: contribution of the surrounding areas

Dust emissions from stockpiles surfaces are often estimated applying mathematical models such as the widely used model proposed by the USEPA. It employs specific emission factors, which are based on the fluid flow patterns over the near surface. But, some of the emitted dust particles settle downstream the pile and can usually be re-emitted which creates a secondary source. The emission from the ground surface around a pile is actually not accounted for by the USEPA model but the method, based on the wind exposure and a reconstruction from different sources defined by the same wind exposure, is relevant. This work aims to quantify the contribution of dust re-emission from the areas surrounding the piles in the total emission of an open storage yard. Three angles of incidence of the incoming wind flow are investigated ( $30^{\circ }, 60^{\circ }$ and $90^{\circ }$ ). Results of friction velocity from numerical modelling of fluid dynamics were used in the USEPA model to determine dust emission. It was found that as the wind velocity increases, the contribution of particles re-emission from the ground area around the pile in the total emission also increases. The dust emission from the pile surface is higher for piles oriented $30^{\circ }$ to the wind direction. On the other hand, considering the ground area around the pile, the $60^{\circ }$ configuration is responsible for higher emission rates (up to 67 %). The global emissions assumed a minimum value for the piles oriented perpendicular to the wind direction for all wind velocity investigated.     

The performance of RAMS in representing the convective boundary layer structure in a very steep valley

Data from a comprehensive field study in the Riviera Valley of Southern Switzerland are used to investigate convective boundary layer structure in a steep valley and to evaluate wind and temperature fields, convective boundary layer height, and surface sensible heat fluxes as predicted by the mesoscale model RAMS. Current parameterizations of surface and boundary layer processes in RAMS, as well as in other mesoscale models, are based on scaling laws strictly valid only for flat topography and uniform land cover. Model evaluation is required to investigate whether this limits the applicability of RAMS in steep, inhomogeneous terrain. One clear-sky day with light synoptic winds is selected from the field study. Observed temperature structure across and along the valley is nearly homogeneous while wind structure is complex with a wind speed maximum on one side of the valley. Upvalley flows are not purely thermally driven and mechanical effects near the valley entrance also affect the wind structure. RAMS captured many of the observed boundary layer characteristics within the steep valley. The wind field, temperature structure, and convective boundary layer height in the valley are qualitatively simulated by RAMS, but the horizontal temperature structure across and along the valley is less homogeneous in the model than in the observations. The model reproduced the observed net radiation, except around sunset and sunrise when RAMS does not take into account the shadows cast by the surrounding topography. The observed sensible heat fluxes fall within the range of simulated values at grid points surrounding the measurement sites. Some of the scatter between observed and simulated turbulent sensible heat fluxes are due to sub-grid scale effects related to local topography.     

Scalar transport from point sources in the flow around a finite-height cylinder

This paper presents a large eddy simulation of mass transfer in the flow around a surface-mounted finite-height circular cylinder. The study was carried out for a cylinder with height-to-diameter ratio of 2.5 and a Reynolds number based on the cylinder diameter of 44000. The approach flow boundary layer had a thickness of about 10% of the cylinder height. A tracer was released at various levels upstream of the cylinder. The effect of the release position in the subsequent spreading and dilution of the plumes is analyzed. It is found that a tracer released at the top or at mid-height is entrained into the recirculation zone behind the cylinder, and therefore presents similar plume evolution in the far wake in both cases. If the tracer is released at around three-quarters of the height of the cylinder, it is not significantly entrained by the main recirculation region, leading to smaller rates of spreading of the plume. Finally, if the tracer is released near the floor, the plume is entrained by the horseshoe vortex that wraps around the cylinder, leading to a large lateral spreading of the plume, remaining always near the floor.     

On thermally forced flows in urban street canyons

During sunny days with periods of low synoptic wind, buoyancy forces can play a critical role on the air flow, and thus on the dispersion of pollutants in the built urban environments. Earlier studies provide evidence that when a surface inside an urban street canyon is at a higher temperature than that of local ambient air, buoyancy forces can modify the mechanically-induced circulation within the canyons (i.e., gaps between buildings). The aspect ratio of the urban canyon is a critical factor in the manifestation of the buoyancy parameter. In this paper, computational fluid dynamics simulations are performed on urban street canyons with six different aspect ratios, focusing on the special case where the leeward wall is at a greater temperature than local ambient air. A non-dimensional measure of the influence of buoyancy is used to predict demarcations between the flow regimes. Simulations are performed under a range of buoyancy conditions, including beyond those of previous studies. Observations from a field experiment and a wind tunnel experiment are used to validate the results.     

Hydrodynamics and turbulence in emergent and sparsely vegetated open channel flow

This present study reports the results of an experimental study characterizing thorough variation of turbulent hydrodynamics and flow distribution in emergent and sparsely vegetated open channel flow. An emergent and rigid sparse vegetation patch with regular spacing between stems along the flow and transverse directions was fixed in the central region of the cross-section of open channel. Experiments were conducted in subcritical flow conditions and velocity measurements were obtained with an acoustic Doppler Velocimetry system. Large variations of the turbulence intensities, Reynolds shear stress, turbulent kinetic energy and vortical motions are found in and around the vegetation patch. At any cross-section through the interior of the vegetation patch, streamwise velocity decreases with increase in streamwise length and the velocity profiles converge from the log-law to a linear profile with increasing slope. Time-averaged lateral and vertical velocities inside the vegetation patch increase with increasing streamwise distance and converge from negative values to positive values. Turbulence intensities interior of the sparse vegetation patch are more than those of without the vegetation patch. Similar to the trend of streamwise velocity profiles inside the vegetation, turbulence intensities and longitudinal-normal Reynolds shear stress profile decreases with streamwise direction. In the interior of the vegetation patch and downstream of the trailing edge, turbulent kinetic energy profiles are exhibiting irregular fluctuations and the maximum values are occurring in the outer layer. Analysis of flow distribution confirms sparse vegetation patch is inducing a serpentine flow pattern in its vicinity. At the leading edge, flow is rushing towards the right hand sidewall, and at the trailing edge, flow is turning to the left hand sidewall. In between the leading and trailing edges, the streamlines are following a zig-zag fashion at varied degree along the streamwise and lateral directions. Immediate upstream of the leading edge and in the interior of the vegetation patch, vortex motion is clearly visible and the vortices are stretched along the width of the channel with streamwise direction.     

Steady Shallow-Water Current and Solute Transport around a Semi-Conical Headland

Laboratory experiments have been carried out to investigate the effects of a sloping wall headland on the flow characteristics and the associated concentration distributions from a point source around the headland. A semi-conical headland with a slope of 1:2 was set up in a flow basin, 4.8 m long and 3.8 m wide. In this paper, the experimental results of a steady shallow-water current are reported. Three dimensional flow velocities in the basin were measured using Sontek-ADV instrument. The dye concentration levels in the basin were measured by two fluorometers. The experimental results showed a large-scale re-circulation region behind the semi-conical headland. The peak turbulence energy, at about 53% of the local kinetic flow energy, coincides with the region of high velocity gradient. Significant vertical flows were observed around the area near the downhill slope of the headland, with a maximum ratio of vertical to horizontal velocities being about 22%. Such relatively significant vertical scouring velocities, coupled with strong turbulence energy and high horizontal velocity gradients in the same region, could cause severe bed erosion. The experimental results have also been compared with the predicted results of a depth-averaged numerical model. The predicted eddy structure and the concentration distribution in the re-circulation area were found to compare favourably with the experimental results. However, the discrepancies in the flow velocities and the concentration levels near the headland were apparent. It was observed that the dye concentration continued to spread in the cross-stream direction after passing the headland, whereas only a limited extent of the lateral spreading was predicted by the numerical model further downstream of the headland.