Search

Prandtl's secondary flows of the second kind generated by laterally-varying roughness are studied in this work using the linearised Reynolds-Averaged Navier-Stokes approach recently proposed. The momentum equations are coupled to the SA transport model for the turbulent viscosity and the surface roughness is captured by adapting established strategies for homogeneous rough walls to spanwise inhomogeneous surfaces. Linearisation of the governing equations and of the roughness model yields a linear framework that elucidates the linear mechanisms that play a role in the generation of secondary flows. In addition, the framework allows exploring the parameter space associated with heterogeneous rough surfaces. Channel flow is considered, with high and low roughness strips arranged symmetrically on the two channel walls. The strip width is varied systematically to ascertain the role played by the length scale characterising the heterogeneity on the size and intensity of secondary flows. It is found that linear mechanisms, i.e. whereby secondary flows are interpreted as the output response of the turbulent mean flow subjected to a forcing localised at the wall, may be sufficient to explain such a role. In fact the model predicts that secondary flows are most intense when the strip width is about 0.7 times the half-channel height, in excellent agreement with available data. Further, a unified framework to analyse combinations of heterogeneous roughness properties and laterally-varying topographies, common in engineering applications, is discussed. This framework yields two separate secondary-flow inducing source mechanisms, i.e. the lateral variation of the virtual origin from which the turbulent structure is allowed to develop and the lateral variation of the streamwise velocity slip, capturing the acceleration/deceleration perceived by the bulk flow over troughs and crests of the topography.

Boundary layer flow over a realistic porous wall might contain both the effects of wall-permeability and wall-roughness. These two effects are typically examined in the context of a rough-wall flow, i.e., by defining a ``roughness'' length or equivalent to capture the effect of the surface on momentum deficit/drag. In this work, we examine the hypothesis of Esteban et al. (2022), that a turbulent boundary layer over a porous wall could be modelled as a superposition of the roughness effects on the permeability effects by using independently obtained information on permeability and roughness. We carry out wind tunnel experiments at high Reynolds number ($14400 \leq Re_{\tau} \leq 33100$) on various combinations of porous walls where different roughnesses are overlaid over a given permeable wall. Measurements are also conducted on the permeable wall as well as the rough walls independently to obtain the corresponding lengthscales. Analysis of mean flow data across all these measurements suggests that an empirical formulation can be obtained where the momentum deficit ($\Delta U^+$) is modelled as a combination of independently obtained roughness and permeability lengthscales. This formulation assumes the presence of outer-layer similarity across these different surfaces, which is shown to be valid at high Reynolds numbers. Finally, this decoupling approach is equivalent to the area-weighted power-mean of the respective permeability and roughness lengthscales, consistent with the approach recently suggested by Hutchins et al. (2023) to capture the effects of heterogeneous rough surfaces., Comment: Under review for publication in JFM Rapids

In this work, we describe the impact of aspect ratio (AR) on the performance of optimally phased, identical flapping flippers in a tandem configuration. Three-dimensional simulations are performed for seven sets of single and tandem finite foils at a moderate Reynolds number, with thrust producing, heave-to-pitch coupled kinematics. Increasing slenderness (or aspect ratio - AR) is found to improve thrust coefficients and thrust augmentation but the benefits level off towards higher values of AR. On the other hand, the propulsive efficiency shows no significant change with increasing AR, while the hind foil outperforms the single by a small margin. Further analysis of the spanwise development and propagation of vortical structures allows us to gain some insights on the mechanisms of these wake interactions and provide valuable information for the design of novel biomimetic propulsion systems.

Flapping flight and swimming are increasingly studied both due to their intrinsic scientific richness and their applicability to novel robotic systems. Strip theory is often applied to flapping wings, but such modeling is only rigorously applicable in the limit of infinite aspect ratio (AR) where the geometry and kinematics are effectively uniform. This work compares the flow features and forces of strip theory and three-dimensional flapping foils, maintaining similitude in the rolling and twisting kinematics while varying the foil AR. We find the key influence of finite AR and spanwise varying kinematics is the generation of a time-periodic spanwise flow which stabilizes the vortex structures and enhances the dynamics at the foil root. An aspect-ratio correction for flapping foils is developed analogous to Prandtl finite wing theory, enabling future use of strip theory in analysis and design of finite aspect ratio flapping foils., Comment: 11 pages, 12 figures

Symmetric flapping foils are known to produce deflected jets at high frequency-amplitude combinations even at a zero mean angle of attack. This reduces the frequency range of useful propulsive configurations without side force. In this study, we numerically analyse the interaction of these deflected jets for tandem flapping foils undergoing coupled heave to pitch motion in a two dimensional domain. The impact of the flapping Strouhal number, foil spacing and phasing on wake interaction is investigated. Our primary finding is that the back foil is capable of cancelling the wake deflection and mean side force of the front foil, even when located up to 5 chord lengths downstream. This is achieved by attracting the incoming dipoles and disturbing their cohesion within the limits of the back foil's range of flapping motion. We also show that the impact on cycle averaged thrust varies from high augmentation to drag generation depending on the wake patterns downstream of the back foil. These findings provide new insights towards the design of biomimetic tandem propulsors, as they expand their working envelope and ability to rapidly increase or decrease the forward speed by manipulating the size of the shed vortices.

Reversed von K\'arm\'an streets are responsible for a velocity surplus in the wake of flapping foils, indicating the onset of thrust generation. However, the wake pattern cannot be predicted based solely on the flapping peak-to-peak amplitude $A$ and frequency $f$ because the transition also depends sensitively on other details of the kinematics. In this work we replace $A$ with the cycle-averaged swept trajectory $\mathcal{T}$ of the foil chord-line. Two dimensional simulations are performed for pure heave, pure pitch and a variety of heave-to-pitch coupling. In a phase space of dimensionless $\mathcal{T}-f$ we show that the drag-to-thrust wake transition of all tested modes occurs for a modified Strouhal $St_{\mathcal{T}}\sim 1$. Physically the product $\mathcal{T}\cdot f$ expresses the induced velocity of the foil and indicates that propulsive jets occur when this velocity exceeds $U_{\infty}$. The new metric offers a unique insight into the thrust producing strategies of biological swimmers and flyers alike as it directly connects the wake development to the chosen kinematics enabling a self similar characterisation of flapping foil propulsion., Comment: Rev 1

A scaling relation for the high frequency regime of streamwise energy spectra (in frequency domain) in the near-wall region is proposed. This is based on the dimensional analysis approach of \cite{perry1977} and \cite{Zamalloa2014} together with the hypothesis that the small-scale fluctuations in the near-wall region should only depend on the viscous scales, analogous to the Prandtl's law-of-the-wall for the mean flow. This allows us to examine the lower bound for the high frequency regime where law-of-the-wall in spectra would hold. Observations in high-Reynolds-number turbulent boundary layer data indicate that a conservative estimate for the start of this high frequency regime is $f^+$ = 0.005 (which corresponds to 200 viscous time-units) across a range of wall-normal positions and Reynolds numbers. This is sufficient to capture the energetic viscous-scaled motions such as the near-wall streaks, which has a time scale of approximately 100 viscous units. This scaling relation and the spectral collapse is consistent with the observations in internal flows at lower Reynolds numbers \cite{Zamalloa2014}., Comment: 12 pages (double spaced), 3 figures