Large-Eddy Simulation of Airﬂow and Pollutant Dispersion in a Model Street Canyon Intersection of Dhaka City

: The atmospheric ﬂow and dispersion of trafﬁc exhaust were numerically studied in this work while considering a model street canyon intersection of a city. The ﬁnite volume method (FVM)-based large-eddy simulation (LES) technique in line with ANSYS Fluent have been used for ﬂow and pollutant dispersion modelling through the consideration of the atmospheric boundary layer (ABL). Hexahedral elements are considered for computational domain discretization in order to numerically solve problems using FVM-LES. The turbulence parameters were superimposed through a spectral synthesizer in the existing LES model through ANSYS Fluent as part of ‘damage control’ due to the unsteady k − (cid:101) simulation. Initially, the code is validated with an experimental study of an urban street canyon where the width and height ratio is in unity. After validation, a model urban street canyon intersection was investigated in this work. The model shows a high pollutant concentration in the intersecting corner areas of the buildings. Additionally, the study of this model intersection shows a high level of pollutant concentration at the leeward wall of downwind building in the case of increased height of an upwind building. Most importantly, it was realized from the street intersection design that three-dimensional interconnection between the dominating canyon vortices and roof level ﬂow plays a pivotal role in pollutant concentration level on the windward walls. The three-dimensional extent of corner eddies and their interconnections with dominating vortices were found to be extremely important as they facilitate enhanced ventilation. Corner eddies only form for the streets towards the freeway and not for the streets towards the intersection. The results and key ﬁndings of this work offer qualitative and quantitative data for the estimation, planning, and implementation of exposure mitigation in an urban environment.


Introduction
The term "street canyon" simply refers to a type of narrow street between buildings lining up in continuous order in both sides, which mostly occurs in urban areas. Street canyons can exhibit a different climate from the nearby areas, specifically when the areas are heavily influenced by the micro-scale meteorological process [1,2]. Another challenge concerning street canyons is ventilation and pollutant dispersion, where the latter indicates the transportation of air pollutants through the outdoor atmosphere. In any case, both processes usually take place through the roof level in a street canyon [1,3,4]. Therefore, air the residents. Additionally, by adhering to fluid dynamical resemblance standards, the effort assured that the obtained results could be transferred to a full-scale design. The results demonstrated the dependability of the procedure that was used. Gromke then constituted the Concentration Data of Street Canyons (CODASC) database to contribute to the wind tunnel experimental data of pollutant concentrations for urban street canyons research with variable aspect ratios and various wind directions. This database aims to contribute wind-tunnel experiment data for the validation of future pollutant dispersion modelling. A detailed study by Gromke et al. [21] demonstrated that a parallel incoming flow lowers concentration by up to 30% on the façades and by up to 60% at pedestrian level. Therefore, the importance of ground or pedestrian level pollutant dispersion needs to be considered in any realistic model. To the authors' knowledge, there is no work available in the literature that considered the presence of an unconventional adjacent street canyon that forms a street intersection. In an intersection design, corner eddies form differently, and their three-dimensional interconnection to dominating vortices inside the canyon is different compared to conventional street canyon designs.
This explains the purpose of this study. The aim of this study is to use the LES turbulence model to study the dispersion of traffic exhaust and flow physics of a threedimensional geometry consisting of four separate ground-surface-mounted buildings with a square intersection. These four buildings form four street canyons in total. Two of these street canyons are perpendicular to the approaching wind flow and other two street canyons are parallel to the approaching wind flow. In numerical modelling, the boundary layer needs to be defined and was reported in the literature [22][23][24]. In this study, a developing boundary layer was considered for the approaching flow to represent the Atmospheric Boundary Layer (ABL). The considered scenario is based on the windtunnel experimental study of Gromke and Ruck [18][19][20]25]. The wind tunnel experiment considered small-scale urban street canyon model. The model studied the velocity profile and concentration measurements of two parallel buildings for an approaching boundary layer flow perpendicular to the canyon centre axis. The current work explores increased H/W (height to width) and W/H (width to height) ratio with canyons entering and leaving the square intersection which makes the flow properties challenging yet more efficient and pragmatic compared to previous studies. This work further focuses on the evolution of pollutant concentrations. Additionally, the origin of the corner vortices is also addressed efficiently.

