Water-Channel Estimation of Eulerian and Lagrangian Time Scales of the Turbulence in Idealized Two-Dimensional Urban Canopies

Lagrangian and Eulerian statistics are obtained from a water-channel experiment of an idealized two-dimensional urban canopy flow in neutral conditions. The objective is to quantify the Eulerian \((T^{\mathrm{E}})\) and Lagrangian \((T^{\mathrm{L}})\) time scales of the turbulence above the canopy layer as well as to investigate their dependence on the aspect ratio of the canopy, AR, as the latter is the ratio of the width (W) to the height (H) of the canyon. Experiments are also conducted for the case of flat terrain, which can be thought of as equivalent to a classical one-directional shear flow. The values found for the Eulerian time scales on flat terrain are in agreement with previous numerical results found in the literature. It is found that both the streamwise and vertical components of the Lagrangian time scale, \(T_\mathrm{u}^\mathrm{L} \) and \(T_\mathrm{w}^\mathrm{L} \), follow Raupach’s linear law within the constant-flux layer. The same holds true for \(T_\mathrm{w}^\mathrm{L} \) in both the canopies analyzed \((AR= 1\) and \(AR= 2\)) and also for \(T_\mathrm{u}^\mathrm{L} \) when \(AR = 1\). In contrast, for \(AR = 2\), \(T_\mathrm{u}^\mathrm{L} \) follows Raupach’s law only above \(z=2H\). Below that level, \(T_\mathrm{u}^\mathrm{L} \) is nearly constant with height, showing at \(z=H\) a value approximately one order of magnitude greater than that found for \(AR = 1\). It is shown that the assumption usually adopted for flat terrain, that \(\beta =T^{\mathrm{L}}/T^{\mathrm{E}}\) is proportional to the inverse of the turbulence intensity, also holds true even for the canopy flow in the constant-flux layer. In particular, \(\gamma /i_\mathrm{u} \) fits well \(\beta _\mathrm{u} =T_\mathrm{u}^\mathrm{L} /T_\mathrm{u}^\mathrm{E} \) in both the configurations by choosing \(\gamma \) to be 0.35 (here, \(i_\mathrm{u} =\sigma _\mathrm{u} / \bar{u} \), where \(\bar{u} \) and \(\sigma _\mathrm{u} \) are the mean and the root-mean-square of the streamwise velocity component, respectively). On the other hand, \(\beta _\mathrm{w} =T_\mathrm{w}^\mathrm{L} /T_\mathrm{w}^\mathrm{E} \) follows approximately \(\gamma /i_\mathrm{w} =0.65/\left( {\sigma _\mathrm{w} /\bar{u} } \right) \) for \(z > 2H\), irrespective of the AR value. The second main objective is to estimate other parameters of interest in dispersion studies, such as the eddy diffusivity of momentum \((K_\mathrm{{T}})\) and the Kolmogorov constant \((C_0)\). It is found that \(C_0\) depends appreciably on the velocity component both for the flat terrain and canopy flow, even though for the latter case it is insensitive to AR values. In all the three experimental configurations analyzed here, \(K_\mathrm{{T}}\) shows an overall linear growth with height in agreement with the linear trend predicted by Prandtl’s theory.     

THE EFFECTS OF VERTICAL VARIATION OF BASIC FLOW AND HORIZONTAL GRADIENT OF TEMPERATURE ON LINEAR AND NONLINEAR GRAVITY WAVES

By using two-dimensional dynamical equations in x-z plane with Boussinesq approximation,the effects of the second-order vertical shear of the basic flow u_zz and the horizontal gradient of temperature (M) on the gravity wave and the isolated gravity wave are discussed.The magnitudes of u_zz and M corresponding to the linear and nonlinear stabilities of the gravity waves are worked out,respectively.The results show that amplitude and width of the isolate dgravity wave are closely related to u_zz and M.It is indicated that the isolated gravity wave with a width of about 10 km can be motivated by the disturbance of sub-synoptic scale in the certain ranges of flow field shear and temperature gradient,while the motivated waves may be associated with the cold surge ahead of a cold front and the other mesoscale synoptic systems.     

Scale Properties of Anisotropic and Isotropic Turbulence in the Urban Surface Layer

The scale properties of anisotropic and isotropic turbulence in the urban surface layer are investigated. A dimensionless anisotropic tensor is introduced and the turbulent tensor anisotropic coefficient, defined as C, where \(C = 3d_{3}\,+\,1 (d_{3}\) is the minimum eigenvalue of the tensor) is used to characterize the turbulence anisotropy or isotropy. Turbulence is isotropic when \(C \approx 1\), and anisotropic when \(C \ll 1\). Three-dimensional velocity data collected using a sonic anemometer are analyzed to obtain the anisotropic characteristics of atmospheric turbulence in the urban surface layer, and the tensor anisotropic coefficient of turbulent eddies at different spatial scales calculated. The analysis shows that C is strongly dependent on atmospheric stability \(\xi = (z-z_{\mathrm{d}})/L_{{\textit{MO}}}\), where z is the measurement height, \(z_{\mathrm{d}}\) is the displacement height, and \(L_{{\textit{MO}}}\) is the Obukhov length. The turbulence at a specific scale in unstable conditions (i.e., \(\xi < 0\)) is closer to isotropic than that at the same scale under stable conditions. The maximum isotropic scale of turbulence is determined based on the characteristics of the power spectrum in three directions. Turbulence does not behave isotropically when the eddy scale is greater than the maximum isotropic scale, whereas it is horizontally isotropic at relatively large scales. The maximum isotropic scale of turbulence is compared to the outer scale of temperature, which is obtained by fitting the temperature fluctuation spectrum using the von Karman turbulent model. The results show that the outer scale of temperature is greater than the maximum isotropic scale of turbulence.     

Optimizing the Determination of Roughness Parameters for Model Urban Canopies

We present an objective optimization procedure to determine the roughness parameters for very rough boundary-layer flow over model urban canopies. For neutral stratification the mean velocity profile above a model urban canopy is described by the logarithmic law together with the set of roughness parameters of displacement height d, roughness length \(z_0\), and friction velocity \(u_*\). Traditionally, values of these roughness parameters are obtained by fitting the logarithmic law through (all) the data points comprising the velocity profile. The new procedure generates unique velocity profiles from subsets or combinations of the data points of the original velocity profile, after which all possible profiles are examined. Each of the generated profiles is fitted to the logarithmic law for a sequence of values of d, with the representative value of d obtained from the minima of the summed least-squares errors for all the generated profiles. The representative values of \(z_0\) and \(u_*\) are identified by the peak in the bivariate histogram of \(z_0\) and \(u_*\). The methodology has been verified against laboratory datasets of flow above model urban canopies.