An adjoint approach is developed to compute the receptivity of the rotating disc boundary layer to surface roughness. The adjoint linearised Navier–Stokes equations, in cylindrical coordinates, are derived and receptivity characteristics are computed for a broad range of azimuthal mode numbers using a fully equivalent velocity–vorticity formulation. For each set of flow conditions (i.e. azimuthal mode number), the adjoint method only requires that the linear and adjoint solutions be computed once. Thus, the adjoint approach offers significant computational and time advantages over alternative receptivity schemes (i.e. direct linearised Navier–Stokes) as they can be used to instantaneously compute the receptivity of boundary layer disturbances to many environmental mechanisms. Stationary cross-flow disturbances are established by randomly distributed surface roughness that is periodic in the azimuthal direction and modelled via a linearisation of the no-slip condition on the disc surface. Each roughness distribution is scaled on its respective root-mean-square. A Monte-Carlo type uncertainty quantification analysis is performed, whereby mean receptivity amplitudes are computed by averaging over many thousands of roughness realisations with variable length and wavelength filters. The amplitude of the cross-flow instability is significantly larger for roughness distributions near the conditions for neutral linear instability, while roughness elements radially outboard have a negligible effect on the receptivity process. Furthermore, receptivity increases sharply for roughness distributions that encompass wavelength scales equivalent to that associated with the cross-flow instability. Finally, mean receptivity characteristics are used to predict the radial range that stationary cross-flow vortices achieve amplitudes sufficient to invalidate the linear stability assumptions.

Lattice Boltzmann simulations were carried out to investigate the noise mitigation mechanisms of a 3-D printed porous trailing-edge insert, elucidating the link between noise reduction and material permeability. The porous insert is based on a unit cell resembling a lattice of diamond atoms. It replaces the last 20 % chord of a NACA 0018 at zero angle-of-attack. A partially blocked insert is considered by adding a solid partition between 84 % and 96 % of the aerofoil chord. The regular porous insert achieves a substantial noise reduction at low frequencies, although a slight noise increase is found at high frequencies. The partially blocked porous insert exhibits a lower noise reduction level, but the noise emission at mid-to-high frequency is slightly affected. The segment of the porous insert near the tip plays a dominant role in promoting noise mitigation, whereas the solid-porous junction contributes, in addition to the rough surface, towards the high-frequency excess noise. The current study demonstrates the existence of an entrance length associated with the porous material geometry, which is linked to the pressure release process that is responsible for promoting noise mitigation. This process is characterised by the aerodynamic interaction between pressure fluctuations across the porous medium, which is found at locations where the porous insert thickness is less than twice the entrance length. Present results also suggest that the noise attenuation level is related to both the chordwise extent of the porous insert and the streamwise turbulent length scale. The porous inserts also cause a slight drag increase compared to their solid counterpart.

A new perspective on the analysis of turbulent boundary layers on streamlined bodies is provided by deriving the axisymmetric Reynolds-averaged Navier–Stokes equations in an orthogonal coordinate system aligned with streamlines, streamline-normal lines and the plane of symmetry. Wall-resolved large-eddy simulation using an unstructured overset method is performed to study flow about the axisymmetric DARPA SUBOFF hull at a Reynolds number of $Re_L = 1.1 \times 10^{6}$ based on the hull length and free-stream velocity. The streamline-normal coordinate is naturally normal to the wall at the hull surface and perpendicular to the free-stream velocity far from the body, which is critical for studying bodies with concave streamwise curvature. The momentum equations naturally reduce to the differential form of Bernoulli's equation and the $s$–$n$ Euler equation for curved streamlines outside of the boundary layer. In the curved laminar boundary layer at the front of the hull, the streamline momentum equation represents a balance of the streamwise advection, streamwise pressure gradient and viscous stress, while the streamline-normal equation is a balance between the streamline-normal pressure gradient and centripetal acceleration. In the turbulent boundary layer on the mid-hull, the curvature terms and streamwise pressure gradient are negligible and the results conform to traditional analysis of flat-plate boundary layers. In the thick stern boundary layer, the curvature and streamwise pressure gradient terms reappear to balance the turbulent and viscous stresses. This balance explains the characteristic variation of static pressure observed for thick boundary layers at the tails of axisymmetric bodies.

Earlier systematic experimental studies of bursting soap films by McEntee & Mysels (J. Phys. Chem., vol. 73, 1969, pp. 3018–3028) revealed the existence of a precursor shock wave preceding the expanding hole in a punctured film, with a disturbed region of shrinking film material in between known as the ‘aureole’. In the present work we report and interpret new phenomena associated with the aureole – the formation of folds on the surface of soap films. In search of the theoretical explanation of the experimentally identified conditions under which the folds appear, we establish that they correspond to catastrophes of collapsing soap films.

