Large eddy simulation (LES) is used to study the evolution and structure of a lock-exchange, Boussinesq gravity current forming in a channel partially blocked by a porous layer. This configuration is used to understand how the characteristics of a surface layer containing floating vegetation affects the generation of thermally driven convective water exchange in a long shallow channel. The porous layer, which represents the roots of the floating vegetation, contains a staggered array of rigid square cylinders of edge length $D$ with solid volume fraction $\unicode[STIX]{x1D719}$. The cylinders extend over a depth $h_{1}<H$ below the free surface, where $H$ is the channel depth. The surface current of lighter fluid splits into two layers, one propagating slowly inside the porous layer and the other flowing beneath the porous layer. The main geometrical parameters of the porous layer, $\unicode[STIX]{x1D719}$ and $h_{1}/H$, have a large effect on the dynamics and structure of the surface current and the temporal variation of the front position. For cases with sufficiently large values of $h_{1}/H$ and $\unicode[STIX]{x1D719}$, the front within the porous layer approaches the triangular shape observed for low Reynolds number lock-exchange currents propagating in a channel containing cylinders over its whole volume ($h_{1}/H=1$), and the surface current transitions to a drag-dominated regime in which the front velocity is proportional to $t^{-1/4}$, where $t$ is the time since the current is initiated. For sufficiently high values of $\unicode[STIX]{x1D719}$, the velocity of the fluid inside the porous layer is close to zero at all locations except for those situated close to the lock gate and for some distance behind the front. Close to the front, lighter fluid from below penetrates into the porous layer due to unstable stratification at the bottom of the porous layer. Simulation results are also used to assess how $\unicode[STIX]{x1D719},h_{1}/H$ and the Reynolds number affect the rate at which the heavier fluid situated initially inside the porous layer is removed by the advancing surface current and the main mixing mechanisms. Based on the estimated time scales for flushing the porous (root) layer, we show that flushing can significantly enhance the overall rate of nutrient removal by the floating vegetation by maintaining a higher concentration of nutrients within the root layer.

In this paper a novel nonlinear feedback control design methodology for incompressible fluid flows aiming at the optimisation of long-time averages of flow quantities is presented. It applies to reduced-order finite-dimensional models of fluid flows, expressed as a set of first-order nonlinear ordinary differential equations with the right-hand side being a polynomial function in the state variables and in the controls. The key idea, first discussed in Chernyshenko et al. (Phil. Trans. R. Soc. Lond. A, vol. 372, 2014, 20130350), is that the difficulties of treating and optimising long-time averages of a cost are relaxed by using the upper/lower bounds of such averages as the objective function. In this setting, control design reduces to finding a feedback controller that optimises the bound, subject to a polynomial inequality constraint involving the cost function, the nonlinear system, the controller itself and a tunable polynomial function. A numerically tractable and efficient approach to the solution of such optimisation problems, based on sum-of-squares techniques and semidefinite programming, is proposed. To showcase the methodology, the mitigation of the fluctuation kinetic energy in the unsteady wake behind a circular cylinder in the laminar regime at $Re=100$, via controlled angular motions of the surface, is numerically investigated. A compact reduced-order model that resolves the long-term behaviour of the fluid flow and the effects of actuation, is first derived using proper orthogonal decomposition and Galerkin projection. In a full-information setting, feedback controllers are then designed to reduce the long-time average of the resolved kinetic energy associated with the limit cycle. These controllers are then implemented in direct numerical simulations of the actuated flow. Control performance, total energy efficiency and the physical control mechanisms identified are analysed in detail. Key elements of the methodology, implications and future work are finally discussed.

