Inspired by natural interfaces with surprising transport properties, innovative modifications of surfaces have been engineered to reduce drag. The common theme across these new developments is the presence of lubricant patches or layers that decrease the direct contact of viscous liquid with non-slippery solid walls. For laminar flow, the traditional assumption regarding the lubricant layer is a constant shear rate or a steady pressure gradient, implying a net flow rate of the lubricant film. By challenging this assumption, Busse et al. (J. Fluid Mech., vol. 727, 2013, pp. 488–508) rigorously found that the hydrodynamic slip is reduced by the presence of a reversal of lubricant flow close to the wall. The analytical results for velocity field and change in drag provide insight into the optimal design of slippery surfaces with lubricant layers for drag reduction.

Recent results published by Gugenberger et al. on surface diffusion (Phys. Rev. E, vol. 78, 2008, 016703), show that the sharp-interface limit of the phase field models often adopted in the literature fails to produce the appropriate boundary conditions. With this knowledge, we consider the sharp-interface limit of phase field models for binary fluids, obtained carefully, where hydrodynamic equations are coupled to phase field evolution based on Cahn–Hilliard or Allen–Cahn theories, in a variety of guises, and unify and contrast their forms and behaviours in the sharp-interface limit. In particular, a tensorial mobility model is analysed, which allows the bulk fluids in the outer region to satisfy classical Navier–Stokes type equations to all orders in the Cahn number.

The thermodynamic effects associated with the growth and collapse of a single cavitation bubble are investigated in the present paper by an experimental approach. The study focuses on the temperature variations in the liquid surrounding the bubble. Experiments are conducted in a cylinder partially filled with water at an ambient temperature and atmospheric pressure. The bubble growth results from the expansion of an initial air bubble, due to the pressure wave generated by a so-called ‘tube-arrest’ method. Several locations of the bubble, at different distances from the bottom wall of the cylinder, are considered. The bottom wall is made of sapphire, which is transparent to both the visible and infrared light spectra which enables temperature measurements by a high-speed thermovision camera at a wavelength of 3–$5~\unicode[.5,0][STIXGeneral,Times]{x03BC} \mathrm{m} $. Water is opaque to the infrared light spectrum, hence only temperatures in the boundary layer and on the liquid vapour interface could be determined. A temperature decrease of ${\sim }3$ K was recorded during the bubble growth while an increase up to 4 K was detected during the collapse. Experimental results are compared to the predictions of the ‘thermal delay’ model based on the assumption that the bubble growth and collapse are due to phase changes only. In this approach, the temperature variations are related to the latent heat exchanges during the vapourization and condensation processes. On the basis of these results, the respective effects of phase change and air dilatation/compression in the bubble dynamics are discussed.

The linear stability of uniform, plane internal wave beams with locally confined spatial profile, in a stratified fluid of constant buoyancy frequency, is discussed. The associated eigenvalue problem is solved asymptotically, assuming perturbations of long wavelength relative to the beam width. In this limit, instability is found only for oblique disturbances which vary in the along-beam and the horizontal transverse directions. The mechanism of instability is a first-harmonic–mean resonant interaction between the underlying wave beam and three-dimensional perturbations that comprise a time-harmonic component, with the beam frequency, and a mean flow. Progressive beams which transport energy in one direction, in particular, are unstable if the beam steepness exceeds a certain threshold value, whereas purely standing beams are unstable even at infinitesimal steepness. A distinguishing feature of this three-dimensional modulational instability is the generation of circulating horizontal mean flows at large distances from the vicinity of the beam.

An investigation into redistribution of vorticity for rapidly accelerating plates with varying kinematics and initial conditions has been performed. Both three-dimensional particle tracking velocimetry and direct force measurements were applied simultaneously. The effective velocity of the feeding shear layer has been identified as the appropriate characteristic velocity rather than the commonly used plunge or free stream velocity. Based on this new normalization for circulation, it has been demonstrated that the existence of initial boundary-layer vorticity on the plunging plate – at least in the near-midplane region – does not contribute to the eventual vortex formation process. In accordance with the literature, however, the tip vortex positioning relative to the plate surface has been identified as an important contributor in the overall force production, particularly once the plate acceleration has ceased.

