AbstractGood representation of turbulence in urban canopy models is necessary for accurate prediction of momentum and scalar distribution in and above urban canopies. To develop and improve turbulence closure schemes for one-dimensional multi-layer urban canopy models, turbulence characteristics are investigated here by analyzing existing large-eddy simulation and direct numerical simulation data. A range of geometries and flow regimes are analyzed that span packing densities of 0.0625 to 0.44, different building array configurations (cubes and cuboids, aligned and staggered arrays, and variable building height), and different incident wind directions ($$0^\circ $$
0
∘
and $$45^\circ $$
45
∘
with regards to the building face). Momentum mixing-length profiles share similar characteristics across the range of geometries, making a first-order momentum mixing-length turbulence closure a promising approach. In vegetation canopies turbulence is dominated by mixing-layer eddies of a scale determined by the canopy-top shear length scale. No relationship was found between the depth-averaged momentum mixing length within the canopy and the canopy-top shear length scale in the present study. By careful specification of the intrinsic averaging operator in the canopy, an often-overlooked term that accounts for changes in plan area density with height is included in a first-order momentum mixing-length turbulence closure model. For an array of variable-height buildings, its omission leads to velocity overestimation of up to $$17\%$$
17
%
. Additionally, we observe that the von Kármán coefficient varies between 0.20 and 0.51 across simulations, which is the first time such a range of values has been documented. When driving flow is oblique to the building faces, the ratio of dispersive to turbulent momentum flux is larger than unity in the lower half of the canopy, and wake production becomes significant compared to shear production of turbulent momentum flux. It is probable that dispersive momentum fluxes are more significant than previously thought in real urban settings, where the wind direction is almost always oblique.
HYPOTHESIS OF PATH INTEGRAL DUALITY: APPLICATIONS TO QED
