We present an investigation of the root-mean-square (r.m.s.) temperature $\unicode[STIX]{x1D70E}_{T}$ and the r.m.s. velocity $\unicode[STIX]{x1D70E}_{w}$ in the bulk of Rayleigh–Bénard turbulence, using new experimental data from the current study and experimental and numerical data from previous studies. We find that, once scaled by the convective temperature $\unicode[STIX]{x1D703}_{\ast }$, the value of $\unicode[STIX]{x1D70E}_{T}$ at the cell centre is a constant ($\unicode[STIX]{x1D70E}_{T,c}/\unicode[STIX]{x1D703}_{\ast }\approx 0.85$) over a wide range of the Rayleigh number ($10^{8}\leqslant Ra\leqslant 10^{15}$) and the Prandtl number ($0.7\leqslant Pr\leqslant 23.34$), and is independent of the surface topographies of the top and bottom plates of the convection cell. A constant close to unity suggests that $\unicode[STIX]{x1D703}_{\ast }$ is a proper measure of the temperature fluctuation in the core region. On the other hand, $\unicode[STIX]{x1D70E}_{w,c}/w_{\ast }$, the vertical r.m.s. velocity at the cell centre scaled by the convective velocity $w_{\ast }$, shows a weak $Ra$-dependence (${\sim}Ra^{0.07\pm 0.02}$) over $10^{8}\leqslant Ra\leqslant 10^{10}$ at $Pr\sim 4.3$ and is independent of plate topography. Similar to a previous finding by He & Xia (Phys. Rev. Lett., vol. 122, 2019, 014503), we find that the r.m.s. temperature profile $\unicode[STIX]{x1D70E}_{T}(z)/\unicode[STIX]{x1D703}_{\ast }$ in the region of the mixing zone with a mean horizontal shear exhibits a power-law dependence on the distance $z$ from the plate, but now the universal profile applies to both smooth and rough surface topographies and over a wider range of $Ra$. The vertical r.m.s. velocity profile $\unicode[STIX]{x1D70E}_{w}(z)/w_{\ast }$ obeys a logarithmic dependence on $z$. The study thus demonstrates that the typical scales for the temperature and the velocity are the convective temperature $\unicode[STIX]{x1D703}_{\ast }$ and the convective velocity $w_{\ast }$, respectively. Finally, we note that $\unicode[STIX]{x1D703}_{\ast }$ may be utilised to study the flow regime transitions in ultrahigh-$Ra$-number turbulent convection.

In this paper we contrast the early impact stage of a highly viscous drop onto a liquid versus a solid substrate. Water drops impacting at low velocities can rebound from a solid surface without contact. This dynamic is mediated through lubrication of a thin air layer between the liquid and solid. Drops can also rebound from a liquid surface, but only for low Weber numbers. Impacts at higher velocities in both cases lead to circular contacts which entrap an air disc under the centre of the drop. Increasing the drop viscosity produces extended air films for impacts on a smooth solid surface even for much larger velocities. These air films eventually break through random wetting contacts with the solid. Herein we use high-speed interferometry to study the extent and thickness profile of the air film for a large-viscosity drop impacting onto a viscous film of the same liquid. We demonstrate a unified scaling of the centreline height of the air film for impacts on both solid and liquid, when using the effective impact velocity. On the other hand, we show that the large-viscosity liquid film promotes air films of larger extent. Furthermore, the rupture behaviour becomes fundamentally different, with the air film between the two compliant surfaces being more stable, lacking the random wetting patches seen on the solid. We map the parameter range where these air films occur and explore the transition from gliding to ring contact at the edge of the drop dimple. After the air film ruptures, the initial contraction occurs very rapidly and for viscosities greater than 100 cSt the retraction velocity of the air film is ${\sim}0.3~\text{m}~\text{s}^{-1}$, independent of the liquid viscosity and impact velocity, in sharp contrast with theoretical predictions.

