The planar laminar flow resulting from the impingement of two gaseous jets of different density issuing into an open space from aligned steadily fed slot nozzles of semi-width $H$ separated by a distance $2L$ is investigated by numerical and analytical methods. Specific consideration is given to the high Reynolds and low Mach number conditions typically present in counterflow-flame experiments, for which the flow is nearly inviscid and incompressible. It is shown that introduction of a density-weighted vorticity–streamfunction formulation effectively reduces the problem to one involving two jets of equal density, thereby removing the vortex-sheet character of the interface separating the two jet streams. Besides the geometric parameter $L/H$, the solution depends only on the shape of the velocity profiles in the feed streams and on the jet momentum-flux ratio. While conformal mapping can be used to determine the potential solution corresponding to uniform velocity profiles, numerical integration is required in general to compute rotational flows, including those arising with Poiseuille velocity profiles, with simplified solutions found in the limits $L/H\ll 1$ and $L/H\gg 1$. The results are used to quantify the near-stagnation-point region, of interest in counterflow-flame studies, including the local value of the strain rate as well as the curvature of the separating interface and the variations of the strain rate away from the stagnation point.

The linear instability mechanisms of horizontal convection in a rectangular cavity forced by a horizontal buoyancy gradient along its surface are investigated using local and global stability analyses for a Prandtl number equal to unity. The results show that the stability of the base flow, a steady circulation characterized by a narrow descending plume and a broad upwelling region, depends on the Rayleigh number, $Ra$. For free-slip boundary conditions at a critical value of $Ra\approx 2\times 10^{7}$, the steady base flow becomes unstable to three-dimensional perturbations, characterized by counter-rotating vortices originating within the plume region. A Wentzel–Kramers–Brillouin (WKB) method applied along closed streamlines demonstrates that this instability is of a Rayleigh–Taylor type and can be used to accurately reconstruct the global instability mode. In the case of no-slip boundary conditions, the base flow also becomes unstable to a self-sustained two-dimensional instability whose critical Rayleigh number is $Ra\approx 1.7\times 10^{8}$. Beyond this critical $Ra$, two-dimensional equilibrium stationary states of the Navier–Stokes equations are computed using the selective frequency damping method. The two-dimensional onset of instability is shown to be characterized by a family of modes also originating within the plume. A local spatio-temporal stability analysis shows that the flow becomes absolutely unstable at the origin of the plume. Taken together, these results illustrate the mechanisms behind the onset of turbulence that has been observed in horizontal convection.

New analytical models are introduced to describe the motion of a Herschel–Bulkley fluid slumping under gravity in a narrow fracture and in a porous medium. A useful self-similar solution can be derived for a fluid injection rate that scales as time $t$; an expansion technique is adopted for a generic injection rate that is power law in time. Experiments in a Hele-Shaw cell and in a narrow channel filled with glass ballotini confirm the theoretical model within the experimental uncertainty.

This paper presents a suite of analytical solutions, for both the free-surface elevation and the flow velocity, for landslide-generated water waves. The one-dimensional (horizontal, 1DH) analytical solutions for water waves generated by a solid landslide moving at a constant speed in constant water depth were obtained for the linear and weakly dispersive wave model as well as the linear and fully dispersive wave model. The area enclosed by the landslide was shown to have stronger lasting effects on the generated water waves than the exact landslide shape. In addition, the resonance solution based on the fully dispersive wave model was examined, and the growth rate was derived. For the 1DH linear shallow water equations (LSWEs) on a constant slope, a closed-form analytical solution, which could serve as a useful benchmark for numerical models, was found for a special landslide forcing function. For the two-dimensional (horizontal, 2DH) LSWEs on a plane beach, we rederived the solutions using the quiescent water initial conditions. The difference between the initial conditions used in the new solutions and those used in previous studies was found to have a permanent effect on the generated waves. We further noted that convergence rate of the 2DH LSWE analytical solutions varies greatly, and advised that case-by-case convergence tests be conducted whenever the modal analytical solutions are numerically evaluated using a finite number of modes.

