An exact solution of a baroclinic three-dimensional vortex dipole in geophysical flows with constant background rotation and constant background stratification is provided under the quasi-geostrophic (QG) approximation. The motion of the dipole is unsteady but the potential vorticity contours move rigidly. The vortex comprises three potential vorticity anomaly modes, with a radial dependence given by the spherical Bessel functions and with azimuthal and polar dependences given by the spherical harmonics. The first mode, or spherical mode, accounts for the horizontal asymmetry of the vortex dipole and curvature of the dipole’s horizontal trajectory. The second mode, or dipolar mode, accounts for the speed of displacement of the vortex dipole. A third mode, or vertical tilting mode, accounts for the dipole’s vertical asymmetry. The QG vertical velocity field has two contributions: the first one is octupolar and depends entirely on the dipolar mode, and the second one is dipolar and depends on the nonlinear interaction between dipolar and vertical tilting modes.

Topology is introducing new tools for the study of fluid waves. The existence of unidirectional Yanai and Kelvin equatorial waves has been related to a topological invariant, the Chern number, that describes the winding of $f$-plane shallow water eigenmodes around band-crossing points in parameter space. In this previous study, the topological invariant was a property of the interface between two hemispheres. Here we ask whether a topological index can be assigned to each hemisphere. We show that this can be done if the shallow water model in the $f$-plane geometry is regularized by an additional odd-viscosity term. We then compute the spectrum of a shallow water model with a sharp equator separating two flat hemispheres, and recover the Kelvin and Yanai waves as two exponentially trapped waves along the equator, with all the other modes delocalized into the bulk. This model provides an exactly solvable example of bulk-interface correspondence in a flow with a sharp interface, and offers a topological interpretation for some of the transition modes described by Iga (J. Fluid Mech., vol. 294, 1995, pp. 367–390). It also paves the way towards a topological interpretation of coastal Kelvin waves along a boundary and, more generally, to an understanding of bulk-boundary correspondence in continuous media.

We determine the smallest instantaneous increase in the strength of an opposing wind that is necessary to permanently reverse the forward displacement flow that is driven by a two-layer thermal stratification. With an interpretation in terms of the flow’s energetics, the results clarify why the ventilation of a confined space with a stably stratified buoyancy field is less susceptible to being permanently reversed by the wind than the ventilation of a space with a uniform buoyancy field. For large opposing wind strengths we derive analytical upper and lower bounds for the system’s marginal stability, which exhibit a good agreement with the exact solution, even for modest opposing wind strengths. The work extends a previous formulation of the problem (Lishman & Woods, Build. Environ., vol. 44 (4), 2009, pp. 666–673) by accounting for the transient dynamics and energetics associated with the homogenisation of the interior, which prove to play a significant role in buffering temporal variations in the wind.

A monochromatic, small amplitude, normally incident standing wave on a sloping beach is unstable to perturbation by subharmonic (half the frequency) edge waves. At equilibrium, edge wave shoreline amplitudes can exceed incident wave amplitudes. Here, the effect of incident wave randomness on subharmonic edge wave excitation is explored following a weakly nonlinear stability analysis under the assumption of narrow-band incident random waves. Edge waves respond to variations in both incident wave phase and amplitude, and the edge wave amplitudes and incident wave groups vary on similar time scales. When bottom friction is included, intermittent subharmonic edge wave excitation is predicted due to the combination of bottom friction and wave phase. Edge wave amplitude can be near zero for long times, but for short periods reaches relatively large values, similar to amplitudes with monochromatic incident waves and no friction.

Several metres below the coastal ocean surface there are areas of high ecological activity that contain thin layers of concentrated motile phytoplankton. Gyrotactic trapping has been proposed as a potential mechanism for layer formation of bottom-heavy swimming algae cells, especially in flows where the vorticity varies linearly with depth (Durham et al., Science, vol. 323(5917), 2009, pp. 1067–1070). Using a continuum model for dilute microswimmer suspensions, we report that an instability of a gyrotactically trapped cell layer can arise in a pressure-driven plane channel flow. The linear stability analysis reveals that the equilibrium cell-layer solution is hydrodynamically unstable due to negative microswimmer buoyancy (i.e. a gravitational instability) over a range of biologically relevant parameter values. The critical cell concentration for this instability is found to be $N_{c}\simeq 10^{4}~\text{cells}~\text{cm}^{-3}$, a value comparable to the typical maximum cell concentration observed in thin layers. This result indicates that the instability may be a potential mechanism for limiting the layer’s maximum cell concentration, especially in regions where turbulence is weak, and motivates the study of its nonlinear evolution, perhaps, in the presence of turbulence.

