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 $AR$) is found to improve thrust coefficients and thrust augmentation but the benefits level off towards higher values of $AR$. However, 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 into the mechanisms of these wake interactions and provide valuable information for the design of novel biomimetic propulsion systems.

In this work, we focus on a multi-hole pressure-probe-based flow measurement system for wind tunnel measurements that provides real-time feedback to a robot probe-manipulator, rendering the system autonomous. The system relies on a novel, computationally efficient flow analysis technique that translates the probe's point measurements of velocity and pressure into an updatable mean flow map that is accompanied by an uncertainty metric. The latter provides guidance to the manipulator when planning the optimal probe path. The probe is then guided by the robot in the flow domain until an available time budget has been exhausted, or until the uncertainty metric falls below a prescribed target threshold in the entire flow domain. We assess the capabilities of our new measurement system using computational fluid dynamics data, for which the ground truth is available in the form of a mean flow field. An application in a real wind tunnel setting is provided as well.

Flagellated bacteria propel themselves by rotating flexible flagella driven by independent motors. Depending on the rotation direction of the motors and the handedness of the helical filaments, the flagella either pull or push the cell body. Motivated by experimental observations of Magnetococcus marinus, we develop an elastohydrodynamic model to study the locomotion of a bi-flagellated bacterium with one puller flagellum and one pusher flagellum. In this model, the boundary integral technique and Kirchhoff rod model are employed respectively to calculate the hydrodynamic forces on the swimmer and model the elastic deformations of the flagella. Our numerical results demonstrate that the model bacterium travels along a double helical trajectory, which is consistent with the experimental observations. Varying the stiffness, orientations or positions of the flagella significantly changes the swimming characteristics. Notably, when either the applied torque is higher than a critical value or the flagellum stiffness is lower than a critical stiffness, the pusher flagellum exhibits overwhirling motion, resulting in a more complicated swimming style and a lower swimming speed. For a moderate flagellum stiffness, the swimming speed is insensitive to the rest configuration orientation over a wide range of orientation angles as the flagella deform to maintain alignment with the swimming direction.

Helical vortex systems, such as those found in the wakes of wind turbines, helicopter rotors and propellers, are subject to instabilities that lead to pairing between adjacent vortex loops. Certain modes of these instabilities can be triggered by an asymmetry in the rotor generating the vortices. In three-vortex systems, like those formed by many industrial rotors, the nonlinear vortex interactions are highly complex, introducing the need for a simple model to predict their dynamics. The current study presents a model for helical vortex systems based on an infinite strip of periodically repeating point vortices, whose motion can be computed using a single equation. This highly simplified model is shown to accurately reproduce the helical vortex dynamics predicted by a more sophisticated filament model and observed in water channel experiments on model rotors. The model is then used to investigate different types of vortex perturbations. Perturbation direction is found to have an important effect on the evolution of the instability, and displacements are observed to induce vortex pairing more quickly than circulation changes. These findings can be used to design asymmetric rotors that induce vortex breakdown more effectively, mitigating detrimental wake effects such as increased fatigue loading on downstream structures.

We demonstrate through the use of a unique acoustically driven microfluidic extensional rheometry platform (ADMiER) that a single measurement – i.e. the time required for a liquid bridge filament comprising a microlitre semen sample to thin and break up under elastocapillary stresses – constitutes an appropriate proxy for quantifying the motile sperm concentration of the sample in place of computer-assisted sperm analysis (CASA) and haemocytometer measurements used in conventional semen assessment – without the need to separately resolve for individual dependencies on each sperm parameter. By benchmarking diagnostic test accuracy results of blind random bull semen samples ($n=35$) against OpenCASA measurements of these parameters, ADMiER is capable of predicting sperm quality to 93.7 % accuracy, 91.4 % sensitivity and 97.5 % specificity, with respect to commonly adopted veterinary industry minimum values for fertility. These results therefore highlight the potential diagnostic capability of the platform as a conceptual first step towards the development of a rapid, low-cost and portable alternative for veterinary male bovine fertility assessment.

