We investigate the motion of a sedimenting drop in the presence of an electric field in an arbitrary direction, otherwise uniform, in the limit of small interface deformation and low-surface-charge convection. We analytically solve the electric potential in and around the leaky dielectric drop, and solve for the Stokesian velocity and pressure fields. We obtain the correction in drop velocity due to shape deformation and surface-charge convection considering small capillary number and small electric Reynolds number which signifies the importance of charge convection at the drop surface. We show that tilt angle, which quantifies the angle of inclination of the applied electric field with respect to the direction of gravity, has a significant effect on the magnitude and direction of the drop velocity. When the electric field is tilted with respect to the direction of gravity, we obtain a non-intuitive lateral motion of the drop in addition to the buoyancy-driven sedimentation. Both the charge convection and shape deformation yield this lateral migration of the drop. Our analysis indicates that depending on the magnitude of the tilt angle, conductivity and permittivity ratios, the direction of the sedimenting drop can be controlled effectively. Our experimental investigation further confirms the presence of lateral migration of the drop in the presence of a tilted electric field, which is in support of the essential findings from the analytical formalism.

The smoothed profile method is extended to study the rheological behaviour of colloidal dispersions under shear flow by using the Lees–Edwards boundary conditions. We start with a reformulation of the smoothed profile method, a direct numerical simulation method for colloidal dispersions, so that it can be used with the Lees–Edwards boundary condition, under steady or oscillatory-shear flow. By this reformulation, all the resultant physical quantities, including local and total shear stresses, become available through direct calculation. Three simple rheological simulations are then performed for (1) a spherical particle, (2) a rigid bead chain and (3) a collision of two spherical particles under shear flow. Quantitative validity of these simulations is examined by comparing the viscosity with that obtained from theory and Stokesian dynamics calculations. Finally, we consider the shear-thinning behaviour of concentrated colloidal dispersions.

The large-scale coherent motions in a realistic swirl fuel-injector geometry are analysed by direct numerical simulations (DNS), proper orthogonal decomposition (POD), and linear global modes. The aim is to identify the origin of instability in this turbulent flow in a complex internal geometry. The flow field in the nonlinear simulation is highly turbulent, but with a distinguishable coherent structure: the precessing vortex core (a spiralling mode). The most energetic POD mode pair is identified as the precessing vortex core. By analysing the fast Fourier transform (FFT) of the time coefficients of the POD modes, we conclude that the first four POD modes contain the coherent fluctuations. The remaining POD modes (incoherent fluctuations) are used to form a turbulent viscosity field, using the Newtonian eddy model. The turbulence sets in from convective shear layer instabilities even before the nonlinear flow reaches the other end of the domain, indicating that equilibrium solutions of the Navier–Stokes are never observed. Linear global modes are computed around the mean flow from DNS, applying the turbulent viscosity extracted from POD modes. A slightly stable discrete $m=1$ eigenmode is found, well separated from the continuous spectrum, in very good agreement with the POD mode shape and frequency. The structural sensitivity of the precessing vortex core is located upstream of the central recirculation zone, identifying it as a spiral vortex breakdown instability in the nozzle. Furthermore, the structural sensitivity indicates that the dominant instability mechanism is the Kelvin–Helmholtz instability at the inflection point forming near vortex breakdown. Adjoint modes are strong in the shear layer along the whole extent of the nozzle, showing that the optimal initial condition for the global mode is localized in the shear layer. We analyse the qualitative influence of turbulent dissipation in the stability problem (eddy viscosity) on the eigenmodes by comparing them to eigenmodes computed without eddy viscosity. The results show that the eddy viscosity improves the complex frequency and shape of global modes around the fuel-injector mean flow, while a qualitative wavemaker position can be obtained with or without turbulent dissipation, in agreement with previous studies. This study shows how sensitivity analysis can identify which parts of the flow in a complex geometry need to be altered in order to change its hydrodynamic stability characteristics.