Filtered Navier-Stokes and Concentration Equations
The time-dependent Navier-Stokes equations are filtered in either Fourier space or configuration space to obtain the governing equations for LES. The eddies with smaller scales, as compared to the filter expanse, or grid spacing adopted for the simulations are efficiently filtered out in the filtering operation. Accordingly, the dynamics of large eddies are governed by the derived equations.
The following forms are obtained after filtering the continuity, momentum, and pollutant concentration equations ∂ ∂t where µ denotes molecular viscosity, c represents the species or pollutant concentration, and D is the diffusivity. The sub-grid-scale stress is denoted by τ ij , Sc t represents the turbulent Schmidt number (Sc t = µ t ρD t , where µ t represents the turbulent viscosity and D t refers to the turbulent diffusivity). Sc t = 0.7 is the default value in ANSYS Fluent™.

Large-Eddy Simulation (LES)
Eddies describe the turbulent motions of a fluid by a broad spectrum of length and time scales. The characteristic expansion of the time-averaged flow (example: shear layer thickness) is a typical comparable method to describe the extent of the largest eddies. The dissipation of TKE (turbulence kinetic energy) is generated by the smallest eddies.

Subgrid-Scale (SGS) Model
Filtered sub-grid-scale stresses are not known and, therefore, need modelling. As the RANS models, ANSYS Fluent™'s sub-grid-scale turbulence model exercises the Boussinesq approximation [26] for calculating sub-grid-scale turbulent stresses from where µ t represents the sub-grid-scale turbulent viscosity. The isotropic component of sub-grid-scale stress, τ kk , is appended to the filtered static pressure expression as it is not generally considered in the numerical models.

Smagorinsky-Lilly Model
Smagorinsky [13] originally introduced this simplistic model. The eddy-viscosity in the Smagorinsky-Lilly model is modelled as: where L s denotes the mixing length of sub-grid-scales, and S ≡ 2S ij S ij . ANSYS Fluent calculates the L s from the term: where κ signifies the von Kármán constant, d signifies the distance to the nearest wall, C s is the Smagorinsky constant, and ∆ = (∆x∆y∆z) 1 3 indicates the filter width. C s = 0.23 was determined by Lilly for homogeneous isotropic turbulence in the inertial sub-area. This value of C s still contributes to the inconsistencies in the appearance of mean shear and transitional flows as the near-solid-boundary and must be lowered in such areas. Briefly, C s does not have a universal constant value; this is considered to be the one of the major shortcomings of this method. C s ≈ 0.1; nevertheless, it gives the most excellent outcomes for an extensive variety of flows, and ANSYS Fluent uses this value as a default.

Setup of Wind-Tunnel Experiment
The experiment was conducted to estimate concentration and flow data in a reduced scale street canyon with an atmospheric boundary layer wind tunnel. An approaching flow with mean velocity (u(y)), profile exponent α = 0.3, and turbulence intensity (I u ) profile exponent α l = 0.36, characterized by the power-law approach, were considered in that work. The expressions of these profiles were: and Air velocity u(y re f = H) = 4.65 m s −1 was used. Here, H presents the building height. The works of Gromke and Ruck [18,25] provide broad details on the computational implementation of atmospheric boundary layer flow as well as the information of the integral length scale profile L ux (y) and spectral arrangements of turbulent kinetic energy S uu (y, f ). A 1:150 scaled design of a separate street canyon of length L = 180 m and street width W = 18 m with two surrounding parallel structures representing buildings of height H = 18 m and width B = 18 m were positioned perpendicular to the approaching flow for the experiment.
Tracer gas emitting from four-line sources blended within the model street and were outlined based on the approach specified in Meroney et al. [27]. They were then used to simulate the discharge of traffic pollutants. The sources passed the street canyon by about 10% on every side to consider the traffic exhaust discharged in sideways crossings. SF 6 (Sulfur Hexafluoride) was the tracer gas. Mean concentration of the SF 6 was estimated on the canyon walls and was then normalized as follows: where c is the measured mean concentration, u H is the developed approaching airflow mean velocity at height H, and Q/l is the tracer gas source intensity per unit length.