In elastic macromolecules, the value of the short-time diffusion coefficient depends on the choice of the point the displacement of which is tracked. On the other hand, the experimentally more relevant long-time diffusion coefficient is independent of the reference point, but its estimation usually requires computationally expensive Brownian dynamics simulations. Here we show how to obtain a precise estimate of the long-time diffusion coefficient of elastic macromolecules in a fast and robust manner, without invoking Brownian dynamics.

Unstable shear layers in environmental and industrial flows roll up into a series of vortices, which often form complex nonlinear merging patterns such as pairs and triplets. These patterns crucially determine the subsequent turbulence, mixing and scalar transport. We show that the late-time, highly nonlinear merging patterns are predictable from the linearized initial state. The initial asymmetry between consecutive wavelengths of the vertical velocity field provides an effective measure of the strength and pattern of vortex merging. The predictions of this measure are substantiated using direct numerical simulations. We also show that this measure has significant implications in determining the route to turbulence and the ensuing turbulence characteristics.

An exact solution of a stable vortex tripole in two-dimensional (2-D) Euler flows is provided. The stable tripole is composed of an inner elliptical vortex and two small-amplitude lateral vortices. The non-vanishing vorticity field of this tripole, referred to as here as an embedded tripole because of the closeness of its vortices, is given in elliptical coordinates $(\unicode[STIX]{x1D707},\unicode[STIX]{x1D708})$ by the even radial and angular order-0 Mathieu functions $\text{Je}_{0}(\unicode[STIX]{x1D707})\text{ce}_{0}(\unicode[STIX]{x1D708})$ truncated at the external branch of the vorticity isoline passing through the two critical points closest to the vortex centre. This tripole mode has a rigid vorticity field which rotates with constant angular velocity equal to $\unicode[STIX]{x1D701}_{0}\text{Je}_{0}(\unicode[STIX]{x1D707}_{1})\text{ce}_{0}(0)/2$, where $\unicode[STIX]{x1D707}_{1}$ is the first zero of $\text{Je}_{0}^{\prime }(\unicode[STIX]{x1D707})$ and $\unicode[STIX]{x1D701}_{0}$ is a constant modal amplitude. It is argued that embedded 2-D tripoles may be conceptually regarded as the superposition of two asymmetric Chaplygin–Lamb dipoles, separated a distance equal to $2R$, as long as their individual trajectory curvature radius $R$ is much shorter than their dipole extent radius.

Lipid membranes are examples of fluid deformable surfaces, which can be viewed as two-dimensional viscous fluids with bending elasticity. With this solid–fluid duality any shape change contributes to tangential flow and vice versa any tangential flow on a curved surface induces shape deformations. This tight coupling between shape and flow makes curvature a natural element of the governing equations. The modelling and numerical tools outlined in Torres-Sánchez et al. (J. Fluid Mech., vol. 872, 2019, pp. 218–271) open a new field of study by enabling the exploration of the role of curvature in this context.

The linear instability of a vortex ring with swirl with Gaussian distributions of azimuthal vorticity and velocity in its core is studied by direct numerical simulation. The numerical study is carried out in two steps: first, an axisymmetric simulation of the Navier–Stokes equations is performed to obtain the quasi-steady state that forms a base flow; then, the equations are linearized around this base flow and integrated for a sufficiently long time to obtain the characteristics of the most unstable mode. It is shown that the vortex rings are subjected to curvature instability as predicted analytically by Blanco-Rodríguez & Le Dizès (J. Fluid Mech., vol. 814, 2017, pp. 397–415). Both the structure and the growth rate of the unstable modes obtained numerically are in good agreement with the analytical results. However, a small overestimation (e.g. 22 % for a curvature instability mode) by the theory of the numerical growth rate is found for some instability modes. This is most likely due to evaluation of the critical layer damping which is performed for the waves on axisymmetric line vortices in the analysis. The actual position of the critical layer is affected by deformation of the core due to the curvature effect; as a result, the damping rate changes since it is sensitive to the position of the critical layer. Competition between the curvature and elliptic instabilities is also investigated. Without swirl, only the elliptic instability is observed in agreement with previous numerical and experimental results. In the presence of swirl, sharp bands of both curvature and elliptic instabilities are obtained for $\unicode[STIX]{x1D700}=a/R=0.1$, where $a$ is the vortex core radius and $R$ the ring radius, while the elliptic instability dominates for $\unicode[STIX]{x1D700}=0.18$. New types of instability mode are also obtained: a special curvature mode composed of three waves is observed and spiral modes that do not seem to be related to any wave resonance. The curvature instability is also confirmed by direct numerical simulation of the full Navier–Stokes equations. Weakly nonlinear saturation and subsequent decay of the curvature instability are also observed.

