Viscous streaming is an efficient mechanism to exploit inertia at the microscale for flow control. While streaming from rigid features has been thoroughly investigated, when body compliance is involved, as in biological settings, little is known. Here, we investigate body elasticity effects on streaming in the minimal case of an immersed soft cylinder. Our study reveals an additional streaming process, available even in Stokes flows. Paving the way for advanced forms of flow manipulation, we illustrate how gained insights may translate to complex geometries beyond circular cylinders.

Compressible mixing layer instabilities are of importance to a wide range of environmental and industrial applications. Past studies have focused on either ideal-gas or real-fluid thermodynamic regimes of single-species mixing layers. However, mixing layers of binary mixtures at supercritical conditions, commonly encountered in fuel injection systems, introduce additional complexities due to the added compositional degree-of-freedom. Moreover, the effect of strong variations in thermodynamic response functions across the Widom line on the binary mixing layer stability remains poorly understood. Thus, the objective of this study is to examine the coupling between the hydrodynamic instability and the real-fluid thermodynamics across the Widom line and its effects on the overall binary mixing layer dynamics. To this end, we develop a linear stability analysis of the full binary-species compressible transport equations coupled with the PC-SAFT equation of state. Analysis shows the existence of a novel instability mechanism that arises from juxtapositioning of the Widom-line transition and the hydrodynamic inflection point. This novel thermodynamically induced instability mechanism has the net effect of destabilizing the binary mixing layer at lowering supercritical conditions towards the critical pressure point. This is in contrast to previous stability analyses of supercritical single-species mixing layers, where increasing pressure destabilizes the flow due to its effect on reducing the density stratification. The discovered thermodynamically induced instability mechanism of binary mixing flows highlights the need for an extension of classical instability criteria to incorporate the effect of strong variations in the thermodynamic response functions across the Widom line on mixing layer instability.

The effect of damping the longest streaks in wall-bounded turbulence is explored using numerical experiments. It is found that long streaks are not required for the self-sustenance of the bursting process, which is relatively little affected by their absence. In particular, there are turbulence states in which the fluctuations of the streamwise velocity have approximately the same length as the bursts, and are thus presumably associated with the bursts themselves, while the burst structure is essentially indistinguishable from flows in which longer velocity fluctuations are present. This suggests that the long streaks found in unmodified flows may be by-products, rather than active parts of the process.

A Janus sphere with a stick–slip pattern can behave quite differently in its hydrodynamics compared with a no-slip or uniform-slip sphere. Here, using the Lorentz reciprocal theorem in conjunction with surface harmonic expansion, we rigorously derive the extended Faxén formula for the stresslet of a weakly stick–slip Janus sphere, capable of describing the anisotropic nature of the stresslet with an arbitrary axisymmetric stick–slip pattern in an arbitrary background flow. We find that slip anisotropy not only causes a variety of additional contributions to the stresslet, but also naturally renders a stresslet–rotation coupling that may turn a suspension of couple-free stick–slip Janus spheres into a dipolar one under the actions of an external couple. Moreover, to correctly account for the impacts of slip anisotropy on the stresslet, it is necessary to include at least the first four surface harmonic contributions. As a result, the anisotropies of both the stresslet and torque on the sphere in a linear flow field are purely reflected by a symmetric quadrupole and hexadecapole. These hydrodynamic quantities can be further mediated by an antisymmetric dipole and octupole due to the gradients of the imposed strain field. The average bulk stress and effective viscosity for a suspension of stick–slip spheres are also determined, showing characteristics quite distinct from those of a suspension of near spheres. If the spheres possess permanent dipole moments, in particular, additional stresslets and couplets can be generated by an applied external couple on each sphere and added into the bulk stress, accompanied by non-Jeffery orientational orbits of such dipolar stick–slip spheres. In addition to the above, the extended Faxén stresslet and torque relations found in this work will also provide the formulae needed for tackling problems involving hydrodynamically interacting stick–slip spheres on which small slip anisotropy may have profound impacts.