Computational Domain and Boundary Conditions
A large-eddy simulation was carried out using the CFD code ANSYS Fluent™ 15.0 to generate the detailed data of experimental work above, focusing on the pollutant concentrations and flow patterns inside a street canyon. Figures 1 and 2 show the computational domain [28][29][30][31][32][33] and boundary conditions considered in this work. A 1:150 scaled design, similar to the experimental setup, was positioned perpendicular to the approaching flow for the simulation.
Hexahedral elements have been used for the computational domain discretization based on the finite volume method (FVM). A total of 3,015,321 nodes and 2,931,392 cells were used for the simulation. Since the appearance of high gradients in the solution variables, i.e., separation, reattachment, and recirculation, occur in the proximity of the buildings and street canyon, these areas are the region of interest and, therefore, most of the cells were arranged in these regions. The generated mesh has a volume-weighted average orthogonal-quality of 0.998 and skewness of 0.03; these values confirm the 'outstanding' [34] quality of the mesh.
User-defined functions (UDFs) were executed in the ANSYS Fluent™ to maintain the WT experiment characteristics. The inlet air velocity profile follows the power-law profile of:   The present boundary conditions guarantee the horizontal homogeneity of the turbulent boundary layer [28][29][30][31][32][33]35,36]. The nonexistence of streamwise gradients in the vertical characterizations of mean air velocity and turbulence quantities are termed 'horizontal homogeneity', i.e., downstream length is preserved in these characterizations (Blocken et al. [35]). A sand grain roughness, K s , represents the surface roughness in ANSYS Fluent™ rather than the aerodynamic roughness, y 0 . Gromke et al. [36] used K s according to aerodynamic roughness height, y 0 = 0.0033 m in the WT experiment, in order to avoid possible difficulties amidst a coarse grid analysis close to the ground caused by a high sand grain roughness value. Gromke et al. [36] realized that K s being equal to y 0 was not accurate in strong means, yet it was justifiable as it yielded desirable results. The turbulence intensity was low close to the ground (z/H < 1), where changes were less than 10%. K s = y 0 = 0.0033 m was also set in this work to address the surface roughness. Salim et al. [28][29][30][31][32][33] similarly examined this issue in an 'empty' domain.
The ground and building surfaces were treated as smooth walls by default in ANSYS Fluent™, according to the above wall roughness parameters, for the present LES.
The converged time-averaged solution, produced by an unsteady realizable k − simulation,was synthesized by applying the ANSYS Fluent™'s inbuilt spectral synthesizer in the current LES. This 'spectral synthesizer' mechanism follows the characteristics of the work of Smirnov et al. [37]. This practice essentially superimposes turbulence parameters in order to interpret the initial and boundary conditions according to the WT experiment to be applied in the LES.

Flow Simulation
The Smagorinsky-Lilly SGS model with 'Bounded Central Differencing' spatial discretization scheme for momentum and 'Second-Order Upwind' for species and energy transport equations were applied for this computation. The pressure and pressure-velocity coupling were solved using the 'Second-Order' spatial discretization and semi-implicit pressure-linked Equation (SIMPLE) scheme in this work. The present LES considers 10 −3 'Absolute Criteria' to address convergence for the scaled residuals.
The time-step size for this simulation is 0.00059. The volume-weighted average cell convective Courant number, also known as CFL, during the simulation is 0.33. At first, 77 flow-through computational time (flow-through time is T = L D /U b , where L D is the streamwise length of the domain and U b is the bulk velocity) was calculated in the simulation. The flow reaches a statistical steady-state at this stage. Then, the calculation was carried out for a further 77 flow-through computational time after resetting the initial flow statistics to assure that the decisive time-averaged outcomes were free of the first conditions. The selected 77 flow-through time was compared to the mean solutions of 90 flow-through time in order to check the convergence of the solution. The mean flow field did not vary between these two sampling sizes, additionally confirming that the solution converged.