Turbulent structures in stably stratified shear layers are studied with direct numerical simulation. Flow visualization confirms the existence of hairpin vortices and highly elongated structures with positive and negative velocity fluctuations, whose streamwise lengths divided by the layer thickness are $O(10^{0})$ and $O(10^{1})$, respectively. The flow at the wavelength related to these structures makes a large contribution to turbulent kinetic energy. These structures become prominent in late time, but with small buoyancy Reynolds numbers indicating suppression of turbulent mixing. Active turbulent mixing associated with the hairpin vortices, however, does occur. The structures and the vertical profile of the integral shear parameter show connections between stable stratified shear layers and wall-bounded shear flows.

An Eulerian–Lagrangian modelling approach was employed in order to investigate the flow field, heat transfer and particle distribution in nanofluid flow in a parallel-plate microchannel, with a focus on relatively low Reynolds numbers ($Re\leqslant 100$). Momentum and thermal interactions between fluid and particle phases were accounted for using a transient two-way coupling algorithm implemented within an in-house code that tracked the simultaneous evolution of the carrier and particulate phases while considering timescale differences between the two phases. The inaccuracy of assuming a homogeneous particle distribution in modelling nanofluid flow in microchannels was established. In particular, shear rate and thermophoresis were found to play a key role in the lateral migration of nanoparticles and in the formation of particle depletion and accumulation regions in the vicinity of the channel walls. At low Reynolds numbers, nanoparticle distribution near the walls was observed to gradually flatten in the streamwise direction. On the other hand, for relatively higher Reynolds numbers, higher particle non-uniformities were observed in the vicinity of the channel walls. Furthermore, it was established that convective heat transfer between channel walls and the bulk fluid can either improve or deteriorate with the addition of nanoparticles, depending on whether the flow exceeded a critical Reynolds number of enhancement. It was also established that Brownian motion and thermophoresis had a major role in nanoparticle deposition on the channel walls. In particular, Brownian motion was the main deposition mechanism for nano-sized particles, whereas due to thermophoresis, nanoparticles were repelled away from channel walls. The result of the competition between the two is that deposition gradually increased along the streamwise direction.

Fragmentation of a vertical buoyant silicone oil jet injected into sugar water is elucidated by refractive index matching and planar laser-induced fluorescence. Compound droplets containing multiple water droplets, some with smaller oil droplets, form regularly at jet Reynolds numbers of $Re=1358$ and 2122 and persist for at least up to 30 nozzle diameters. In contrast, they rarely appear at $Re=594$. The origin of some of the encapsulated water droplets can be traced back to the entrained water ligaments during the initial roll-up of Kelvin–Helmholtz vortices. Analysis using random forest-based procedures shows that the fraction of compound droplets does not vary significantly with $Re$, but increases rapidly with droplet diameter, reaching 78 % for 2 mm droplets. Consequently, the size distributions of compound droplets have peaks that increase in magnitude and shift to a lower diameter with increasing $Re$. On average, the interior pockets raise the oil–water interfacial area by 15 %, increasing with diameter and axial location. Also, while the oil droplets are deformed by the jet’s shear field, the interior interfaces remain nearly spherical, consistent with prior studies of the deformation of isolated compound droplets for relevant capillary numbers and viscosity ratio.

