Search

This paper investigates flow past a backward-facing step (BFS) in a duct at Reynolds number $Re = 5080$ based on step height, mean inflow velocity and fluid kinematic viscosity. The flow configuration matches a combustion experiment conducted by Pitz and Daily in 1983. High-resolution velocity fields are obtained in OpenFOAM by direct numerical simulation, and the flow field analysed by Lagrangian approaches. Trajectories of fluid particles in areas of interest are obtained by high-order numerical integration, and used to compute finite-time Lyapunov exponents (FTLEs) and polar rotation angles. Lagrangian coherent structures (LCSs) are extracted using geodesic theory, including hyperbolic LCSs and elliptic LCSs. We use complementary qualitative and quantitative LCS analyses to uncover the underlying flow structures. Notably, we find that a flow pathway in which fluid particles rarely diverge from adjacent particles is opened and closed by FTLE ridges determined by the periodic shedding of vortices from the BFS. Two dominant vortices with significant Lagrangian coherence, generated respectively by the separated boundary layer and shear layer, are self-sustaining and of comparable strength. Hyperbolic repelling LCSs act as transport barriers between the pathway and cores of the coherent vortices, thus playing a major part in the fluid entrainment process. Interactions between these different geometric regions partitioned by LCSs lead to intrinsic complexity in the BFS flow.

This paper presents the results of the time resolved flow field measurements around a realistic rowing oar blade that moves along a realistic path through water. To the authors' knowledge no prior account of this complex flow field has been given. Simultaneously with the flow field measurements, the hydrodynamic forces acting on the blade were measured. These combined measurements allow us to identify the relevant flow physics that governs rowing propulsion, and subsequently use this information to adjust the oar blade configuration to improve rowing propulsion. Analysis of the instationary flow field around the oar blade during the drive phase indicated how the initial formation, and subsequent development, of leading-edge and trailing-edge vortices are related to the generation of instationary lift and drag forces, and how these forces contribute to rowing propulsion. It is shown that the observed individual flow mechanisms are similar to the flow mechanisms observed in bird flight, but that the overall propulsive mechanism for rowing propulsion is fundamentally different. To quantify the rowing propulsion efficiency, we introduced the energetic efficiency and the impulse efficiency, where the latter can be interpreted as the alignment of the generated impulse with the propulsive direction. It is found that in the conventional oar blade configuration, the generated impulse is not aligned with the propulsive direction, indicating that the propulsion is suboptimal. By adjusting the angle at which the blade is attached to the oar, the generation of leading- A nd trailing-edge vortices is altered such that the generated impulse better aligns with the propulsive direction, thus increasing the efficiency.

The theory of Kolmogorov, enunciated for very large Reynolds numbers, has progressively been shown to be inoperative for characterizing flows of practical relevance. Yet, in a recent study by Alves Portela et al. (J. Fluid Mech., vol. 896, 2020, A16), the turbulence statistics in the very near wake of a square prism at modest Reynolds numbers, reveal a significant portion of scales complying with a cascade of Kolmogorov type. By resorting to a generalized version of the Karman–Howarth–Kolmogorov equation, this intriguing observation is shown to be an illusion, hiding a measurable influence of coherent structures and statistical inhomogeneity. This striking conclusion highlights that a complete statistical theory of turbulence cannot dispense with the influence of large scales, possibly coherent, motions.

