This paper quantifies the instantaneous form of large-scale turbulent structures in canonical smooth-wall boundary layers, demonstrating that they adhere to a form that is consistent with the self-sustaining streak instability model suggested by Flores & Jiménez (Phys. Fluids, vol. 22, 2010, 071704) and Hwang & Cossu (Phys. Fluids, vol. 23, 2011, 061702). Our motivation for this study stems from previous observations of large-scale streaks that have been spatially locked in position within spanwise-heterogeneous boundary layers. Here, using similar tools, we demonstrate that the randomly occurring large-scale structures in canonical layers show similar behaviour. Statistically, we show that the signature of large-scale coherent structures exhibits increasing meandering behaviour with distance from the wall. At the upper edge of the boundary layer, where these structures are severely misaligned from the main-flow direction, the induced velocities associated with the strongly yawed vortex packets/clusters yield a significant spanwise-velocity component leading to an apparent oblique coherence of spanwise-velocity fluctuations. This pronounced meandering behaviour also gives rise to a dominant streamwise periodicity at a wavelength of approximately $6\unicode[STIX]{x1D6FF}$. We further statistically show that the quasi-streamwise roll-modes formed adjacent to these very large wavy motions are often one-sided (spanwise asymmetric), in stark contrast to the counter-rotating form suggested by conventional conditionally averaged representations. To summarise, we sketch a representative picture of the typical large-scale structures based on the evidence gathered in this study.

We investigate the stability of the Prandtl model for katabatic slope flows using both linear stability theory and direct numerical simulations. Starting from Prandtl’s analytical solution for uniformly cooled laminar slope flows, we use linear stability theory to identify the onset of instability and features of the most unstable modes. Our results show that the Prandtl model for parallel katabatic slope flows is prone to transverse and longitudinal modes of instability. The transverse mode of instability manifests itself as stationary vortical flow structures aligned in the along-slope direction, whereas the longitudinal mode of instability emerges as waves propagating in the base-flow direction. Beyond the stability limits, these two modes of instability coexist and form a complex flow structure crisscrossing the plane of flow. The emergence of a particular form of these instabilities depends strongly on three dimensionless parameters, which are the slope angle, the Prandtl number and a newly introduced stratification perturbation parameter, which is proportional to the relative importance of the disturbance to the background stratification due to the imposed surface buoyancy flux. We demonstrate that when this parameter is sufficiently large, then the stabilising effect of the background stratification can be overcome. For shallow slopes, the transverse mode of instability emerges despite meeting the Miles–Howard stability criterion of $Ri>0.25$. At steep slope angles, slope flow can remain linearly stable despite attaining Richardson numbers as low as $3\times 10^{-3}$.

This study investigates the stability and transition of Görtler vortices in a hypersonic boundary layer using linear stability theory and direct numerical simulations. In the simulations, Görtler vortices are separately excited by wall blowing and suction with spanwise wavelengths of 3, 6 and 9 mm. In addition to primary streaks with the same wavelength as the blowing and suction, secondary streaks with half the wavelength also emerge in the 6 and 9 mm cases. The streaks develop into mushroom structures before breaking down. The breakdown processes of the three cases are dominated by a sinuous-mode instability, a varicose-mode instability and a combination of the two, respectively. Both fundamental and subharmonic instabilities are relevant in all cases. Multiple modes are identified in the secondary-instability stage, some of which originate from the primary instabilities (first and second Mack modes). We demonstrate that the first Mack mode can be destabilized to either a varicose-mode or sinuous-mode streak instability depending on its frequency and wavelength, whereas the second Mack mode undergoes a stabilizing stage before turning into a varicose mode in the 6 and 9 mm cases. An energy analysis reveals the stabilizing and destabilizing mechanisms of the primary instabilities under the influence of Görtler vortices, highlighting the role played by the spanwise production based on the spanwise gradient of the streamwise velocity in both varicose and sinuous modes. The effects introduced by the secondary streaks are examined by filtering the secondary streaks in two new simulations with nominally identical conditions to those of the 6 and 9 mm cases. Remarkably, the secondary streaks can destabilize the Görtler vortices, therefore advancing the transition. The stability theory results are in good agreement with those from direct numerical simulations.