We present a numerical method for simulating the flow induced by bubbles rising at large Reynolds number. This method is useful to simulate configurations of large dimensions involving a great number of bubbles. The action that each bubble exerts on the liquid is modelled as a volume source of momentum distributed over a few mesh-grid elements. The flow in the vicinity of the bubbles is thus not finely resolved. The bubbles are treated as Lagrangian particles that move under the influence of the hydrodynamic force exerted by the liquid. The determination of this force on a given bubble requires knowledge of the liquid flow that is undisturbed by this bubble. A model is developed to accurately estimate this disturbance for large-Reynolds-number objects and get rid of any spurious self-induced effect. Thanks to that, a homogeneous swarm of rising bubbles is simulated. Comparisons with experiments show a good agreement with the flow scales larger than the bubbles, which turn out to be controlled by the interactions between bubble wakes and rather independent of unresolved smaller scales. This method can be used to study the coupling between bubble-induced agitation and large-scale motions, such as those produced in industrial bubble columns.

A viscoelastic thread in vacuum is known to evolve into a beads-on-a-string structure at large deformations. If the thread is immersed in another fluid, the surrounding medium may influence the topological structure of it, which remains unexplored. To get some insights into the nonlinear behaviour of a viscoelastic thread in a two-phase flow system, a one-dimensional model is developed under the slender body approximation, in which the thread of a highly viscoelastic fluid described by the Oldroyd-B or Giesekus constitutive equation is immersed in a Newtonian viscous fluid of much smaller density and viscosity inside a cylindrical tube. The effect of the outer viscous fluid layer and the confinement of the tube is examined. It is found that the outer fluid may change substantially the beads-on-a-string structure of the viscoelastic thread. Particularly, it may induce the formation of secondary droplets on the filament between adjacent primary droplets, even for large wavenumbers. The outer fluid exerts a resistance force on the extensional flow in the filament, but the necking of the thread is slightly accelerated, due to the redistribution of capillary and elastic forces along the filament accompanied by the formation of secondary droplets. Decreasing the tube radius leads to an increase in secondary droplet size but affects little the morphology of the thread. The non-uniformity of the filament between droplets is more pronounced for a Giesekus viscoelastic thread, and pinch-off of a Giesekus thread always occurs in the neck region connecting the filament to the primary droplet in the presence of the outer viscous fluid.

The motion of a pair of bubbles rising side by side under the influence of external magnetic fields is numerically examined. Through solving the fully three-dimensional Navier–Stokes equations, the results reveal that the bubble interactions are rather sensitive to the field direction and strength. At first, we identify that, in a hydrodynamic flow, whether the two bubbles will bounce or coalesce depends on the developments of the counter-rotating streamwise vortices during the collision. In particular, for an originally bouncing bubble pair, a streamwise magnetic field tends to promote their coalescence by weakening the strengths of the standing streamwise vortices, and such a weakening effect is caused by the asymmetric distribution of the Lorentz force in the presence of another bubble such that a torque is induced to offset the original streamwise vortices. Under a horizontal magnetic field, on the other hand, the influences are highly dependent on the angle between the bubble centroid line and the field: a transverse field or a moderate spanwise field always leads the bubble pair to coalescence while a strong spanwise field has the opposite effect. This anisotropic effect comes from the Lorentz force induced flow diffusion along the magnetic field, which not only produces two pairs of streamwise vortices at the bubble rear, but also homogenizes the pressure along the magnetic lines. As the competition between the two mechanisms varies with the magnetic direction and strength, the interaction between the bubble pair also changes. We show that the external magnetic fields control the bubble interaction through reconstructing the vortex structures, and hence the core mechanisms are identified.

We investigate the dynamics of a dilute suspension of small, heavy particles superposed on a reservoir of still, pure fluid. The study is performed by means of numerical simulations of the Saffman model for a dilute particle suspension (Saffman, J. Fluid Mech., vol. 13, issue 1, 1962, pp. 120–128). In the presence of gravity forces, the interface between the two phases is unstable and evolves in a turbulent mixing layer which broadens in time. In the case of negligible particle inertia, the particle-laden phase behaves as a denser fluid, and the dynamics of the system recovers to that of the incompressible Rayleigh–Taylor set-up. Conversely, particles with large inertia affect the evolution of turbulent flow, delaying the development of turbulent mixing and breaking the up–down symmetry within the mixing layer. The inertial dynamics also leads to particle clustering, characterised by regions with higher particle density than the initial uniform density, and by the increase of the local Atwood number.

