Abstract
This study expands the study conducted by Zhang et al. (Boundary-Layer Meteorol, 2022, Vol. 183, 97–123) to elucidate turbulent structures within an ideal two-dimensional street canyon, and determine their contribution to pollutant removal. In response to the limitations of the reference study wherein spanwise turbulent structures were not elucidated, an advanced technique called spectral proper orthogonal decomposition (SPOD) with two-dimensional Fourier transformation was applied in the current study. This approach combined proper orthogonal decomposition along the spatial streamwise and vertical directions and Fourier decomposition along the time and spatial spanwise direction. Using this technique, various fluctuation patterns were reasonably decomposed and visualised according to their scales. In addition, their intensities and contributions to pollutant removal were quantitatively analysed using the SPOD spectrum and cospectra. The time scales of most energetic modes were found to be proportional to their spanwise length scales, regardless of their scales and flow mechanisms. The ranges of the spanwise wavenumber and frequency, at which pollutant removal events occurred at the roof level, were specified. These ranges coincided with those of small-scale structures caused by Kelvin–Helmholtz instabilities. Further complex analysis of the correlation between large- and small-scale structures showed that large-scale structures have a high possibility of indirectly enhancing or weakening pollutant removal by amplifying or suppressing small-scale structures. This occurred stochastically along both time and spanwise directions. Quantitatively, the amplitude of small-scale structures was amplified or suppressed by 13–16% on average by high- or low-momentum fluids.
Similar content being viewed by others
References
Araya DB, Colonius T, Dabiri JO (2017) Transition to bluff-body dynamics in the wake of vertical-axis wind turbines. J Fluid Mech 813:346–381. https://doi.org/10.1017/jfm.2016.862
Aubry N, Holmestil P, John L, Stone E (1988) The dynamics of coherent structures in the wall region of a turbulent boundary layer. J Fluid Mech 192:115–173. https://doi.org/10.1017/S0022112088001818
Barbano F, Brattich E, Di Sabatino S (2021) Characteristic scales for turbulent exchange processes in a real urban canopy. BLM 178:119–142. https://doi.org/10.1007/s10546-020-00554-5
Berkooz G, Holmes P, Lumley JL (1993) The proper orthogonal decomposition in the analysis of turbulent flows. Annu Rev Fluid Mech 25:539–575. https://doi.org/10.1146/annurev.fl.25.010193.002543
Blackman K, Perret L (2016) Non-linear interactions in a boundary layer developing over an array of cubes using stochastic estimation. Phys Fluids 28:095108. https://doi.org/10.1063/1.4962938
Blackman K, Perret L, Savory E (2018) Effects of the upstream-flow regime and canyon aspect ratio on non-linear interactions between a street-canyon flow and the overlying noundary layer. BLM 169:537–558. https://doi.org/10.1007/s10546-018-0378-y
Bright VB, Bloss WJ, Cai X (2013) Urban street canyons: Coupling dynamics, chemistry and within-canyon chemical processing of emissions. Atmos Environ 68:127–142. https://doi.org/10.1016/j.atmosenv.2012.10.056
Carassale L, Solari G, Tubino F (2007) Proper orthogonal decomposition in wind engineering - Part 2 Theoretical aspects and some applications. Wind and Struct 10:177–208
Castro IP, Cheng H, Reynolds R (2006) Turbulence over urban-type roughness: deductions from wind-tunnel measurements. BLM 118:109–131. https://doi.org/10.1007/s10546-005-5747-7
Chandran D, Baidya R, Monty JP, Marusic I (2017) Two-dimensional energy spectra in high-Reynolds-number turbulent boundary layers. J Fluid Mech 826:R1. https://doi.org/10.1017/jfm.2017.359
Cheng WC, Liu CH (2011) Large-eddy simulation of flow and pollutant transports in and above two-dimensional idealized street canyons. BLM 139:411–437. https://doi.org/10.1007/s10546-010-9584-y
Cheng WC, Liu CH, Leung DYC (2008) Computational formulation for the evaluation of street canyon ventilation and pollutant removal performance. Atmos Environ 42:9041–9051. https://doi.org/10.1016/j.atmosenv.2008.09.045
Cintolesi C, Pulvirenti B, Di SS (2021) Large-eddy simulations of pollutant removal enhancement from urban canyons. BLM. https://doi.org/10.1007/s10546-021-00610-8
Coceal O, Dobre A, Thomas TG, Belcher SE (2007) Structure of turbulent flow over regular arrays of cubical roughness. J Fluid Mech 589:375–409. https://doi.org/10.1017/S002211200700794X
Cui Z, Cai X, Baker C J (2004) Large-eddy simulation of turbulent flow in a street canyon. Q J R Meteorol Soc 130:1373–1394. https://doi.org/10.1256/qj.02.150
Derebail Muralidhar S, Podvin B, Mathelin L, Fraigneau Y (2019) Spatio-temporal proper orthogonal decomposition of turbulent channel flow. J Fluid Mech 864:614–639. https://doi.org/10.1017/jfm.2019.48
Di Bernardino A, Monti P, Leuzzi G, Querzoli G (2018) Pollutant fluxes in two-dimensional street canyons. Urban Climate 24:80–93. https://doi.org/10.1016/j.uclim.2018.02.002
Gromke C, Ruck B (2007) Influence of trees on the dispersion of pollutants in an urban street canyon—Experimental investigation of the flow and concentration field. Atmos Environ 41:3287–3302. https://doi.org/10.1016/J.ATMOSENV.2006.12.043
Han BS, Baik JJ, Kwak KH, Park SB (2018) Large-eddy simulation of reactive pollutant exchange at the top of a street canyon. Atmos Environ 187:381–389. https://doi.org/10.1016/j.atmosenv.2018.06.012
Inagaki A, Kanda M (2010) Organized structure of active turbulence over an array of cubes within the logarithmic layer of atmospheric flow. BLM 135:209–228. https://doi.org/10.1007/s10546-010-9477-0
Jaroslawski T, Perret L, Blackman K, Savory E (2019) The spanwise variation of roof-level turbulence in a street-canyon flow. BLM 170:373–394. https://doi.org/10.1007/s10546-018-0405-z
Jaroslawski T, Savory E, Perret L (2020) Roof-level large- and small-scale coherent structures in a street canyon flow. Environ Fluid Mech 20:739–763. https://doi.org/10.1007/s10652-019-09721-w
Kanda M (2006) Large-eddy simulations on the effects of surface geometry of building arrays on turbulent organized structures. BLM 118:151–168. https://doi.org/10.1007/s10546-005-5294-2
Kastner-Klein P, Berkowicz R, Britter R (2004) The influence of street architecture on flow and dispersion in street canyons. Meteorol Atmos Phys 87:121–131. https://doi.org/10.1007/s00703-003-0065-4
Kikumoto H, Ooka R (2012) A numerical study of air pollutant dispersion with bimolecular chemical reactions in an urban street canyon using large-eddy simulation. Atmos Environ 54:456–464. https://doi.org/10.1016/j.atmosenv.2012.02.039
Lee JH, Sung HJ, Krogstad PÅ (2011) Direct numerical simulation of the turbulent boundary layer over a cube-roughened wall. J Fluid Mech 669:397–431. https://doi.org/10.1017/S0022112010005082
Letzel MO, Krane M, Raasch S (2008) High resolution urban large-eddy simulation studies from street canyon to neighbourhood scale. Atmos Environ 42:8770–8784. https://doi.org/10.1016/j.atmosenv.2008.08.001
Li X-X, Liu C-H, Leung DYC (2008) Large-eddy simulation of flow and pollutant dispersion in high-aspect-ratio urban street canyons with wall model. BLM 129:249–268
Liu CH, Leung DYC, Barth MC (2005) On the prediction of air and pollutant exchange rates in street canyons of different aspect ratios using large-eddy simulation. Atmos Environ 39:1567–1574. https://doi.org/10.1016/j.atmosenv.2004.08.036
Liu CH, Wong CCC (2014) On the pollutant removal, dispersion, and entrainment over two-dimensional idealized street canyons. Atmos Res 135–136:128–142. https://doi.org/10.1016/j.atmosres.2013.08.006
Louka P, Belcher SE, Harrison RG (2000) Coupling between air flow in streets and the well-developed boundary layer aloft. Atmos Environ 34:2613–2621. https://doi.org/10.1016/S1352-2310(99)00477-X
Lumley JL (1970) Stochastic tools in turbulence. Academic Press, Cambridge
Lumley JL (1967) The structure of inhomogeneous turbulent flows. Atmospheric turbulence and wave propagation pp. 166–178
Michioka T, Sato A (2012) Effect of incoming turbulent structure on pollutant removal from two-dimensional street canyon. BLM 145:469–484. https://doi.org/10.1007/s10546-012-9733-6
Michioka T, Sato A, Takimoto H, Kanda M (2011) Large-eddy simulation for the mechanism of pollutant removal from a two-dimensional street canyon. BLM 138:195–213. https://doi.org/10.1007/s10546-010-9556-2
Michioka T, Takimoto H, Sato A (2014) Large-eddy simulation of pollutant removal from a three-dimensional street canyon. BLM 150:259–275. https://doi.org/10.1007/s10546-013-9870-6
Morris J, Henry M (1955) Flow in rough conduits. Trans Am Soc Civ Eng 120:373–398. https://doi.org/10.1061/TACEAT.0007206
Oke TR (2002) Boundary layer climates, 2nd edn. Routledge, London
Park SB, Baik JJ, Han BS (2021) Coherent flow structures and pollutant dispersion in a street canyon. BLM 182:363–378. https://doi.org/10.1007/S10546-021-00669-3
Park SJ, Kim JJ, Kim MJ et al (2015) Characteristics of flow and reactive pollutant dispersion in urban street canyons. Atmos Environ 108:20–31. https://doi.org/10.1016/j.atmosenv.2015.02.065
Perret L, Rivet C (2013) Dynamics of a turbulent boundary layer over cubical roughness elements: Insight from PIV measurements and POD analysis. In: Proceedings of Eighth International Symposium on Turbulence and Shear Flow Phenomena. Poitiers, France
Salim SM, Cheah SC, Chan A (2011) Numerical simulation of dispersion in urban street canyons with avenue-like tree plantings: Comparison between RANS and LES. Build Environ 46:1735–1746. https://doi.org/10.1016/J.BUILDENV.2011.01.032
Salizzoni P, Marro M, Soulhac L et al (2011) Turbulent transfer between street canyons and the overlying atmospheric boundary layer. BLM 141:393–414. https://doi.org/10.1007/s10546-011-9641-1
Schmidt OT, Colonius T (2020) Guide to spectral proper orthogonal decomposition. AIAA J 58:1023–1033. https://doi.org/10.2514/1.J058809
Schmidt OT, Towne A, Colonius T et al (2017) Wavepackets and trapped acoustic modes in a turbulent jet: Coherent structure eduction and global stability. J Fluid Mech 825:1153–1181. https://doi.org/10.1017/jfm.2017.407
Sirovich L (1987) Turbulence and the dynamics of coherent structures. I Coherent Struct Q Appl Math 45:561–571
Solazzo E, Britter RE (2007) Transfer processes in a simulated urban street canyon. BLM 124:43–60. https://doi.org/10.1007/s10546-007-9176-7
Takimoto H, Inagaki A, Kanda M et al (2013) Length-scale similarity of turbulent organized structures over surfaces with different roughness types. BLM 147:217–236. https://doi.org/10.1007/s10546-012-9790-x
Towne A, Schmidt OT, Colonius T (2018) Spectral proper orthogonal decomposition and its relationship to dynamic mode decomposition and resolvent analysis. J Fluid Mech 847:821–867. https://doi.org/10.1017/jfm.2018.283
Tutkun M, George WK (2017) Lumley decomposition of turbulent boundary layer at high Reynolds numbers. Phys Fluids. https://doi.org/10.1063/1.4974746
Uehara K, Murakami S, Oikawa S, Wakamatsu S (2000) Wind tunnel experiments on how thermal stratification affects flow in and above urban street canyons. Atmos Environ 34:1553–1562. https://doi.org/10.1016/S1352-2310(99)00410-0
Volino RJ, Schultz MP, Flack KA (2009) Turbulence structure in a boundary layer with two-dimensional roughness. J Fluid Mech 635:75–101. https://doi.org/10.1017/S0022112009007617
Welch PD (1967) The use of fast Fourier transform for the estimation of power spectra: a method based on time averaging over short, modified periodograms. IEEE Trans Audio Electroacoust 15:70–73. https://doi.org/10.1109/TAU.1967.1161901
Woll JW (1959) Homogeneous stochastic processes. Pac J Math. https://doi.org/10.2140/PJM.1959.9.293
Zhang B, Ooka R, Kikumoto H (2022) Spectral proper orthogonal decomposition analysis of turbulent flow in a two-dimensional street canyon and its role in pollutant removal. BLM 183:97–123. https://doi.org/10.1007/s10546-021-00676-4
Zhang B, Ooka R, Kikumoto H (2021) Identification of three-dimensional flow features around a square-section building model via spectral proper orthogonal decomposition. Phys Fluids 33:35151. https://doi.org/10.1063/5.0041395
Zhang YW, Gu ZL, Cheng Y, Lee SC (2011) Effect of real-time boundary wind conditions on the air flow and pollutant dispersion in an urban street canyon-Large eddy simulations. Atmos Environ 45:3352–3359. https://doi.org/10.1016/j.atmosenv.2011.03.055
Zheng X, Montazeri H, Blocken B (2021) Large-eddy simulation of pollutant dispersion in generic urban street canyons: guidelines for domain size. J Wind Eng Ind Aerodyn 211:104527. https://doi.org/10.1016/J.JWEIA.2021.104527
Zhong J, Cai XM, Bloss WJ (2016) Coupling dynamics and chemistry in the air pollution modelling of street canyons: a review. Environ Pollut 214:690–704
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Appendices
Appendix 1: An Algorithm for Spectral Proper Orthogonal Decomposition with Two-Dimensional Fourier Transformation
An algorithm for the original SPOD was proposed by Towne et al. (2018) to enhance the convergence and enforceability of the estimation of modes and mode energies from the discrete velocity data. This algorithm combines Welch’s method for spectral estimation and the snapshot POD algorithm (Sirovich 1987). This algorithm is applicable for estimating modes and mode energies in SPOD with 2DFT, although some complements need to be made.
When the Fourier transform is applied over the y-coordinate on both sides of (8), it becomes:
The right side shows the form of the original SPOD such that \({{\varvec{\uppsi}}}_{n}\) and \(\lambda_{n}\) can be directly estimated using the algorithm proposed by Towne et al. (2018), whereas the input of the algorithm (the left side of (18)) is no longer the fluctuating velocity \({\mathbf{u^{\prime}}}\) — it is a one-dimensional Fourier transformation. Therefore, one can simply discretise the y-coordinate into yj (j = 1, 2, …, Ny), and a discrete Fourier transform can be applied over the y-coordinate before the implementation of the original SPOD algorithm.
Based on this idea, the algorithm for SPOD with 2DFT is expressed as follows:
Step (i): Calculate the one-dimensional discrete Fourier transform of the fluctuation velocity \({\mathbf{u^{\prime}}}\) over the y-coordinate:
where the fast Fourier transform (FFT) algorithm is recommended.
Step (ii): Replace \({\mathbf{u^{\prime}}}\) in the original SPOD algorithm by \( {\tilde{\mathbf{u}}}^{\prime } \left( {x,z,t;\eta _{j} } \right) \), and obtain the mode \({{\varvec{\uppsi}}}_{n} \left( {x,z;\eta_{j} ,f} \right)\) and mode energies \(\lambda_{n} \left( {\eta_{j} ,f} \right)\) with discrete frequency indices.
Step (iii): Repeat (ii) over j = 1, 2, …, Ny until the modes and mode energies at all discrete wavenumbers are obtained.
The original SPOD algorithm can be processed separately in steps (ii) and (iii) for each discrete wavenumber. Therefore, parallel processing is recommended to accelerate the computation. This does not contradict the energy cascade because statistics can be performed separately on non-independent variables. In addition, because \( {\tilde{\mathbf{u}}}^{\prime } \left( {x,z,t;\eta _{j} } \right) \) consists of complex numbers, the modes are no longer symmetric with respect to f = 0. Both positive and negative f might be needed to be computed with the same wavenumber. However, because the mode energy is symmetric, and the mode and mode coefficients are conjugate symmetric with respect to the original point on the f-η plane, only two quadrants need to be considered.
Appendix 2: Cospectra in the Framework of Spectral Proper Orthogonal Decomposition with Two-Dimensional Fourier Transformation
The SPOD cospectrum between the streamwise and vertical velocity is defined as:
where \(\psi_{n}^{(u)}\) and \(\psi_{n}^{(w)}\) denote the u and w components of vector \({{\varvec{\uppsi}}}_{n}\), respectively. This cospectrum shows the planar density of the Reynolds shear stress on the f-η plane, and its integral equals the Reynolds shear stress at the spatial position (x, z); thus,
The SPOD cospectrum between the vertical velocity and concentration \(c\left( {x,y,z,t} \right)\) is defined as:
where \(\hat{c^{\prime}}\) denotes the 2DFT of \(c^{\prime}\), whose spatial and temporal discretisation should be the same as the velocity. This cospectrum shows the planar density of the vertical turbulent mass flux on the f-η plane, and its integral equals the vertical turbulent mass flux at the spatial position (x, z); therefore,
These SPOD cospectra were estimated using Welch's method (Welch 1967), wherein the time signal was first divided into blocks, and the energy of the entire signal was estimated by averaging the energy values over all blocks. When computing \(S_{n}^{(wc)}\), the block dividing and window function (if used) in the 2DFT of the concentration should be retained the same as the velocity when performing the SPOD algorithm. In the application, results in only two quadrants need to be computed because these cospectra are symmetric with respect to the original point on the f–η plane.
Appendix 3: Additional Remark in Plotting Spectra and Cospectra
In Figs. 2 and 5, the spectra and cospectra are shown by plotting the nondimensionalised value \(f\eta S\) (where S can be either spectra or cospectra, see (13)–(15)) in the linear scale against the frequency/wavenumber in the logarithmic scale. This way of plotting is considered proper for quantitatively analysing the energy distributions. Note that the integrals of spectra or cospectra over frequency and wavenumber become energy values (see (12), (21), and (23)). For any integration range \(\Omega_{f\eta }\) on the f-η plane, we can write the following transformation:
where \(\left( {\ln 10} \right)^{2}\) is a constant. This transformation indicates that plotting \(f\eta S\) against \(\log_{10} f\) and \(\log_{10} \eta\) keeps the shapes of energy distributions undistorted in the provided figures of spectra and cospectra, and this is recommended in spectral analysis.
Data Availability Statements
The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.
Rights and permissions
Springer Nature or its licensor holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Zhang, B., Ooka, R. & Kikumoto, H. Spatiotemporal Spectral Analysis of Turbulent Structures and Pollutant Removal in Two-Dimensional Street Canyon. Boundary-Layer Meteorol 185, 63–91 (2022). https://doi.org/10.1007/s10546-022-00724-7
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10546-022-00724-7