Rotorcraft can encounter highly unsteady flow when descending at a steep angle, leading to a flow condition called vortex ring state, which is associated with strong oscillatory airloads and substantial losses in mean rotor thrust. This study examines the aerodynamic coupling between closely arranged rotors in vertical flight and assesses the extent to which rotor–rotor interactions affect the rotor performance in this flight stage. Wind tunnel experiments were performed on a small-scale, dual-rotor set-up with adjustable rotor spacing, and the effect of rotor separation on thrust generation was quantified. Pairs of 4 in., 5 in. and 6 in. rotors ($3.0 \times 10^4< Re<8.1 \times 10^4$) were investigated, with load cell measurements showing significant thrust losses and concomitantly increased thrust oscillations as descent rate increased. Peak losses and fluctuations were consistently recorded at descent rates of 1.2–1.3 times the hover induced velocity for all rotor sizes and separations. While tests showed that the mean aerodynamic performance of dual-rotor systems is generally similar to that of single rotors, appreciable changes to the descent characteristics could be observed at low rotor separations. Particle image velocimetry flow visualization suggests considerable changes to the flow field as rotor separation decreases, where individual vortex ring systems merge into a single vortex ring structure.

This paper sets out to investigate the vortex flow of spinnaker yacht sails, which are low-aspect-ratio highly cambered wings used to sail downwind. We tested three model-scale sails with the same sections but different twists over a range of angles of attack in a water tunnel at a Reynolds number of 21 000. We measured the forces with a balance and the velocity field with particle image velocimetry. The sails experience massively separated three-dimensional flow and leading-edge vortices convect at half of the free-stream velocity in a turbulent shear layer. Despite the massive flow separation, the twist of the sail does not change the lift curve slope, in agreement with strip theory. As the angle of attack and the twist vary, flow reattachment might occur in the time-average sense, but this does not necessarily result in a higher lift to drag ratio as the vorticity field is marginally affected. Finally, we investigated the effect of secondary vorticity, vortex stretching and diffusion on the vorticity fluxes. Overall, these results provide new insights into the vortex flow and associated force generation mechanism of wings with massively separated flow.

The physics and mechanism of sheet/cloud cavitation in a convergent–divergent channel are investigated using synchronized dynamic surface pressure measurement and high-speed imaging in a water tunnel to probe the cavity shedding mechanism. Experiments are conducted at a fixed Reynolds number of Re = 7.8 × 105 for different values of the cavitation number σ between 1.20 and 0.65, ranging from intermittent inception cavitation, sheet cavitation to quasi-periodic cloud cavitation. Two distinct cloud cavitation regimes, i.e. the re-entrant jet and shockwave shedding mechanism, are observed, accompanied by complex flow phenomenon and dynamics, and are examined in detail. An increase in pressure fluctuation intensity at the numbers 3 and 4 transducer locations are captured during the transition from re-entrant jet to shockwave shedding mechanism. The spectral content analysis shows that, in cloud cavitation, several frequency peaks are identified with the dominant frequency caused by the large-scale cavity shedding process and the secondary frequency related to re-entrant jet/shockwave dynamics. Statistical analysis based on defined grey level profiles reveals that, in cloud cavitation, the double-peak behaviours of the probability density functions with negative skewness values are found to be owing to the interactions of the re-entrant jet/shockwave with cavities in the region of 0.25 ~ 0.65 mean cavity length (Lc). In addition, multi-scale proper orthogonal decomposition analysis with an emphasis on the flow structures in the region of 0.25 ~ 0.65 Lc reveals that, under the shockwave shedding mechanism, both the re-entrant jet and shockwave are captured and their interactions are responsible for the dynamics and statistics of cloud shedding process.

Natural ventilation can play an important role towards preventing the spread of airborne infections in indoor environments. However, quantifying natural ventilation flow rates is a challenging task due to significant variability in the boundary conditions that drive the flow. In the current study, we propose and validate an efficient strategy for using computational fluid dynamics to assess natural ventilation flow rates under variable conditions, considering the test case of a single-room home in a dense urban slum. The method characterizes the dimensionless ventilation rate as a function of the dimensionless ventilation Richardson number and the wind direction. First, the high-fidelity large-eddy simulation (LES) predictions are validated against full-scale ventilation rate measurements. Next, simulations with identical Richardson numbers, but varying dimensional wind speeds and temperatures, are compared to verify the proposed similarity relationship. Last, the functional form of the similarity relationship is determined based on 32 LES. Validation of the surrogate model against full-scale measurements demonstrates that the proposed strategy can efficiently inform accurate building-specific similarity relationships for natural ventilation flow rates in complex urban environments.