The separating and reattaching flows and the wake of a finite rectangular plate are studied by means of direct numerical simulation data. The large amount of information provided by the numerical approach is exploited here to address the multi-scale features of the flow and to assess the self-sustaining mechanisms that form the basis of the main unsteadinesses of the flows. We first analyse the statistically dominant flow structures by means of three-dimensional spatial correlation functions. The developed flow is found to be statistically dominated by quasi-streamwise vortices and streamwise velocity streaks as a result of flow motions induced by hairpin-like structures. On the other hand, the reverse flow within the separated region is found to be characterized by spanwise vortices. We then study the spectral properties of the flow. Given the strongly inhomogeneous nature of the flow, the spectral analysis has been conducted along two selected streamtraces of the mean velocity field. This approach allows us to study the spectral evolution of the flow along its paths. Two well-separated characteristic scales are identified in the near-wall reverse flow and in the leading-edge shear layer. The first is recognized to represent trains of small-scale structures triggering the leading-edge shear layer, whereas the second is found to be related to a very large-scale phenomenon that embraces the entire flow field. A picture of the self-sustaining mechanisms of the flow is then derived. It is shown that very-large-scale fluctuations of the pressure field alternate between promoting and suppressing the reverse flow within the separation region. Driven by these large-scale dynamics, packages of small-scale motions trigger the leading-edge shear layers, which in turn created them, alternating in the top and bottom sides of the rectangular plate with a relatively long period of inversion, thus closing the self-sustaining cycle.

Flow separation and its control have been the subject of intensive research for decades. Flow separation occurs when the boundary layer loses contact with the associated confining wall, which is usually caused by a pressure gradient acting against the local flow direction. Numerous strategies exist to control flow separation, and in this study we demonstrate experimentally that vertical flow separation over steep slopes in shallow free-surface flows may be suppressed by contracting the flow horizontally upstream of the slope. We found that, unexpectedly, introducing lateral non-uniformity in the upstream flow field could suppress vertical flow separation for steep slopes up to 1 in 2. This study reveals the possibility of two different flow states over steep slopes; (i) a vertically attached flow combined with horizontal convergence, and (ii) a vertically detached flow combined with horizontal divergence. A detailed analysis of the dynamics of the two different flow states is presented. Although a predictive relation determining the transition point between the two flow states was not found in the current study, the observed phenomena were shown to be strongly related to the magnitude of the lateral gradient at the upstream edge of the slope. The results demonstrate a significant influence of the vertical flow state - separated or attached - on the shear stress at the confining boundaries of the flow.

International audience; The caudal fin of swimming animals can be modelled as a thrust-producing flapping foil. When considered alone, such a foil produces on average a jet wake with a positive net momentum. It has been argued that the instability characteristics of these averaged wakes are linked to the propulsion efficiency of swimming animals. Here, we reconsider this question by taking into account both the thrust and the drag exerted on a self-propelled swimming body. To do so, we study the stability of a family of momentumless wakes, constructed as the Oseen approximation of a force doublet moving at constant velocity. By performing a local stability analysis, we first show that these wakes undergo a transition from absolute to convective instability. Then, using the time stepper approach by integrating the linearised Navier–Stokes system, we investigate the global stability and reveal the influence of a non-parallel base flow as well as the role of the locally absolutely unstable upstream region in the wake. Finally, to complete the global scenario, we address the nonlinear evolution of the wake disturbance. These results are then discussed in the context of aquatic locomotion. According to the present stability results, and assuming the Oseen approximation whose validity has been assessed only for moderate Reynolds number, the momentumless wake of aquatic animals is generally stable, whereas the corresponding thrust part is unstable. It is therefore essential to consider all forces exerted on a self-propelled animal when discussing its wake stability and its propulsion efficiency.

A direct numerical simulation (DNS) of a plane wall jet is performed at a Reynolds number of $Re_{j}=7500$. The streamwise length of the domain is long enough to achieve self-similarity for the mean flow and the Reynolds shear stress. This is the highest Reynolds number wall jet DNS for a large domain achieved to date. The high resolution simulation reveals the unsteady flow field in great detail and shows the transition process in the outer shear layer and inner boundary layer. Mean flow parameters of maximum velocity decay, wall shear stress, friction coefficient and jet spreading rate are consistent with several other studies reported in the literature. Mean flow, Reynolds normal and shear stress profiles are presented with various scalings, revealing the self-similar behaviour of the wall jet. The Reynolds normal stresses do not show complete similarity for the given Reynolds number and domain length. Previously published inner layer budgets based on LES are inaccurate and those that have been measured are only available in the outer layer. The current DNS provides fully balanced, explicitly calculated budgets for the turbulence kinetic energy, Reynolds normal stresses and Reynolds shear stress in both the inner and outer layers. The budgets are scaled with inner and outer variables. The inner-scaled budgets in the near wall region show great similarity with turbulent boundary layers. The only remarkable difference is for the turbulent diffusion in the wall-normal Reynolds stress and Reynolds shear stress budgets. The outer layer interacts with the inner layer through turbulent diffusion and the excess energy from the wall-normal direction is transferred to the spanwise direction.