We consider thin liquid films on the inside (rimming flows) or outside (coating flows) of a cylinder with horizontal axis, rotating about this axis. If the liquid’s net volume is small, such films are known to evolve towards a steady state with a smooth surface, whereas, for larger amounts, the flow develops a ‘shock’ similar to a tidal bore. In this work, smooth films are shown to be unstable. Since the strongest instability occurs at wavelengths comparable to the film’s thickness, our analysis is based on the full Navier–Stokes equations, not on the lubrication approximation (which has been traditionally used in this problem). It is also shown that, for cylinders of sufficiently small radii, the instability can be suppressed by surface tension.

Characterizing scalar dispersion is a key concern in a wide variety of applications, including both steady-state and time-dependent studies of wastewater outfalls, salinity distribution in estuaries, and the spreading of pollutants from industrial spills. As the size of a scalar plume grows with respect to the size of the containing water body, the effective dispersion varies, from the well-known ${ \sigma }_{x}^{2} \sim {t}^{3} $ behaviour for a plume enveloped in a region of linear shear, to the ${ \sigma }_{x}^{2} \sim t$ behaviour at the limit of a laterally well-mixed plume. We introduce an additional regime in which the plume extends across the full range of the available shear, but is not significantly affected by the lateral bounds of the water body. Through an analytic treatment we show that this regime exhibits a ${ \sigma }_{x}^{2} \sim {t}^{2} $ behaviour, independent of lateral mixing coefficient. Particle tracking results in an idealized, tidal channel–shoal basin demonstrate this regime as particle clouds straddle the channel–shoal interface. Quantitative analysis of spatial moments as plumes transition between regimes show good correlation between the observed parameters and parameters predicted by the analytic framework.

This article investigates the appearance of instabilities in two planar coflowing fluid sheets with different densities and viscosities via experiments, numerical simulation and linear stability analysis. At low dynamic pressure ratios a convective instability is shown to appear for which the frequency of the waves in the primary atomization region is influenced by both liquid and gas velocities. For large dynamic pressure ratios an asymptotic regime is obtained in which frequency is solely controlled by gas velocity and the instability becomes absolute. The transition from convective to absolute is shown to be influenced by the velocity defect induced by the presence of the separator plate. We show that in this regime the splitter plate thickness can also affect the nature of the instability if it is larger than the gas vorticity thickness. Computational and experimental results are in agreement with the predictions of a spatio-temporal stability analysis.

Two-dimensional and three-dimensional Rayleigh–Bénard convection is compared using results from direct numerical simulations and previous experiments. The phase diagrams for both cases are reviewed. The differences and similarities between two- and three-dimensional convection are studied using $Nu(Ra)$ for $\mathit{Pr}= 4. 38$ and $\mathit{Pr}= 0. 7$ and $Nu(Pr)$ for $Ra$ up to $1{0}^{8} $. In the $Nu(Ra)$ scaling at higher $Pr$, two- and three-dimensional convection is very similar, differing only by a constant factor up to $\mathit{Ra}= 1{0}^{10} $. In contrast, the difference is large at lower $Pr$, due to the strong roll state dependence of $Nu$ in two dimensions. The behaviour of $Nu(Pr)$ is similar in two and three dimensions at large $Pr$. However, it differs significantly around $\mathit{Pr}= 1$. The Reynolds number values are consistently higher in two dimensions and additionally converge at large $Pr$. Finally, the thermal boundary layer profiles are compared in two and three dimensions.

We study the conditions under which discontinuous mechanical properties of a collapsible tube can induce transcritical flows, i.e. the transition through the critical state where the speed index (analogous to the Mach or the Froude numbers for compressible and free surface flows, respectively) is one. Such a critical transition may strongly modify the flow properties, cause a significant reduction in the cross-sectional area of the tube, and limit the flow. General relationships are obtained for a short segment using a one-dimensional model under steady flow conditions. Marginal curves delimiting the transcritical regions are identified in terms of the speed index and the cross-sectional area ratio. Since there are many examples of such flows in physiology and medicine, we also analyse the specific case of prosthesis (graft or stent) implantation in blood vessels. We then compute transcritical conditions for the case of stiffness and reference area variations, considering a collapsible tube characterized by physiological parameters representative of both arteries and veins. The results suggest that variations in mechanical properties may induce transcritical flow in veins but is unrealistic in arteries.