We carry out an experimental study of the granular surface flow of nearly monodisperse glass beads on a conical heap formed on a rough circular disc by a narrow stream of the particles from a hopper, with the pouring point displaced from the centre of the disc. During the growth phase, an axisymmetric heap is formed, which grows either by periodic avalanches or by non-periodic avalanches that occur randomly over the azimuthal location of the heap, depending on the operating conditions and system properties. The dynamics of heap growth is characterized by the variation of the heap height, angle of repose and the angular velocity of the periodic avalanche with time, for different mass flow rates from the hopper. When the base of the heap reaches the edge of the disc closest to the pouring point, the heap stops growing and a steady surface flow of particles is developed on the heap surface, with particles flowing over the edge of the disc into a collection tray. The geometry is a unique example of a granular flow on an erodible bed without any bounding side walls. The corresponding steady state geometry of the asymmetric heap is characterized by means of surface contours and angles of repose. The streamwise and transverse surface velocities are measured using high-speed video photography and image analysis for different mass flow rates. The flowing layer thickness is measured by immersing a coated needle in the flow at different positions on the mid-line of the flow. The surface angle of the flowing layer is found to be significantly smaller than the angle of repose and to be independent of the mass flow rate. The velocity profiles at different streamwise positions for different mass flow rates are found to be geometrically similar and are well described by Gaussian functions. The flowing layer thickness is calculated from a model using the measured surface velocities. The predicted values match the measured values quite well.

Supersonic turbulent flows at Mach 2.7 over concave surfaces for two different radii of curvature were investigated and compared with a flat plate turbulent boundary layer using direct numerical simulations. The streamwise velocity reduces in the outer part of the boundary layer due to compression, while it increases near the wall due to curvature, with a higher shape factor for the concave cases. The near-wall spanwise streak spacing reduces compared to the flat plate, with large-scale streaks and turbulence amplification also observed. Streamwise velocity iso-surfaces and streamlines show the generation of Görtler-like vortices, consistent with significant centrifugal effects. Abundant small vortices are shown to be associated with large baroclinic production of vorticity that is caused by the density and pressure gradients that are associated with concave compression. Profiles of turbulent kinetic energy and turbulent Mach number exhibit a characteristic two-layer structure in the concave boundary layer cases. In the outer layer, turbulence is greatly amplified, whereas a local balance exists in the inner layer. Turbulent energy budget analysis shows that both production and dissipation increase near the concave wall, whereas in the outer part of the boundary layer, the production is increased and ultimately balanced by convection and turbulent transport.

The Green–Naghdi equations are an extension of the shallow-water equations that capture the effects of finite fluid depth at arbitrary order in the characteristic height to width aspect ratio $H/L$. The shallow-water equations capture these effects to first order only, resulting in a relatively simple two-dimensional fluid-dynamical model for the layer horizontal velocity and depth. The Green–Naghdi equations, like the shallow-water equations, are two-dimensional fluid equations expressing momentum and mass conservation. There are different ‘levels’ of the Green–Naghdi equations of rapidly increasing complexity. In the present paper we focus on the behaviour of the lowest-level Green–Naghdi equations for a rotating shallow fluid layer, paying close attention to the flow structure at small spatial scales. We compare directly with the shallow-water equations and study the differences arising in their solutions. By recasting the equations into a form which both explicitly conserves Rossby–Ertel potential vorticity and represents the leading-order departure from geostrophic–hydrostatic balance, we are able to accurately describe both the ‘slow’ predominantly sub-inertial balanced dynamics and the ‘fast’ residual imbalanced dynamics. This decomposition has proved fruitful in studies of shallow-water dynamics but appears not to have been used before in studies of Green–Naghdi dynamics. Importantly, we find that this decomposition exposes a fundamental inconsistency in the Green–Naghdi equations for horizontal scales less than the mean fluid depth, scales for which the Green–Naghdi equations are supposed to more accurately model. Such scales exhibit pronounced activity compared to the shallow-water equations, and in particular spectra for certain fields like the divergence are flat or rising at high wavenumbers. This indicates a lack of convergence at small scales, and is also consistent with the poor convergence of total energy with resolution compared to the shallow-water equations. We suggest a mathematical reformulation of the Green–Naghdi equations which may improve convergence at small scales.

