Investigation of turbulence and interfacial exchange features of the gap area within the fully developed Shallow-Submerged canopy flow

The flow patterns within a longitudinal gap area formed by discontinuous distributions of submerged canopy, as well as the momentum and mass exchange characteristics between the gap area and the overlying free-flow, were studied using high-resolution Large Eddy Simulation (LES). The gap area is located within the fully developed region of submerged canopy flow. The simulations considered four aspect ratios of the gap area, L/h (ranging from 1 to 4), where L and h represent the span of the gap area and the height of the canopy, respectively, and two canopy densities, ϕ = 0.08 and 0.15. Results indicate that recirculation vortices appear only at ϕ = 0.15 and penetrate the downstream canopy patches to varying extents, with the distance of invasion decreasing as L/h increases. Influenced by the recirculation vortices, the bed shear stress in the downstream part of the gap area and the initial section of the downstream canopy patch is significantly increased compared to the fully developed region of the upstream canopy patch. Within the parameter range covered, vertical penetration of the mixing layer into the gap area always occurs, consistently falling into the shear layer growth regime, with significant enhancement of turbulence within the gap area even at L/h = 1. The high momentum entrainment elevates longitudinal velocity and turbulence intensity in the upper corner regions of the downstream canopy patch's leading edge, combined with localized increases in bed shear stress, potentially destabilizing plants at the forefront of the downstream canopy patch. Although the total momentum exchange across the interface between the gap area and the upper free-flow layer is approximately independent of L/h, larger canopy densities lead to stronger momentum transport, and turbulent transport always dominates. Mass exchange generally increases with L/h, with more efficient vertical mass exchange in ϕ = 0.15 cases at L/h = 3 and 4.