We investigate the stability of an elastic plate in supersonic gas flow. This problem has been studied in many papers regarding panel flutter, where uniform flow is usually considered. In this paper, we take the boundary layer on the plate into account and investigate its influence on plate stability. Three problem formulations are studied. First, we investigate the stability of travelling waves in an infinite-length plate. Second, the nature of the instability (absolute or convective instability) is examined. Finally, by using solutions of the first two problems, instability of a long finite-length plate is studied by using Kulikovskii’s global instability criterion. The following results are obtained. All the eigenmodes of a finite-length plate are split into two types, which we call subsonic and supersonic. The influence of the boundary layer on these eigenmodes can be of two kinds. First, for a generalized convex boundary layer profile (typical for accelerating flow), supersonic eigenmodes are stabilized by the boundary layer, whereas subsonic disturbances are destabilized. Second, for a profile with a generalized inflection point (typical for constant and decelerating flows), supersonic eigenmodes are destabilized in a thin boundary layer and stabilized in a thick layer; subsonic eigenmodes are damped. The correspondence between the influence of the boundary layer on panel flutter and the stability of the boundary layer over a rigid wall is established. Examples of stable boundary layer profiles of both types are given.

The effect of different levels of background turbulence on the dynamics and mixing of an axisymmetric turbulent jet at different Reynolds numbers has been investigated. Approximately homogeneous and isotropic background turbulence was generated by a random jet array and had a negligible mean flow (${\langle {U}_{\alpha } \rangle }/ {u}_{\alpha \mathit{rms}} \ll 1$). Velocity measurements of a jet issuing into two different levels of background turbulence were conducted for three different jet Reynolds numbers. The results showed that the mean axial velocities decay faster with increasing level of background turbulence (compared with a jet in quiescent surroundings), while the mean radial velocities increase, especially close to the edges of the jet. Furthermore, the axial root-mean-square velocities of the jet increased in the presence of background turbulence, as did the jet’s width. However, the mass flow rate of the jet decreased, from which it can be inferred that the entrainment into the jet is reduced in a turbulent background. The effect of background turbulence on the entrainment mechanisms is discussed.

The geologic sequestration of carbon dioxide ( CO2) in structural and stratigraphic traps is a viable option to reduce anthropogenic emissions. While dissolution of the CO2 stored in these traps reduces the long-term leakage risk, the dissolution process remains poorly understood in systems that reflect the appropriate subsurface geometry. Here, we study dissolution in a porous layer that exhibits a feature relevant for CO2 storage in structural and stratigraphic traps: a finite CO2 source along the top boundary that extends only part way into the layer. This feature represents the finite extent of the interface between free-phase CO2 pooled in a trap and the underlying brine. Using theory and simulations, we describe the dissolution mechanisms in this system for a wide range of times and Rayleigh numbers, and classify the behaviour into seven regimes. For each regime, we quantify the dissolution flux numerically and model it analytically, with the goal of providing simple expressions to estimate the dissolution rate in real systems. We find that, at late times, the dissolution flux decreases relative to early times as the flow of unsaturated water to the CO2 source becomes constrained by a lateral exchange flow though the reservoir. Application of the models to several representative reservoirs indicates that dissolution is strongly affected by the reservoir properties; however, we find that reservoirs with high permeabilities ($k\geq 1$ Darcy) that are tens of metres thick and several kilometres wide could potentially dissolve hundreds of megatons of CO2 in tens of years.

We experimentally investigate the three-dimensional wake behind a finite wall-mounted square cylinder at $\mathit{Re}= 12\hspace{0.167em} 000$ and aspect ratio of 4. Focus is placed on the base flow and oscillatory fluctuation. Time-resolved three-dimensional velocity fields are constructed from high-frame-rate particle image velocimetry (PIV) and simultaneously recorded surface pressure measurements. All three velocity components are resolved in a rectangular near-wake region by two orthogonal dense arrays of parallel PIV planes. A key enabler is a generalized phase average incorporating a slowly varying base flow, a variable oscillation amplitude and higher harmonics. These generalizations reduce the instantaneous residual 30 % below those of a traditional phase average. Moreover, the resolved variations reveal analytical constraints of the mean flow and oscillation levels, such as the mean-field paraboloid. The proposed methodology for generalized phase averaging and for construction of three-dimensional velocity fields from two-dimensional PIV data is applicable to a large class of turbulent flows with oscillatory dynamics.