The fundamental challenge to characterize and quantify thermal transport in the strongly nonlinear regime of Rayleigh–Bénard convection – the buoyancy-driven flow of a horizontal layer of fluid heated from below – has perplexed the fluid dynamics community for decades. Rayleigh proposed controlling the temperature of thermally conducting boundaries in order to study the onset of convection, in which case vertical heat transport gauges the system response. Conflicting experimental results for ostensibly similar set-ups have confounded efforts to discriminate between two competing theories for how boundary layers and interior flows interact to determine transport through the convecting layer asymptotically far beyond onset. In a conceptually new approach, Bouillaut, Lepot, Aumaître and Gallet (J. Fluid Mech., vol. 861, 2019, R5) devised a procedure to radiatively heat a portion of the fluid domain bypassing rigid conductive boundaries and allowing for dissociation of thermal and viscous boundary layers. Their experiments reveal a new level of complexity in the problem suggesting that heat transport scaling predictions of both theories may be realized depending on details of the thermal forcing.

The shear-induced collective diffusivity down a concentration gradient in a viscous emulsion is computed using direct numerical simulation. A layer of randomly packed drops subjected to a shear flow, shows the layer width to increase with the $1/3$ power of time, consistent with a semi-dilute theory that assumes a diffusivity linear with concentration. This characteristic scaling and the underlying theory are used to compute the collective diffusivity coefficient. This is the first ever computation of this quantity for a system of deformable particles using fully resolved numerical simulation. The results match very well with previous experimental observations. The coefficient of collective diffusivity varies non-monotonically with the capillary number, due to the competing effects of increasing deformation and drop orientation. A phenomenological correlation for the collective diffusivity coefficient as a function of capillary number is presented. We also apply an alternative approach to compute collective diffusivity, developed originally for a statistically homogeneous rigid sphere suspension – computing the dynamic structure factor from the simulated droplet positions and examining its time variation at small wavenumber. We show that the results from this alternative approach qualitatively agree with our computation of collective diffusivity including the prediction of the non-monotonic variation of diffusivity with the capillary number.

This work proposes a feedback-loop strategy to suppress intrinsic oscillations of resonating flows in the fully nonlinear regime. The frequency response of the flow is obtained from the resolvent operator about the mean flow, extending the framework initially introduced by McKeon & Sharma (J. Fluid Mech., vol. 658, 2010, pp. 336–382) to study receptivity mechanisms in turbulent flows. Using this linear time-invariant model of the nonlinear flow, modern control methods such as structured ${\mathcal{H}}_{\infty }$-synthesis can be used to design a controller. The approach is successful in damping self-sustained oscillations associated with specific eigenmodes of the mean-flow spectrum. Despite excellent performance, the linear controller is however unable to completely suppress flow oscillations, and the controlled flow is effectively attracted towards a new dynamical equilibrium. This new attractor is characterized by a different mean flow, which can in turn be used to design a second controller. The method can then be iterated on subsequent mean flows, until the coupled system eventually converges to the base flow. An intuitive parallel can be drawn with Newton’s iteration: at each step, a linearized model of the flow response to a perturbation of the input is sought, and a new linear controller is designed, aiming at further reducing the fluctuations. The method is illustrated on the well-known case of two-dimensional incompressible open-cavity flow at Reynolds number $Re=7500$, where the fully developed flow is initially quasiperiodic (2-torus state). The base flow is reached after five iterations. The present work demonstrates that nonlinear control problems may be solved without resorting to nonlinear reduced-order models. It also shows that physically relevant linear models can be systematically derived for nonlinear flows, without resorting to black-box identification from input–output data; the key ingredient being frequency-domain models based on the linearized Navier–Stokes equations about the mean flow. Applicability to amplifier flows and turbulent dynamics has, however, yet to be investigated.