We study the motion of a buoyant or a nearly neutrally buoyant nano-sized spheroid in a fluid filled tube without or with an imposed pressure gradient (weak Poiseuille flow). The fluctuating hydrodynamics approach and the deterministic method are both employed. We ensure that the fluctuation–dissipation relation and the principle of thermal equipartition of energy are both satisfied. The major focus is on the effect of the confining boundary. Results for the velocity and the angular velocity autocorrelations (VACF and AVACF), the diffusivities and the drag and the lift forces as functions of the shape, the aspect ratio, the inclination angle and the proximity to the wall are presented. For the parameters considered, the boundary modifies the VACF and AVACF such that three distinct regimes are discernible – an initial exponential decay followed by an algebraic decay culminating in a second exponential decay. The first is due to the thermal noise, the algebraic regime is due both to the thermal noise and the hydrodynamic correlations, while the second exponential decay shows the effect of momentum reflection from the confining wall. Our predictions display excellent comparison with published results for the algebraic regime (the only regime for which earlier results exist). We also discuss the role of the off-diagonal elements of the mobility and the diffusivity tensors that enable the quantifications of the degree of lift and margination of the nanocarrier. Our study covers a range of parameters that are of wide applicability in nanotechnology, microrheology and in targeted drug delivery.

The influence of dense-gas effects on compressible wall-bounded turbulence is investigated by means of direct numerical simulations of supersonic turbulent channel flows. Results are obtained for PP11, a heavy fluorocarbon representative of dense gases, the thermophysics properties of which are described by using a fifth-order virial equation of state and advanced models for the transport properties. In the dense-gas regime, the speed of sound varies non-monotonically in small perturbations and the dependency of the transport properties on the fluid density (in addition to the temperature) is no longer negligible. A parametric study is carried out by varying the bulk Mach and Reynolds numbers, and results are compared to those obtained for a perfect gas, namely air. Dense-gas flow exhibits almost negligible friction heating effects, since the high specific heat of the fluids leads to a loose coupling between thermal and kinetic fields, even at high Mach numbers. Despite negligible temperature variations across the channel, the mean viscosity tends to decrease from the channel walls to the centreline (liquid-like behaviour), due to its complex dependency on fluid density. On the other hand, strong density fluctuations are present, but due to the non-standard sound speed variation (opposite to the mean density evolution across the channel), the amplitude is maximal close to the channel wall, i.e. in the viscous sublayer instead of the buffer layer like in perfect gases. As a consequence, these fluctuations do not alter the turbulence structure significantly, and Morkovin’s hypothesis is well respected at any Mach number considered in the study. The preceding features make high Mach wall-bounded flows of dense gases similar to incompressible flows with variable properties, despite the significant fluctuations of density and speed of sound. Indeed, the semi-local scaling of Patel et al. (Phys. Fluids, vol. 27 (9), 2015, 095101) or Trettel & Larsson (Phys. Fluids, vol. 28 (2), 2016, 026102) is shown to be well adapted to compare results from existing surveys and with the well-documented incompressible limit. Additionally, for a dense gas the isothermal channel flow is also almost adiabatic, and the Van Driest transformation also performs reasonably well. The present observations open the way to the development of suitable models for dense-gas turbulent flows.

Flapping and revolving wings can produce attached leading-edge vortices when the angle of attack is large. In this work, a low-order model is proposed for the edge vortices that develop on a revolving plate at $90^{\circ }$ angle of attack, which is the simplest limiting case, yet shows remarkable similarity with the generally known leading-edge vortices. The problem is solved analytically, providing short closed-form expressions for the circulation and the position of the vortex. The good agreement with the numerical solution of the Navier–Stokes equations suggests that, for the conditions examined, the vorticity production at the sharp edge and its subsequent three-dimensional transport are the main effects that shape the edge vortex.

The unsteady global dynamics of a gravitational liquid sheet interacting with a one-sided adjacent air enclosure (commonly referred to as nappe oscillation configuration) is addressed under the assumptions of potential flow and the presence of surface tension effects. From a theoretical viewpoint the problem is challenging, because from previous literature it is known that the equation governing the evolution of small disturbances exhibits a singularity at the vertical station where the local flow velocity equals the capillary wave velocity (local critical condition), although the solution to the problem has not yet been found. The equation governing the local dynamics resembles one featuring the forced vibrations of a string of finite length, formulated in the reference frame moving with the flow velocity, and exhibits both slow and fast characteristic curves. From the global system perspective the nappe behaves as a driven damped spring–mass oscillator, where the inertial effects are linked to the liquid sheet mass and the spring is represented by the equivalent stiffness of the air enclosure acting on the displacement of the compliant nappe centreline. A suited procedure is developed to remove the singularity of the integro-differential operator for Weber numbers less than unity. The investigation is carried out by means of a modal (i.e. time asymptotic) linear approach, which is corroborated by numerical simulations of the governing equation and supported by systematic comparisons with experimental data from the literature, available in the supercritical regime only. As regards the critical regime for the unit Weber number, the major theoretical result is a sharp increase in oscillation frequency as the flow Weber number is gradually reduced from supercritical to subcritical values due to the shift of the prevailing mode from the slow one to the fast one.