We numerically investigate the passive flight of a flexible $\unicode[STIX]{x039B}$-flyer in a vertically oscillating airflow with zero mean stream. The flexibility of the flyer is introduced by a torsion spring installed at the hinged joint. We study the effects of spring stiffness, density, resting angle and actuation efforts on the hovering performance. The results suggest that the occurrence of resonance in flexible flyers can result in significantly different performances in flexible and rigid flyers. It is found that flexibility can have two opposing effects, reducing or increasing the actuation efforts for hovering, depending on the range of driving frequency. This result is explained by the modulation of relative motion between the flyer and the imposed background flow due to the involvement of passive angular oscillation. The angular oscillation patterns, the wake symmetry properties and the postural stability behaviours under different driving conditions are also explored. Based on the findings of the present study, the ideal parameter values for stable hovering are suggested. The results of this study offer novel insight into the mechanism by which the flexibility of the flyer affects the passive hovering performance.

An asymmetric body with a sharp leading edge and a rounded trailing edge produces a smaller wave disturbance moving forwards than backwards, and this is reflected in the wave drag coefficient. This experimental fact is not captured by Michell’s theory for wave drag (Michell Lond. Edinb. Dubl. Phil. Mag. J. Sci., vol. 45 (272), 1898, pp. 106–123). In this study, we use a tow-tank experiment to investigate the effects of asymmetry on wave drag, and show that these effects can be replicated by modifying Michell’s theory to include the growth of a symmetry-breaking boundary layer. We show that asymmetry can have either a positive or a negative effect on drag, depending on the depth of motion and the Froude number.

Thin viscous liquid films sitting on a solid substrate support nonlinear capillary waves, driven by surface shear stresses at a liquid–gas interface. When surface tension is spatially dependent other mechanisms, such as the thermocapillary effect, influence the dynamics of thin films. In this article we show that in liquids with broken time-reversal symmetry the character of the aforementioned waves and of the thermocapillary effect are significantly modified due to the presence of odd or Hall viscosity in the liquid. This is because odd viscosity gives rise to new terms in the pressure gradient of the flow thus modifying the evolution equation of the liquid–gas interface accordingly. These terms in turn break the reflection symmetry of the evolution equation leading the system to evolve from a pitchfork to a Hopf bifurcation. The odd-viscosity incipient waves can stabilize unstable thin liquid films. For instance, we show that they can suppress the thermocapillary instability. We establish the parameter ranges that odd viscosity has to satisfy in order to initiate those waves that will lead to stability.

This study analyses the interplay between classical acoustic modes and intrinsic thermoacoustic (ITA) modes in a simple thermoacoustic system. The analysis is performed using a frequency-domain low-order network model as well as a time-domain spatially discretised model. Anti-correlated modal sensitivities are found to arise due to a pairwise interplay between acoustic and ITA modes. The magnitude of the sensitivities increases as the interplay between the modes grows stronger. The results show a global behaviour of the modes linked to the presence of exceptional points in the spectrum. The time-domain analysis results in a delay-differential equation and allows the investigation of non-normal behaviour and its consequences. Pseudospectral analysis reveals that energy amplification is crucially linked to an interplay between acoustic and ITA modes. While higher non-orthogonality between two modes is correlated with peaks in modal sensitivity, transient energy growth does not necessarily involve the most sensitive modes. In particular, growth estimates based on the Kreiss constant demonstrate that transient amplification relies critically on the proximity of the non-normal modes to the imaginary axis. The time scale for transient amplification is identified as the flame time delay, which is further corroborated by determining the optimal initial conditions responsible for the bulk of the non-modal energy growth. The flame is identified as an active and dominant contributor to energy gain. The frequency of the optimal perturbation matches the acoustic time scale, once more confirming an interplay between acoustic and ITA structures. Flame-based amplification factors of two to five are found, which are significant when feeding into the acoustic dynamics and eventually triggering nonlinear limit-cycle behaviour.

