We investigate a hierarchy of eddy-viscosity terms in proper orthogonal decomposition (POD) Galerkin models to account for a large fraction of unresolved fluctuation energy. These Galerkin methods are applied to large eddy simulation (LES) data for a flow around a vehicle-like bluff body called an Ahmed body. This flow has three challenges for any reduced-order model: a high Reynolds number, coherent structures with broadband frequency dynamics, and meta-stable asymmetric base flow states. The Galerkin models are found to be most accurate with modal eddy viscosities as proposed by Rempfer & Fasel (J. Fluid Mech., vol. 260, 1994a, pp. 351–375; J. Fluid Mech. vol. 275, 1994b, pp. 257–283). Robustness of the model solution with respect to initial conditions, eddy-viscosity values and model order is achieved only for state-dependent eddy viscosities as proposed by Noack, Morzyński & Tadmor (Reduced-Order Modelling for Flow Control, CISM Courses and Lectures, vol. 528, 2011). Only the POD system with state-dependent modal eddy viscosities can address all challenges of the flow characteristics. All parameters are analytically derived from the Navier–Stokes-based balance equations with the available data. We arrive at simple general guidelines for robust and accurate POD models which can be expected to hold for a large class of turbulent flows.

In this work we study computationally the dynamics of a liquid bridge formed between a two-dimensional trapezoidal cavity, which represents an axisymmetric cell or a plane groove engraved in a roll, and a moving plate. The flow is a model of the liquid transfer process in gravure printing systems. The considered plate kinematics represents the actual motion of a roll-to-roll system, which includes extension, shear and rotation relative to the cavity. The fluid flow is modelled by solving the Stokes equations, discretized with the finite element method; the evolving free surfaces are accommodated by employing a pseudosolid mesh deforming algorithm. The results show that as the roll radius is reduced, thus increasing the lateral and rotational motions of the top plate relative to the cavity, a larger volume of liquid is transferred to the plate. However, due to lateral displacement of the contact lines, special care must be taken concerning the wettability properties of the substrate to avoid errors in the pattern fidelity. The predictions also show a strong nonlinear behaviour of the liquid fraction extracted from a cavity as a function of the capillary number. At high capillary numbers the fluid dynamics is mainly controlled by the extensional motion due to the strong contact line pinning. However, at low values of the capillary number, the contact lines have higher mobility and the liquid fraction primarily depends on the lateral and rotational plate velocity. These mechanisms tend to drag the fluid outside the cavity and increase the liquid fraction transferred to the plate, as has been observed in experiments.

Phoretic self-propulsion is a unique example of force- and torque-free motion on small scales. The classical framework describing the flow field around a particle swimming by self-diffusiophoresis neglects the advection of the solute field by the flow and assumes that the chemical interaction layer is thin compared to the particle size. In this paper we quantify and characterize the effect of solute advection on the phoretic swimming of a sphere. We first rigorously derive the regime of validity of the thin-interaction-layer assumption at finite values of the Péclet number (${Pe}$). Under this assumption, we solve computationally the flow around Janus phoretic particles and examine the impact of solute advection on propulsion and the flow created by the particle. We demonstrate that although advection always leads to a decrease of the swimming speed and flow stresslet at high values of the Péclet number, an increase can be obtained at intermediate values of ${Pe}$. This possible enhancement of swimming depends critically on the nature of the chemical interactions between the solute and the surface. We then derive an asymptotic analysis of the problem at small ${Pe}$ which allows us to rationalize our computational results. Our computational and theoretical analysis is accompanied by a parallel study of the influence of reactive effects at the surface of the particle (Damköhler number) on swimming.