Dispersion Method
The needed traffic exhaust locations were earmarked in the regions of volume in the geometry and defined as separate fluid zones in order to simulate the line sources. A discharge rate of Q = 10 gs −1 with fluid material SF 6 as the source was then installed for these zones to follow the WT settings. Figure 3 shows more details.  A total of 20, 000 time-steps were calculated for the present LES time-averaged results following an initial 20, 000 time-step calculation to reach a steady state. A sampling period of 105 s was considered for the wind tunnel concentration data [18]. It is seen from a qualitative overview of Figure 4 that the present LES can offer efficient results of the spatial pattern of the pollutant. The region of the centre line (z/H = 0) at both walls should be most articulated regarding the concentration pattern, since the highest concentration occurs here and is therefore regarded as the most crucial area.  [29] contained some flow-field and pollutant concentration data along the z/H = 0 line within the street canyon. Salim et al. [29] and several other studies confirmed that LES is the most reliable prediction tool for validating flow fields and turbulence structures in corresponding ex-periments [28,29,38,39]. It can be noticed from Figures 5 and 6 that the LES prediction of pollutant concentrations at various points of spanwise locations is indeed good.

Area of Interest: Dhaka
Dhaka is the capital and the biggest city of Bangladesh. It is the main location for government administration, business, and other industries. Due to the rapid growth and lack of urban planning, environmental and socio-economic problems are increasing [40]. Street canyon intersections in Dhaka see the highest number of traffic compared to any other streets. Most of the street intersections are highly unconventional in design when considering building height and width ratios. That is why an unconventional street intersection design was studied in this work, simplifying the most commonly seen intersection designs in Dhaka, Bangladesh (see Figure 2).
One of the dominant contributors of pollutant dispersion in Dhaka is unfit motor vehicles within narrow streets causing immense traffic congestion. Black carbon (BC) is one of toxic components reported in studies of unfit vehicles in Dhaka. According to Salam et al. [40], the averaged fine particulate matter (FPM) in a typical Dhaka street can go 15 units above the standard value, which raises the temperature into a hot condition even in winter. Therefore, relying on weather conditions is not a realistic approach in terms of street canyons in Dhaka due to the unpredictability of temperatures. While the impact of temperature was reported in the literature, most studies focused on well-developed and planned urban areas in Europe and Australia, for example [41,42]. However, in a city such as Dhaka with almost 21 million people in a 306.4 km 2 area, both roof level and ground level pollutant concentrations need to be monitored concurrently due to the narrow streets and large population.

