We identify ‘minimal seeds’ for turbulence, i.e. initial conditions of the smallest possible total perturbation energy density $E_{c}$ that trigger turbulence from the laminar state, in stratified plane Couette flow, the flow between two horizontal plates of separation $2H$, moving with relative velocity $2{\rm\Delta}U$, across which a constant density difference $2{\rm\Delta}{\it\rho}$ from a reference density ${\it\rho}_{r}$ is maintained. To find minimal seeds, we use the ‘direct-adjoint-looping’ (DAL) method for finding nonlinear optimal perturbations that optimise the time-averaged total dissipation of energy in the flow. These minimal seeds are located adjacent to the edge manifold, the manifold in state space that separates trajectories which transition to turbulence from those which eventually decay to the laminar state. The edge manifold is also the stable manifold of the system’s ‘edge state’. Therefore, the trajectories from the minimal seed initial conditions spend a large amount of time in the vicinity of some states: the edge state; another state contained within the edge manifold; or even in dynamically slowly varying regions of the edge manifold, allowing us to investigate the effects of a stable stratification on any coherent structures associated with such states. In unstratified plane Couette flow, these coherent structures are manifestations of the self-sustaining process (SSP) deduced on physical grounds by Waleffe (Phys. Fluids, vol. 9, 1997, pp. 883–900), or equivalently finite Reynolds number solutions of the vortex–wave interaction (VWI) asymptotic equations initially derived mathematically by Hall & Smith (J. Fluid Mech., vol. 227, 1991, pp. 641–666). The stratified coherent states we identify at moderate Reynolds number display an altered form from their unstratified counterparts for bulk Richardson numbers $\mathit{Ri}_{B}=g{\rm\Delta}{\it\rho}H/({\it\rho}_{r}{\rm\Delta}U^{2})=O(\mathit{Re}^{-1})$, and exhibit chaotic motion for larger $\mathit{Ri}_{B}$. We demonstrate that at hith Reynolds number the suppression of vertical motions by stratification strongly disrupts input from the waves to the roll velocity structures, thus preventing the waves from reinforcing the viscously decaying roll structures adequately, when $\mathit{Ri}_{B}=O(\mathit{Re}^{-2})$.

Steady flows that optimize heat transport are obtained for two-dimensional Rayleigh–Bénard convection with no-slip horizontal walls for a variety of Prandtl numbers $\mathit{Pr}$ and Rayleigh number up to $\mathit{Ra}\sim 10^{9}$. Power-law scalings of $\mathit{Nu}\sim \mathit{Ra}^{{\it\gamma}}$ are observed with ${\it\gamma}\approx 0.31$, where the Nusselt number $\mathit{Nu}$ is a non-dimensional measure of the vertical heat transport. Any dependence of the scaling exponent on $\mathit{Pr}$ is found to be extremely weak. On the other hand, the presence of two local maxima of $\mathit{Nu}$ with different horizontal wavenumbers at the same $\mathit{Ra}$ leads to the emergence of two different flow structures as candidates for optimizing the heat transport. For $\mathit{Pr}\lesssim 7$, optimal transport is achieved at the smaller maximal wavenumber. In these fluids, the optimal structure is a plume of warm rising fluid, which spawns left/right horizontal arms near the top of the channel, leading to downdraughts adjacent to the central updraught. For $\mathit{Pr}>7$ at high enough $\mathit{Ra}$, the optimal structure is a single updraught lacking significant horizontal structure, and characterized by the larger maximal wavenumber.

It is currently known that Tollmien–Schlichting (TS) waves can be attenuated by the introduction of spanwise mean velocity gradients in an otherwise two-dimensional boundary layer (BL). The stabilizing effect, associated with an extra turbulence production term, is strong enough to obtain a delay in transition to turbulence induced by TS waves, with the implication of reducing skin-friction drag. Miniature vortex generators (MVGs), mounted in an array, have successfully been used to obtain velocity modulations by the generation of alternating high- and low-speed streaks in the spanwise direction to control the BL. Experimentally, an initial amplification of the TS waves has been reported, which takes place in the near-wake region of the MVG array. The higher the streak amplitude, the stronger the downstream stabilizing effect becomes, but with the drawback of experiencing an even stronger initial amplification. This can lead to a sub-critical transitional Reynolds number, which would not only mean that the control has failed but, even worse, also lead to an advancement of the transition location. Here, direct numerical simulations and a local spatial stability analysis have been performed in order to reach a deeper understanding of this behaviour. The results agree well with experiments and we propose an explanation of the described behavior in terms of stability properties of the controlled BL. This important knowledge can be used in future designs of BL modulators, which can lead to improved stability of the control and to an extended region of laminar flow.