The influence of the nozzle aspect ratio ($AR=1$, 2 and 4), stroke length ($L_{0}=1.85$, 3.7 and 5.55) and Reynolds number ($Re=79$, 158, 316 and 632) on the behaviour of elliptic synthetic jets is studied experimentally. Laser-induced fluorescence and two-dimensional and stereoscopic particle image velocimetry are used to analyse the vortex dynamics and evolution mechanism. It is found that the fluid elements around the major axis of an elliptic vortex ring move downstream faster and tend to approach the centreline, while the fluid elements around the minor axis move downstream at a slower speed and away from the centreline, thereby resulting in the occurrence of the well-known axis-switching phenomenon for elliptic synthetic jets. During this process, a pair of arc-like vortices forms ahead of the primary vortex ring, and they are constituted by streamwise vortices in the leg part and spanwise vortices in the head part; two pairs of streamwise vortices form from the inside of the primary vortex ring and develop in the tails. The streamwise vortices are pushed away progressively from the centreline by the synthetic jet vortex rings that are formed during the subsequent periods. These additional vortical structures for non-circular synthetic jets show regular and periodic characteristics, which are quite different from the previous findings for non-circular jets. Their mutual interaction with the vortex ring causes significant changes in the topology of elliptic synthetic jets, which further results in the variation of the statistical characteristics. Increasing the aspect ratio, stroke length and Reynolds number will make the evolution of the synthetic jet become more unstable and complex. In addition, the entrainment rate of an elliptical synthetic jet is larger than that of a circular synthetic jet and it increases with the nozzle aspect ratio ($AR\leqslant 4$) and Reynolds number. It is indicated that the formation of streamwise vortices could enhance the entrainment rate. This finding provides substantial evidence for the potential application of elliptic synthetic jets for effective flow control.

In this paper we investigate the effect of stable stratification on plane Couette flow when gravity is oriented in the spanwise direction. When the flow is turbulent we observe near-wall layering and associated new mean flows in the form of large-scale spanwise-flattened streamwise rolls. The layers exhibit the expected buoyancy scaling $l_{z}\sim U/N$ where $U$ is a typical horizontal velocity scale and $N$ the buoyancy frequency. We associate the new coherent structures with a stratified modification of the well-known large-scale secondary circulation in plane Couette flow. We find that the possibility of the transition to sustained turbulence is controlled by the relative size of this buoyancy scale to the spanwise spacing of the streaks. In parts of parameter space transition can also be initiated by a newly discovered linear instability in this system (Facchini et al., J. Fluid Mech., vol. 853, 2018, pp. 205–234). When wall turbulence can be sustained the linear instability opens up new routes in phase space for infinitesimal disturbances to initiate turbulence. When the buoyancy scale suppresses turbulence the linear instability leads to more ordered nonlinear states, with the possibility for intermittent bursts of secondary shear instability.

We consider peeling of an elastic sheet away from an elastic substrate through propagation of a fluid-filled crack along the interface between the two. The peeling is driven by a bending moment applied to the sheet and is resisted by viscous flow towards the crack tip and by the toughness of any bonding between the sheet and the substrate. Travelling-wave solutions are determined using lubrication theory coupled to the full equations of elasticity and fracture. The propagation speed $v$ scales like $M^{3}/\unicode[STIX]{x1D707}\bar{E}^{2}d^{5}=Bd\unicode[STIX]{x1D705}^{3}/144\unicode[STIX]{x1D707}$, where $d$ is the sheet’s thickness, $B=\bar{E}d^{3}/12$ its stiffness, $\bar{E}=E/(1-\unicode[STIX]{x1D708}^{2})$ its plane-strain modulus, $\unicode[STIX]{x1D707}$ the fluid viscosity, $M$ the applied bending moment and $\unicode[STIX]{x1D705}=M/B$ the sheet’s curvature due to bending; and the prefactor depends on the dimensionless toughness. If the toughness is small then there is a region of dry shear failure ahead of the fluid-filled region. The expressions for the propagation speed have been used to derive new similarity solutions for the spread of an axisymmetric fluid-filled blister in a variety of regimes: constant-flux injection resisted by elastohydrodynamics in the tip leads to spread proportional to $t^{4/13}$, $t^{4/17}$ and $t^{7/19}$ for peeling-by-bending, gravitational spreading and peeling-by-pulling, respectively.