Understanding particle drift in suspension flows is of the highest importance in numerous engineering applications where particles need to be separated and filtered out from the suspending fluid. Commonly known drift mechanisms such as the Magnus force, Saffman force and Segré–Silberberg effect all arise only due to inertia of the fluid, with similar effects on all non-spherical particle shapes. In this work, we present a new shape-selective lateral drift mechanism, arising from particle inertia rather than fluid inertia, for ellipsoidal particles in a parabolic velocity profile. We show that the new drift is caused by an intermittent tumbling rotational motion in the local shear flow together with translational inertia of the particle, while rotational inertia is negligible. We find that the drift is maximal when particle inertial forces are of approximately the same order of magnitude as viscous forces, and that both extremely light and extremely heavy particles have negligible drift. Furthermore, since tumbling motion is not a stable rotational state for inertial oblate spheroids (nor for spheres), this new drift only applies to prolate spheroids or tri-axial ellipsoids. Finally, the drift is compared with the effect of gravity acting in the directions parallel and normal to the flow. The new drift mechanism is stronger than gravitational effects as long as gravity is less than a critical value. The critical gravity is highest (i.e. the new drift mechanism dominates over gravitationally induced drift mechanisms) when gravity acts parallel to the flow and the particles are small.

Recently, the authors considered a thin steady developed viscous free-surface flow passing the sharp trailing edge of a horizontally aligned flat plate under surface tension and the weak action of gravity, acting vertically, in the asymptotic slender-layer limit (J. Fluid Mech., vol. 850, 2018, pp. 924–953). We revisit the capillarity-driven short-scale viscous–inviscid interaction, on account of the inherent upstream influence, immediately downstream of the edge and scrutinise flow detachment on all smaller scales. We adhere to the assumption of a Froude number so large that choking at the plate edge is insignificant but envisage the variation of the relevant Weber number of $O(1)$. The main focus, tackled essentially analytically, is the continuation of the structure of the flow towards scales much smaller than the interactive ones and where it no longer can be treated as slender. As a remarkable phenomenon, this analysis predicts harmonic capillary ripples of Rayleigh type, prevalent on the free surface upstream of the trailing edge. They exhibit an increase of both the wavelength and amplitude as the characteristic Weber number decreases. Finally, the theory clarifies the actual detachment process, within a rational description of flow separation. At this stage, the wetting properties of the fluid and the microscopically wedge-shaped edge, viewed as infinitely thin on the larger scales, come into play. As this geometry typically models the exit of a spout, the predicted wetting of the wedge is related to what in the literature is referred to as the teapot effect.

Four different applications of spectral proper orthogonal decomposition (SPOD) are demonstrated on large-eddy simulation data of a turbulent jet. These are: low-rank reconstruction, denoising, frequency–time analysis and prewhitening. We demonstrate SPOD-based flow-field reconstruction using direct inversion of the SPOD algorithm (frequency-domain approach) and propose an alternative approach based on projection of the time series data onto the modes (time-domain approach). We further present a SPOD-based denoising strategy that is based on hard thresholding of the SPOD eigenvalues. The proposed strategy achieves significant noise reduction while facilitating drastic data compression. In contrast to standard methods of frequency–time analysis such as wavelet transform, a proposed SPOD-based approach yields a spectrogram that characterises the temporal evolution of spatially coherent flow structures. A convolution-based strategy is proposed to compute the time-continuous expansion coefficients. When applied to the turbulent jet data, SPOD-based frequency–time analysis reveals that the intermittent occurrence of large-scale coherent structures is directly associated with high-energy events. This work suggests that the time-domain approach is preferable for low-rank reconstruction of individual snapshots, and the frequency-domain approach for denoising and frequency–time analysis.

We consider the translation of a rigid, axisymmetric, tightly fitting object through a cylindrical elastic tube filled with viscous fluid, using a combination of theory and direct numerical simulations. The intruding object is assumed to be wider than the undeformed tube radius, forcing solid–solid contact in the absence of relative motion. The motion of the object establishes a thin fluid film that lubricates this contact. Our theory couples lubrication theory to a geometrically nonlinear membrane description of the tube's elasticity, and applies to a slender intruding object and a thin tube with negligible bending rigidity. We show using asymptotic and numerical solutions of the theory, that the thickness of the thin fluid film scales with the square root of the relative speed for small speeds, set by a balance of hoop stresses, membrane tension and fluid pressure. While membrane tension is relatively small at the entrance of the film, it dominates near the exit and produces undulations of the film thickness, even in the limit of vanishing speeds and slender objects. We find that the drag force on the intruding object depends on the slope of its surface at the entrance to the thin fluid film, and scales as the square root of the relative speed. The predictions of the lubricated membrane theory for the shape of the film and the force on the intruder are in quantitative agreement with three-dimensional direct numerical simulations of the coupled fluid–elastic problem.