This paper presents an analytic solution for gust–aerofoil interaction noise for flat plates with spanwise-varying periodic leading edges in uniform mean flow. The solution is obtained by solving the linear inviscid equations via separation of variables and the Wiener–Hopf technique, and is suitable for calculating the far-field noise generated by any leading edge with a single-valued piecewise linear periodic spanwise geometry. Acoustic results for homogeneous isotropic turbulent flow are calculated by integrating the single-gust solution over a wavenumber spectrum. The far-sound pressure level is calculated for five test-case geometries; sawtooth serration, slitted $v$-root, slitted $u$-root, chopped peak and square wave, and compared to experimental measurements. Good agreement is seen over a range of frequencies and tip-to-root ratios (varying the sharpness of the serration). The analytic solution is then used to calculate the propagating pressure along the leading edge of the serration for fixed spanwise wavenumbers, i.e. only the contribution to the surface pressure which propagates to the far field. Using these results, two primary mechanisms for noise reduction are discussed; tip and root interference, and a redistribution of energy from cuton modes to cutoff modes. A secondary noise-reduction mechanism due to nonlinear features is also discussed and seen to be particularly important for leading edges with very narrow slits.

Application of an electric field across the pressure-driven stratified flow of a pair of miscible fluids inside a microchannel manifests interesting electrohydrodynamic (EHD) instabilities. Experiments uncover distinctive instability regimes with an increase in electric field Rayleigh number ($Ra^{\unicode[STIX]{x1D713}}$) – a linear-onset regime, a time-periodic nonlinear regime analogous to the von Kármán vortex street in the downstream and a regime with coherent flow patterns. The experiments also reveal that such linear and nonlinear instabilities can be stimulated non-invasively in a microchannel to mix or de-mix fluids simply by turning the electric field on or off, indicating the suitability of the process for on-demand micromixing. The characteristics of these instabilities have been theoretically investigated with the help of an Orr–Sommerfeld framework, which discloses the possibility of five distinctive finite-wavenumber modes for the instability. The EHD stresses originating due to the application of electric field stimulate a pair of shorter-wavelength electric field modes beyond a critical value of $Ra^{\unicode[STIX]{x1D713}}$. Increase in the levels of charge injection and EHD stresses lower the critical $Ra^{\unicode[STIX]{x1D713}}$ of these modes. The relatively longer-wavelength viscous mode is found to appear when the viscosity stratification between the fluid layers is high. Beyond a threshold Schmidt number ($Sc$), a diffusive mode is also found to appear near the mixed interfacial region. A thinner interface between the fluids at a higher $Sc$ helps this mode to behave as the interfacial mode of immiscible fluids. Contrast of ionic mobility in the fluids leads to the appearance of the K-mode of instability at much shorter wavelengths. The reported phenomena can be of significance in the domains of microscale mixing, pumping, heat exchange, mass transfer and reaction engineering.