The application of the Kutta condition to unsteady flows has been controversial over the years, with increased research activities over the 1970s and 1980s. This dissatisfaction with the Kutta condition has been recently rejuvenated with the increased interest in low-Reynolds-number, high-frequency bio-inspired flight. However, there is no convincing alternative to the Kutta condition, even though it is not mathematically derived. Realizing that the lift generation and vorticity production are essentially viscous processes, we provide a viscous extension of the classical theory of unsteady aerodynamics by relaxing the Kutta condition. We introduce a trailing-edge singularity term in the pressure distribution and determine its strength by using the triple-deck viscous boundary layer theory. Based on the extended theory, we develop (for the first time) a theoretical viscous (Reynolds-number-dependent) extension of the Theodorsen lift frequency response function. It is found that viscosity induces more phase lag to the Theodorsen function particularly at high frequencies and low Reynolds numbers. The obtained theoretical results are validated against numerical laminar simulations of Navier–Stokes equations over a sinusoidally pitching NACA 0012 at low Reynolds numbers and using Reynolds-averaged Navier–Stokes equations at relatively high Reynolds numbers. The physics behind the observed viscosity-induced lag is discussed in relation to wake viscous damping, circulation development and the Kutta condition. Also, the viscous contribution to the lift is shown to significantly decrease the virtual mass, particularly at high frequencies and Reynolds numbers.

Rational large Reynolds number matched asymptotic expansions of three-dimensional nonlinear magneto-hydrodynamic (MHD) states are the concern of this contribution. The nonlinear MHD states, assumed to be predominantly driven by a unidirectional shear, can be sustained without any linear instability of the base flow and hence are responsible for subcritical transition to turbulence. Two classes of nonlinear MHD states are found. The first class of nonlinear states emerged out of a nice combination of the purely hydrodynamic vortex/wave interaction theory by Hall & Smith (J. Fluid Mech., vol. 227, 1991, pp. 641–666) and the resonant absorption theories on Alfvén waves, developed in the solar physics community (e.g. Sakurai et al. Solar Phys., vol. 133, 1991, pp. 227–245; Goossens et al. Solar Phys., vol. 157, 1995, pp. 75–102). Similar to the hydrodynamic theory, the mechanism of the MHD states can be explained by the successive interaction of the roll, streak and wave fields, which are now defined both for the hydrodynamic and magnetic fields. The derivation of this ‘vortex/Alfvén wave interaction’ state is rather straightforward as the scalings for both of the hydrodynamic and magnetic fields are identical. It turns out that the leading-order magnetic field of the asymptotic states appears only when a small external magnetic field is present. However, it does not mean that purely shear-driven dynamos are not possible. In fact, the second class of ‘self-sustained shear-driven dynamo theory’ shows a magnetic generation that is slightly smaller in size in the absence of any external field. Despite its small size, the magnetic field causes the novel feedback mechanism in the velocity field through resonant absorption, wherein the magnetic wave becomes more strongly amplified than the hydrodynamic counterpart.

An important class of canonical problems that is employed in quantifying the slipperiness of microstructured superhydrophobic surfaces is concerned with the calculation of the hydrodynamic loads on adjacent solid bodies whose size is large relative to the microstructure period. The effect of superhydrophobicity is most pronounced when the latter period is comparable to the separation between the solid probe and the superhydrophobic surface. We address the above distinguished limit, considering a simple configuration where the superhydrophobic surface is formed by a periodically grooved array, in which air bubbles are trapped in a Cassie state, and the solid body is an infinite cylinder. In the present part, we consider the case where the grooves are aligned perpendicular to the cylinder and allow for three modes of rigid-body motion: rectilinear motion perpendicular to the surface; rectilinear motion parallel to the surface, in the groove direction; and angular rotation about the cylinder axis. In this scenario, the flow is periodic in the direction parallel to the axis. Averaging over the small-scale periodicity yields a modified lubrication description where the small-scale details are encapsulated in two auxiliary two-dimensional cell problems which respectively describe pressure- and boundary-driven longitudinal flow through an asymmetric rectangular domain, bounded by a compound surface from the bottom and a no-slip surface from the top. Once the integral flux and averaged shear stress associated with each of these cell problems are calculated as a function of the slowly varying cell geometry, the hydrodynamic loads experienced by the cylinder are provided as quadratures of nonlinear functions of the latter distributions over a continuous sequence of cells.