Pollutant Concentration Data Overview
The developed numerical method in Section 3 was used to run a simulation considering an urban intersection to analyse pollutant concentrations on the street canyon walls. Figure 2 shows the new building distribution and walls considered for pollutant concentration analysis. A 1:150 scaled design was positioned perpendicular to the approaching flow for the simulation. The validated computational domain [28][29][30][31][32][33] and boundary conditions, described in Figure 1, were considered for this simulation as well.
The meshing method used in the developed numerical strategy, previously discussed in Section 3, was maintained in this simulation. The meshing method generates 2,802,615 nodes and 2,734,066 cells with hexahedral elements. Again, most nodes and cells were generated in the proximity of the buildings and street canyons since these are the areas of interest. Table 1 shows further information of meshes from the ANSYS design module. Meanwhile, the converged time-averaged solutions, referred in Section 3.2.1 are presented in Figure 7 demonstrating for both velocity( Figure 7a) and TKE (Figure 7b) char-acterisations in the considered inlet boundary conditions. Overall, the correlations between the WTE and UDF for both attributes were found to be in great agreement. Figure 8 shows the grid independence test in terms of the velocity by considering two mesh sets Grid 1 and Grid 2 and found good agreement. Here, Grid 1 was utilized for the whole simulation. The time-step size for this simulation was 0.0018. The volume-weighted average cell convective Courant number, also known as CFL, during the simulation was 0.89. Initially, 31 flowthrough computational time was calculated to generate time-averaged data. The statistical steady-state was realized at this state. Then, the calculation was sampled for 1, 3, and 5 flow-through computational time after resetting the initial time-averaged flow statistics to demonstrate the time-averaged pollutant concentration evolution for the considered problem. Two more line sources, parallel to the incoming wind flow, were added to the considered model for the current study; more details can be seen in Figure 2. The newly added line sources also maintain the validated dispersion method (see Section 3.2.3) with a discharge rate of Q = 10 gs −1 and SF 6 as the pollutant fluid material.    Figures 9 and 10 present the time-averaged normalized pollutant concentration data at the building walls that were considered for analysis in the modelled urban environmentlabelled as A, B, C, and D in Figure 2. A total of 15,000 time-steps were computed for this simulation, maintaining the developed LES method (in Section 3). Figures 11-14 show the evolution of time-averaged normalized concentration data on the canyon walls based on the normalized flow time that was considered for this analysis. These data show that pollutant concentration was initially higher on the walls. Most of the pollutant concentration was close to the ground surface and in the midsection of the walls after the first few normalized flow time points. The wind flow started to come inside the canyons at this stage, and the dispersion of pollutants began. Unchanging dispersion processes, which is the final stage of the pollutant concentration, were observed at the later stages of the flow; the numerical solution converged at this point. By analysing the converged data, it could be seen that the concentration of pollutants was mostly in the area that was closer to the intersection. The pollutant concentration was relatively higher on the section of the street canyon walls closest to the intersection, rather than the walls' section further from the intersection. This phenomenon is also true for the elevation of the building walls. Therefore, the higher elevation showed lower pollutant concentrations.

Pollutant Concentration on Upwind Buildings
Amongst the canyon walls perpendicular to the incoming wind flow, wall A1 was positioned at the upwind side of the domain and had a relatively higher concentration. Building A had a height of 1.5H, which is 0.5H (higher than the reference height (H)). This elevated height prevents the wind from entering or leaving the street canyon (the street canyon between Building A and C) effectively, as can be seen from the streamlines of the velocity magnitude U /u H in Figures 15 and 16. Both figures further showed that there was not any big circulating vortex in the midsection of the street canyon used to carry the pollutants out of the street canyon, but just a small circulating vortex that formed close to Building A. According to Figure 16, the vertical flow had a relatively low velocity inside the street canyon, thus explaining the higher level of pollutant concentration seen in Figure 11 (Wall A1).  Pollutant concentrations on wall B1, also located at the upwind side of the building, were relatively lower compared to those on wall A1. The wind flow continuously entered and exited the street canyon (i.e., the street canyon between Buildings B and D) effectively, as seen from the streamlines of the velocity magnitude U /u H in Figures 15 and 17, since Building B retains the standard reference height of H. Figures 15 and 17 also confirmed the presence of a large circulating vortex inside the street canyon, which is the main factor for pollutant dispersion, along with a smaller vortex formation. The vertical flow velocity had a relatively higher magnitude inside the canyon as seen in Figure 17. The velocity magnitude and vertical velocity data justified the level of pollutant concentration shown in Figure 12 (wall B1).
In the street canyon, two features of the concentration pattern at the inner canyon walls could be observed, as depicted in Figure 4. First, we observed the deterioration of pollutant concentration at the street ending closest to the freeway at both walls A and B. Second, we observed significantly higher concentrations on wall A. However, in the intersection design, the decline in pollutant concentration was not seen for the street ending closest to the intersection; instead the pollutant concentration increased in this area. However, the reduction in concentration was only noticed at the street ending closest to the freeway. Walls A1 and B1, which were perpendicular to the incoming flow, showed a reduced concentration in the intersection design, as described in Figure 9, when compared to wall A in the street canyon. Peak values of concentrations on walls A1 and B1 were found to be highly dependent on the heights of Buildings A and B in the intersection design. The drop in the concentrations at the streets ending towards the freeway was the result of improved ventilation for both street canyon and intersection. Nevertheless, in the street intersection design, higher pollutant concentration at the street ending towards the intersection was still evident due to ineffective ventilation.
Walls A2 and B2 were situated in the upwind side of the domain, and they were parallel to the incoming wind flow. These two walls had a very low level of pollutant concentrations. Figure 15 showed that the incoming wind mostly flowed through the parallel street canyon (i.e., the small street canyon between Buildings A and B) without any interruption. This flow characteristics explain the level of pollutant concentrations on wall A2, as seen in Figure 11 (wall A2).