Interactions between an oblique shock wave generated by a sharp fin placed on a cylindrical surface and the incoming boundary layer are investigated to unravel the mean features of the resulting shock/boundary layer interaction (SBLI) unit. This fin-on-cylinder SBLI unit has several unique features caused by the three-dimensional (3-D) relief offered by the cylindrical surface that noticeably alter the shock structure. Complementary experimental and computational studies are made to delineate both the surface and off-body flow features of the fin-on-cylinder SBLI unit and to obtain a detailed understanding of the mechanisms that dictate the mean flow and wall pressure features of the SBLI unit. Results show that the fin-on-cylinder SBLI exhibits substantial deviation from quasi-conical symmetry that is observed in planar fin SBLI. Furthermore, the separated flow growth rate appears to decrease with downstream distance and the separation size is consistently smaller than the planar fin SBLI with the same inflow and fin configurations. The causes for the observed diminution of the separated flow and its downstream growth rate were investigated in the light of changes caused by the cylinder curvature on the inviscid as well as separation shock. It was found that the inviscid shock gets progressively weakened in the region close to the triple point with downstream distance due to the 3-D relief effect from cylinder curvature. This weakening of the inviscid shock feeds into the separation shock, which is also independently impacted by the 3-D relief, to result in the observed modifications in the fin-on-cylinder SBLI unit.

Parametric subharmonic instability arising in a narrow-band wave spectrum is investigated. Using a statistical equation that describes weakly nonlinear interactions in a random wave field, we perform analytical and numerical stability analyses for a modulating wave train. The analytically obtained growth rate $\unicode[STIX]{x1D706}=(-\unicode[STIX]{x1D707}+\sqrt{\unicode[STIX]{x1D707}^{2}+4CE_{B}})/2$ agrees favourably with the results from direct numerical experiments, where $\unicode[STIX]{x1D707}$ is the half-value width of the background wave frequency spectrum, $E_{B}$ is the background wave energy density, and $C$ is a constant. This expression has two asymptotic limits: $\unicode[STIX]{x1D706}\sim \sqrt{CE_{B}}$ for $\unicode[STIX]{x1D707}\ll \sqrt{CE_{B}}$ and $\unicode[STIX]{x1D706}\sim CE_{B}/\unicode[STIX]{x1D707}$ for $\unicode[STIX]{x1D707}\gg \sqrt{CE_{B}}$. In the terms of weak turbulence, these two growth rates correspond to the ones occurring in the dynamic and kinetic time scales. In this way, our formulation successfully unifies the two conventional types of parametric subharmonic instability and offers a new criterion to determine the applicability of the classical kinetic equation in three-wave systems.

We present the first application of an artificial neural network trained through a deep reinforcement learning agent to perform active flow control. It is shown that, in a two-dimensional simulation of the Kármán vortex street at moderate Reynolds number ($Re=100$), our artificial neural network is able to learn an active control strategy from experimenting with the mass flow rates of two jets on the sides of a cylinder. By interacting with the unsteady wake, the artificial neural network successfully stabilizes the vortex alley and reduces drag by approximately 8 %. This is performed while using small mass flow rates for the actuation, of the order of 0.5 % of the mass flow rate intersecting the cylinder cross-section once a new pseudo-periodic shedding regime is found. This opens the way to a new class of methods for performing active flow control.

Recent experiments on thin films flowing down a vertical fibre with varying nozzle diameters present a wealth of new dynamics that illustrate the need for more advanced theory. We present a detailed analysis using a full lubrication model that includes slip boundary conditions, nonlinear curvature terms and a film stabilization term. This study brings to focus the presence of a stable liquid layer playing an important role in the full dynamics. We propose a combination of these physical effects to explain the observed velocity and stability of travelling droplets in the experiments and their transition to isolated droplets. This is also supported by stability analysis of the travelling wave solution of the model.