We investigate the two-dimensional (2D) instability recently discussed by Gallet et al. (Phys. Fluids, vol. 22, 2010, 034105) and Ilin & Morgulis (J. Fluid Mech., vol. 730, 2013, pp. 364–378) which arises when a radial cross-flow is imposed on a centrifugally stable swirling flow. By finding a simpler rectilinear example of the instability – a sheared half-plane, the minimal ingredients for the instability are identified and the destabilising/stabilising effect of inflow/outflow boundaries clarified. The instability – christened ‘boundary inflow instability’ here – is of critical layer type where this layer is either at the inflow wall and the growth rate is $O(\sqrt{{\it\eta}})$ (as found by Ilin & Morgulis (J. Fluid Mech., vol. 730, 2013, pp. 364–378)), or in the interior of the flow and the growth rate is $O({\it\eta}\log 1/{\it\eta})$, where ${\it\eta}$ measures the (small) inflow-to-tangential-flow ratio. The instability is robust to changes in the rotation profile, even to those which are very Rayleigh-stable, and the addition of further physics such as viscosity, three-dimensionality and compressibility, but is sensitive to the boundary condition imposed on the tangential velocity field at the inflow boundary. Providing the vorticity is not fixed at the inflow boundary, the instability seems generic and operates by the inflow advecting vorticity present at the boundary across the interior shear. Both the primary bifurcation to 2D states and secondary bifurcations to 3D states are found to be supercritical. Assuming an accretion flow driven by molecular viscosity only, so ${\it\eta}=O(Re^{-1})$, the instability is not immediately relevant for accretion disks since the critical threshold is $O(Re^{-2/3})$ and the inflow boundary conditions are more likely to be stress-free than non-slip. However, the analysis presented here does highlight the potential for mass entering a disk to disrupt the orbiting flow if this mass flux possesses vorticity.

The present study describes, first, an efficient algorithm for computing solutions in terms of capillary–gravity solitary waves of the irrotational Euler equations with a free surface and, second, provides numerical evidences of the existence of an infinite number of generalised solitary waves (solitary waves with undamped oscillatory wings). Using conformal mapping, the unknown fluid domain, which is to be determined, is mapped into a uniform strip of the complex plane. In the transformed domain, a Babenko-like equation is then derived and solved numerically.

Direct numerical simulations of the fully developed turbulent flow through a porous square duct are performed to study the effect of the permeable wall on the secondary cross-stream flow. The volume-averaged Navier–Stokes equations are used to describe the flow in the porous phase, a packed bed with porosity ${\it\varepsilon}_{c}=0.95$. The porous square duct is computed at $\mathit{Re}_{b}\simeq 5000$ and compared with the numerical simulations of a turbulent duct with four solid walls. The two boundary layers on the top wall and porous interface merge close to the centre of the duct, as opposed to the channel, because the sidewall boundary layers inhibit the growth of the shear layer over the porous interface. The most relevant feature in the porous duct is the enhanced magnitude of the secondary flow, which exceeds that of a regular duct by a factor of four. This is related to the increased vertical velocity, and the different interaction between the ejections from the sidewalls and the porous medium. We also report a significant decrease in the streamwise turbulence intensity over the porous wall of the duct (which is also observed in a porous channel), and the appearance of short spanwise rollers in the buffer layer, replacing the streaky structures of wall-bounded turbulence. These spanwise rollers most probably result from a Kelvin–Helmholtz type of instability, and their width is limited by the presence of the sidewalls.

Many micro-organisms find themselves immersed in fluids displaying non-Newtonian rheological properties such as viscoelasticity and shear-thinning viscosity. The effects of viscoelasticity on swimming at low Reynolds numbers have already received considerable attention, but much less is known about swimming in shear-thinning fluids. A general understanding of the fundamental question of how shear-thinning rheology influences swimming still remains elusive. To probe this question further, we study a spherical squirmer in a shear-thinning fluid using a combination of asymptotic analysis and numerical simulations. Shear-thinning rheology is found to affect a squirming swimmer in non-trivial and surprising ways; we predict and show instances of both faster and slower swimming depending on the surface actuation of the squirmer. We also illustrate that while a drag and thrust decomposition can provide insights into swimming in Newtonian fluids, extending this intuition to problems in complex media can prove problematic.