We examine nonlinear geostrophic adjustment in a rapidly rotating (small Rossby number Ro) stably neutrally stratified (SNS) fluid consisting of a stratified upper layer with $N>f$ ($N$ is the buoyancy frequency, $f$ the Coriolis parameter) and a homogeneous lower layer, the density and other fields being continuous at the interface between the layers. The angular speed of rotation is non-parallel to gravity; the traditional and hydrostatic approximations are not used. The wave spectrum in the model consists of internal and gyroscopic waves. During the adjustment an arbitrary long-wave perturbation is split in a unique way into slow quasi-geostrophic (QG) and fast ageostrophic components with typical time scales $(Ro\, f)^{-1}$ and $f^{-1}$, respectively. The QG flow is governed by two coupled nonlinear equations of conservation of QG potential vorticity (PV) in the layers. The fast component is a sum of internal waves and inertial oscillations (long gyroscopic waves) confined to the homogeneous layer and modulated by an amplitude depending on coordinates and slow time. On times $t\sim (\, f\, Ro)^{-1}$ the slow component is not influenced by the fast one but the inertial oscillations amplitude is coupled to the QG flow and obeys an equation practically coinciding with that in the barotropic case (Reznik, J. Fluid Mech., vol. 743, 2014, pp. 585–605). A non-stationary boundary layer with large vertical gradients of horizontal velocity develops in the stratified layer near the interface to prevent penetration of the inertial oscillations into the stratified fluid; an analogous weaker boundary layer arises near the upper rigid lid. At large times the internal waves gradually decay because of dispersion and the resulting motion consists of the slow QG component and inertial oscillations confined to the barotropic lower layer.

The transition in dynamics from low-amplitude, aperiodic, combustion noise to high-amplitude, periodic, combustion instability in confined, combustion environments was studied experimentally in a laboratory-scale combustor with two different flameholding devices in a turbulent flow field. We show that the low-amplitude, irregular pressure fluctuations acquired during stable regimes, termed ‘combustion noise’, display scale invariance and have a multifractal signature that disappears at the onset of combustion instability. Traditional analysis often treats combustion noise and combustion instability as acoustic problems wherein the irregular fluctuations observed in experiments are often considered as a stochastic background to the dynamics. We demonstrate that the irregular fluctuations contain useful information of prognostic value by defining representative measures such as Hurst exponents that can act as early warning signals to impending instability in fielded combustors.

The present study investigates the mechanisms associated with tonal noise emission from a NACA 0012 aerofoil at moderate incidence ($0^{\circ },1^{\circ },2^{\circ }$ and $4^{\circ }$ angle of attack) and with Reynolds numbers ranging from 100 000 to 270 000. Simultaneous time-resolved particle image velocimetry (PIV) of the aeroacoustic source region near the trailing edge and acoustic measurements in the far field are performed in order to establish the correspondence between the flow structure and acoustic emissions. Results of these experiments are presented and analysed in view of past research for a number of selected cases. Characteristics of the acoustic emission and principal features of the average flow field agree with data presented in previous studies on the topic. Time-resolved analysis shows that downstream convecting vortical structures, resulting from growing shear layer instabilities, coherently pass the trailing edge at a frequency equal to that of the dominant tone. Therefore, the scattering of the vortical structures and their associated wall pressure fluctuations are identified as tone generating mechanisms for the cases investigated here. Moreover, wavelet analysis of the acoustic pressure and velocity signals near the trailing edge show a similar periodic amplitude modulation which is associated with multiple tonal peaks in the acoustic spectrum. Periodic amplitude modulation of the acoustic pressure and velocity fluctuations on the pressure side are also observed when transition is forced on the suction side, showing that pressure-side events alone can be the cause.