We quantify experimentally the dispersal characteristics of dense particle clouds in high-speed interactions with an atmosphere. Focused on the fundamentals, the experiments, conducted in a large-scale shock tube, involve a well-characterized ‘curtain’ of (falling) particles that fully occupies the cross-sectional area of the expansion section. The particle material (glass) and size (${\sim}$1 mm) are fixed, as is the curtain thickness (${\sim}$30 mm) and the particle volume fractions in it, varying from ${\sim}$58 % at the top of the curtain to ${\sim}$24 % near the bottom. Thus, the principal experimental variable is the impacting shock strength, with Mach numbers varying from 1.2 to 2.6, and flow speeds that cover from subsonic ($M_{IS}\sim 0.3$) to transonic and supersonic ($M_{IS}\sim 1.2$). The peak shock pressure ratio, 7.6, yields a flow speed of ${\sim}\!630~\text{m}~\text{s}^{-1}$, and a curtain expansion rate at ${\sim}$20 000 g. We record visually (high-speed, particle-resolving shadowgraphic method) the reflected/transmitted pressure waves and the transmitted contact wave, as well as the curtain displacements, and we measure the reflected/transmitted pressure transients. Data analysis yields simple rules for the amplitudes of the reflected pressure waves and the rapid cloud expansions observed, and we discover a time scaling that hints at a universal regime for cloud expansion. The data and these data-analysis results can provide the validation basis for numerical simulations meant to enable a deeper understanding of the key physics that drive this rather complex dispersal process.

The leading-edge boundary layer (LEBL) in the front part of swept airplane wings is prone to three-dimensional subcritical instability, which may lead to bypass transition. The resulting increase of airplane drag and fuel consumption implies a negative environmental impact. In the present paper, we present a temporal biglobal secondary stability analysis (SSA) and direct numerical simulations (DNS) of this flow to investigate a subcritical transition mechanism. The LEBL is modelled by the swept Hiemenz boundary layer (SHBL), with and without wall suction. We introduce a pair of steady, counter-rotating, streamwise vortices next to the attachment line as a generic primary disturbance. This generates a high-speed streak, which evolves slowly in the streamwise direction. The SSA predicts that this flow is unstable to secondary, time-dependent perturbations. We report the upper branch of the secondary neutral curve and describe numerous eigenmodes located inside the shear layers surrounding the primary high-speed streak and the vortices. We find secondary flow instability at Reynolds numbers as low as $Re\approx 175$, i.e. far below the linear critical Reynolds number $Re_{crit}\approx 583$ of the SHBL. This secondary modal instability is confirmed by our three-dimensional DNS. Furthermore, these simulations show that the modes may grow until nonlinear processes lead to breakdown to turbulent flow for Reynolds numbers above $Re_{tr}\approx 250$. The three-dimensional mode shapes, growth rates, and the frequency dependence of the secondary eigenmodes found by SSA and the DNS results are in close agreement with each other. The transition Reynolds number $Re_{tr}\approx 250$ at zero suction and its increase with wall suction closely coincide with experimental and numerical results from the literature. We conclude that the secondary instability and the transition scenario presented in this paper may serve as a possible explanation for the well-known subcritical transition observed in the leading-edge boundary layer.

We calculate here the force on a probe at a viscous, compressible interface, laden with soluble surfactant that equilibrates on a finite time scale. The motion of the probe through the interface drives variations in the surfactant concentration at the interface that in turn leads to a Marangoni flow that contributes to the force on the probe. We demonstrate that the Marangoni force on the probe depends non-trivially on the surface shear and dilatational viscosities of the interface indicating the difficulty in extracting these material properties from force measurements at compressible interfaces.

In this study we give a characterization of semi-geostrophic turbulence by performing freely decaying simulations for the case of constant uniform potential vorticity, a set of equations known as the surface semi-geostrophic approximation. The equations are formulated as conservation laws for potential temperature and potential vorticity, with a nonlinear Monge–Ampère type inversion equation for the streamfunction, expressed in a transformed coordinate system that follows the geostrophic flow. We perform model studies of turbulent surface semi-geostrophic flows in a domain doubly periodic in the horizontal and limited in the vertical by two rigid lids, allowing for variations of potential temperature at one of the boundaries, and we compare the results with those obtained in the corresponding surface quasi-geostrophic case. The results show that, while the surface quasi-geostrophic dynamics is dominated by a symmetric population of cyclones and anticyclones, the surface semi-geostrophic dynamics features a more prominent role of fronts and filaments. The resulting distribution of potential temperature is strongly skewed and peaked at non-zero values at and close to the active boundary, while symmetry is restored in the interior of the domain, where small-scale frontal structures do not penetrate. In surface semi-geostrophic turbulence, energy spectra are less steep than in the surface quasi-geostrophic case, with more energy concentrated at small scales for increasing Rossby number. The energy related to frontal structures, the lateral strain rate and the vertical velocities are largest close to the active boundary. These results show that the semi-geostrophic model could be of interest for studying the lateral mixing of properties in geophysical flows.