Vertical-axis wind turbines are excellent candidates to diversify wind energy technology, but their aerodynamic complexity limits industrial deployment. To improve the efficiency and lifespan of vertical-axis wind turbines, we desire data-driven models and control strategies that take into account the timing and duration of subsequent events in the unsteady flow development. Here, we aim to characterise the chain of events that leads to dynamic stall on a vertical-axis wind turbine blade and to quantify the influence of the turbine operation conditions on the duration of the individual flow development stages. We present time-resolved flow and unsteady load measurements of a wind turbine model undergoing dynamic stall for a wide range of tip-speed ratios. Proper orthogonal decomposition is used to identify dominant flow structures and to distinguish six characteristic stall stages: the attached flow, shear-layer growth, vortex formation, upwind stall, downwind stall and flow reattachment stage. The timing and duration of the individual stages are best characterised by the non-dimensional convective time. Dynamic stall stages are also identified based on aerodynamic force measurements. Most of the aerodynamic work is done during the shear-layer growth and the vortex formation stage which underlines the importance of managing dynamic stall on vertical-axis wind turbines.

Strongly coupled sequences of shock waves, known as shock trains, are present in high-speed propulsion systems, where the presence of sidewalls substantially modifies the boundary layer thickness, skin friction and streamwise pressure distribution. In the present contribution, scale-resolved numerical simulations are performed on supersonic channel (infinite span) and square duct flows to evaluate the effect of sidewall confinement with and without shock trains. Comparable secondary flow vortices are observed in the duct case with and without the presence of the shock train. The absence of a separation region at the leading shock of the duct case results in lower flow deflection compared with the channel case, leading to a reduced shock strength. The principal effect of the sidewalls is to cause a shock train that is approximately twice as long and composed of a larger number of shocks. A modification of previous models, based on a momentum thickness-based blockage parameter, leads to an improved collapse of the channel and duct cases.

We perform a two-dimensional numerical study on the thermal effect of porous media on global heat transport and flow structure in Rayleigh–Bénard (RB) convection, focusing on the role of thermal conductivity $\lambda$ of porous media, which ranges from $0.1$ to $50$ relative to the fluid. The simulation is carried out in a square RB cell with the Rayleigh number $Ra$ ranging from $10^7$ to $10^9$ and the Prandtl number $Pr$ fixed at $4.3$. The porosity of the system is fixed at $\phi =0.812$, with the porous media modelled by a set of randomly displayed circular obstacles. For a fixed $Ra$, the increase of conductivity shows a small effect on the total heat transfer, slightly depressing the Nusselt number. The limited influence comes from the small number of obstacles contacting with thermal plumes in the system as well as the counteraction of the increased plume area and the depressed plume strength. The study shows that the global heat transfer is insensitive to the conduction effect of separated porous media in the bulk region, which may have implications for industrial designs.

Shock waves interacting with turbulent boundary layers on wings can result first in self-sustained flow unsteadiness and eventually in structural vibration. Due to its importance to modern wing design and aircraft certification, the transonic flow physics continue to be investigated intensively. Herein we focus the discussion on three main aspects. First, we assess a practical implementation of an iterative resolvent algorithm in the linear harmonic incarnation of an industrial computational fluid dynamics code for computing optimal forcing and response modes. This heavily relies on the efficient solution of large sparse linear systems of equations. Second, we showcase its application as a predictive tool to detect transonic buffet flow unsteadiness, well before a global stability analysis can first identify its dynamics through weakly damped eigenmodes, using the NASA common research model at wind-tunnel conditions. Third, we discuss its ability to uncover modal physics, not identifiable through global stability analysis, revealing higher-frequency wake and wingtip vortex modes while shedding some light on the elusive finite wing equivalent of the aerofoil buffet mode. We demonstrate that earlier computational limitations of resolvent analysis, when solving the truncated singular value decomposition using matrix-forming methods with direct matrix factorisation, have been overcome ready for industrial use.

We investigate the ventilation conditions required to control the propagation of smoke, produced by a tunnel fire, in the presence of two inertial forcings: a transverse extraction system and a longitudinal flow. For that purpose, we performed a series of experiments in a reduced-scale tunnel, using a mixture of air and helium to simulate the release of hot smoke during a fire. Experiments were designed to focus on the ventilation flows that allow the buoyant release to be confined between two adjacent extraction vents. Different source conditions, in terms of density and velocity of the buoyant release, were analysed along with different vent configurations. Experiments allowed us to quantify the increase of the extraction velocity needed to confine the buoyant smoke, overcoming the effect of an imposed longitudinal velocity. Vents with a rectangular shape, and spanning over the whole tunnel width, provide the best performance. Finally, we studied the stratification conditions of the flow, individuating four regimes. Interestingly, when the stratification conditions fade out, as both the longitudinal flow and vertical extraction flows increase, the flow dynamics becomes almost independent of the forcing induced by the presence of buoyant smoke, which eventually acts as a passive scalar transported by the flow.