As opposed to the case of shallow mixing layers forming between parallel streams of unequal velocities and equal densities, the spatial development of the mixing interface (MI) between two parallel streams of unequal velocities and sufficiently large density contrast is controlled by the formation of a spatially developing, lock-exchange-like flow in transverse planes. Buoyancy effects are driven by the presence of a vertical density interface near the splitter plate. Eddy-resolving numerical simulations conducted in a wide and very long channel are used to investigate the mean flow structure and the effects of the lock-exchange-like flow and the associated coherent structures on mixing and the capacity of the flow to entrain sediment from the channel bed. When the two streams have unequal densities, quasi-two-dimensional Kelvin–Helmholtz (KH) vortices still form near the MI origin, but their coherence is lost over a much shorter distance compared with the case with no density contrast, and their cores are severely stretched in the transverse direction. A main cell of recirculating (cross-) flow forms in the mean flow, in between the fronts of the near-bed intrusion of heavier fluid and the free-surface intrusion of lighter fluid. The instantaneous flow fields contain streamwise-oriented vortical cells along the interface separating the regions containing heavier and lighter fluids. These vortical cells play an important role in enhancing mixing, similar to the KH billows forming in a classical lock-exchange flow. The regions of highest turbulence amplification are situated next to the boundaries of the main recirculating cell. Once the oscillating fronts of the intrusions get sufficiently close to the channel sidewalls, streamwise cells of secondary flow form near the channel boundaries. For cases with a strong density contrast, the free-surface mixing pattern is not a good indicator of mixing between the two streams. For the same flow velocities in the incoming channels, the two streams mix faster with increasing density difference between the two streams. This is because, as opposed to the KH vortices, coherent structures induced by the formation of the lock-exchange-like flow maintain their coherence and capacity to induce mixing at very large distances from the splitter plate. Meanwhile, the redistribution of the streamwise momentum leading to uniform, fully developed flow over the whole cross-section is delayed by buoyancy effects. Away from the splitter plate, the region with the highest sediment entrainment potential is situated next to the edge of the main recirculation cell on the high-speed side of the channel. For the same flow velocities in the incoming channels, the sediment entrainment capacity of the flow is much larger in the simulations conducted with density contrast between the incoming streams and peaks for the case when the faster stream contains the denser fluid.

Transport phenomena in high Reynolds number wall-bounded stratified flows are dominated by the interplay between the turbulence structures generated at the wall and the buoyancy-induced large-scale waves populating the channel core. In this study, we want to investigate the flow physics of wall-bounded stratified turbulence at relatively high shear Reynolds number $Re_\tau$ and for mild to moderate stratification level – quantified here by the shear Richardson number varying in the range $0\leqslant Ri_{\tau } \leqslant 300$. By increasing stratification, active turbulence is sustained only in the near-wall region, whereas intermittent turbulence, modulated by the presence of non-turbulent wavy structures (internal gravity waves), is observed at the channel core. In such conditions, the wall-normal transport of momentum and heat is considerably reduced compared with the case of non-stratified turbulence. A careful characterization of the flow-field statistics shows that, despite temperature and wall-normal velocity fluctuations being very large at the channel centre, the mean value of their product – the buoyancy flux – vanishes for $Ri_{\tau } \geqslant 200$. We show that this behaviour is due to the presence of a $\sim {\rm \pi}/2$ phase delay between the temperature and the wall-normal velocity signals: when wall-normal velocity fluctuations are large (in magnitude), temperature fluctuations are almost zero, and vice versa. This constitutes a blockage effect to the wall-normal exchange of energy. In addition, we show that the friction factor scales as $C_f \sim Ri_{\tau }^{-1/3}$, and we propose a new scaling for the Nusselt number, $Nu \cdot Re_{\tau }^{-2/3} \sim Ri_{\tau }^{-1/3}$. These scaling laws, which seem to be robust over the explored range of parameters, complement and extend previous experimental and numerical data, and are expected to help the development of improved models and parametrizations of stratified flows at large $Re_{\tau }$.

The supersonic wake of a circular cylinder in Mach 3 flow was studied through high-speed, focussing schlieren photography. The mean and unsteady behaviour of the separated shear layers, the reattachment process, the recompression wave and the early wake are analysed, and discussed in detail. The fluctuations in the wake are stronger and more coherent than those within the approaching shear layers and the recirculation region. The recompression of the shear layers energises the finer scales in the flow which leads to a departure from a $-$1 spectral roll-off observed in the schlieren spectra further upstream. The recompression wave exhibits low-frequency unsteadiness and a ripple-type motion which occurs as it is perturbed by shocklets radiating from the coherent structures in the wake. The wake consists of coherent disturbances with the same characteristic frequency as that for an incompressible flow over a cylinder; however, this instability is suppressed as the wake accelerates, presumably due to increasing compressibility. The primary instability of the wake flow has a characteristic frequency nearly twice that of its incompressible counterpart and it is shown to be driven by the presence of aeroacoustic resonance in the wake. It is also shown that the resonance, which leads to the formation of broadband standing waves in the wake, is the result of an interaction between the wake instabilities and upstream propagating acoustic waves in the wake. The acoustic waves originate upstream of the reattachment region and are believed to be generated by the unsteady separation on the cylinder surface.