In many physiological settings, microorganisms must swim through viscous fluids with suspended polymeric networks whose length scales are comparable to that of the organism. Here we present a model of a flagellar swimmer moving through a compliant viscoelastic network immersed in a three-dimensional viscous fluid. The swimmer moves with a prescribed gait, exerting forces on the fluid and the heterogeneous network. The viscoelastic structural links of this network are stretched or compressed in response to the fluid flow caused by these forces, and these elastic deformations also generate forces on the viscous fluid. Here we track the swimmer as it leaves a region of Newtonian fluid, enters and moves through a heterogeneous network and finally enters a Newtonian region again. We find that stiffer networks give a boost to the velocity of the swimmer. In addition, we find that the efficiency of swimming is dependent upon the evolution of the compliant network as the swimmer progresses through it.

The identification of coherent structures from experimental or numerical data is an essential task when conducting research in fluid dynamics. This typically involves the construction of an empirical mode base that appropriately captures the dominant flow structures. The most prominent candidates are the energy-ranked proper orthogonal decomposition (POD) and the frequency-ranked Fourier decomposition and dynamic mode decomposition (DMD). However, these methods are not suitable when the relevant coherent structures occur at low energies or at multiple frequencies, which is often the case. To overcome the deficit of these ‘rigid’ approaches, we propose a new method termed spectral proper orthogonal decomposition (SPOD). It is based on classical POD and it can be applied to spatially and temporally resolved data. The new method involves an additional temporal constraint that enables a clear separation of phenomena that occur at multiple frequencies and energies. SPOD allows for a continuous shifting from the energetically optimal POD to the spectrally pure Fourier decomposition by changing a single parameter. In this article, SPOD is motivated from phenomenological considerations of the POD autocorrelation matrix and justified from dynamical systems theory. The new method is further applied to three sets of PIV measurements of flows from very different engineering problems. We consider the flow of a swirl-stabilized combustor, the wake of an airfoil with a Gurney flap and the flow field of the sweeping jet behind a fluidic oscillator. For these examples, the commonly used methods fail to assign the relevant coherent structures to single modes. The SPOD, however, achieves a proper separation of spatially and temporally coherent structures, which are either hidden in stochastic turbulent fluctuations or spread over a wide frequency range. The SPOD requires only one additional parameter, which can be estimated from the basic time scales of the flow. In spite of all these benefits, the algorithmic complexity and computational cost of the SPOD are only marginally greater than those of the snapshot POD.

The effect of an external pressure disturbance, being displaced with a constant speed along the free surface of a viscous thin film, is studied theoretically in the lubrication approximation in one- and two-dimensional geometries. In the comoving frame, the imposed pressure field creates a stationary deformation of the interface – a wake – that spatially vanishes in the far region. The shape of the wake and the way it vanishes depend on both the speed and size of the external source and the properties of the film. The wave resistance, namely the force that has to be externally furnished in order to maintain the wake, is analysed in detail. For finite-size pressure disturbances, it increases with the speed, up to a certain transition value, above which a monotonic decrease occurs. The role of the horizontal extent of the pressure field is studied as well, revealing that for a smaller disturbance the latter transition occurs at a higher speed. Eventually, for a Dirac pressure source, the wave resistance either saturates for a one-dimensional geometry, or diverges for a two-dimensional geometry.

A train of high-speed microdrops impacting on a liquid pool can create a very deep and narrow cavity, reaching depths more than 1000 times the size of the individual drops. The impact of such a droplet train is studied numerically using boundary integral simulations. In these simulations, we solve the potential flow in the pool and in the impacting drops, taking into account the influence of liquid inertia, gravity and surface tension. We show that for microdrops the cavity shape and maximum depth primarily depend on the balance of inertia and surface tension and discuss how these are influenced by the spacing between the drops in the train. Finally, we derive simple scaling laws for the cavity depth and width.