Pollutant Concentration on Downwind Buildings
As previously described in Section 4.3, upwind Building A had a relatively higher height than the reference height (H). The elevated height causes a poor ventilation facility inside this canyon (i.e., the canyon between Buildings A and C), thus preventing big circulating vortices from developing in this area. Only small vortices were formed inside this canyon, which failed to disperse the traffic-exhaust pollutants effectively. These characteristics of the flow caused a high pollutant concentration on wall C1, as shown in Figure 13 (wall C1), even though Building C was situated in the downwind side of the domain. Meanwhile, wall D1, also located at the downwind side of the domain, showed the lowest level of pollutant concentration, as shown in Figure 14 (wall D1), compared to other perpendicular street canyon walls even though Building D had a height of 1.5H, which is 0.5H higher than the reference height of H. However, this elevation of building D did not prevent the incoming flow from effectively entering or exiting the street canyon (i.e., the canyon between Buildings A and C).
The deterioration in concentration at wall B was seen for the street ending closest to the freeway, as shown in Figure 4, in the street canyon design. However, the downwind side of building wall B showed lower pollutant concentrations when compared to wall A. Furthermore, in the street intersection design, pollutant concentration increased at walls C1 and D1 for the street ending closest to the intersection, as seen in Figure 10. The reduction in concentration was also seen on walls C1 and D1 for the street ending closest to the freeway. As previously described, this observation proved that the streets ending closest to freeways have much better ventilation facilities when compared to streets ending closest to the intersection. Nonetheless, walls C1 and D1, which are both perpendicular to the incoming flow, had a big difference in pollutant concentrations. As mentioned earlier, in the street canyon design the downwind wall B experienced comparatively lower pollutant concentration. However, in the intersection design wall C1 experienced a much higher concentration compared to wall D1. This behaviour could be attributed to the elevated height of upwind Building A. Nevertheless, compared to the street canyon design, downwind side buildings still see much higher pollutant concentration for the street intersection design due to ineffective ventilation at the street ending closest to the intersection.
The pollutant concentration level on walls C2 and D2 was comparatively higher. Building C and D sit on the downwind side of the domain, and have a width of 1.5H, which is 0.5H higher than the reference width (H). Before entering the street canyon (i.e., the street canyon between Building B and D) the wind flow was interrupted and separated in several places, as seen in Figure 15. This interruption and separation of the incoming wind flow together with the increased width of theses two buildings have a significant influence on the pollutant concentration level, which explains the relatively higher pollutant concentration level on wall C2, as seen in Figure 13 (wall C2), and on wall D2, as seen in Figure 14 (wall D2).

Sampling Overview
The calculation was sampled at 1, 3, and 5 normalized flow-through computational time points after resetting the initial time-averaged flow statistics in order to demonstrate the time-averaged pollutant concentration evolutions for the considered problem. The sampling was conducted once the statistical steady-state was optimised to illustrate statistical data based on normalised flow instances. This process of capturing pollutant concentration evolution was previously considered for a street canyon [28][29][30]32]. Figures 11-14 show these normalized concentration data on the walls considered for analysis.