Red blood cells are the subject of diverse studies. One branch is the observation and theoretical modelling of their behaviour in a shear flow. This work deals with the flow of single red cells suspended in solutions much more viscous than blood plasma. Below a critical shear rate (${\dot {\gamma } }_{t} $) the red cells rotate with little change of their resting shape. Above that value they become elongated and aligned in the shear field. We measured ${\dot {\gamma } }_{t} $ at viscosities (${\eta }_{0} $) ranging from 10.7 to 104 mPa s via observation along the vorticity of a Poiseuille flow in a glass capillary; ${\eta }_{0} {\dot {\gamma } }_{t} $ decreased steeply with increasing ${\eta }_{0} $ up to a value of 25 mPa s and remained constant for higher values. Present theoretical models are not in keeping with the measured data. Modifications of basic model assumptions are suggested.

Fronts, or regions with strong horizontal density gradients, are ubiquitous and dynamically important features in the ocean and atmosphere. In the atmosphere, fronts are associated with some of the most severe weather events, while in the ocean, fronts are associated with enhanced turbulence, water mass transformation and biological activity. Here, we examine the dynamics involved in the formation of fronts, or frontogenesis, in detail using a generalized mathematical framework. This extends previous work which has generally revolved around two limiting cases: fronts generated through forcing due to a convergent large-scale flow, and fronts generated spontaneously during the geostrophic adjustment of an initially unbalanced flow. Here, we introduce a new generalized momentum coordinate to simultaneously describe forced and spontaneous frontogenesis. The nonlinear, inviscid, Boussinesq, hydrostatic governing equations for uniform PV flow are solved for arbitrary Rossby and Froude number. The solution is then examined in three distinct cases. Firstly, for a zero potential vorticity (PV) flow bounded by rigid lids, a general solution is derived for the transient response of the fluid to an arbitrary initial mass imbalance and deformation field. The deformation frontogenesis solution of Hoskins & Bretherton (J. Atmos. Sci., vol. 29, 1972, pp. 11–37) and the mass imbalance solution of Blumen (J. Phys. Oceanogr., vol. 30, 2000, pp. 31–39) emerge as two limits of this general solution. Secondly, the problem of geostrophic adjustment of an initial mass imbalance (no deformation field) is considered for uniform PV flow bounded by rigid lids. The general solution is derived, composed of an adjusted state and a transient component describing the propagation of inertia–gravity waves. The criteria for the occurrence of a frontal discontinuity is determined in terms of the Rossby and Froude numbers. The uniform PV solution reduces identically to the zero PV solution of Blumen in the limit of vanishing background stratification. Thirdly, we examine the more general case of uniform PV flow with a deformation field and either balanced or unbalanced initial conditions. In this case the solution is composed of a time-varying mean state – matching the Hoskins & Bretherton solution in the limit of small strain – and an inertia gravity wave field, the dynamics of which are examined in detail. Our analysis provides a unifying framework capable of describing frontal formation and geostrophic adjustment in a wide variety of settings.

The long-time viscous flow about two identical rotating circular cylinders in a side-by-side arrangement is investigated using an adaptive numerical scheme based on the vortex method. The Stokes solution of the steady flow about the two-cylinder cluster produces a uniform stream in the far field, which is the so-called Jeffery’s paradox. The present work first addresses the validation of the vortex method for a low-Reynolds-number computation. The unsteady flow past an abruptly started purely rotating circular cylinder is therefore computed and compared with an exact solution to the Navier–Stokes equations. The steady state is then found to be obtained for $t\gg 1$ with ${\mathit{Re}}_{\omega } {r}^{2} \ll t$, where the characteristic length and velocity are respectively normalized with the radius ${a}_{1} $ of the circular cylinder and the circumferential velocity ${\Omega }_{1} {a}_{1} $. Then, the influence of the Reynolds number ${\mathit{Re}}_{\omega } = { a}_{1}^{2} {\Omega }_{1} / \nu $ about the two-cylinder cluster is investigated in the range $0. 125\leqslant {\mathit{Re}}_{\omega } \leqslant 40$. The convection influence forms a pair of circulations (called self-induced closed streamlines) ahead of the cylinders to alter the symmetry of the streamline whereas the low-Reynolds-number computation (${\mathit{Re}}_{\omega } = 0. 125$) reaches the steady regime in a proper inner domain. The self-induced closed streamline is formed at far field due to the boundary condition being zero at infinity. When the two-cylinder cluster is immersed in a uniform flow, which is equivalent to Jeffery’s solution, the streamline behaves like excellent Jeffery’s flow at ${\mathit{Re}}_{\omega } = 1. 25$ (although the drag force is almost zero). On the other hand, the influence of the gap spacing between the cylinders is also investigated and it is shown that there are two kinds of flow regimes including Jeffery’s flow. At a proper distance from the cylinders, the self-induced far-field velocity, which is almost equivalent to Jeffery’s solution, is successfully observed in a two-cylinder arrangement.