A numerical study on the response of a two-dimensional bluff body wake subjected to harmonic forcing imposed by two rear pitching flaps is performed. The wake is generated by a rectangle at a height-based Reynolds number $Re=100$, characterised by laminar vortex shedding. Two forcing strategies are examined corresponding to in-phase ‘snaking’ and out-of-phase ‘clapping.’ The effects of the bluff body aspect ratio ($AR=1,2,4$), flapping frequency, flapping amplitude, flap length and Reynolds number are investigated. For the snaking motion, a strong fundamental resonance of the root mean square (r.m.s.) drag is observed when the wake is forced near the vortex shedding frequency. For the clapping motion, a weak subharmonic resonance is observed when the forcing is applied near twice the vortex shedding frequency resulting in an increase of the lift r.m.s. whereas the drag r.m.s. remains unaffected. Both resonances intensify the vortex shedding and a concomitant mean drag increase is observed for the snaking motion. Forcing away from the resonant regimes, both motions result in considerable drag reduction through wake symmetrisation and propulsion mechanisms. The formation of two vortex dipoles per oscillation period due to the flapping motion, which weaken the natural vortex shedding, has been identified as the main symmetrisation mechanism. A single scaling parameter is proposed to collapse the mean drag reduction of the forced flow for both motions over a wide range of flapping frequencies, amplitudes and flap lengths. Finally, the assessment of the performance of the forcing strategies has revealed that clapping is more effective than snaking.

This study presents direct numerical simulations of turbulent Rayleigh–Bénard convection in non-colloidal suspensions, with special focus on the heat transfer modifications in the flow. Adopting a Rayleigh number of $10^8$ and Prandtl number of 7, parametric investigations of the particle volume fraction $0\leq \varPhi \leq 40\,\%$ and particle diameter $1/20\leq d^*_p\leq 1/10$ with respect to the cavity height, are carried out. The particles are neutrally buoyant, rigid spheres with physical properties that match the fluid phase. Up to $\varPhi =25\,\%$, the Nusselt number increases weakly but steadily, mainly due to the increased thermal agitation that overcomes the decreased kinetic energy of the flow. Beyond $\varPhi =30\,\%$, the Nusselt number exhibits a substantial drop, down to approximately 1/3 of the single-phase value. This decrease is attributed to the dense particle layering in the near-wall region, confirmed by the time-averaged local volume fraction. The dense particle layer reduces the convection in the near-wall region and negates the formation of any coherent structures within one particle diameter from the wall. Significant differences between $\varPhi \leq 30\,\%$ and 40 % are observed in all statistical quantities, including heat transfer and turbulent kinetic energy budgets, and two-point correlations. Special attention is also given to the role of particle rotation, which is shown to contribute to maintaining high heat transfer rates in moderate volume fractions. Furthermore, decreasing the particle size promotes the particle layering next to the wall, inducing a similar heat transfer reduction as in the highest particle volume fraction case.

Wake evolution of an oscillating foil with combined heaving and pitching motion is evaluated numerically for a range of phase offsets ($\phi$), chord-based Strouhal numbers ($St_c$) and Reynolds numbers ($Re$). The increase in $\phi$ from $90^\circ$ to $180^\circ$ at a given $St_c$ and $Re$ coincides with a transition of pitch- to heave-dominated kinematics that further reveals novel transitions in wake topology characterized by bifurcated vortex streets. At $Re= 1000$, each of the dual streets constitutes a dipole-like paired configuration of counter-rotating coherent structures that resemble qualitatively the formation of $2P$ mode. A new mathematical relation between the relative circulation of coherent dipole-like paired structures and kinematic parameters is proposed, including heave-based ($St_h$), pitch-based ($St_{\theta }$) and combined motion ($St_A$) Strouhal numbers, as well as $\phi$. This model can predict accurately the wake transition towards $2P$ mode characterized by a bifurcation, at low $Re= 1000$. At $Re= 4000$, however, the relationship was inaccurate in predicting the wake transition. A shear splitting process is observed at $Re= 4000$, which leads to the formation of reverse Bénard–von Kármán mode in conjunction with $2P$ mode. Increasing $\phi$ further depicts a consistent prolongation of the splitting process, which coincides with a unique transition in terms of absence and reappearance of bifurcated dipole-like pairs at $\phi = 120^\circ$ and $180^\circ$, respectively. Changes in the spatial arrangement of $2P$ pairs observed consistently for oscillating foils with the combined motion constitute a novel wake transition that becomes more dominant at higher Reynolds numbers.