We report on direct numerical simulations to examine the spectral behaviour of turbulence close to and at a flat, stress-free surface. We find, consistent with field measurements near such a free surface, that an inertial-range type of behaviour is obtained for the horizontal components of the velocity at and near the stress-free surface, at horizontal wavelengths for which the vertical velocity is much smaller than the horizontal components. At approximately an integral length scale from the stress-free surface, the flow has adjusted back to more classical isotropic turbulence. The behaviour of the turbulence near the stress-free surface is similar to that observed recently for strongly stratified flows, and we argue that the causes of that behaviour are the same in both flows: the suppression of the large-scale vertical velocity and the allowance of strong vertical shearing of the horizontal velocity leading to a downscale transfer of energy and to the development of the $-5/3$ spectra for the horizontal velocities.

Magnetic droplets obtained by induced phase separation in a magnetic colloid show a large variety of shapes when exposed to an external field. However, the description of the shapes is often limited. Here, we formulate an algorithm based on three-dimensional boundary-integral equations for strongly magnetic droplets in a high-frequency rotating magnetic field, allowing us to find their figures of equilibrium in three dimensions. The algorithm is justified by a series of comparisons with known analytical results. We compare the calculated equilibrium shapes with experimental observations and find a good agreement. The main features of these observations are the oblate–prolate transition, the flattening of prolate shapes with the increase of magnetic field strength and the formation of starfish-like equilibrium shapes. We show both numerically and in experiments that the magnetic droplet behaviour may be described with a triaxial ellipsoid approximation. Directions for further research are mentioned, including the dipolar interaction contribution to the surface tension of the magnetic droplets, accounting for the large viscosity contrast between the magnetic droplet and the surrounding fluid.

A droplet may ‘walk’ across the surface of a vertically vibrating bath of the same fluid, due to the propulsive interaction with its wave field. This hydrodynamic pilot-wave system exhibits many dynamics previously believed to exist only in the quantum realm. Starting from first principles, we derive a discrete-time fluid model, whereby the bath–droplet interactions are modelled as instantaneous. By analysing the stability of the fixed points of the system, we explain the dynamics of a walking droplet and capture the quantisations for multiple-droplet interactions. Circular orbits in a harmonic potential are studied, and a double quantisation of chaotic trajectories is obtained through systematic statistical analysis.

The equilibration time $\unicode[STIX]{x1D70F}$ in response to a change in flux from $Q$ to $\unicode[STIX]{x1D6EC}Q$ after an injection period $T$ applied to either a low-Reynolds-number gravity current or one propagating through a porous medium, in both axisymmetric and one-dimensional geometries, is shown to be of the form $\unicode[STIX]{x1D70F}=Tf(\unicode[STIX]{x1D6EC})$, independent of all the remaining physical parameters. Numerical solutions are used to investigate $f(\unicode[STIX]{x1D6EC})$ for each of these situations and compare very well with experimental results in the case of an axisymmetric current propagating over a rigid horizontal boundary. Analysis of the relaxation towards self-similarity provides an illuminating connection between the excess (deficit) volume from early times and an asymptotically equivalent shift in time origin, and hence a good quantitative estimate of $\unicode[STIX]{x1D70F}$. The case $\unicode[STIX]{x1D6EC}=0$ of equilibration after ceasing injection at time $T$ is a singular limit. Extensions to high-Reynolds-number currents and to the case of a constant-volume release followed by constant-flux injection are discussed briefly.

For the constant-stress layer of wall turbulence, two-point correlations of velocity fluctuations are studied theoretically by using the attached-eddy hypothesis, i.e. a phenomenological model of a random superposition of energy-containing eddies that are attached to the wall. While previous studies had invoked additional assumptions, we focus on the minimum assumptions of the hypothesis to derive the most general forms of the correlation functions. They would allow us to use or assess the hypothesis without any effect of those additional assumptions. We also study the energy spectra and the two-point correlations of the rate of momentum transfer and of the rate of energy dissipation.