The influence of fixed temperature and fixed heat flux thermal boundary conditions on rapidly rotating convection in the plane layer geometry is investigated for the case of stress-free mechanical boundary conditions. It is shown that whereas the leading-order system satisfies fixed temperature boundary conditions implicitly, a double boundary layer structure is necessary to satisfy the fixed heat flux thermal boundary conditions. The boundary layers consist of a classical Ekman layer adjacent to the solid boundaries that adjust viscous stresses to zero, and a layer in thermal wind balance just outside the Ekman layers that adjusts the normal derivative of the temperature fluctuation to zero. The influence of these boundary layers on the interior geostrophically balanced convection is shown to be asymptotically weak, however. Upon defining a simple rescaling of the thermal variables, the leading-order reduced system of governing equations is therefore equivalent for both boundary conditions. These results imply that any horizontal thermal variation along the boundaries that varies on the scale of the convection has no leading-order influence on the interior convection, thus providing insight into geophysical and astrophysical flows where stress-free mechanical boundary conditions are often assumed.

We present experimental evidence that a minute amount of polymer additives can significantly enhance heat transport in the bulk region of turbulent thermal convection. The effects of polymer additives are found to be the enhancement of coherent heat fluxes and suppression of incoherent heat fluxes. The enhanced heat transport is associated with the increased coherency of thermal plumes, as a result of the suppression of small-scale turbulent fluctuations by polymers. The incoherent heat flux, arising from turbulent background fluctuations, makes no net contribution to heat transport. The fact that polymer additives can increase the coherency of thermal plumes is supported by the measurements of a number of local quantities, such as the extracted plume amplitude and width, the velocity autocorrelation functions and the velocity–temperature cross-correlation coefficient. The results from local measurements also suggest the existence of a threshold value for the polymer concentration, only above which significant modification of the plume coherent properties and enhancement of the local heat flux can be observed. Estimation of the plume emission rate suggests a stabilization of the thermal boundary layer by polymer additives.

We numerically investigate the effects of non-Newtonian fluid properties, including shear thinning and elasticity, on the locomotion of Taylor’s swimming sheet with arbitrary amplitude. Our results show that elasticity hinders the swimming speed, but a shear-thinning viscosity in the absence of elasticity enhances the speed. The combination of the two effects, modelled using a Giesekus constitutive equation, hinders the swimming speed. We find that the swimming speed of an infinitely long waving sheet in an inelastic shear-thinning fluid has a maximum, whose value depends on the sheet undulation amplitude and the fluid rheological properties. The power consumption, on the other hand, follows a universal scaling law.

Enhancement of heat transport across a fluid layer is of fundamental interest as well as great technological importance. For decades, Rayleigh–Bénard convection has been a paradigm for the study of convective heat transport, and how to improve its overall heat-transfer efficiency is still an open question. Here, we report an experimental and numerical study that reveals a novel mechanism that leads to much enhanced heat transport. When vertical partitions are inserted into a convection cell with thin gaps left open between the partition walls and the cooling/heating plates, it is found that the convective flow becomes self-organized and more coherent, leading to an unprecedented heat-transport enhancement. In particular, our experiments show that with six partition walls inserted, the heat flux can be increased by approximately 30 %. Numerical simulations show a remarkable heat-flux enhancement of up to 2.3 times (with 28 partition walls) that without any partitions.

Laboratory experiments involving downward air–helium fountains are presented. The large density differences between these releases and the ambient allow us to investigate how non-Boussinesq effects modify fountain heights and fountain fluctuations in comparison with the Boussinesq case (i.e. marginal density differences). In these experiments, the source Froude number is varied over a wide range covering (i) the very weak, (ii) the weak and (iii) the forced fountain regimes. It is shown that the classical Boussinesq correlations can be extended to the non-Boussinesq case provided that the Froude number is multiplied by the square root of the ratio between the released fluid density and that of the ambient. In the range investigated, no influence of the source Reynolds number is observed.