We study the case of inertial particles heated by thermal radiation while settling by gravity through a turbulent transparent gas. We consider dilute and optically thin regimes in which each particle receives the same heat flux. Numerical simulations of forced homogeneous turbulence are performed taking into account the two-way coupling of both momentum and temperature between the dispersed and continuous phases. Particles much smaller than the smallest flow scales are considered and the point-particle approximation is adopted. The particle Stokes number (based on the Kolmogorov time scale) is of order unity, while the nominal settling velocity is up to an order of magnitude larger than the Kolmogorov velocity, marking a critical difference with previous two-way coupled simulations. It is found that non-heated particles enhance turbulence when their settling velocity is sufficiently high compared to the Kolmogorov velocity. Energy spectra show that the non-heated particle settling impacts both the very small and very large flow scales, while the intermediate scales are weakly affected. When heated, particles shed plumes of buoyant gas, further modifying the turbulence structure. At the considered radiation intensities, clustering is strong but the classic mechanism of preferential concentration is modified, while preferential sweeping is eliminated or even reversed. Particle heating also causes a significant reduction of the mean settling velocity, which is caused by rising buoyant plumes in the vicinity of particle clusters. The turbulent kinetic energy is affected non-monotonically as the radiation intensity is increased due to the competing effects of the downward gravitational force and the upward buoyancy force. The thermal radiation influences all scales of the turbulence. The effects of settling and buoyancy on the turbulence anisotropy are also discussed.

The modal instability encountered by the incompressible flow inside a toroidal pipe is studied, for the first time, by means of linear stability analysis and direct numerical simulation (DNS). In addition to the unquestionable aesthetic appeal, the torus represents the smallest departure from the canonical straight pipe flow, at least for low curvatures. The flow is governed by only two parameters: the Reynolds number $\mathit{Re}$ and the curvature of the torus ${\it\delta}$, i.e. the ratio between pipe radius and torus radius. The absence of additional features, such as torsion in the case of a helical pipe, allows us to isolate the effect that the curvature has on the onset of the instability. Results show that the flow is linearly unstable for all curvatures investigated between 0.002 and unity, and undergoes a Hopf bifurcation at $\mathit{Re}$ of about 4000. The bifurcation is followed by the onset of a periodic regime, characterised by travelling waves with wavelength $\mathit{O}(1)$ pipe diameters. The neutral curve associated with the instability is traced in parameter space by means of a novel continuation algorithm. Tracking the bifurcation provides a complete description of the modal onset of instability as a function of the two governing parameters, and allows a precise calculation of the critical values of $\mathit{Re}$ and ${\it\delta}$. Several different modes are found, with differing properties and eigenfunction shapes. Some eigenmodes are observed to belong to groups with a set of common characteristics, deemed ‘families’, while others appear as ‘isolated’. Comparison with nonlinear DNS shows excellent agreement, confirming every aspect of the linear analysis, its accuracy, and proving its significance for the nonlinear flow. Experimental data from the literature are also shown to be in considerable agreement with the present results.

The paper presents a newly developed method for the analysis of the flow around open rotors characterised by hubs of general shape. The exact and implicit solution of the axysimmetric, inviscid and incompressible flow is represented as the superposition of infinite ring vortices properly arranged along the hub surface and the rotor wake. The solution is made explicit through a semi-analytical and iterative procedure. The proposed semi-analytical approach can deal with hubs of arbitrary shape and with quite general rotor load distributions. The method strongly couples the flow induced by the rotor and the hub. Moreover, the contraction/divergence and the rotation of the wake can be fully taken into account. The results of the semi-analytical method are also compared with those obtained with a widely diffused actuator disk model based on computational fluid dynamics (CFD) techniques. Finally, in comparison with more advanced methods, such as those relying on a CFD approach, this method is characterised by an extremely reduced computational cost. The computer code is freely available on contacting the authors.

Direct numerical simulations are used to study the interaction of a cavity-backed circular orifice with grazing laminar and turbulent boundary layers and incident sound waves. The flow conditions and geometry are representative of single degree-of-freedom acoustic liners applied in the inlet and exhaust ducts of aircraft engines and are the same as those from experiments conducted at NASA Langley. The simulations identify the fluid mechanics of how the sound field and state of the grazing boundary layer impact the in-orifice flow and suggest a simple flow analogy that enables scaling estimates. From the scaling estimates the simulations are then used to develop reduced-order models for the in-orifice flow and a time-domain impedance model is constructed. The liner is found to increase drag at all conditions studied by an amount that increases with the incident sound pressure amplitude.