The knuckleball is considered to be one of the hardest pitches to hit in baseball due to its seemingly unpredictable motion. It has gained popularity in cricket in recent times. It is shown that the delivery referred to as knuckleball in cricket, at present, does not exhibit a zigzag motion and is, therefore, a misnomer. We propose a delivery in cricket that is associated with an erratic trajectory similar to the knuckleball pitch in baseball. Force measurement experiments in a wind tunnel on a new cricket ball in various orientations of the seam to the incoming flow and at different Reynolds number are carried out. The results are utilized to estimate the trajectory of knuckleball deliveries. The key parameters are the seam angle, speed and spin rate of the ball at the time of its release. Their effect on the trajectory is studied in detail. The optimal combination of these parameters that result in a knuckleball, which is likely hard for the batter to play, is identified.

The flow dynamics of small-width wall-attached jets generated by a Coand-effect nozzle is investigated by unsteady Reynolds-averaged Navier–Stokes simulations. The data are validated by comparison with hot-wire velocity measurements performed on the same flow configurations. The jets exhibit a complex topology strongly influenced not only by the spanwise vorticity (as usually observed in wall jets) but also by a vorticity component normal to the wall and induced by the shear layer developing on the jet sides. This results in an original U-shaped jet whose characteristics are studied in detail for three different mass flow rates. The robustness of the flow topology on a larger range of injected mass flow rates is finally presented and discussed in terms of the injected momentum near the wall. The resulting flow profiles point out that our injector is expected to be a promising candidate for active flow control in gas-turbine compressors for aeronautical and energy applications.

Displacement ventilation, where cool external air enters a room through low-level vents and warmer air leaves through high-level vents, is characterised by vertical gradients in pressure arising from the warmer indoor temperatures. Models usually assume that horizontal variations of temperature difference are small in comparison and are, therefore, unimportant. Small-scale laboratory experiments and computational fluid dynamics were used to examine these flows, driven by a uniformly heated floor. These experiments and simulations show that the horizontal variations of temperature difference can be neglected for predictions of the bulk ventilation rate; however, they also evidence that these horizontal variations can be significant and play a critical role in establishing the pattern of flow within the room – this renders the horizontal position of the low- and high-level vents (relative to one another) important. We consider two cases: single-ended (where inlet and outlet are at the same end of the room) and opposite-ended. In both cases the ventilation flow rate is the same. However, in the opposite-ended case a dead zone is established in the upper part of the room which results in significant horizontal variations. We consider the formation of this dead zone by examining the streamline patterns and the age of air within the room. We discuss the implications for occupant exposure to pollutants and airborne disease.

The present paper simplifies the naturally formed dunes (riverbeds) as large-scale three-dimensional staggered wavy walls to investigate the features of the accompanying secondary flows and streamwise vortices via large-eddy simulation. A comparison between the swirling strength and the mean velocities suggests where a secondary flow induces upwash or downwash motions. Moreover, we propose a pseudo-convex wall mechanism to interpret the directionality of the secondary flow. The centrifugal instability criterion is then used to reveal the generation of the streamwise vortices. Based on these analytical results, we found that the streamwise vortices are generated in the separation and reattachment points on both characteristic longitudinal–vertical and horizontal cross-sections, which is related to the curvature effect of the turbulent shear layer. Furthermore, the maximum Görtler number characterized by the ratio of centrifugal to viscous effects suggests that, for fixed ratio of spanwise- to streamwise-wavelength cases, the strongest centrifugal instability occurring on the longitudinal–vertical cross-section gradually dominates with the increases in amplitude. A similar trend for the cases with varied spanwise wavelength can also be found. It is also found that the streamwise vortices are generated more readily via transverse flow around the crest near the separation and reattachment points when the ratio of spanwise- to streamwise-wavelength equals 1.

Predicting the laminar to turbulent transition is an important aspect of computational fluid dynamics because of its impact on skin friction. Traditional transition prediction methods such as local stability theory or the parabolized stability equation method do not allow for the consideration of strongly non-parallel boundary layer flows, as in the presence of surface defects (bumps, steps, gaps, etc.). A neural network approach, based on an extensive database of two-dimensional incompressible boundary layer stability studies in the presence of gap-like surface defects, is used. These studies consist of linearized Navier–Stokes calculations and provide information on the effect of surface irregularity geometry and aerodynamic conditions on the transition to turbulence. The physical and geometrical parameters characterizing the defect and the flow are then provided to a neural network whose outputs inform about the effect of a given gap on the transition through the ${\rm \Delta} N$ method (where N represents the amplification of the boundary layer instabilities).