In the present paper we focus on deriving the modified Smoluchowski slip velocity of second-order fluids, for electroosmotic flows over plane surfaces with arbitrary non-uniform surface potential in the presence of thin electric double layers (EDLs). We employ matched asymptotic expansion to stretch the electric double layer and subsequently apply regular asymptotic expansions taking the Deborah number ($De$) as the gauge function. Modified slip velocities correct up to $O(De^{2})$ are presented. Two sample cases are considered to demonstrate the effects of viscoelasticity on slip velocity: (i) an axially periodic patterned potential and (ii) a step-change-like variation in the surface potential. The central result of our analysis is that, unlike Newtonian fluids, the electroosmotic slip velocity for second-order fluids does not, in general, align with the direction of the applied external electric field. Proceeding further forward, we show that the slip velocity in a given direction may, in fact, depend on the applied electric field strength in a mutually orthogonal direction, considering three dimensionality of the flow structure. In addition, we demonstrate that the modified slip velocity is not proportional to the zeta potential, as in the cases of Newtonian fluids; rather it depends strongly on the gradients of the interfacial potential as well. Our results are likely to have potential implications so far as the design of charge modulated microfluidic devices transporting rheologically complex fluids is concerned, such as for mixing and bio-reactive system analysis in lab-on-a-chip-based micro-total-analysis systems handling bio-fluids.

This paper presents measurements of mixing efficiency of the two counter-flowing gravity currents created by symmetric lock exchange in a channel. The novel feature of this work is that the buoyancy Reynolds number of the currents is higher than in previous experiments, so that the mixing is not significantly affected by viscosity. We find that the mixing efficiency asymptotes to 0.08 at high Reynolds numbers. We present a model of the mixing based on the evolution of idealized mean profiles of velocity and density at the interface between the two currents, the results of which are in good agreement with the measurements of mixing efficiency.

An analytical method for determining the shape of hollow vortices in shear flows is presented in detail. In a non-dimensional formulation, it is shown that the problem has one degree of freedom represented by the free choice of the non-dimensionalized speed $\unicode[STIX]{x1D705}$ at the boundary of the vortex. The solutions form two families of shapes corresponding to vortex circulation and shear-flow vorticity having the opposite or same sign. When the signs are opposite, the shape family resembles that described by Llewellyn Smith & Crowdy (J. Fluid Mech., vol. 691, 2012, pp. 178–200) for hollow vortices in a potential flow with strain. As for that flow, there is a minimum value of $\unicode[STIX]{x1D705}$ below which there is no solution as the boundary of the vortex self-intersects, while, when the signs are the same, solutions exist for $0<\unicode[STIX]{x1D705}$.

From rain drops landing on the ocean to inkjet printing, the impact of droplets onto moving liquid surfaces is a ubiquitous process in nature and in industry. A rich range of phenomena can arise. The behaviour depends on the inertia, the properties of the drops and the relative speeds in the impact zone. While the result ranges from tranquil coalescence to violent splashing, intermediate regimes also occur, including partial and complete bouncing and even ‘surfing’ of the droplet. These regimes are determined by the ratio of the drop and surface velocities and the liquid properties. A regime diagram can be constructed in which distinct dynamical regimes are clearly identified.

The effect of a uniform electric field on the motion of a drop in an unbounded plane Poiseuille flow is studied analytically. The drop and suspending media are considered to be Newtonian and leaky dielectric. We solve for the two-way coupled electric and flow fields analytically by using a double asymptotic expansion for small charge convection and small shape deformation. We obtain two important mechanisms of cross-stream migration of the drop: (i) shape deformation and (ii) charge convection. The second one is a new source of cross-stream migration of the drop in plane Poiseuille flow which is due to an asymmetric charge distribution on the drop surface. Our study reveals that charge convection can cause a spherical non-deformable drop to migrate in the cross-stream direction. The combined effect of charge convection and shape deformation significantly alters the drop velocity, drop trajectory and steady state transverse position of the drop. We predict that, depending on the orientation of the applied uniform electric field and the relevant drop/medium electrohydrodynamic parameters, the drop may migrate either towards the centreline of the flow or away from it. We obtain that the final steady state transverse position of the drop is independent of its initial transverse position in the flow field. Most interestingly, we show that the drop can settle in an off-centreline steady state transverse position. Two-dimensional numerical simulations are also performed to study the drop motion in the combined presence of plane Poiseuille flow and a tilted electric field. The drop trajectory and steady state transverse position of the drop obtained from numerical simulations are in qualitative agreement with the analytical results.