We discuss energetic restrictions on the entrainment coefficient${\it\alpha}$for axisymmetric jets and plumes. The resulting entrainment relation includes contributions from the mean flow, turbulence and pressure, fundamentally linking${\it\alpha}$to the production of turbulence kinetic energy, the plume Richardson number$\mathit{Ri}$and the profile coefficients associated with the shape of the buoyancy and velocity profiles. This entrainment relation generalises the work by Kaminskiet al. (J. Fluid Mech., vol. 526, 2005, pp. 361–376) and Fox (J. Geophys. Res., vol. 75, 1970, pp. 6818–6835). The energetic viewpoint provides a unified framework with which to analyse the classical entrainment models implied by the plume theories of Mortonet al.(Proc. R. Soc. Lond.A, vol. 234, 1955, pp. 1–23) and Priestley & Ball (Q. J. R. Meteorol. Soc., vol. 81, 1954, pp. 144–157). Data for pure jets and plumes in unstratified environments indicate that to first order the physics is captured by the Priestley and Ball entrainment model, implying that (1) the profile coefficient associated with the production of turbulence kinetic energy has approximately the same value for pure plumes and jets, (2) the value of${\it\alpha}$for a pure plume is roughly a factor of$5/3$larger than for a jet and (3) the enhanced entrainment coefficient in plumes is primarily associated with the behaviour of the mean flow and not with buoyancy-enhanced turbulence. Theoretical suggestions are made on how entrainment can be systematically studied by creating constant-$\mathit{Ri}$flows in a numerical simulation or laboratory experiment.

Gas injection into a liquid cross-flow is examined for the case where the gas is injected beneath a horizontal flat surface. For moderate Froude numbers, the gas pocket that is formed will rise toward the flow boundary under the action of buoyancy, a condition that is conducive to the formation of gas layers for friction-drag reduction on the surface. At the location of gas injection, a plume whose geometry is related to the mass and momentum flux of the injected gas and liquid cross-flow is formed, and the influence of buoyancy is minimal. However, as the gas pocket convects downstream, buoyancy brings the gas back upward to the flow boundary, and leads to the bifurcation of the pocket into two distinct branches, forming a stable ‘V’-shape. Under some conditions, the flow between the two gas branches is almost entirely liquid, while for others there exists a bubbly flow or a continuous sheet of gas between the branches. The sweep angle and cross-sectional geometry of the gas branches are related to free-stream speed and boundary-layer thickness of the liquid cross-flow, the mass-injection rate of the gas, the diameter of the injection orifice and the gas outlet mean velocity and gas–jet angle. Data for a range of experimental conditions are used to scale the flow and results are compared to numerical computations of the flow, and these data are used to illustrate the underlying flow processes responsible leading to the formation the stable and straight gas branches. A simple model based on the balance of forces around a stable gas branch is presented and used to scale the observed data, and we use the results of this analysis and the computations to discuss how the process of gas injection may interact with the formation of the stable gas pockets farther downstream.