This paper describes the mechanism underpinning modal staging behaviour in screeching jets. An upstream-propagating subsonic guided-jet mode is shown to be active in all stages of screech. Axial variation of shock-cell spacing manifests in the spectral domain as a series of suboptimal peaks. It is demonstrated that the guided-jet mode may be energized by interactions of the Kelvin–Helmholtz wavepacket with not only the primary shock wavenumber peak, but also suboptimal peaks; interaction with suboptimals is shown to be responsible for closing the resonance loop in multiple stages of jet screech. A consideration of the full spectral representation of the shocks reconciles several of the classical models and results for jet screech that had heretofore been paradoxical. It is demonstrated that there are multiple standing waves present in the near field of screeching jets, corresponding to the superposition of the various waves active in these jets. Multimodal behaviour is explored for jets in a range of conditions, demonstrating that multiple peaks in the frequency spectra can be due to either changes in which peak of the shock spectra the Kelvin–Helmholtz wavepacket is interacting with, or a change in azimuthal mode, or both. The absence of modal staging in high-aspect-ratio non-axisymmetric jets is also explained in the context of the aforementioned mechanism. The paper closes with a new proposed theory for frequency selection in screeching jets, based on the observation that these triadic interactions appear to underpin selection of the guided-jet mode wavelength in all measured cases.

The two-dimensional acoustic wave equation for inviscid compressible boundary layer flows, i.e. the Pridmore-Brown equation with an exponential velocity profile for homentropic flows, is studied for the reflection and over-reflection of acoustic waves based on the exact solution in terms of the confluent Heun function. The reflection coefficient $R$, which is the ratio of the amplitude of the reflected to that of the incoming acoustic wave, is determined as a function of the streamwise wavenumber $\alpha$, the Mach number $M$ and the incident angle $\phi$ of the acoustic waves. Over-reflection refers to $R>1$, i.e. the reflected wave has a larger amplitude than the incident wave. We prove that, in the supersonic context, energy is always transferred from the base flow to the reflected wave, i.e. $R<1$ does not exist. Meanwhile, this fact is intimately linked to the critical layer. We show that the presence of the critical layer leads to an optimal energy exchange from the base flow into the acoustic wave, i.e. the critical layer ensures $R>1$. In our analysis, we observe a special phenomenon, resonant over-reflection, which is proven to be closely related to resonant frequencies $\omega _r$ of unstable modes of the temporal stability of the base flow. At resonant frequencies of the first unstable mode, the over-reflection coefficient exhibits an unusual peak in an extremely narrow frequency interval. The maximum values of these peaks are largely synchronized with the variation of the growth rate $\omega _i$ of the unstable modes. In addition, resonant over-reflection appears also at resonant frequencies of other higher unstable modes, but their peaks of the over-reflection coefficient are always smaller than that induced by the first unstable mode.

A spatially developing turbulent boundary layer subject to a space- and time-dependent pressure gradient is analysed via large-eddy simulation. The unsteadiness is prescribed by imposing an oscillating suction–blowing velocity profile at the top boundary of the computational domain. The alternating favourable and adverse pressure gradients cause the flow to separate and reattach to the wall periodically. A range of reduced frequencies $k$ was investigated, spanning from a very rapid flutter-like motion to a slow, quasi-steady flapping. The Reynolds number based on the boundary-layer displacement thickness $\delta _o^{*}$ at the inflow plane is $Re_*=1000$. Both time- and phase-averaged fields are analysed and results are compared with steady conditions. The reduced frequency $k$ has a significant effect on the transient flow-separation process. For high $k$ the separation bubble does not grow as thick as in the corresponding steady case, but the length of the bubble remains comparable; hysteresis is observed in the near-wall region. As $k$ is reduced, a threshold is met at which the separation bubble grows in the wall-normal direction. However, the length of the bubble is significantly reduced again when compared with the steady case. At this frequency, the region of slow-moving fluid generated by the flow reversal is advected downstream, causing a decorrelation between the forcing (the imposed free-stream velocity) and the velocity and pressure downstream of the separation bubble. Moreover, hysteresis effects are shifted away from the wall. At the lowest frequency a quasi-steady solution is approached; however, transient effects are still present in the backflow region.