The selection and placement of actuators and sensors to control compressible viscous flows is addressed by developing a novel methodology based upon the eigensystem structural sensitivity of the linearized evolution operator appropriate for linear feedback control. Forward and adjoint global modes are used to construct a space of possible perturbations to the linearized operator, which yields a small optimization problem for selecting the parameters that best achieve the control objective, including where they should be placed. The method is demonstrated by informing actuation to suppress amplification of the instabilities in boundary layer separation in a high-subsonic diffuser. Complete stabilization is observed in the separated shear layer for short downstream distances at modest Reynolds number. Higher Reynolds numbers and longer distances are expected to be more challenging to stabilize; here the control informed by the procedure still substantively suppresses amplification of instabilities. It is also demonstrated that more complex actuator–sensor selections may not yield superior controllers.

The influence of near-wall density and viscosity gradients on near-wall turbulence in a channel is studied by means of direct numerical simulation of the low-Mach-number approximation of the Navier–Stokes equations. Different constitutive relations for density $\unicode[STIX]{x1D70C}$ and viscosity $\unicode[STIX]{x1D707}$ as a function of temperature are used in order to mimic a wide range of fluid behaviours and to develop a generalised framework for studying turbulence modulations in variable-property flows. Instead of scaling the velocity solely based on local density, as done for the van Driest transformation, we derive an extension of the scaling that is based on gradients of the semilocal Reynolds number, defined as $Re_{\unicode[STIX]{x1D70F}}^{\star }\equiv Re_{\unicode[STIX]{x1D70F}}\sqrt{(\overline{\unicode[STIX]{x1D70C}}/\overline{\unicode[STIX]{x1D70C}}_{w})}/(\overline{\unicode[STIX]{x1D707}}/\overline{\unicode[STIX]{x1D707}}_{w})$ (the bar and subscript $w$ denote Reynolds averaging and wall value respectively, while $Re_{\unicode[STIX]{x1D70F}}$ is the friction Reynolds number based on wall values). This extension of the van Driest transformation is able to collapse velocity profiles for flows with near-wall property gradients as a function of the semilocal wall coordinate. However, flow quantities like mixing length, turbulence anisotropy and turbulent vorticity fluctuations do not show a universal scaling very close to the wall. This is attributed to turbulence modulations, which play a crucial role in the evolution of turbulent structures and turbulence energy transfer. We therefore investigate the characteristics of streamwise velocity streaks and quasistreamwise vortices and find that, similarly to turbulence statistics, the turbulent structures are also strongly governed by $Re_{\unicode[STIX]{x1D70F}}^{\star }$ profiles and that their dependence on individual density and viscosity profiles is minor. Flows with near-wall gradients in $Re_{\unicode[STIX]{x1D70F}}^{\star }$ ($\text{d}Re_{\unicode[STIX]{x1D70F}}^{\star }/\text{d}y\neq 0$) show significant changes in inclination and tilting angles of quasistreamwise vortices. These structural changes are responsible for the observed modulation of the Reynolds stress generation mechanism and the inter-component energy transfer in flows with strong near-wall $Re_{\unicode[STIX]{x1D70F}}^{\star }$ gradients.

Nonlinear triadic interactions are at the heart of our understanding of turbulence. In flows where waves are present, modes must not only be in a triad to interact, but their frequencies must also satisfy an extra condition: the interactions that dominate the energy transfer are expected to be resonant. We derive equations that allow direct measurement of the actual degree of resonance of each triad in a turbulent flow. We then apply the method to the case of rotating turbulence, where eddies coexist with inertial waves. We show that for a range of wavenumbers, resonant and near-resonant triads are dominant, the latter allowing a transfer of net energy towards two-dimensional modes that would be inaccessible otherwise. The results are in good agreement with approximations often done in theories of rotating turbulence, and with the mechanism of parametric instability proposed to explain the development of anisotropy in such flows. We also observe that, at least for the moderate Rossby numbers studied here, marginally near-resonant and non-resonant triads play a non-negligible role in the coupling of modes.