The statistical properties of turbulent fluids depend on how local the energy transfers among scales are, i.e. whether the energy transfer at some given scale is due to the eddies at that particular scale, or to eddies at larger (non-local) scale. This locality in the energy transfers may have consequences for the relative dispersion of passive particles. In this paper, we consider a class of generalized two-dimensional flows (produced by the so-called $\unicode[STIX]{x1D6FC}$ -turbulence models), theoretically possessing different properties in terms of locality of energy transfers. It encompasses the standard barotropic quasi-geostrophic (QG) and the surface quasi-geostrophic (SQG) models as limiting cases. The relative dispersion statistics are examined, both as a function of time and as a function of scale, and compared to predictions based on phenomenological arguments assuming the locality of the cascade. We find that the dispersion statistics follow the predicted values from local theories, as long as the parameter $\unicode[STIX]{x1D6FC}$ is small enough (dynamics close to that of the SQG model), for sufficiently small initial pair separations. In contrast, non-local dispersion is observed for the QG model, a robust result when looking at relative displacement probability distributions. However, we point out that spectral energy transfers do have a non-local contribution for models with different values of $\unicode[STIX]{x1D6FC}$ , including the SQG case. This indicates that locality/non-locality of the turbulent cascade may not always imply locality/non-locality in the relative dispersion of particles and that the self-similar nature of the turbulent cascade is more appropriate for determining the relative dispersion locality.

When an acoustic wave travels in a lossy medium such as a liquid, it progressively transfers its pseudo-momentum to the fluid, which results in a steady flow called acoustic streaming. This phenomenon involves a balance between sound attenuation and shear, such that the streaming flow does not vanish in the limit of vanishing viscosity. Hence, the effect of viscosity has long been ignored in acoustic streaming experiments. Here, we investigate the acoustic streaming in sessile droplets exposed to surface acoustic waves. According to experimental data, the flow structure and velocity magnitude are both strongly influenced by the fluid viscosity. We compute the sound wave propagation and hydrodynamic flow motion using a numerical method that reduces memory requirements via a spatial filtering of the acoustic streaming momentum source terms. These calculations agree qualitatively well with experiments and reveal how the acoustic field in the droplet, which is dominated by a few caustics, controls the flow pattern. We evidence that chaotic acoustic fields in droplets are dominated by a few caustics. It appears that the caustics drive the flow, which allows for qualitative prediction of the flow structure. Finally, we apply our numerical method to a broader span of fluids and frequencies. We show that the canonical case of the acoustic streaming in a hemispherical sessile droplet resting on a lithium niobate substrate only depends on two dimensionless numbers related to the surface and bulk wave attenuation. Even in such a baseline configuration, we observe and characterize four distinct flow regimes.

Heat flux and drag correlations are developed for high-speed flow over spherical geometries that are accurate for any Knudsen number ranging from continuum to free-molecular conditions. A stagnation point heat flux correlation is derived as a correction to the continuum (Fourier model) heat flux and also reproduces the correct heat flux in the free-molecular limit by use of a bridging function. In this manner, the correlation can be combined with existing continuum correlations based on computational fluid dynamics simulations, yet it can now be used accurately in the transitional and free-molecular regimes. The functional form of the stagnation point heat flux correlation is physics based, and was derived via the Burnett and super-Burnett equations in a recent article, Singh & Schwartzentruber (J. Fluid Mech., vol. 792, 2016, pp. 981–996). In addition, correlation parameters from the literature are used to construct simple expressions for the local heat flux around the sphere as well as the integrated drag coefficient. A large number of direct simulation Monte Carlo calculations are performed over a wide range of conditions. The computed heat flux and drag data are used to validate the correlations and also to fit the correlation parameters. Compared to existing continuum-based correlations, the new correlations will enable engineering analysis of flight conditions at higher altitudes and/or smaller geometry radii, useful for a variety of applications including blunt body planetary entry, sharp leading edges, low orbiting satellites, meteorites and space debris.