Spectral degeneracies where eigenvalues and eigenvectors simultaneously coalesce, also known as exceptional points, are a natural consequence of the strong non-normality of the Orr–Sommerfeld operator describing the evolution of infinitesimal disturbances in parallel shear flows. While the resonances associated with these points give rise to algebraic growth, the development of non-modal stability theory exploiting specific perturbation structures with much larger potential for transient energy growth has led to waning interest in spectral degeneracies. The appearance of subharmonic eigenvalue orbits, recently discovered in the periodic spectrum of pulsating Poiseuille flow, can be traced back to the coalescence of eigenvalues at exceptional points. We present a thorough analysis of the spectral properties of the linear operator to identify exceptional points and accurately map the prevalence of subharmonic eigenvalue orbits for a large range of pulsation amplitudes and frequencies. This information is then combined with solutions of the linear initial value problem to analyse the impact of the appearance of these orbits on the temporal evolution of linear disturbances in pulsating Poiseuille flow. The periodic amplification phases are shown to be heralded by repeated non-normal growth bursts that are intensified by the formation of subharmonic orbits involving the leading eigenvalues. These bursts are associated with the change of alignment of the perturbation from the decaying towards the amplified branch of the subharmonic eigenvalue orbits in a so-called branch transition process.

The propagation speed, shape and stability of the rim generated by a liquid-curtain breakup are studied. In the experiment, a liquid curtain surrounded by a slot die, edge guides and the surface of a roller breaks at the contact point between the edge guide and roller in a low-Weber-number range, and the rim propagates in the horizontal direction. Except for the initial time, the rim is almost straight and has a nearly constant propagation speed. For an Ohnesorge number much smaller than 1, unevenness occurs on the rim and the droplets separate from it. When the Ohnesorge number is of the order of unity, the rim becomes convex vertically downward, and the liquid lump flows down. The shape, propagation speed and surface stability of the rim are discussed by analysing the equation proposed by Entov & Yarin (J. Fluid Mech., vol. 140, 1984, pp. 91–111). It is shown that the volume flow rate condition at the slot die exit is important to explain the propagation of the rim. Additionally, in the initial stage of the curtain breakup, the Plateau–Rayleigh instability causes unevenness on the rim surface, and after the rim reaches the slot die exit, the Rayleigh–Taylor instability generates a liquid lump on the rim, which grows into droplets when the Ohnesorge number is much less than 1.

The impact of a liquid drop on a solid surface involves many intertwined physical effects, and is influenced by drop velocity, surface tension, ambient pressure and liquid viscosity, among others. Experiments by Kolinski et al. (Phys. Rev. Lett., vol. 112, no. 13, 2014b, p. 134501) show that the liquid–air interface begins to deviate away from the solid surface even before contact. They found that the lift-off of the interface starts at a critical time that scales with the square root of the kinematic viscosity of the liquid. To understand this, we study the approach of a liquid drop towards a solid surface in the presence of an intervening gas layer. We take a numerical approach to solve the Navier–Stokes equations for the liquid, coupled to the compressible lubrication equations for the gas, in two dimensions. With this approach, we recover the experimentally captured early time effect of liquid viscosity on the drop impact, but our results show that lift-off time and liquid kinematic viscosity have a more complex dependence than the square-root scaling relationship. We also predict the effect of interfacial tension at the liquid–gas interface on the drop impact, showing that it mediates the lift-off behaviour.