The primary and secondary laminar flows in annuli with longitudinal grooves and driven by pressure gradients have been analysed. There exist geometric configurations reducing pressure losses in primary flows in spite of an increase of the wall wetted area. The parameter ranges when such flows exist have been determined using linear stability theory. Two types of secondary flows have been identified. The first type has the form of the classical travelling waves driven by shear and modified by the grooves. The axisymmetric waves dominate for sufficiently large radii of the annuli while different spiral waves dominate for small radii. The secondary flow topology is unique in the former case and has the form of axisymmetric rings propagating in the axial direction. Topologies in the latter case are not unique, as spiral waves with left and right twists can emerge under the same conditions, resulting in flow structures varying from spatial rings to rhombic forms. The most intense motion of this type occurs near the walls. The second type of secondary flow has the form of travelling waves driven by inertial effects with characteristics very distinct from the shear waves. Its critical Reynolds number increases proportionally to $S^{-2}$, where $S$ denotes the groove amplitude, while the amplification rates increase proportionally to $S^{2}$. These waves exist only if $S$ is above a well-defined minimum and their axisymmetric forms dominate, with the most intense motion occurring near the annulus mid-section. Geometries that give preference to the latter waves have been identified. It is shown that the drag-reducing topographies stabilize the classical travelling waves; these waves are driven by viscous shear, so reduction of this shear decreases their amplification. The same topographies destabilize the new waves; these waves are driven by an inviscid mechanism associated with the formation of circumferential inflection points, and an increase of the groove amplitude increases their amplification. The flow conditions when the presence of grooves can be ignored, i.e. the annuli can be treated as being hydraulically smooth, have been determined.

This work demonstrates that the popular arithmetic mean conformation tensor frequently used in the analysis of turbulent viscoelastic flows is not a good representative of the ensemble. Alternative means based on recent developments in the literature are proposed, namely, the geometric and log-Euclidean means. These means are mathematically consistent with the Riemannian structure of the manifold of positive-definite tensors, on which the conformation tensor lives, and have useful properties that make them attractive alternatives to the arithmetic mean. Using a turbulent FENE-P channel flow dataset, it is shown that these two alternatives are physically representative of the ensemble. By definition, these means minimize the geodesic distance to realizations and exactly preserve the scalar geometric mean of the volume and of the principal stretches. The proposed geometric and log-Euclidean means have clear physical interpretations and are attractive quantities for turbulence modelling.

We investigate the dynamics of the buoyancy-driven rise of a bubble inside a viscoplastic material when it is subjected to an acoustic pressure field. To this end, we develop a simplified model based on the Lagrangian formalism assuming a pulsating bubble with a spherical shape. Moreover, to account for the effects of a deformable bubble, we also perform detailed two-dimensional axisymmetric simulations. Qualitative agreement is found between the simplified approach and the detailed numerical simulations. Our results reveal that the acoustic excitation enhances the mobility of the bubble, by increasing the size of the yielded region that surrounds the bubble, thereby decreasing the effective viscosity of the liquid and accelerating the motion of the bubble. This effect is significantly more pronounced at the resonance frequency, and it is shown that bubble motion takes place even for Bingham numbers (Bn) that can be orders of magnitude higher than the critical Bn for bubble entrapment in the case of a static pressure field.

We consider a potential flow model of axisymmetric waves travelling on a ferrofluid jet. The ferrofluid coats a copper wire, through which an electric current is run. The induced azimuthal magnetic field magnetises the ferrofluid, which in turn stabilises the well known Plateau–Rayleigh instability seen in axisymmetric capillary jets. This model is of interest because the stabilising mechanism allows for axisymmetric magnetohydrodynamical solitary waves. A numerical scheme capable of computing steady periodic, solitary and generalised solitary wave solutions is presented. It is found that the solution space for the model is very similar to that of the classical problem of two-dimensional gravity–capillary waves.