We study turbulent flows in a smooth straight pipe of circular cross-section up to friction Reynolds number $({\textit {Re}}_{\tau }) \approx 6000$ using direct numerical simulation (DNS) of the Navier–Stokes equations. The DNS results highlight systematic deviations from Prandtl friction law, amounting to approximately $2\,\%$, which would extrapolate to approximately $4\,\%$ at extreme Reynolds numbers. Data fitting of the DNS friction coefficient yields an estimated von Kármán constant $k \approx 0.387$, which nicely fits the mean velocity profile, and which supports universality of canonical wall-bounded flows. The same constant also applies to the pipe centreline velocity, thus providing support for the claim that the asymptotic state of pipe flow at extreme Reynolds numbers should be plug flow. At the Reynolds numbers under scrutiny, no evidence for saturation of the logarithmic growth of the inner peak of the axial velocity variance is found. Although no outer peak of the velocity variance directly emerges in our DNS, we provide strong evidence that it should appear at ${\textit {Re}}_{\tau } \gtrsim 10^4$, as a result of turbulence production exceeding dissipation over a large part of the outer wall layer, thus invalidating the classical equilibrium hypothesis.

The current study uses large eddy simulations to investigate the transient response of a utility-scale wind turbine wake to dynamic changes in atmospheric and operational conditions, as observed in previous field-scale measurements. Most wind turbine wake investigations assume quasi-steady conditions, but real wind turbines operate in a highly stochastic atmosphere, and their operation (e.g. blade pitch, yaw angle) changes constantly in response. Furthermore, dynamic control strategies have been recently proposed to optimize wind farm power generation and longevity. Therefore, improved understanding of dynamic wake behaviours is essential. First, changes in blade pitch are investigated and the wake expansion response is found to display hysteresis as a result of flow inertia. The time scales of the wake response to different pitch rates are quantified. Next, changes in wind direction with different time scales are explored. Under short time scales, the wake deflection is in the opposite direction of that observed under quasi-steady conditions. Finally, yaw changes are implemented at different rates, and the maximum inverse wake deflection and time scale are quantified, showing a clear dependence on yaw rate. To gain further physical understanding of the mechanism behind the inverse wake deflection, the streamwise vorticity in different parts of the wake is quantified. The results of this study provide guidance for the design of advanced wake flow control algorithms. The lag in wake response observed for both blade pitch and yaw changes shows that proposed dynamic control strategies must implement turbine operational changes with a time scale of the order of the rotor time scale or slower.

We present direct numerical simulations of developing turbulent channel flows subjected to thermal expansion or contraction downstream of a heated or cooled wall. Using different constitutive relations for viscosity we analyse the response of variable property flows to streamwise acceleration/deceleration by separating the effect of streamwise acceleration/deceleration from the effect of wall-normal property variations. We demonstrate that, beyond a certain streamwise location, the flow can be considered in a state of ‘quasi-equilibrium’ regarding semilocally scaled variables. As such, we claim that the development of turbulent quantities due to streamwise acceleration/deceleration is localized to the region of impulsive heating/cooling, while changes in turbulence occurring farther downstream can be attributed solely to property variations. This finding allows us to study turbulence modulation in accelerating/decelerating flows using the semilocal scaling framework. By investigating the energy redistribution among the turbulent velocity fluctuations, we conclude that a change in bulk streamwise velocity has a non-local effect which originates from the change in mean shear and modifies the energy pathways through velocity-pressure-gradient correlations. On the other hand, the wall-normal property gradients have a local effect and act through the modification of the viscous dissipation. We show that it is possible to superimpose and compare the two different effects when using the semilocal scaling framework.

A new scaling is derived that yields a Reynolds-number-independent profile for all components of the Reynolds stress in the near-wall region of wall-bounded flows, including channel, pipe and boundary layer flows. The scaling demonstrates the important role played by the wall shear stress fluctuations and how the large eddies determine the Reynolds number dependence of the near-wall turbulence behaviour.