For wave loads on cylinders constituting a long but finite array in the presence of incident waves, variations in the magnitude of the load with the non-dimensional wavenumber exhibit interesting features. Towering spikes and nearby secondary peaks (troughs) associated with trapped modes have been studied extensively. Larger non-trapped regions other than these two are termed Region III in this study. Studies of Region III are rare. We find that fluctuations in Region III are regular; the horizontal distance between two adjacent local maximum/minimum points, termed fluctuation spacing, is constant and does not change with non-dimensional wavenumbers. Fluctuation spacing is related only to the total number of cylinders in the array, identification serial number of the cylinder concerned and wave incidence angle. Based on the interaction theory and constructive/destructive interference, we demonstrate that the fluctuation characteristics can be predicted using simple analytical formulae. The formulae for predicting fluctuation spacing and the abscissae of every peak and trough in Region III are proposed. We reveal the intrinsic mechanism of the fluctuation phenomenon. When the diffraction waves emitted from the cylinders at the ends of the array and the cylinder concerned interfere constructively/destructively, peaks/troughs are formed. The fluctuation phenomenon in Region III is related to solutions of inhomogeneous equations. By contrast, spikes and secondary peaks are associated with solutions of the eigenvalue problem. This study of Region III complements existing understanding of the characteristics of the magnitude of wave load. The engineering significances of the results are discussed as well.

Continents exert a thermal blanket effect to the mantle underneath by locally accumulating heat and modifying the flow structures, which in turn affects continent motion. This dynamic feedback is studied numerically with a simplified model of an insulating plate over a thermally convecting fluid with infinite Prandtl number at Rayleigh numbers of the order of $10^{6}$. Several plate–fluid coupling modes are revealed as the plate size varies. In particular, small plates show long durations of stagnancy over cold downwellings. Between long stagnancies, bursts of velocity are observed when the plate rides on a single convection cell. As plate size increases, the coupled system transitions to another type of short-lived stagnancy, in which case hot plumes develop underneath. For an even larger plate, a unidirectional moving mode emerges as the plate modifies impeding flow structures it encounters while maintaining a single convection cell underneath. These identified modes are reminiscent of some real cases of continent–mantle coupling. Results show that the capability of a plate to overcome impeding flow structures increases with plate size, Rayleigh number and intensity of internal heating. Compared to cases with a fixed plate, cases with a freely drifting plate are associated with higher Nusselt number and more convection cells within the flow domain.

In this paper, the lateral migration of a neutrally buoyant spherical particle in the pressure-driven rectangular channel flow of an Oldroyd-B fluid is numerically investigated with a fictitious domain method. The aspect ratio of the channel cross-section considered is 1 and 2, respectively. The particle lateral motion trajectories are shown for the bulk Reynolds number ranging from 1 to 100, the ratio of the solvent viscosity to the total viscosity being 0.5, and a Weissenberg number up to 1.5. Our results indicate that the lateral equilibrium positions located on the cross-section midline, diagonal line, corner and channel centreline occur successively as the fluid elasticity is increased, for particle migration in square channel flow with finite fluid inertia. The transition of the equilibrium position depends strongly on the elasticity number (the ratio of the Weissenberg number to the Reynolds number) and weakly on the Reynolds number. The diagonal-line equilibrium position occurs at an elasticity number ranging from roughly 0.001 to 0.02, and can coexist with the midline and corner equilibrium positions. When the fluid inertia is negligibly small, particles migrate towards the channel centreline, or the closest corner, depending on their initial positions and the Weissenberg number, and the corner attractive area first increases and then decreases as the Weissenberg number increases. For particle migration in a rectangular channel with an aspect ratio of 2, the transition of the equilibrium position from the midline, ‘diagonal line’ (the line where two lateral shear rates are equal to each other), off-centre long midline and channel centreline takes place as the Weissenberg number increases at moderate Reynolds numbers. An off-centre equilibrium position on the long midline is observed for a large blockage ratio of 0.3 (i.e. the ratio of the particle diameter to the channel height is 0.3) at a low Reynolds number. This off-centre migration is driven by shear forces, unlike the elasticity-induced rapid inward migration, which is driven by the normal force (pressure or first normal stress difference).