A novel data-driven modal decomposition of fluid flow is proposed, comprising key features of proper orthogonal decomposition (POD) and dynamic mode decomposition (DMD). The first mode is the normalized real or imaginary part of the DMD mode that minimizes the time-averaged residual. The $N$th mode is defined recursively in an analogous manner based on the residual of an expansion using the first $N-1$ modes. The resulting recursive DMD (RDMD) modes are orthogonal by construction, retain pure frequency content and aim at low residual. Recursive DMD is applied to transient cylinder wake data and is benchmarked against POD and optimized DMD (Chen et al., J. Nonlinear Sci., vol. 22, 2012, pp. 887–915) for the same snapshot sequence. Unlike POD modes, RDMD structures are shown to have purer frequency content while retaining a residual of comparable order to POD. In contrast to DMD, with exponentially growing or decaying oscillatory amplitudes, RDMD clearly identifies initial, maximum and final fluctuation levels. Intriguingly, RDMD outperforms both POD and DMD in the limit-cycle resolution from the same snapshots. Robustness of these observations is demonstrated for other parameters of the cylinder wake and for a more complex wake behind three rotating cylinders. Recursive DMD is proposed as an attractive alternative to POD and DMD for empirical Galerkin models, in particular for nonlinear transient dynamics.

Cantilevered elastic sheets and rods immersed in a steady uniform flow are known to undergo instabilities that give rise to complex dynamics, including limit cycle behaviour and chaotic motion. Recent work has examined their stability in an inverted configuration where the flow impinges on the free end of the cantilever with its clamped edge downstream: this is commonly referred to as an ‘inverted flag’. Theory has thus far accurately captured the stability of wide inverted flags only, i.e. where the dimension of the clamped edge exceeds the cantilever length; the latter is aligned in the flow direction. Here, we theoretically examine the stability of slender inverted flags and rods under steady uniform flow. In contrast to wide inverted flags, we show that slender inverted flags are never globally unstable. Instead, they exhibit bifurcation from a state that is globally stable to multiple equilibria of varying stability, as flow speed increases. This theory is compared with new and existing measurements on slender inverted flags and rods, where excellent agreement is observed. The findings of this study have significant implications to investigations of biological phenomena such as the motion of leaves and hairs, which can naturally exhibit a slender geometry with an inverted configuration.

For thermal mixed-convection flows, the Nusselt number is a function of Reynolds number, Grashof number and the angle between the forced- and natural-convection directions. We consider flow over a heated cylinder for which there is no universal correlation that accurately predicts Nusselt number as a function of these parameters, especially in opposing-convection flows, where the natural convection is against the forced convection. Here, we revisit this classical problem by employing modern tools from machine learning to develop a general multi-fidelity framework for constructing a stochastic response surface for the Nusselt number. In particular, we combine previously developed experimental correlations (low-fidelity model) with direct numerical simulations (high-fidelity model) using Gaussian process regression and autoregressive stochastic schemes. In this framework the high-fidelity model is sampled only a few times, while the inexpensive empirical correlation is sampled at a very high rate. We obtain the mean Nusselt number directly from the stochastic multi-fidelity response surface, and we also propose an improved correlation. This new correlation seems to be consistent with the physics of this problem as we correct the vectorial addition of forced and natural convection with a pre-factor that weighs differently the forced convection. This, in turn, results in a new definition of the effective Reynolds number, hence accounting for the ‘incomplete similarity’ between mixed convection and forced convection. In addition, due to the probabilistic construction, we can quantify the uncertainty associated with the predictions. This information-fusion framework is useful for elucidating the physics of the flow, especially in cases where anomalous transport or interesting dynamics may be revealed by contrasting the variable fidelity across the models. While in this paper we focus on the thermal mixed convection, the multi-fidelity framework provides a new paradigm that could be used in many different contexts in fluid mechanics including heat and mass transport, but also in combining various levels of fidelity of models of turbulent flows.