The scalings of the local entrainment velocity $v_{n}$ of the turbulent/non-turbulent interface and of the turbulence dissipation rate are closely related to each other in an axisymmetric and self-similar turbulent wake. The turbulence dissipation scaling implied by the Kolmogorov equilibrium cascade phenomenology is consistent with a Kolmogorov scaling of $v_{n}$ whereas the non-equilibrium dissipation scaling reported for various turbulent flows in Vassilicos (Annu. Rev. Fluid Mech., vol. 47, 2015, pp. 95–114), Dairay et al. (J. Fluid Mech., vol. 781, 2015, pp. 166–195), Goto & Vassilicos (Phys. Lett. A, vol. 379 (16), 2015, pp. 1144–1148) and Obligado et al. (Phys. Rev. Fluids, vol. 1 (4), 2016, 044409) is consistent with a different scaling of $v_{n}$. We present results from a direct numerical simulation of a spatially developing axisymmetric and self-similar turbulent wake which supports this conclusion and the assumptions that it is based on.

Three methods are considered for estimating the downstream evolution of wavepackets in turbulent jets based on upstream measurements. The parabolised stability equations are used to compute a transfer function between axially and radially separated points in the flow, and the performance of this theoretical model is compared with that of two empirical approaches, direct transfer function calculation and autoregressive moving-average exogenous system identification, both of which require unsteady experimental data. The three approaches, which perform equally well, prove suitable for estimation of the downstream evolution of wavepackets using pressure data measured in the near-nozzle region. Over distances of the order of a couple of jet diameters, correlations of up to 80 % are observed between estimation and measurement. The performance deteriorates as axial separation between input and output is increased. While the two empirical approaches are limited in terms of both the number of input–output pairs and the number of flow variables that can be reasonably considered, the parabolised stability equations-based approach has no such limitation and can be used to perform full-field estimates comprising all of the dependent variables; in this it constitutes a potentially formidable means by which to perform single-input–multiple-output estimation. It has the further advantage of not requiring unsteady data for its construction, the only necessary ingredients being the mean flow and the linearised equations of motion.

We study turbulent flows in pressure-driven planar channels with imposed unstable thermal stratification, using direct numerical simulations in a wide range of Reynolds and Rayleigh numbers and reaching flow conditions which are representative of fully developed turbulence. The combined effect of forced and free convection produces a peculiar pattern of quasi-streamwise rollers occupying the full channel thickness, with aspect ratio considerably higher than unity; it has been observed that they have an important redistributing effect on temperature and momentum, providing for a substantial fraction of the heat and momentum flux at bulk Richardson numbers larger than $0.01$. The mean values and the variances of the flow variables do not appear to follow Prandtl’s scaling in the free-convection regime, except for the temperature and vertical velocity fluctuations, which are more directly affected by wall-attached turbulent plumes. We find that the Monin–Obukhov theory nevertheless yields a useful representation of the main flow features. In particular, the widely used Businger–Dyer flux-profile relationships are found to provide a convenient way of approximately accounting for the bulk effects of friction and buoyancy, although the individual profiles may have wide scatter from the alleged trends. Significant deviations are found in direct numerical simulations with respect to the commonly used parametrization of the momentum flux in the light-wind regime, which may have important practical impact in wall models of atmospheric dynamics. Finally, for modelling purposes, we devise a set of empirical predictive formulae for the heat flux and friction coefficients, which are within approximately $10\,\%$ standard deviation from the numerical results in a wide range of flow parameters.

The statistics of rotational motion of small, inertialess triaxial ellipsoids are computed along Lagrangian trajectories extracted from direct numerical simulations of homogeneous isotropic turbulence. The total particle angular velocity and its components along the three principal axes of the particle are considered, expanding on the results of Chevillard & Meneveau (J. Fluid Mech., vol. 737, 2013, pp. 571–596) who showed results of the rotation rate of the particle’s principal axes. The variance of the particle angular velocity, referred to as the particle enstrophy, is found to increase as particles become elongated, regardless of whether they are axisymmetric. This trend is explained by considering the contributions of vorticity and strain rate to particle rotation. It is found that the majority of particle enstrophy is due to fluid vorticity. Strain-rate-induced rotations, which are sensitive to shape, are mostly cancelled by strain–vorticity interactions. The remainder of the strain-rate-induced rotations are responsible for weak variations in particle enstrophy. For particles of all shapes, the majority of the enstrophy is in rotations about the longest axis, which is due to alignment between the longest axis and fluid vorticity. The integral time scale for particle angular velocities about different axes reveals that rotations are most persistent about the longest axis, but that a full revolution is rare.