The onset of two-dimensional instabilities in the flow past two side-by-side circular cylinders is numerically investigated in the ranges $0.1\leq g\leq 3$ and $\mathit{Re}, with $g$ being the non-dimensional gap spacing between the surfaces of the two cylinders and $\mathit{Re}$ the Reynolds number. A comprehensive, global stability analysis of the symmetric base flow is carried out, indicating that three harmonic modes and one steady antisymmetric mode become unstable at different values of $g$ and $\mathit{Re}$. These modes are known to promote distinct flow regimes at increasing values of $g$: single bluff-body, asymmetric, in-phase and antiphase synchronized vortex shedding. For each mode, the inherent structural sensitivity is examined in order to identify the core region of the related instability mechanism. In addition, by exploiting the structural sensitivity analysis to base flow modifications, a passive control strategy is proposed for the simultaneous suppression of the two synchronized shedding modes using two small secondary cylinders. Its effectiveness is then validated a posteriori by means of direct numerical simulations.

Wind tunnel measurements were performed for the wake produced by a three-bladed wind turbine immersed in uniform flow. These tests show the presence of a vorticity structure in the near-wake region mainly oriented along the streamwise direction, which is denoted as the hub vortex. The hub vortex is characterized by oscillations with frequencies lower than that connected to the rotational velocity of the rotor, which previous works have ascribed to wake meandering. This phenomenon consists of transversal oscillations of the wind turbine wake, which might be excited by the vortex shedding from the rotor disc acting as a bluff body. In this work, temporal and spatial linear stability analyses of a wind turbine wake are performed on a base flow obtained with time-averaged wind tunnel velocity measurements. This study shows that the low-frequency spectral component detected experimentally matches the most amplified frequency of the counter-winding single-helix mode downstream of the wind turbine. Then, simultaneous hot-wire measurements confirm the presence of a helicoidal unstable mode of the hub vortex with a streamwise wavenumber roughly equal to that predicted from the linear stability analysis.

We examine the dynamics of a semi-infinite vortex sheet attached not to a semi-infinite plate but instead to a rigid right-angled wedge, with the sheet aligned along one of its edges. Our approach to this problem, which was suggested by David Crighton, accords well with the fundamental ethos of Crighton’s work, which was characterized by ‘the application of rigorous mathematical approximations to fluid mechanical idealizations of practically relevant problems’ (Ffowcs Williams, Annu. Rev. Fluid Mech., vol.34, 2002, pp.37–49). The resulting linearised unsteady potential flow is forced by an oscillatory dipole in the uniform stream passing along the top of the wedge, while there is stagnant fluid in the remaining quadrant. Spatial instability is considered according to well-established methods: causality is enforced by allowing the frequency to become temporarily complex. The essentially quadrant-type geometry replaces the usual Wiener–Hopf technique by the Mellin transform. The core difficulty is that a first-order difference equation of period 4 requires a solution of period unity. As a result, the complex fourth roots $(\pm 1\pm \text{i})$ of $-4$ appear in the complementary function. The Helmholtz instability wave is excited and requires careful handling to obtain explicit results for the amplitude of the instability wave.

Dieser Beitrag ist mit Zustimmung des Rechteinhabers aufgrund einer (DFG geförderten) Allianz- bzw. Nationallizenz frei zugänglich. This publication is with permission of the rights owner freely accessible due to an Alliance licence and a national licence (funded by the DFG, German Research Foundation) respectively. We propose a generalization of proper orthogonal decomposition (POD) for optimal flow resolution of linearly related observables. This Galerkin expansion, termed ‘observable inferred decomposition’ (OID), addresses a need in aerodynamic and aeroacoustic applications by identifying the modes contributing most to these observables. Thus, OID constitutes a building block for physical understanding, least-biased conditional sampling, state estimation and control design. From a continuum of OID versions, two variants are tailored for purposes of observer and control design, respectively. Firstly, the most probable flow state consistent with the observable is constructed by a ‘least-residual’ variant. This version constitutes a simple, easily generalizable reconstruction of the most probable hydrodynamic state to preprocess efficient observer design. Secondly, the ‘least-energetic’ variant identifies modes with the largest gain for the observable. This version is a building block for Lyapunov control design. The efficient dimension reduction of OID as compared to POD is demonstrated for several shear flows. In particular, three aerodynamic and aeroacoustic goal functionals are studied: (i) lift and drag fluctuation of a two-dimensional cylinder wake flow; (ii) aeroacoustic density fluctuations measured by a sensor array and emitted from a two-dimensional compressible mixing layer; and (iii) aeroacoustic pressure monitored by a sensor array and emitted from a three-dimensional compressible jet. The most ‘drag-related’, ‘lift-related’ and ‘loud’ structures are distilled and interpreted in terms of known physical processes.