Heterogeneous roughness in the form of streamwise aligned strips is known to generate large scale secondary motions under turbulent flow conditions that can induce the intriguing feature of larger flow rates above rough than smooth surface parts. The hydrodynamical definition of a surface roughness includes a large scale separation between the roughness height and the boundary layer thickness which is directly related to the fact that the drag of a laminar flow is not altered by the presence of roughness. Existing simplified approaches for direct numerical simulation of roughness strips do not fulfil this requirement of an unmodified laminar base flow compared with a smooth wall reference. It is shown that disturbances induced in a modified laminar base flow can trigger large-scale motions with resemblance to turbulent secondary flow. We propose a simple roughness model that allows us to capture the particular features of turbulent secondary flow without impacting the laminar base flow. The roughness model is based on the prescription of a spanwise slip length, a quantity that can directly be translated into the Hama roughness function for a homogeneous rough surface. The heterogeneous application of the slip-length boundary condition results in very good agreement with existing experimental data in terms of the secondary flow topology. In addition, the proposed modelling approach allows us to quantitatively evaluate the drag increasing contribution of the secondary flow. Both the secondary flow itself and the related drag increase reveal a very small dependence on the gradient of the transition between rough and smooth surface parts only. Interestingly, the observed drag increase due to secondary flows above the modelled roughness is significantly smaller than the one previously reported for roughness resolving simulations. We hypothesise that this difference arises from the fact that roughness resolving simulations cannot truly fulfil the requirement of large scale separation.

We consider high-order strongly nonlinear long wave models expanded in a single small parameter measuring the ratio of the water depth to the characteristic wavelength. By examining its dispersion relation, the high-order system for the bottom velocity is found stable to all disturbances at any order of approximation. On the other hand, systems for other velocities can be unstable and even ill-posed, as signified by the unbounded maximum growth. Under the steady assumption, a new third-order solitary wave solution of the Euler equations is obtained using the high-order strongly nonlinear system and is expanded in an amplitude parameter, which is different from that used in weakly nonlinear theory. The third-order solution is shown to well describe various physical quantities induced by a finite-amplitude solitary wave, including the wave profile, horizontal velocity profile, particle velocity at the crest and bottom pressure. For numerical computations, the first- and second-order strongly nonlinear systems for the bottom velocity are considered. It is shown that finite difference schemes are unstable due to truncation errors introduced in approximating high-order spatial derivatives and, therefore, a more accurate spatial discretization scheme is necessary. Using a pseudo-spectral method based on finite Fourier series combined with an iterative scheme for the inversion of a non-local operator, the strongly nonlinear systems are solved numerically for the propagation of a single solitary wave and the head-on collision of two counter-propagating solitary waves of finite amplitudes, and the results are compared with previous laboratory measurements.

Recently, the existence of so-called granular bubbles and droplets has been demonstrated experimentally. Granular bubbles and droplets are clusters of particles that respectively rise and sink if submerged in an aerated and vibrated bed of another granular material of different size and/or density. However, currently, there is no model that explains the coherent motion of these clusters and predicts the transition between a rising and sinking motion. Here, we propose an analytical model predicting accurately the neutral buoyancy limit of a granular bubble/droplet. This model allows the compilation of a regime map identifying five distinct regimes of granular bubble/droplet motion.

We have performed high-fidelity simulations of turbulent open-channel flows over submerged rigid canopies made of cylindrical filaments of fixed length $l=0.25H$ ($H$ being the domain depth) mounted on the wall with angle of inclination $\theta$. The inclination is the free parameter that sets the density of the canopy by varying its frontal area. The density of the canopy, based on the solidity parameter $\lambda$, is a widely accepted criterion defining the ongoing canopy flow regime, with low values ($\lambda \ll 0.1$) indicating the sparse regime, and higher values ($\lambda > 0.1$) the dense regime. All the numerical predictions have been obtained considering the same nominal bulk Reynolds number (i.e. $Re_b=U_b H/\nu = 6000$). We consider nine configurations of canopies, with $\theta$ varying symmetrically around $0^{\circ }$ in the range $\theta \in [\pm 78.5^{\circ }$], where positive angles define canopies inclined in the flow direction (with the grain) and $\theta =0^{\circ }$ corresponds to the wall-normally mounted canopy. The study compares canopies with identical solidity obtained inclining the filaments in opposite angles, and assesses the efficacy of the solidity as a representative parameter. It is found that when the canopy is inclined, the actual flow regime differs substantially from the one of a straight canopy that shares the same solidity, indicating that criteria based solely on this parameter are not robust. Finally, a new phenomenological model describing the interaction between the coherent structures populating the canopy region and the outer flow is given.