We perform direct numerical simulations (DNS) of emulsions in homogeneous isotropic turbulence using a pseudopotential lattice-Boltzmann (PP-LB) method. Improving on previous literature by minimizing droplet dissolution and spurious currents, we show that the PP-LB technique is capable of long stable simulations in certain parameter regions. Varying the dispersed-phase volume fraction $\unicode[STIX]{x1D719}$, we demonstrate that droplet breakup extracts kinetic energy from the larger scales while injecting energy into the smaller scales, increasingly with higher $\unicode[STIX]{x1D719}$, with approximately the Hinze scale (Hinze, AIChE J., vol. 1 (3), 1955, pp. 289–295) separating the two effects. A generalization of the Hinze scale is proposed, which applies both to dense and dilute suspensions, including cases where there is a deviation from the $k^{-5/3}$ inertial range scaling and where coalescence becomes dominant. This is done using the Weber number spectrum $We(k)$, constructed from the multiphase kinetic energy spectrum $E(k)$, which indicates the critical droplet scale at which $We\approx 1$. This scale roughly separates coalescence and breakup dynamics as it closely corresponds to the transition of the droplet size ($d$) distribution into a $d^{-10/3}$ scaling (Garrett et al., J. Phys. Oceanogr., vol. 30 (9), 2000, pp. 2163–2171; Deane & Stokes, Nature, vol. 418 (6900), 2002, p. 839). We show the need to maintain a separation of the turbulence forcing scale and domain size to prevent the formation of large connected regions of the dispersed phase. For the first time, we show that turbulent emulsions evolve into a quasi-equilibrium cycle of alternating coalescence and breakup dominated processes. Studying the system in its state-space comprising kinetic energy $E_{k}$, enstrophy $\unicode[STIX]{x1D714}^{2}$ and the droplet number density $N_{d}$, we find that their dynamics resemble limit cycles with a time delay. Extreme values in the evolution of $E_{k}$ are manifested in the evolution of $\unicode[STIX]{x1D714}^{2}$ and $N_{d}$ with a delay of ${\sim}0.3{\mathcal{T}}$ and ${\sim}0.9{\mathcal{T}}$ respectively (with ${\mathcal{T}}$ the large eddy timescale). Lastly, we also show that flow topology of turbulence in an emulsion is significantly more different from single-phase turbulence than previously thought. In particular, vortex compression and axial straining mechanisms increase in the droplet phase.

The evolution of an $\text{SF}_{6}$ layer surrounded by air is experimentally studied in a semi-annular convergent shock tube by high-speed schlieren photography. The gas layer with a sinusoidal outer interface and a circular inner interface is realized by the soap-film technique such that the initial condition is well controlled. Results show that the thicker the gas layer, the weaker the interface–coupling effect and the slower the evolution of the outer interface. Induced by the distorted transmitted shock and the interface coupling, the inner interface exhibits a slow perturbation growth which can be largely suppressed by reducing the layer thickness. After the reshock, the inner perturbation increases linearly at a growth rate independent of the initial layer thickness as well as of the outer perturbation amplitude and wavelength, and the growth rate can be well predicted by the model of Mikaelian (Physica D, vol. 36, 1989, pp. 343–357) with an empirical coefficient of 0.31. After the linear stage, the growth rate decreases continuously, and finally the perturbation freezes at a constant amplitude caused by the successive stagnation of spikes and bubbles. The convergent geometry constraint as well as the very weak compressibility at late stages are responsible for this instability freeze-out.

This paper presents a numerical performance evaluation of the leading edge tubercles hydrofoil with particular focus on a fully turbulent flow regime. Efforts were focused on the setting up of an appropriate numerical approach required for an in-depth analysis of this phenomenon, being able to predict the main flow features and the hydrodynamic performance of the foil when operating at high Reynolds numbers. The numerical analysis was conducted using an improved delayed detached eddy simulation for Reynolds numbers corresponding to the transitional and fully turbulent flow regimes at different angles of attack for the pre-stall and post-stall regimes. The results show that tubercles operating in turbulent flow improve the hydrodynamic performance of the foil when compared to a transitional flow regime. Flow separation was identified behind the tubercle troughs, but was significantly reduced when operating in a turbulent regime and for which we have identified the main flow mechanisms. This finding confirms that the tubercle effect identified in a transitional regime is not lost in a turbulent flow. Furthermore, when the hydrofoil operates in the turbulent flow regime, the transition to a turbulent regime takes place further upstream. This phenomenon suppresses a formation of a laminar separation bubble and therefore the hydrofoil exhibits a superior hydrodynamic performance when compared to the same foil in the transitional regime.