The coherent structures of a turbulent wake generated behind a bluff three-dimensional axisymmetric body are investigated experimentally at a diameter-based Reynolds number of $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}{\sim }2\times 10^5$. Proper orthogonal decomposition of base pressure measurements indicates that the most energetic coherent structures retain the structure of the symmetry-breaking laminar instabilities and are manifested as unsteady vortex shedding with azimuthal wavenumber $m={\pm }1$. In a rotating reference frame, the shedding preserves the reflectional symmetry and is linked with a reflectionally symmetric mean pressure distribution on the base. Due to a slow rotation of the symmetry plane of the turbulent wake around the axis of the body, statistical axisymmetry is recovered in the time average. The ratio of the time scales associated with the slow rotation of the symmetry plane and the vortex shedding is of order 100.

International audience; We compute the spatial optimal energy amplification of steady inflow perturbations in a non-parallel wake and analyse their stabilizing action on the global mode instability. The optimal inflow perturbations, which are assumed spanwise periodic and varicose, consist in streamwise vortices that induce the downstream spatial transient growth of streamwise streaks. The maximum energy amplification of the streaks increases with the spanwise wavelength of the perturbations, in accordance with previous results obtained for the temporal energy growth supported by parallel wakes. A family of increasingly streaky wakes is obtained by forcing optimal inflow perturbations of increasing amplitude and then solving the nonlinear Navier-Stokes equations. We show that the linear global instability of the wake can be completely suppressed by forcing optimal perturbations of sufficiently large amplitude. The attenuation and suppression of self-sustained oscillations in the wake by optimal 3D perturbations is confirmed by fully nonlinear numerical simulations. We also show that the amplitude of optimal spanwise periodic (3D) perturbations of the basic flow required to stabilize the global instability is much smaller than the one required by spanwise uniform (2D) perturbations despite the fact that the first order sensitivity of the global eigenvalue to basic flow modifications is zero for 3D spanwise periodic modifications and non-zero for 2D modifications. We therefore conclude that first-order sensitivity analyses can be misleading if used far from the instability threshold, where higher order terms are the most relevant

This work is dedicated to systematically studying and predicting the wake characteristics of a yawed wind turbine immersed in a turbulent boundary layer. To achieve this goal, wind tunnel experiments were performed to characterize the wake of a horizontal-axis wind turbine model. A high-resolution stereoscopic particle image velocimetry system was used to measure the three velocity components in the turbine wake under different yaw angles and tip-speed ratios. Moreover, power and thrust measurements were carried out to analyse the performance of the wind turbine. These detailed wind tunnel measurements were then used to perform a budget study of the continuity and Reynolds-averaged Navier–Stokes equations for the wake of a yawed turbine. This theoretical analysis revealed some notable features of the wakes of yawed turbines, such as the asymmetric distribution of the wake skew angle with respect to the wake centre. Under highly yawed conditions, the formation of a counter-rotating vortex pair in the wake cross-section as well as the vertical displacement of the wake centre were shown and analysed. Finally, this study enabled us to develop general governing equations upon which a simple and computationally inexpensive analytical model was built. The proposed model aims at predicting the wake deflection and the far-wake velocity distribution for yawed turbines. Comparisons of model predictions with the wind tunnel measurements show that this simple model can acceptably predict the velocity distribution in the far wake of a yawed turbine. Apart from the ability of the model to predict wake flows in yawed conditions, it can provide valuable physical insight on the behaviour of turbine wakes in this complex situation.