Processes 09 01754
Processes 09 01754
Processes 09 01754
Article
Study of Heavy Gas Pollutants’ Dispersion in Street
Canyon Terrain
Xiaoye Nie 1 , Yuyan Fan 1 , Hong Gao 1,2, * and Qiang Lin 1,3, *
1 School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, China;
[email protected] (X.N.); [email protected] (Y.F.)
2 Tianjin Key Laboratory of Chemical Process Safety and Equipment Technology, Tianjin 300072, China
3 Zhejiang Institute of Tianjin University, Ningbo 315201, China
* Correspondence: [email protected] (H.G.); [email protected] (Q.L.)
Abstract: This study focused on heavy gas dispersion under the terrain conditions of street canyons.
The effects of street aspect ratio and height ratio were investigated, and the influence of environmental
wind speed in the typical ideal street canyon terrain was explored. The results indicated that the
surrounding flow field distributions in street terrains were dominated by higher buildings. In
addition, when the building height was held constant, the flow field was affected by the joint
influence of the two isolated buildings. The interception effect of the street canyon on upstream
pollutants declined with the decrease in the street canyon’s aspect ratio. In addition, when the height
ratios were different, a large quantity of upstream pollutants accumulated on the windward side
of higher buildings. The relative concentration per unit area inside the canyon was affected by the
air circulation inside and outside the canyon and the size of the dispersion space. The increase in
the environmental wind speed promotes the entry of pollutants into the street while aggravating
the overall dispersion of the pollutants. Therefore, the emergence of the most unsafe wind speeds
caused most of the pollutants to gather in the street canyons.
Citation: Nie, X.; Fan, Y.; Gao, H.;
Lin, Q. Study of Heavy Gas Keywords: wind tunnel; CFD simulation; pollutant dispersion; heavy gas; safety
Pollutants’ Dispersion in Street
Canyon Terrain. Processes 2021, 9,
1754. https://doi.org/10.3390/
pr9101754 1. Introduction
The fluorine chemical industry has expanded with the development of the military
Academic Editor: Li Xi
industry since the 1930s, and fluorine chemical products are now widely used in the
Received: 22 August 2021
military industry, aerospace, electronics, light industry, and agriculture [1]. With the
Accepted: 25 September 2021 increase in demand for fluorine chemical products, the production process and scale of
Published: 30 September 2021 fluorine chemicals have also been continuously developed and expanded, and the related
safety issues have also attracted widespread attention.
Publisher’s Note: MDPI stays neutral The leakage of toxic gases is the focus of attention in numerous fields, including that
with regard to jurisdictional claims in of fluorine chemicals. Chlorine and hydrogen fluoride, as raw materials and products in
published maps and institutional affil- the fluorine chemical process, exist in multiple links of the production process. Due to
iations. the extremely high toxicity of the two gases, leakage has the potential to cause significant
harm [2–4]. Most toxic gases readily combine and react with water molecules in the air to
form heavy gases, due to their excessive molecular mass or relatively active physical and
chemical properties. Therefore, the study of heavy gas dispersion is a popular research
Copyright: © 2021 by the authors.
topic [5–7]. In addition, researchers tend to use safe and reliable gases in experiments
Licensee MDPI, Basel, Switzerland. and simulations of the dispersion of heavy gases. Tan et al. [8,9] conducted a series
This article is an open access article of wind tunnel experiments using CO2 gas to simulate the cross block, and observed
distributed under the terms and the accumulation and collapse of heavy gas inside the block. Xin et al. [10] constructed a
conditions of the Creative Commons numerical simulation model for chlorine dispersion, and analyzed the individual and social
Attribution (CC BY) license (https:// risks under various leakage conditions. Mao et al. [11] built an experimental system to
creativecommons.org/licenses/by/ explore the changes in the concentration field of indoor heavy gas leakage and dispersion,
4.0/).
and influencing factors, and found that the gas concentration fluctuations at different
locations are different, which makes it difficult to estimate the risk of pollutants.
The street canyon terrain is a kind of architectural layout widely used in production
and living areas. The complex wind field around it can easily intercept nearby pollutants
and pose a threat to internal personnel. Marucci et al. [12] studied the influence of heating
different walls in a two-dimensional street on the surrounding flow field and the dispersion
of internal pollutants. The results indicated that the heating of the windward walls had
a greater impact on various parameters. Yang et al. [13] used a fire dynamics simulator
to simulate the dispersion of chlorine horizontally leaking outside the two-dimensional
street canyon, and used the LES method to restore the vortex shape and pollutant flow
state in several street canyon configurations. Xu et al. [14] proposed a computational fluid
dynamics (CFD) method for simulating gas dispersion in urban environments, explored the
emission and dispersion process of line-source heavy gas pollutants in extremely long street
canyons, and compared concentrations of pollutants near the canyon walls. Most research
on street canyons has adopted the infinite two-dimensional canyon model, which can
represent some topographical features of street canyons, but ignores the horizontal airflow
exchange between street canyons and the external airflow. In many studies on the leakage
of pollutants from cities or factory buildings [15–17], short street canyons are common, and
pollutants spread significantly along street canyons. Therefore, investigating the terrain of
three-dimensional short street canyons is crucial. Many scholars have undertaken efforts
to improve the accuracy of numerical simulation. Ferreira et al. [18] proposed a novel jet
model, DESQr (Diameter of Equivalent Simulation for Quick Run), which was implemented
in the framework of FDS (Fire Dynamics Simulator) and able to capture turbulent eddies
in a simulation. Vianna et al. [19] modified the classic porosity distributed resistance
formulation, and suggested a new model that represented the part of the filtered geometry
from the original model through the porosity values of the unstructured tetrahedral mesh.
Ferreira and Tatiele et al. [20] combined the Gilbert–Johnson–Keerthi distance algorithm
and the finite volume to obtain refined porosity values and processed complex geometries
within a feasible calculation time.
In the present work, the distributions of the flow field around street canyon terrains
and the effects on the dispersion of heavy gas pollutants were investigated. This paper
is divided into five sections. Section 2 describes the wind tunnel experimental devices,
concentration display system, and three experimental programs. Section 3 introduces
the conditions of three-dimensional numerical simulation. Wind tunnel experiments
and numerical simulations were carried out on different street canyons. The impacts of
different aspect ratios and height ratios on streets are discussed in Section 4. In addition,
the existence of the most unsafe wind speed under the typical ideal street canyon terrain
was explored. In this study, the dispersion of pollutants was restored by adopting the laser
particle tracking technology, which realized the overall display of the concentration rather
than providing only the results of individual points and assisted with CFD simulation to
explore the manner in which pollutants enter the street canyon.
the entrance of the experimental section, of which the cross-section was 0.35 × 0.45 m. The
size ratio of the wind tunnel experimental model to the actual street was 1:400. The front
end of the experimental section was equipped with wedges and rough elements to simulate
the real urban wind profile, and the IFA300 constant temperature hot wire anemometer
was used to measure the wind speed. The final measurement of the set ambient wind
conformed to the power law with the reference point height, z0 , being 0.1 m, and the
characteristic wind speed, u0 , being 1.3731 m/s:
α
z
u = u0 (1)
z0
where α is the wind profile index, and its fitted value was 0.243. The wind profile index is
related to atmospheric stability and terrain conditions. “Technical Methods for Establishing
Local Air Pollutant Emission Standards (GB/T 3840-1991)” [21] provides a series of index
reference values under different atmospheric stability and terrain conditions, as shown in
Table 1. The fitting result accorded with the distribution law of the urban boundary layer,
and the atmospheric stability was between C and D. Atmospheric stability varies from A
to F from strong to stable, and when the stability is between C and D, it is close to neutral.
The experimental wind speed and fitting curve are shown in Figure 2. In this work, the
density Froude number, Fr, and the building Reynolds number, ReH, were selected as the
similarity criterion:
u0 2
Fr = 0 (2)
g0 L
uH H
ReH = ( ) (3)
ν
where:
u0 represents the characteristic velocity,
g0 ’ represents the relative acceleration (=g(ρs − ρa )/ρa , ρs is the density of the leaking
gas, ρa is the air density),
L represents the characteristic length,
ν represents the air kinematic viscosity,
uH represents the wind speed at the height of the building,
H represents the building height.
A B C D E, F
City 0.10 0.15 0.20 0.25 0.30
Rural 0.07 0.07 0.10 0.15 0.25
Processes 2021, 9, 1754 4 of 20
During the experiment, the Froude numbers remained consistent and the Reynolds
number was independent. The building Reynolds number was 2112.16, which was greater
than its critical value of 2.1 × 103 [22,23].
C I
C= = (4)
C0 I0
where:
Figure 3. (a) Grayscale picture of instantaneous concentration field, (b) grayscale picture of average
concentration field, (c) comparison of numerical simulations and experimental results.
Processes 2021, 9, 1754 5 of 20
3. Numerical Simulation
A simulation model was established according to the size of the wind tunnel experi-
ments. The typical ideal street canyon configuration with street aspect ratio Rw and height
ratio RH of 1.0 is shown in Figure 4a. The building height was H, and the upstream building
model was located at the midline 2 H downstream of the leakage source. The width of the
street formed by the two buildings was H. The ICEM non-structural tetrahedral grid was
used to mesh the model. In addition to refining the grids of ground, buildings, and the
vicinity of the leak source, the independence of the mesh was also verified. Three kinds of
grids were set up for comparison, and the number of grid cells was 1,603,409, 2,528,821,
and 6,740,591, as shown in Figure 4b. Finally, the grid with 2,528,821 cells was selected.
Under the same conditions, the number of grid cells in the model under each configuration
condition in the simulation process exceeded 2.5 million. The numerical simulations were
performed with all parameters set according to the wind tunnel experiment, and the user-
defined function (UDF) was used to input the inlet wind speed to fit the environmental
result. The ambient wind and the leakage source inlet were set as the velocity inlet, and the
wind tunnel outlet was set as the outflow outlet. In addition, the buildings and the ground
Processes 2021, 9, 1754 6 of 20
were prescribed non-slip walls, and other boundary conditions were kept symmetrical. The
roughness was the default value on the standard wall condition. In addition, the solution
method was the Simple mode, and the residual error control was 10−4 .
Figure 4. (a) Typical ideal street canyon model; (b) grid independent verification.
Figure 5. Time-averaged streamlines and pressure contours on the horizontal plane (z/H = 0.01) of
streets with different aspect ratios: (a) Rw = 2.0; (b) Rw = 1.0; (c) Rw = 0.5.
Figure 6 shows the time-averaged streamlines and flow velocity contours of street
buildings with different aspect ratios in the middle horizontal plane (z/H = 0.5). None of
the main separation curve a exists, curve b remains, and the entire streamlines are in a state
of contraction. The distribution of vortices in the flow field in the middle of the buildings
is consistent with the flow field on the ground. By comparing the pressure contours in
Figure 5, the extremely low pressure around the buildings is concentrated on the center
of each vortex, whereas the extremely low velocity is focused on the position where the
vortex intersects.
Processes 2021, 9, 1754 8 of 20
Figure 6. Time-averaged streamlines and flow velocity contours on the horizontal plane (z/H = 0.5)
of streets with different aspect ratios: (a) Rw = 2.0; (b) Rw = 1.0; (c) Rw = 0.5.
Figure 7 shows the time-averaged streamlines and flow velocity contours of the
vertical plane in the middle of street buildings (y = 0) with different aspect ratios. The
longitudinal streamlines of the three street canyons are similar, and the vortex centers
are close to the top of the downstream buildings. As the street aspect ratio decreases, the
vortex shape on the vertical section tends to be complete. The air generally flows in near
the windward side of the downstream buildings and flows out near the leeward side of the
upstream buildings.
Processes 2021, 9, 1754 9 of 20
Figure 7. Time-averaged streamlines and flow velocity contours on the vertical plane (y = 0) of streets
with different aspect ratios: (a) Rw = 2.0; (b) Rw = 1.0; (c) Rw = 0.5.
Figure 8 shows the instantaneous concentration fields of the wind tunnel experiment
on the middle vertical plane of street buildings (y = 0) with different aspect ratios. Results
show that the gas released during the experiment is dispersed in agglomerates. During
the dispersion process, the air mass expands, and the diameter of the pollutant plume
continues to increase. Buildings aggravate the turbulence of the surrounding airflow,
thereby intensifying the dispersion of air masses. After passing through the buildings,
the thick air masses disperse into flocs and continue to spread. As shown in Figure 8,
when pollutants enter street canyons from the top of the windward side of the downstream
buildings, it is consistent with the longitudinal vortex distributions in Figure 7. Due to the
different widths of street canyons, the dispersion spaces of pollutants after entering each
street canyon are also different. During the experiment, it was observed that most of the
pollutants were concentrated on the top of the deep street canyon, and flowed in a circular
shape, and were not easily dispersed to the bottom, whereas pollutants easily spread to the
ground in the typical ideal street canyon, and the degree of air mass agglomeration was
also lower than in the deep street canyon.
Processes 2021, 9, 1754 10 of 20
Figure 8. Instantaneous concentration fields on the vertical plane (y = 0) of streets with different aspect ratio: (a) Rw = 2.0;
(b) Rw = 1.0; (c) Rw = 0.5.
Figure 10. Contours of street ventilation surfaces with different aspect ratios: (a) Rw = 2.0; (b) Rw = 1.0; (c) Rw = 0.5.
Processes 2021, 9, 1754 11 of 20
Figure 11 shows the percentage of pollutant flux per unit area of street ventilation
surfaces with different aspect ratios. As the street aspect ratio Rw declines, the proportion of
pollutant flux entering the street canyon through both ends of the street increases, whereas
that from the top decreases, and the gap between the two declines as Rw decreases. Results
show pollutants entering the canyons through the top account for a larger proportion of the
streets in different configurations. As shown in Figures 5 and 6, as the street aspect ratio Rw
decreases, the development space of the vortex formed at the ends of the street increases,
and the vortex shape is complete on the horizontal plane. When Rw =0.5, the distribution
range of these two vortices almost occupies the entire street canyon. In addition, this part
of the flow carries more pollutants into the street canyon, increasing the proportion of the
street canyon lateral flux.
Figure 11. Flux proportion of pollutant per unit area of street ventilation surfaces with different
aspect ratios.
Figure 12 shows the time-averaged concentration fields of the wind tunnel experiment
on the vertical plane in the middle of street buildings (y = 0) with different aspect ratios.
Most of the pollutants are intercepted by upstream buildings or continue to spread above,
and the pollutant concentrations in the street canyons are very low. The relative concen-
tration per unit area C* of the street canyon was analyzed, as shown in Figure 13. It can
be seen that, as the aspect ratio of street canyon increases, the relative concentration per
unit area C* of the vertical plane in the middle of the street canyon also increases, which
is possibly related to the size of the dispersion space. Due to the blocking effect of the
buildings, the quantity of pollutants that can enter the canyons is limited. Therefore, the
larger the dispersion space in the canyon, the more the pollutants are dispersed, and C* is
relatively low.
Figure 12. Time-averaged concentration fields on the vertical plane (y = 0) with different aspect ratio streets: (a) Rw = 2.0;
(b) Rw = 1.0; (c) Rw = 0.5.
Processes 2021, 9, 1754 12 of 20
Figure 13. Relative concentration per unit area C* on the vertical plane (y = 0) with different aspect
ratio streets.
Figure 14. Time-averaged streamlines and flow velocity contours on the horizontal plane (z = 0) of
streets with different height ratios: (a) RH = 0.5; (b) RH = 1.0; (c) RH = 2.0.
Comparing Figure 15a and c of the internal flow fields of the street canyons, it can be
seen that the streamlines are emitted from the side of the higher buildings, which proves
that the directions of airflow flowing into the street canyons are not consistent, as shown in
Figure 16. Figure 16 shows the time-averaged streamlines and flow velocity contours of the
middle vertical plane of street buildings (y = 0) with different height ratios. Clearly, high
buildings dominate the distribution of the surrounding flow fields. In the step-up street
canyon, displayed in Figure 16a, the vortex formed in the street canyon is more complete
than the vortex in Figure 16b due to the greater interception of the high building, and the
center of the vortex (see Figure 16a) moves upstream. In the step-down street canyon, the
main vortex formed by the high building is located at the top of the wake recirculation
cavity, as shown in Figure 16c. The deformation vortex having a direction opposite to
that of the normal wake vortex inside the street canyon is formed by the interaction of
the buildings on both sides. According to Figures 15 and 16, the step-up street canyon is
located inside the recirculation area of the upstream low building, and the downstream
high building more strongly intercepts the wind flow. Therefore, the internal airflow
velocity of the step-up street is lower than that of the other two situations, which is likely
to affect the circulation and dispersion of pollutants.
Processes 2021, 9, 1754 14 of 20
Figure 15. Time-averaged streamlines and flow velocity contours on the horizontal plane (z/H = 0.5)
of streets with different height ratios: (a) RH = 0.5; (b) RH = 1.0; (c) RH = 2.0.
Figure 16. Time-averaged streamlines and pressure contours on the vertical plane (y = 0) of streets
with different height ratios: (a) RH = 0.5; (b) RH = 1.0; (c) RH = 2.0.
Processes 2021, 9, 1754 15 of 20
Figure 17 shows the concentration contours of street ventilation surfaces with different
height ratios. The larger the windward area, the stronger the interception effect on pollu-
tants. High buildings intercept most of the pollutants because the pollutants accumulate
on their windward side. As shown in Figure 17a, the concentration of pollutants in the
street canyon is relatively high. However, as shown in Figure 17c, only a small quantity
of pollutants enters the flow from both ends of the street canyon through the bypass, and
very few pollutants can enter from the top of the street canyon over the high building. In
addition, the high concentration of pollutants at the top of the street canyon shifts from
near the downstream building to the vicinity of the upstream building, which is consistent
with the movement of the inner vortex center, as shown in Figure 17a.
Figure 17. Contours of street ventilation surfaces with different height ratios: (a) RH = 0.5; (b) RH = 1.0; (c) RH = 2.0.
Figure 18 shows the flux proportion of pollutants per unit area of street ventilation
surfaces with different height ratios. In contrast to the typical ideal street canyon, the
pollutant flux ratios of each ventilation surface in the step-up or step-down street canyon
are similar. However, as shown in Figure 17, a significant difference exists between the net
pollutant flux of different terrains. Regarding the step-up street canyon, the pollutant flux
ratio still conforms to the rule that the top is greater than the two ends, but the difference
between the flux ratios is very small, indicating that the downstream higher building plays
a leading role in intercepting pollutants. In addition, pollutants accumulate and circulate
throughout the entire street canyon. On the contrary, for the step-down street canyon,
pollutants are intercepted outside the upstream high building, and the net flux from each
ventilation surface into the street canyon is very low, resulting in a similar phenomenon.
Figure 18. Flux proportion of pollutants per unit area of street ventilation surfaces with different
height ratios.
Processes 2021, 9, 1754 16 of 20
Figure 19 shows the time-averaged concentration fields of the wind tunnel experiment
on the vertical plane in the middle of street buildings (y = 0) with different height ratios.
Figure 20 shows the relative concentration per unit area C* of the average concentration in
the street canyon. The step-up street canyon can easily retain upstream pollutants because
the internal concentration is significantly higher than that of the other two street canyons.
In addition, the step-down street can effectively prevent the spread of pollutants, and its
internal pollutant concentration is the lowest among the three configured street canyons.
Figure 19. Time-averaged concentration fields on the vertical plane (y = 0) with different height ratios streets: (a) RH = 2.0;
(b) RH = 1.0; (c) RH = 0.5.
Figure 20. Relative concentration per unit area C* on the vertical plane (y = 0) of streets with different
height ratios.
As a result, the proportion of pollutants entering from the top of the street canyon is
correspondingly reduced.
Figure 21. Flux proportion of pollutants per unit area of street ventilation surfaces with different
characteristic wind speeds.
Figure 22 shows the time-averaged concentration fields of the wind tunnel experiment
on the middle vertical plane of the typical ideal street canyon (y = 0) under different wind
speeds. The distribution of pollutants under the four configurations is similar overall. The
difference is that, as the environmental wind speed increases, the pollutant plume height
declines, the pollutant dispersion speed increases, and the concentration at the plume tail
decreases. The relative concentration per unit area C* in the street canyon first rises and
then decreases with the increase in environmental wind speed, as shown in Figure 23. The
impacts of increased environmental wind speed on the dispersion of pollutants are mainly
reflected in two aspects. The increase in wind speed improves the pressure gap between
the inside and outside of the street canyon, which intensifies air exchange and promotes
the entry of pollutants. However, the dispersion of pollutants is also aggravated, making
it difficult for pollutants to accumulate. The overall concentration of the plume within a
certain volume is reduced, and the pollutants that can enter the canyon also decrease. Thus,
the most unsafe wind speed exists when the leakage source is outside the street canyon,
thereby maximizing the possible accumulation of pollutants in the street canyon. In this
experiment, this unsafe wind speed was about 1.1 m/s.
Figure 22. Time-averaged concentration fields on the vertical plane (y = 0) of a typical ideal street
with different characteristic wind speeds: (a) u0 = 0.9 m/s; (b) u0 = 1.1 m/s; (c) u0 = 1.3 m/s;
(d) u0 = 1.5 m/s.
Processes 2021, 9, 1754 18 of 20
Figure 23. Relative concentration per unit area C* on the vertical plane (y = 0) of the typical ideal
street with different characteristic wind speeds.
5. Conclusions
This paper focused on the effect of street canyon terrain on the surrounding flow
field and the dispersion of heavy gas pollutants, focusing on the impacts of various aspect
and height ratios and different wind speeds in the typical ideal street canyon. Wind
tunnel experiments and laser particle tracking technology were adopted to visualize the
dispersion of pollutants, in combination with CFD simulation, to explore the manner in
which pollutants enter street canyons. The following conclusions can be drawn:
(1) When the heights of the street buildings are equal, and only the aspect ratio
changes, the impacts of buildings on both sides of the street on the wind field can be clearly
expressed. However, when the building heights of the street are different, the effect of
higher buildings on the surrounding flow field dominates, and determines the distribution
of the vortex in the street canyon.
(2) Pollutants disperse mainly through the top of the street into the canyon inside.
However, this proportion decreases with a lower street aspect ratio or a higher height ratio.
In addition, the interception effect of the street canyon on the upstream pollutants then
weakens, and the relative concentration per unit area C* also decreases.
(3) The growth of the environmental wind speed inhibits the initial upward mo-
mentum of pollutants, resulting in a decrease in the proportion of pollutants entering
the street canyon from the top. Furthermore, the increase in environmental wind speed
increases the pressure difference between the street canyon and outside, promoting the
airflow circulation. However, it also aggravates the dispersion of pollutants and reduces
the overall concentration of the plume. Therefore, the combined influence of positive and
negative effects generates an unsafe wind speed that results in the highest concentration of
pollutants inside the street canyon.
References
1. Hu, W. Production Technology of Fluorine Chemical Industry(I); Science Press: Beijing, China, 2010; p. 1. (In Chinese)
2. Law, W.P.; Erain, N.; Ramli, N.I.; Gimbun, J. Assessment of chlorine leak dispersion around Gebeng industrial area and potential
evacuation route. Atmos. Res. 2019, 216, 117–129. [CrossRef]
3. Li, J.; Zhang, B.; Tang, S.; Tong, R. Application of FLUENT in simulating outcomes from chlorine leakage accidents in a typical
chemical factory. Toxicol. Ind. Health 2013, 32, 919–935. [CrossRef] [PubMed]
4. Yang, S.; Jeon, K.; Kang, D.; Han, C. Accident analysis of the Gumi hydrogen fluoride gas leak using CFD and comparison with
post-accidental environmental impacts. J. Loss Prevent. Proc. 2017, 48, 207–215. [CrossRef]
5. Deaves, D.M. Dense gas dispersion modelling. J. Loss Prevent. Proc. 1992, 5, 219–227. [CrossRef]
6. Rauchegger, C.; Bayley, S.; Schröder, V.; Thevenin, D. Dispersion of heavy gases-experimental results and numerical simulations.
Process Saf. Prog. 2015, 34, 280–285. [CrossRef]
7. Sklavounos, S.; Rigas, F. Validation of turbulence models in heavy gas dispersion over obstacles. J. Hazard. Mater. 2004, 108, 9–20.
[CrossRef] [PubMed]
8. Tan, W.; Wang, K.; Li, C.; Liu, L.; Wang, Y.; Zhu, G. Experimental and numerical study on the dispersion of heavy gases in urban
environments. Process Saf. Environ. 2018, 116, 640–653. [CrossRef]
9. Tan, W.; Li, C.; Wang, K.; Zhu, G.; Wang, Y.; Liu, L. Dispersion of carbon dioxide plume in street canyons. Process Saf. Environ.
2018, 116, 235–242. [CrossRef]
10. Xin, B.; Wan, L.; Dang, W. Dynamic Simulation and Quantitative Risk Assessment of Indoor Heavy Gas Diffusion. J. Phys. Conf.
Ser. 2020, 1549, 042080. [CrossRef]
11. Mao, N.; Wang, Z.; Zhou, C.; Tong, X. Qualitative and Quantitative Investigation on Concentration of Indoor Dense Gas Caused
by Leakage and Diffusion. Environ. Prog. Sustain. 2019, 38, 424–434. [CrossRef]
12. Marucci, D.; Carpentieri, M. Effect of local and upwind stratification on flow and dispersion inside and above a bi-dimensional
street canyon. Build. Environ. 2019, 156, 74–88. [CrossRef]
13. Yang, R.; Zhang, J.; Shen, S.; Li, X.; Chen, J. Numerical Investigation of the Impact of Different Configurations and Aspect Ratios
on Dense Gas Dispersion in Urban Street Canyons. Tsinghua Sci. Technol. 2007, 12, 345–351. [CrossRef]
14. George, X.; Lim, A.; Gopalan, H.; Jing, L.; Poh, H.J. CFD Simulation of Chemical Gas Dispersion Under Atmospheric Boundary
Conditions. Int. J. Comp. Meth. 2019, 16, 1940011.
15. Fiates, J.; Santos, R.; Neto, F.F.; Francesconi, A.Z.; Simoes, V.; Vianna, S.S.V. An alternative CFD tool for gas dispersion modelling
of heavy gas. J. Loss Prevent. Proc. 2016, 44, 583–593. [CrossRef]
16. Garcia-Sanchez, C.; Tendeloo, G.V.; Gorle, C. Quantifying inflow uncertainties in RANS simulations of urban pollutant dispersion.
Atmos. Environ. 2017, 161, 263–273. [CrossRef]
17. Pontiggia, M.; Derudi, M.; Alba, M.; Scaioni, M.; Rota, R. Hazardous gas releases in urban areas: Assessment of consequences
through CFD modelling. J. Hazard. Mater. 2010, 176, 589–596. [CrossRef]
18. Ferreira, E.S.; Vianna, S. Large Eddy Simulation Combined with Equivalent Diameter for Turbulent Jet Modelling and Gas
Dispersion. Braz. J. Chem. Eng. 2016, 33, 525–540. [CrossRef]
19. Sávio, S.V.; Vianna, R. Stewart Cant. Modified porosity approach and laminar flamelet modelling for advanced simulation of
accidental explosions. J. Loss Prev. Process. Ind. 2010, 23, 3–14.
20. Ferreira, T.D.; Santos, R.G.; Vianna, S.S. A coupled finite volume method and Gilbert–Johnson–Keerthi distance algorithm for
computational fluid dynamics modelling. Comput. Methods Appl. Mech. Eng. 2019, 352, 417–436. [CrossRef]
21. GB/T 3840-1991. Technical Methods for Establishing Local Air Pollutant Emission Standards, National Standard of China. 1991.
Available online: http://openstd.samr.gov.cn/bzgk/gb/newGbInfo?hcno=D308374668604004F828200717A3DB5C (accessed on
24 September 2021).
22. Cui, P.Y.; Li, Z.; Tao, W.Q. Numerical investigations on Re-independence for the turbulent flow and pollutant dispersion under
the urban boundary layer with some experimental validations. Int. J. Heat Mass Transf. 2017, 106, 422–436. [CrossRef]
23. Ohba, M. Experimental studies for effects of separated flow on gaseous diffusion around two model buildings: Effects of flow
conditions of oncoming flow. J. Archit. Plan. Environ. Eng. (Trans. AIJ) 1989, 406, 21–30.
24. Luo, Q. Wind Tunnel Experimental Study on Air Pollutant Transportation and Dispersion for Hilly Terrain. Ph.D. Thesis, Tianjin
University, Tianjin, China, 2008.
25. Chan, T.L.; Dong, G.; Leung, C.W. Validation of a two-dimensional pollutant dispersion model in an isolated street canyon. Atmos.
Environ. 2002, 36, 861–872. [CrossRef]
26. Shi, T.; Ming, T.; Wu, Y.; Peng, C.; Richter, R.D. The effect of exhaust emissions from a group of moving vehicles on pollutant
dispersion in the street canyons. Build. Environ. 2020, 181, 107120. [CrossRef]
27. Chang, J.C.; Hanna, S.R. Technical Description and User’s Guide for the BOOT Statistical Model Evaluation Software Package, Version
2.0.; George Mason University and Harvard School of Public Health: Fairfax, VA, USA, 2005.
28. Jiang, G.; Yoshie, R. Side ratio effects on flow and pollutant dispersion around an isolated high-rise building in a turbulent
boundary layer. Build. Environ. 2020, 180, 107078. [CrossRef]
Processes 2021, 9, 1754 20 of 20
29. Barros, D.; Borée, J.; Noack, B.R.; Spohn, A.; Ruiz, T. Bluff Body Drag Manipulation Using Pulsed Jets and Coanda Effect. J. Fluid
Mech. 2016, 805, 422. [CrossRef]
30. Posa, A.; Oresta, P.; Lippolis, A. Analysis of a directional hydraulic valve by a Direct Numerical Simulation using an immersed-
boundary method. Energy Convers. Manag. 2013, 65, 497–506. [CrossRef]