We consider free-surface flow over two localised obstacles using the framework of the forced Korteweg–de Vries equation in a suite of numerical simulations. Our main concern is with the transcritical regime when the oncoming flow has a Froude number close to unity. The flow behaviour can be characterised by the Froude number and the maximum heights of the obstacles. In the transcritical regime at early times, undular bores are produced upstream and downstream of each obstacle. Our main aim is to describe the interaction of these undular bores between the obstacles, and to find the outcome at very large times. We find that the flow development can be defined in three stages. The first stage is described by the well-known development of undular bores upstream and downstream of each obstacle. The second stage is the interaction between the undular bore moving downstream from the first obstacle and the undular bore moving upstream from the second obstacle. The third stage is the very large time evolution of this interaction, when one of the obstacles controls criticality. For equal obstacle heights, our analytical and numerical results indicate that either one of the obstacles can control flow criticality, that being the first obstacle when the flow is slightly subcritical and the second obstacle otherwise. For unequal obstacle heights the larger obstacle controls criticality. The results obtained here complement a recent numerical study using the fully nonlinear, but non-dispersive, shallow water equations.

Simultaneous measurements of the phase-averaged velocity distribution and the underlying discharge quantities of a dielectric barrier discharge plasma actuator (PA) are performed at $10~\text{k}\text{Hz}$ discharge frequency to investigate the interplay of the discharge and the surrounding flow. The underlying velocity information for the force estimation is obtained by means of phase-averaged particle image velocimetry; the discharge quantities are determined from a Lissajous-figure analysis. The results uncover a clear cause–effect relation between the phase-dependent effective discharge capacitance of the PA and the resulting spatiotemporal volume-force distributions. From this novel insight, it must be concluded that the instantaneous effective discharge intensity dominates the momentum-transfer rate rather than the formerly assumed operating voltage.

We analyse velocity fluctuations inside coherent vortices generated as a result of the inverse cascade in the two-dimensional (2-D) turbulence in a finite box. As we demonstrated in Kolokolov & Lebedev (Phys. Rev. E, vol. 93, 2016, 033104), the universal velocity profile, established in Laurie et al. (Phys. Rev. Lett., vol. 113, 2014, 254503), corresponds to the passive regime of the flow fluctuations. This property enables one to calculate correlation functions of the velocity fluctuations in the universal region. We present the results of the calculations that demonstrate a non-trivial scaling of the structure function. In addition the calculations reveal strong anisotropy of the structure function.

A theoretical model is developed to predict the maximum spreading of liquid metal drops when impacting onto dry surfaces under the influence of a vertical magnetic field. This model, which is constructed based on the energy conversion principle, agrees very well with the numerical results, covering a wide range of impact speeds, contact angles and magnetic strengths. When there is no magnetic field, we found that the maximum spreading factor can be predicted well by an interpolating scheme between the viscous and capillary effects, as proposed by Laan et al. (Phys. Rev. Appl., vol. 2 (4), 2014, 044018). However, when gradually increasing the magnetic field strength, the induced Lorentz forces are dominant over the viscous and capillary forces, taking the spreading behaviour into the ‘Joule regime’, where the Joule dissipation is significant. For most situations of practical interest, namely when the strength of the magnetic field is less than 3 T, all three energy conversion routes are important. Therefore, we determine the correct scaling behaviours for the magnetic influence by first equating the loss of kinetic energy to the Joule dissipation in the Joule regime, then by interpolating it with the viscous dissipation and the capillary effects, which allows for a universal rescaling. By plotting the numerical results against the theoretical model, all the results can be rescaled onto a single curve regardless of the materials of the liquid metals or the contact angles of the surfaces, proving that our theoretical model is correct in predicting the maximum spreading factor by constructing a balanced formula between kinetic energy, capillary energy, viscous dissipation energy and Joule dissipation energy.