In this paper, we present computational results of some two-fold azimuthally symmetric travelling waves and their stability. Calculations over a range of Reynolds numbers ($Re$) reveal connections between a class of solutions computed by Wedin & Kerswell (J. Fluid Mech., vol. 508, 2004, pp. 333–371) (henceforth called the WK solution) and the $Re\rightarrow \infty$ vortex–wave interaction theory of Hall & Smith (J. Fluid Mech., vol. 227, 1991, pp. 641–666) and Hall & Sherwin (J. Fluid Mech., vol. 661, 2010, pp. 178–205). In particular, the continuation of the WK solutions to larger values of $Re$ shows that the WK solution bifurcates from a shift-and-rotate symmetric solution, which we call the WK2 state. The WK2 solution computed for $Re\leqslant 1.19\times 10^{6}$ shows excellent agreement with the theoretical $Re^{-5/6}$, $Re^{-1}$ and $O(1)$ scalings of the waves, rolls and streaks respectively. Furthermore, these states are found to have only two unstable modes in the large $Re$ regime, with growth rates estimated to be $O(Re^{-0.42})$ and $O(Re^{-0.92})$, close to the theoretical $O(Re^{-1/2})$ and $O(Re^{-1})$ asymptotic results for edge and sinuous instability modes of vortex–wave interaction states (Deguchi & Hall, J. Fluid Mech., vol. 802, 2016, pp. 634–666) in plane Couette flow. For the nonlinear viscous core states (Ozcakir et al., J. Fluid Mech., vol. 791, 2016, pp. 284–328), characterized by spatial a shrinking of the wave and roll structure towards the pipe centre with increasing $Re$, we continued the solution to $Re\leqslant 8\times 10^{6}$ and we find only one unstable mode in the large Reynolds number regime, with growth rate scaling as $Re^{-0.46}$ within the class of symmetry-preserving disturbances.

Wing shape is an important factor affecting the aerodynamic performance of wings of monocopters and flapping-wing micro air vehicles. Here, an evolutionary structural optimisation method is adapted to optimise wing shape to enhance the lift force due to aerodynamic pressure on the wing surfaces. The pressure distribution is observed to vary with the span-based Reynolds number over a range covering most insects and samaras. Accordingly, the optimised wing shapes derived using this evolutionary approach are shown to adjust with Reynolds number. Moreover, these optimised shapes exhibit significantly higher lift coefficients (${\sim}50\,\%$) than the initial rectangular wing forebear. Interestingly, the optimised shapes are found to have a large area outboard, broadly in line with the features of high-lift forewings of multi-winged insects. According to specific aerodynamic performance requirements, this novel method could be employed in the optimisation of improved wing shapes for micro air vehicles.

An analytical model is presented for the wake interaction tones produced by a contra-rotating propeller. We re-cast the usual far-field radiation formulae as a double integral over a nominal propeller source annulus. Assuming that the number of blades on both propellers is large, we evaluate the integral asymptotically in terms of its leading-order contributions from interior stationary or boundary critical points which represent the specific locations on the propeller annulus that dominate the sound radiation. The asymptotic approach is powerful producing results in the form of one-line algebraic formulae that contain no integrals or special functions yet remain accurate. The asymptotics show that sweep is not necessarily beneficial and can cause the blade design to become critical for particular tones and directions in terms of a continuum of interior points distributed along a line on the propeller source annulus producing a higher-order result and thus an enhanced radiated sound field. The paper also shows how the interior points are completely consistent with the sub- or super-critical gust response of a swept blade. Tones with low and zero azimuthal mode order are treated as special cases and the asymptotics show that, as the mode order reduces, the radiated sound becomes concentrated around the flight axis where even higher-order solutions are possible, including rings and annuli of stationary points around the propeller annulus. Full numerical calculations confirm the accuracy of the asymptotic approach.