Turbulent flow evolution and energy cascades are significantly different in two-dimensional (2-D) and three-dimensional (3-D) flows. Studies have investigated these differences in obstacle-free turbulent flows, but solid boundaries have an important impact on the cross-over from 3-D to 2-D turbulence dynamics. In this work, we investigate the span effect on the turbulence nature of flow past a circular cylinder at $Re=10\,000$. It is found that even for highly anisotropic geometries, 3-D small-scale structures detach from the walls. Additionally, the natural large-scale rotation of the Kármán vortices rapidly two-dimensionalise those structures if the span is 50 % of the diameter or less. We show this is linked to the span being shorter than the Mode B instability wavelength. The conflicting 3-D small-scale structures and 2-D Kármán vortices result in 2-D and 3-D turbulence dynamics which can coexist at certain locations of the wake depending on the domain geometric anisotropy.

The interfacial behaviour of surfactant-laden drops squeezing through tight constrictions in a uniform far-field flow is modelled with respect to capillary number, drop-to-medium viscosity ratio and surfactant contamination. The surfactant is treated as insoluble and non-diffusive, and drop surface tension is related to surfactant concentration by a linear equation of state. The constriction is formed by three solid spheres held rigidly in space. A characteristic aspect of this confined and contaminated multiphase system is the rapid development of steep surfactant-concentration gradients during the onset of drop squeezing. The interplay between two physical effects of surfactant, namely the greater interface deformability due to decreased surface tension and interface immobilization due to Marangoni stresses, results in particularly rich drop-squeezing dynamics. A three-dimensional boundary-integral algorithm is used to describe drop hydrodynamics, and accurate treatment of close squeezing and trapped states is enabled by advanced singularity subtraction techniques. Surfactant transport and hydrodynamics are coupled via the surface convection equation (or convection–diffusion equation, if artificial diffusion is included), the interfacial stress balance and a solid-particle contribution based on the Hebeker representation. For extreme conditions, such as drop-to-medium viscosity ratios significantly less than unity, it is found that upwind-biased methods are the only stable approaches for modelling surfactant transport. Two distinct schemes, upwind finite volume and flow-biased least squares, are found to provide results in close agreement, indicating negligible numerical diffusion. Surfactant transport is enhanced by low drop-to-medium viscosity ratios, at which extremely sharp concentration gradients form during various stages of the squeezing process. The presence of surfactant, even at low degrees of contamination, significantly decreases the critical capillary number for droplet trapping, due to the accumulation of surfactant at the downwind pole of the drop and its subsequent elongation. Increasing the degree of contamination significantly affects surface mobility and further decreases the critical capillary number as well as drop squeezing times, up to a threshold above which the addition of surfactant negligibly affects squeezing dynamics.

Turbulent flows in a thin layer can develop an inverse energy cascade leading to spectral condensation of energy when the layer height is smaller than a certain threshold. These spectral condensates take the form of large-scale vortices in physical space. Recently, evidence for bistability was found in this system close to the critical height: depending on the initial conditions, the flow is either in a condensate state with most of the energy in the two-dimensional (2-D) large-scale modes, or in a three-dimensional (3-D) flow state with most of the energy in the small-scale modes. This bistable regime is characterised by the statistical properties of random and rare transitions between these two locally stable states. Here, we examine these statistical properties in thin-layer turbulent flows, where the energy is injected by either stochastic or deterministic forcing. To this end, by using a large number of direct numerical simulations (DNS), we measure the decay time $\unicode[STIX]{x1D70F}_{d}$ of the 2-D condensate to 3-D flow state and the build-up time $\unicode[STIX]{x1D70F}_{b}$ of the 2-D condensate. We show that both of these times $\unicode[STIX]{x1D70F}_{d},\unicode[STIX]{x1D70F}_{b}$ follow an exponential distribution with mean values increasing faster than exponentially as the layer height approaches the threshold. We further show that the dynamics of large-scale kinetic energy may be modelled by a stochastic Langevin equation. From time-series analysis of DNS data, we determine the effective potential that shows two minima corresponding to the 2-D and 3-D states when the layer height is close to the threshold.