This study investigates the free transverse flow-induced vibration (FIV) of an elastically mounted low-mass-ratio square cylinder in a free stream, at three different incidence angles: ${{\alpha }}=0^\circ $, $20^\circ $ and $45^\circ $. This geometric setup presents a body with an angle of attack, sharp corners and some afterbody, and therefore is a generic body that can be used to investigate a wide range of FIV phenomena. A recent study by Nemes et al. (J. Fluid Mech., vol. 710, 2012, pp. 102–130) provided a broad overview of the flow regimes present as a function of both the angle of attack ${{\alpha }}$ and reduced flow velocity ${U^{*}}$. Here, the focus is on the three aforementioned representative angles of attack: ${{\alpha }}=0^\circ $, where the FIV is dominated by transverse galloping; ${{\alpha }}=45^\circ $, where the FIV is dominated by vortex-induced vibration (VIV); and an intermediate value of ${{\alpha }}=20^\circ $, where the underlying FIV phenomenon has previously been difficult to determine. For the ${{\alpha }}=0^\circ $ case, the amplitude of oscillation increases linearly with the flow speed except for a series of regimes that occur when the vortex shedding frequency is in the vicinity of an odd-integer multiple of the galloping oscillation frequency, and the vortex shedding synchronizes to this multiple of the oscillation frequency. It is shown that only odd-integer multiple synchronizations should occur. These synchronizations explain the ‘kinks’ in the galloping amplitude response for light bodies first observed by Bearman et al. (J. Fluids Struct., vol. 1, 1987, pp. 19–34). For the ${{\alpha }}=45^\circ $ case, the VIV response consists of a number of subtle, but distinctly different regimes, with five regimes of high-amplitude oscillations, compared to two found in the classic VIV studies of a circular cylinder. For the intermediate ${{\alpha }}=20^\circ $ case, a typical VIV ‘upper branch’ occurs followed by a ‘higher branch’ of very large-amplitude response. The higher branch is caused by a subharmonic synchronization between the vortex shedding and the body oscillation frequency, where two cycles of vortex shedding occur over one cycle of oscillation. It appears that this subharmonic synchronization is a direct result of the asymmetric body. Overall, the FIV of the square cylinder is shown to be very rich, with a number of distinct regimes, controlled by both ${{\alpha }}$ and ${U^{*}}$. Importantly, ${{\alpha }}$ controls the underlying FIV phenomenon, as well as controlling the types of possible synchronization between the oscillation and vortex shedding.

Despite the confirmation of slip flows and successful drag reduction (DR) in small-scaled laminar flows, the full impact of superhydrophobic (SHPo) DR remained questionable because of the sporadic and inconsistent experimental results in turbulent flows. Here we report a systematic set of bias-free reduction data obtained by measuring the skin-friction drags on a SHPo surface and a smooth surface at the same time and location in a turbulent boundary layer (TBL) flow. Each monolithic sample consists of a SHPo surface and a smooth surface suspended by flexure springs, all carved out from a $2.7 \times 2.7 {\mathrm{mm}}^{2}$ silicon chip by photolithographic microfabrication. The flow tests allow continuous monitoring of the plastron on the SHPo surfaces, so that the DR data are genuine and consistent. A family of SHPo samples with precise profiles reveals the effects of grating parameters on turbulent DR, which was measured to be as much as ${\sim }75\, \%$.

The equilibrium shape of liquid drops on elastic substrates is determined by minimizing elastic and capillary free energies, focusing on thick incompressible substrates. The problem is governed by three length scales: the size of the drop $R$, the molecular size $a$ and the ratio of surface tension to elastic modulus $\gamma /E$. We show that the contact angles undergo two transitions upon changing the substrate from rigid to soft. The microscopic wetting angles deviate from Young’s law when $\gamma /(Ea)\gg 1$, while the apparent macroscopic angle only changes in the very soft limit $\gamma /(ER)\gg 1$. The elastic deformations are worked out for the simplifying case where the solid surface energy is assumed to be constant. The total free energy turns out to be lower on softer substrates, consistent with recent experiments.

We present a computational study of cluster-induced turbulence (CIT), where the production of fluid-phase kinetic energy results entirely from momentum coupling with finite-size inertial particles. A separation of length scales must be established when evaluating the particle dynamics in order to distinguish between the continuous mesoscopic velocity field and the uncorrelated particle motion. To accomplish this, an adaptive spatial filter is employed on the Lagrangian data with an averaging volume that varies with the local particle-phase volume fraction. This filtering approach ensures sufficient particle sample sizes in order to obtain meaningful statistics while remaining small enough to avoid capturing variations in the mesoscopic particle field. Two-point spatial correlations are computed to assess the validity of the filter in extracting meaningful statistics. The method is used to investigate, for the first time, the properties of a statistically stationary gravity-driven particle-laden flow, where particle–particle and fluid–particle interactions control the multiphase dynamics. Results from fully developed CIT show a strong correlation between the local volume fraction and the granular temperature, with maximum values located at the upstream boundary of clusters (i.e. where maximum compressibility of the particle velocity field exists), while negligible particle agitation is observed within clusters.