This study presents a buoyancy-driven stability analysis in a three-dimensional inclined porous medium with a capillary transition zone that is formed between a non-wetting and an underlying wetting phase. In this two-phase, two-component, partially miscible system, a solute from a non-wetting phase diffuses into a porous layer saturated with a wetting-phase fluid, creating a dense diffusive boundary layer beneath an established capillary transition zone. Transient concentration and gravity-driven velocity fields are derived for the wetting phase while the saturation field remains fixed. Linear stability analysis with the quasi-steady-state approximation is employed to determine the onset of solutal convective instability for buoyancy-dominant, in-transition and capillary-dominant systems. The analysis of the problem leads to a differential eigenvalue problem composed of a system of three complex-valued equations that are numerically solved to determine the critical times, critical wavenumbers and neutral stability curves as a function of inclination angle for different Bond numbers. The layer inclination is shown to play an essential role in the stability of the problem, where the gravity-driven flow removes solute concentrations in the diffusive boundary layer. The results indicate that the horizontal porous layer exhibits the fastest onset of instability, and longitudinal rolls are always more unstable than oblique and transverse rolls. The inclination angle has a more substantial impact on stabilizing the diffusive boundary layer in the buoyancy-dominant than in the capillary-dominant systems. Furthermore, for both buoyancy-dominant and capillary-dominant systems, the critical times and wavenumbers vary exponentially with inclination angle ≤ 60° and follow the Stirling model.

Motivated by the mushy zones of sea ice, volcanoes and icy moons of the outer solar system, we perform a theoretical and numerical study of boundary-layer convection along a vertical heated wall in a bounded ideal mushy region. The mush is comprised of a porous and reactive binary alloy with a mixture of saline liquid in a solid matrix, and is studied in the near-eutectic approximation. Here, we demonstrate the existence of four regions and study their behaviour asymptotically. Starting from the bottom of the wall, the four regions are (i) an isotropic corner region; (ii) a buoyancy dominated vertical boundary layer; (iii) an isotropic connection region; and (iv) a horizontal boundary layer at the top boundary with strong gradients of pressure and buoyancy. Scalings from numerical simulations are consistent with the theoretical predictions. Close to the heated wall, the convection in the mushy layer is similar to a rising buoyant plume abruptly stopped at the top, leading to increased pressure and temperature in the upper region, whose impact is discussed as an efficient melting mechanism.

A new nonlinear Schrödinger equation (NLSE) is presented for ocean surface waves. Earlier derivations of NLSEs that describe the evolution of deep-water waves have been limited to a narrow bandwidth, for which the bound waves at second order in wave steepness are described in leading-order approximations. This work generalizes these earlier works to allow for deep-water waves of a broad bandwidth with large directional spreading. The new NLSE permits simple numerical implementations and can be extended in a straightforward manner in order to account for waves on water of finite depth. For the description of second-order waves, this paper proposes a semianalytical approach that can provide accurate and computationally efficient predictions. With a leading-order approximation to the new NLSE, the instability region and energy growth rate of Stokes waves are investigated. Compared with the exact results based on McLean (J. Fluid Mech., vol. 511, 1982, p. 135), predictions by the new NLSE show better agreement than by Trulsen et al. (Phys. Fluids, vol. 12, 2000, pp. 2432–2437). With numerical implementations of the new NLSE, the effects of wave directionality are investigated by examining the evolution of a directionally spread focused wave group. A downward shift of the spectral peak is observed, owing to the asymmetry in the change rate of energy in a more complex manner than that for uniform Stokes waves. Rapid oblique energy transfers near the group at linear focus are observed, likely arising from the instability of uniform Stokes waves appearing in a narrow spectrum subject to oblique sideband disturbances.

Starting from the Boltzmann–Enskog kinetic equations, the charge transport equation for bidisperse granular flows with contact electrification is derived with separate mean velocities, total kinetic energies, charges and charge variances for each solid phase. To close locally averaged transport equations, a Maxwellian distribution is presumed for both particle velocity and charge. The hydrodynamic equations for bidisperse solid mixtures are first revisited and the resulting model consisting of the transport equations of mass, momentum, total kinetic energy, which is the sum of the granular temperature and the trace of fluctuating kinetic tensor, and charge is then presented. The charge transfer between phases and the charge build-up within a phase are modelled with local charge and effective work function differences between phases and the local electric field. The revisited hydrodynamic equations and the derived charge transport equation with constitutive relations are assessed through hard-sphere simulations of three-dimensional spatially homogeneous, quasi-one-dimensional spatially inhomogeneous bidisperse granular gases and a three-dimensional segregating bidisperse granular flow with conducting walls.