Key Findings on Each Building Block
The considered geometry was considered based on four building blocks. Since the position and dimension of the blocks varied in the computation, the key findings on each building block are discussed below.
Building Block A: The peak concentration was larger at initial instances of the normalized flow-time on wall A1-seen in Figure 11. After just 1 flow-through time (first instance), the peak concentration was approximately 26.67% higher compared to the final time-averaged concentration. Wall A2 showed approximately 33.33% higher peak concentration after the first instance, compared to the final time-averaged data. The peak concentration was reduced in vertical, streamwise, and spanwise expanses on both walls A1 and A2. The changes in vertical, streamwise, and spanwise extents were significant and noticeable. After the third and fifth instances, the concentration values kept changing, as seen in Figure 11. Every flow-through time (instance) caused the concentration to reduce by dispersing the pollutant, and then the concentration reached the final time-averaged stage when changes were no longer possible.
Building Block B: Wall B1 showed a high concentration value as anticipated since it was on the leeward side, as seen in Figure 12. This wall also showed a higher concentration after the first instance, and an approximately 66.67% higher peak concentration compared to the final mean concentration. The third and fifth instances demonstrated approximately 44.44% and 22.22% higher peak concentration, respectively. However, vertical and spanwise expanses of the concentration had significantly changed at these points. Wall B2 had an approximately 33.33% higher peak concentration, similar to wall A2. The changes in vertical and streamwise expanses of the concentration were noticeable for wall B2 after each instance (third and fifth) compared to the first instance.
Building Block C: Wall C1, even though it was on the windward side, experienced a high pollutant concentration, as displayed in Figure 13. In initial instances, and at the first, third, and fifth instances, it showed an approximately 63.64%, 59.09%, and 54.5% higher peak concentration compared to the final mean concentration, respectively. Furthermore, wall C2 showed up to an approximately 42.86% higher peak concentration in initial instances. In addition, similar every other wall, these two walls exhibited changes in pollutant concentration in vertical, streamwise, and spanwise expanses.
Building Block D: Wall D1, which was on the windward side, demonstrated a low pollutant concentration overall, as shown in Figure 14. This relatively lower concentration was in line with the initial hypothesis and similar to the windward wall of the street canyon. The lower concentration values were seen from the first instance. However, the peak concentration was still approximately 50% higher compared to the final mean concentration. Although the distribution of the concentration changes after the third instance, the peak concentration still remained approximately 50% higher compared to the final mean value. After the fifth instance, the peak concentration stayed approximately 25% higher. However, wall D2 shows a different pattern compared to other walls. The maximum difference of peak concentration was 58.33% higher compared to final mean concentration. This was seen after the third instance before the peak concentration started decreasing again. Similar to the other walls, walls D1 and D2 showed changes in pollutant concentration in vertical, streamwise, and spanwise expanses in every instance.

Significance of the Approach and Future Directions
The demand for understanding the pollutant dispersion in the presence of an unconventional adjacent street canyon should serve as a realistic yardstick rather than a simplified urban street design with fixed aspect ratios. The profitability of the template pollutant dispersion design presented in this study should benefit inhabitants of denselypopulated urban area such as Dhaka. As a result, if the pollutant concentration can be monitored through an automated model based on concentration level monitoring, a proper urbanization of a city can be implemented for the benefit of human beings. The pollutant concentration values including the peaks and the rate of plummeting reported in this study should be useful for analysing the 'building geometry' and its impact on airflow. Having access to pollutant concentration data on a regular basis, air quality control professionals could liaise with councils or property owners to understand the separation of buildings or walls including an in-depth understanding of height and width. As previously mentioned, most of the published works were based on geometries of urban buildings in terms of roof-level concentration monitoring and were based on geographic location. This finding makes sense to some extent. However, in a place such as Dhaka where traffic congestion is quite familiar in any commercial area, the street-level pollutant concentration needs to be considered. In addition, this could help air quality researchers in Dhaka create a map (similar to ESRI ArcGIS or Google Earth) to record and create a database for Dhaka's pollutant concentrations at both roof and street levels. Bangladesh working on becoming fully digital to benefit citizens through the programme "Bangladesh Vision 2041", where a smartphone application could be developed and made available for everyone for live air quality monitoring [43,44]. Therefore, a range of parametric modelling was implemented while developing this model, such as varying the aspect ratio, timescale, and grid/node size of the numerical simulations, and the conclusion on the peak concentration was reached. This will help provide significantly more information to urban dwellers or any other concerned citizen.
Based on the findings of the current study and the literature review, adequate evidence now exists to definitely emphasize the potential of the adjacent street canyon model into commercial applications or improvement of existing applications. Even in a developed country such as Australia, the Environment Protection Authority (EPA) Victoria often finds it difficult to track air pollutant concentration at a large scale in the state of Victoria [45,46]. The process would also replace the manual labour required to monitor the concentration levels at a short timescale because the changes in pollutant concentration would be recorded automatically through the model. While a specific area of interest was considered in this study, there is a room for improvement across a larger region with the possible implications of geographical coordinates. This process is already in practice in developed countries, but it should find future applications in Bangladesh.

