In this study, macro-rough flows over beds with different permeability values are simulated using the large-eddy simulation, and the results are analysed by applying the double-averaging (DA) methodology. Spheres of different sizes and arrangements were used to form the beds, which are deemed to be permeable granular beds. The influence of bed permeability on the turbulence dynamics and structure is investigated. It was observed that the scales of the spanwise vortical structures over more permeable beds are larger than those over less permeable beds. This is attributed to large-scale spanwise-alternate strips of varying Reynolds shear stress (RSS), emerging from the surface of macro-rough elements for the permeable beds. The DA stress balance suggests that the time-averaged spanwise vortical structure leads to a damping in DA RSS and an unusual peak of the form-induced stress in the main flow. In the streamwise direction, both large turbulent structures that originate from the Kelvin–Helmholtz-type instability and small turbulent structures that are associated with the turbulent transport across the gaps of the roughness elements are more prevalent over highly permeable beds. Near the bed, the relative magnitude of turbulent events shows a transition from a ejections-dominating to sweeps-dominating zone with vertical distance. Further, several hydrodynamic characteristics normalized by inner scales (kinematic viscosity to shear velocity ratio) show a greater dependency on permeability Reynolds number than those normalized by sediment size. The study provides an insight into the mechanism of mass transfer near the fluid–particle interface, which is vital to benthic and aquatic ecology.

In this study, we present a three-dimensional numerical model for the interaction of flow with submerged flexible vegetation, based on a large-eddy simulation and the immersed boundary method. The model innovatively realises the interaction between the flow and highly flexible vegetation with clustered leaves. Besides being a three-dimensional model of motion with full degrees of freedom, this study improves the consideration of the motion of the vegetation in all directions, and in addition the energy and momentum transfer in the spanwise direction. Furthermore, we perform a flume experiment for the flow with submerged flexible vegetation, the results of which are used to validate the simulation effects of the numerical model. It is found that the numerical model can effectively simulate the velocity profiles and the movement of vegetation induced by the flow. Using the model to analyse the flow–vegetation interaction, we find that the movement of vegetation is closely related to the flow velocity. As the flow velocity increases, both the offset angle and the vegetation swaying amplitude increase. Compared to vertical rigid vegetation, the tilting of flexible vegetation does not significantly change the velocity difference and the magnitude of the turbulent kinetic energy between the inside and the outside of the vegetation canopy, but it does weaken the disturbance to flow, thus reducing the resistance to flow. However, the swaying of vegetation dose significantly increase the velocity difference between the inside and the outside of the canopy. It forms Kelvin–Helmholtz hairpin vortices intensifying the turbulence production, and enhancing the disturbance and resistance to flow.