A class of unsteady free convection flows over a differentially cooled horizontal surface is considered. The cooling, specified in terms of an imposed negative buoyancy or buoyancy flux, varies laterally as a step function with a single step change. As thermal boundary layers develop on either side of the step change, an intrinsically unsteady, boundary-layer-like flow arises in the transition zone between them. Self-similarity model solutions of the Boussinesq equations of motion, thermal energy, and mass conservation, within a boundary-layer approximation, are obtained for flows of unstratified fluids driven by a surface buoyancy or buoyancy flux, and flows of stably stratified fluids driven by a surface buoyancy flux. The motion is characterized by a shallow, primarily horizontal flow capped by a weak return flow. Stratification weakens the primary flow and strengthens the return flow. The flows intensify as the step change in surface forcing increases or as the Prandtl number decreases. Simple formulas are obtained for the propagation speeds, trajectories and the evolution of velocity maxima and other local extrema. Similarity-model predictions are verified through numerical simulations in which no boundary-layer approximations are made.

We investigate the Rayleigh–Bénard convection problem within the context of a diffusion flame formed in a horizontal channel where the fuel and oxidizer concentrations are prescribed at the porous walls. This problem seems to have received no attention in the literature. When formulated in the low-Mach-number approximation the model depends on two main non-dimensional parameters, the Rayleigh number and the Damköhler number, which govern gravitational strength and reaction speed respectively. In the steady state the system admits a planar diffusion flame solution; the aim is to find the critical Rayleigh number at which this solution becomes unstable to infinitesimal perturbations. In the Boussinesq approximation, a linear stability analysis reduces the system to a matrix equation with a solution comparable to that of the well-studied non-reactive case of Rayleigh–Bénard convection with a hot lower boundary. The planar Burke–Schumann diffusion flame, which has been previously considered unconditionally stable in studies disregarding gravity, is shown to become unstable when the Rayleigh number exceeds a critical value. A numerical treatment is performed to test the effects of compressibility and finite chemistry on the stability of the system. For weak values of the thermal expansion coefficient $\alpha $, the numerical results show strong agreement with those of the linear stability analysis. It is found that as $\alpha $ increases to a more realistic value the system becomes considerably more stable, and also exhibits hysteresis at the onset of instability. Finally, a reduction in the Damköhler number is found to decrease the stability of the system.

Transcritical flow of a stratified fluid past a broad localised topographic obstacle is studied analytically in the framework of the forced extended Korteweg–de Vries, or Gardner, equation. We consider both possible signs for the cubic nonlinear term in the Gardner equation corresponding to different fluid density stratification profiles. We identify the range of the input parameters: the oncoming flow speed (the Froude number) and the topographic amplitude, for which the obstacle supports a stationary localised hydraulic transition from the subcritical flow upstream to the supercritical flow downstream. Such a localised transcritical flow is resolved back into the equilibrium flow state away from the obstacle with the aid of unsteady coherent nonlinear wave structures propagating upstream and downstream. Along with the regular, cnoidal undular bores occurring in the analogous problem for the single-layer flow modelled by the forced Korteweg–de Vries equation, the transcritical internal wave flows support a diverse family of upstream and downstream wave structures, including kinks, rarefaction waves, classical undular bores, reversed and trigonometric undular bores, which we describe using the recent development of the nonlinear modulation theory for the (unforced) Gardner equation. The predictions of the developed analytic construction are confirmed by direct numerical simulations of the forced Gardner equation for a broad range of input parameters.