Classical diffusiophoresis describes the motion of particles in an electrolyte or non-electrolyte solution with an imposed concentration gradient. We investigate the autophoresis of two particles in an electrolyte solution where the concentration gradient is produced by either adsorption or desorption of ions at the particle surfaces. We find that when the sorption fluxes are large, the ion concentration near the particle surfaces, and consequently the Debye length, is strongly modified, resulting in a nonlinear dependence of the phoretic speed on the sorption flux. In particular, we show that the phoretic velocity saturates at a finite value for large desorption fluxes, but depends superlinearly on the flux for adsorption fluxes, where both conclusions are in contrast with previous results that predict a linear relationship between autophoretic velocity and sorption flux. Our theory can also be applied to precipitation/dissolution and other surface chemical processes.

We formulate the Oldroyd-B model for viscoelastic fluids in terms of magnetic-like fields obeying a set of equations analogous to Maxwell’s equations. In the limit of infinite relaxation time for the polymer, the polymeric stress tensor can be identified with the Maxwell stress tensor of a magnetic field. Away from this asymptotic case, the stress tensor of the polymer cannot be decomposed in terms of a tensor product of a magnetic field any more and several theoretical issues arise. We show that the analogy between the Oldroyd-B model and Maxwell’s equations can still be rigorously extended provided that one defines three magnetic-like fields obeying Maxwell’s equations with magnetic currents and charges. This solves the theoretical caveats and leads to a better understanding of the viscoelastic instability. In particular, we evidence a gauge symmetry which unifies some previous works, and we investigate several gauge choices. As an illustration we apply our method to viscoelastic Taylor–Couette flow but this theory of ‘viscoelastic fields’ is general and may be useful in a large variety of viscoelastic flows. The present study may also be of interest from the electromagnetic point of view, as it provides real systems possessing magnetic-like charges (monopoles) and currents.

The large Reynolds number asymptotic approximations of the neutral curve of Taylor–Couette flow subject to an axial uniform magnetic field are analysed. The flow has been extensively studied since the early 1990s as the magneto-rotational instability (MRI) occurring in the flow may explain the origin of the instability observed in some astrophysical objects. Elsewhere, the ideal approximation has been used to study high-speed flows, which sometimes produces paradoxical results. For example, ideal flows must be completely stabilised for a sufficiently strong applied magnetic field. On the other hand, the vanishing magnetic Prandtl number limit of the stability should be purely hydrodynamic, so instability must occur when Rayleigh’s stability condition is violated. Our first discovery is that this apparent contradiction can be resolved by showing the abrupt appearance of the hydrodynamic instability at a certain critical value of the magnetic Prandtl number. This is found using the asymptotically large Reynolds number limit but with a sufficiently long wavelength to retain some diffusive effects. Our second finding concerns the so-called Velikhov–Chandrasekhar paradox, namely the mismatch of the zero external magnetic field limit of the Velikhov–Chandrasekhar stability criterion and Rayleigh’s stability criterion. We show for fully wide-gap cases that the high Reynolds number asymptotic analysis of the MRI naturally yields the simple stability condition that describes smooth transition from Rayleigh to Velikhov–Chandrasekhar stability criteria with increasing Lundquist number.

New solutions describing the interaction of helical pairs of counter-rotating vortices are obtained using a vortex filament approach. The vortices are assumed to have a small core size allowing the calculation of the self-induced velocities from the Biot–Savart law using the cutoff theory. These new vortex structures do not possess any helical symmetry but they exhibit a spatial periodicity and are stationary in a rotating and translating frame. Their properties, such as radial deformation, frame velocity and induced flow, are provided as a function of the four geometric parameters characterizing each solution. Approximate solutions are also obtained when the mutual interaction is weak. This allows us to provide explicit expressions for the rotation and translation velocities of the structure in this limit. First-order corrections describing helix deformation are also calculated and used for comparison with the numerical results. The variation of the vortex core size induced by the helix deformation is also analysed. We show that these variations have a weak effect on the shape and characteristics of the solutions, for the range of parameters that we have considered. The results are finally applied to rotor wakes. It is explained how these solutions could possibly describe the far wake of an helicopter rotor in vertical flight.