Conclusions
Atmospheric flow and pollutant dispersion processes inside a conventional urban street canyon and a novel model urban intersection were presented in this study. At first, numerical computations were performed for the urban street canyon, and the obtained results were compared with previous wind tunnel experiments. The validated numerical method was applied later in order to carry out a numerical study for a model urban intersection. The presence of an unconventional adjacent street canyon forming a street intersection is one of the unique features of this study. The numerical calculation was sampled at one, three, and five normalized flow through time points. Prior to that, the initial time-averaged flow statistics were reset to explain the evolution of pollutant concentration within the considered geometry. For the considered urban intersection, the parts close to the urban-intersection always showed the maximum concentrations. Increasing the height of a building in the upwind side of the domain prevents flows to come inside of the A-C perpendicular street canyon and form strong canyon vortices, which essentially led to ineffective pollutant dispersion inside this street canyon. It was also observed that the highest pollutant concentration by wall A1 and A2 exhibited approximately 26.67% and 33.33% increases at the first instance compared to the final time-averaged concentration. Furthermore, wall B2 demonstrated an identical increase to wall A2 (33.33%), whereas the peak concentration at wall B1 showed a sharp increase of 66.67% after the first instance. The rate of increase plummeted to 44.44% and 22.22% after the third and fifth instances. Increased width of buildings in the downwind side of the domain also causes a high pollutant concentration inside of C-D parallel street canyon. Interestingly, the maximum pollutant concentrations at wall C1 after first, third, and fifth instances remained quite close to 63.64%, 59.09%, and 54.5%, respectively, whereas wall C2 initially recorded an almost 42.86% increase in the concentration, which is much greater than that of walls A2 and B2. Unlike the other walls, D demonstrated a lower pollutant concentration after the fifth instance (only 25% higher than the initial). Furthermore, it was also proven that unbalanced natural ventilation inside the street canyons worked as a catalyst for largely increasing the pollutant concentrations at the pedestrian level and on the street canyon walls. The importance of the height and width of buildings in the upwind side of the domain, considering the reference height H and width W, was cemented through numerical work that showed how the canyon vortex or canyon vortex-like structures were likely to develop. Over-height buildings, in the upwind side of the domain, block the flow entering into the canyon and thus fail to remove the polluted air. The findings of this research were found to be in accordance with this conceptual understanding as well.
The key findings of this study suggest that the 3D interconnection between dominating canyon vortices and roof level flow determine the pollutant concentration level on windward walls. Therefore, the extent of corner eddies and associated correlations with dominating vortices are responsible for enhanced ventilation. The flow structures and concentration profiles observed in the present study can be further analyzed and improved, particularly those identifying changes such as the presence of trees in the canyons. Ordinarily, numerical modelling considers porous structures to simulate trees in pollutant dispersion studies. More investigations on the flow and concentration profiles of the considered model urban intersection with tree plating inside the canyon's together with the appropriate numerical method, would be of interest in order to broaden this knowledge.  Acknowledgments: All the authors gratefully acknowledge the fifth author and the University of Leeds for providing the facilities needed to use ANSYS Fluent.

Conflicts of Interest:
The authors declare no conflict of interest.