An analytical correlation is developed for stagnation-point heat flux on spherical objects travelling at high velocity which is accurate for conditions ranging from the continuum to the free-molecular flow regime. Theoretical analysis of the Burnett and super-Burnett equations is performed using simplifications from shock-wave and boundary-layer theory to determine the relative contribution of higher-order heat flux terms compared with the Fourier heat flux (assumed in the Navier–Stokes equations). A rarefaction parameter ($W_{r}\equiv M_{\infty }^{2{\it\omega}}/Re_{\infty }$), based on the free-stream Mach number ($M_{\infty }$), the Reynolds number ($Re_{\infty }$) and the viscosity–temperature index (${\it\omega}$), is identified as a better correlating parameter than the Knudsen number in the transition regime. By studying both the Burnett and super-Burnett equations, a general form for the entire series of higher-order heat flux contributions is obtained. The resulting heat flux expression includes terms with dependence on gas properties, stagnation to wall-temperature ratio and a main dependence on powers of the rarefaction parameter $W_{r}$. The expression is applied as a correction to the Fourier heat flux and therefore can be combined with any continuum-based correlation of choice. In the free-molecular limit, a bridging function is used to ensure consistency with well-established free-molecular flow theory. The correlation is then fitted to direct simulation Monte Carlo (DSMC) solutions for stagnation-point heat flux in high-speed nitrogen flows. The correlation is shown to accurately capture the variation in heat flux predicted by the DSMC method in the transition flow regime, while limiting to both continuum and free-molecular values.

Blocked, continuously stratified, crest-controlled flows have hydraulically supercritical downslope flow in the lee of a ridge-like obstacle. The downslope flow separates from the obstacle and, depending on conditions further downstream, transitions to a subcritical state. A controlled, stratified overflow and its transition to a subcritical state are investigated here in a set of three-dimensional numerical experiments in which the height of a second, downstream ridge is varied. The downslope flow is associated with an isopycnal and streamline bifurcation, which acts to form a nearly-uniform-density isolating layer and a sharp pycnocline that separates deeper blocked and stratified fluid between the ridges from the flow above. The height of the downstream obstacle is communicated upstream via gravity waves that propagate along the density interface and set the separation depth of the downslope flow. The penetration depth of the downslope flow, its susceptibility to shear instabilities, and the amount of energy dissipated in the turbulent outflow all increase as the height of a downstream ridge, which effectively sets the downstream boundary conditions, is reduced.

We develop a generalised unsteady plume theory and compare it with a new direct numerical simulation (DNS) dataset for an ensemble of statistically unsteady turbulent plumes. The theoretical framework described in this paper generalises previous models and exposes several fundamental aspects of the physics of unsteady plumes. The framework allows one to understand how the structure of the governing integral equations depends on the assumptions one makes about the radial dependence of the longitudinal velocity, turbulence and pressure. Consequently, the ill-posed models identified by Scase & Hewitt (J. Fluid Mech., vol. 697, 2012, pp. 455–480) are shown to be the result of a non-physical assumption regarding the velocity profile. The framework reveals that these ill-posed unsteady plume models are degenerate cases amongst a comparatively large set of well-posed models that can be derived from the generalised unsteady plume equations that we obtain. Drawing on the results of DNS of a plume subjected to an instantaneous step change in its source buoyancy flux, we use the framework in a diagnostic capacity to investigate the properties of the resulting travelling wave. In general, the governing integral equations are hyperbolic, becoming parabolic in the limiting case of a ‘top-hat’ model, and the travelling wave can be classified as lazy, pure or forced according to the particular assumptions that are invoked to close the integral equations. Guided by observations from the DNS data, we use the framework in a prognostic capacity to develop a relatively simple, accurate and well-posed model of unsteady plumes that is based on the assumption of a Gaussian velocity profile. An analytical solution is presented for a pure straight-sided plume that is consistent with the key features observed from the DNS.

Townsend’s attached eddy hypothesis assumes the existence of a set of energetic and geometrically self-similar eddies in the logarithmic layer in wall-bounded turbulent flows, which can be scaled with their distance to the wall. To examine the possible self-similarity of the energetic eddies in fully developed turbulent pipe flow, we performed stereo particle image velocimetry measurements together with a proper orthogonal decomposition analysis. For two Reynolds numbers, $Re_{{\it\tau}}=1330$ and 2460, the resulting modes/eddies were shown to exhibit self-similar behaviour for eddies with wall-normal length scales spanning a decade. This single length scale provides a complete description of the cross-sectional shape of the self-similar eddies.