The two-dimensional free-surface problem of an ideal jet impinging on an uneven wall is studied using complex-variable and transform techniques. A relation between the flow angle on the free surface and the wall angle is first obtained. Then, by using a Hilbert transform and the generalized Schwarz–Christoffel transformation technique, a system of nonlinear integro-differential equations for the flow angle and the wall angle is formulated. For the case of symmetric flow, a compatibility condition for the system is automatically satisfied. In some special cases, for instance when the wall is a wedge, the problem reduces to the evaluation of several integrals. Moreover, in the case of a jet impinging normally on a flat wall, the classical result is recovered. For the asymmetric case, a relation is obtained between the point in the reference ζ-plane which corresponds to the position of the stagnation point in the physical plane, the flow speed and the shape of the wall. The solution to a linearized problem is given, for comparison. Some numerical solutions are presented, showing the shape of the free surface corresponding to a number of different wall shapes.

At the body-diameter Reynolds number 2000, the axisymmetric wake behind a slender streamlined body of revolution remained laminar without any indication of breakup far beyond the five-body-lengths test section. The mean-velocity profiles were measured with close spacings in order to obtain experimental information on the transient process from boundary layer to fully developed wake. A semi-empirical theory was developed to describe the flow-field characteristics.

Numerical simulations and laboratory experiments were jointly conducted to investigate a bathtub vortex under the influence of a protruding cylinder in a rotating tank. In the set-up, a central drain hole is placed at the bottom of the tank and a top-down cylinder is suspended from the rigid top lid, with fluid supplied from the sidewall for mass conservation. The cylinder is protruded to produce the Taylor column effect. The flow pattern depends on the Rossby number ($\mathit{Ro}= U/ fR$), the Ekman number ($\mathit{Ek}= \nu / f{R}^{2} )$ and the height ratio, $h/ H$, where $R$ is the radius of the cylinder, $f$ is the Coriolis parameter, $\nu $ is the kinematic viscosity of the fluid, $h$ is the vertical length of the cylinder and $H$ is the height of the tank. It is found appropriate to choose $U$ to be the average inflow velocity of fluid entering the column beneath the cylinder. Steady-state solutions obtained by numerically solving the Navier–Stokes equations in the rotating frame are shown to have a good agreement with flow visualizations and particle tracking velocimetry (PTV) measurements. It is known that at $\mathit{Ro}\sim 1{0}^{- 2} $, the central downward flow surrounded by the neighbouring Ekman pumping forms a classic one-celled bathtub vortex structure when there is no protruding cylinder ($h/ H= 0$). The influence of a suspended cylinder ($h/ H\not = 0$) leads to several findings. The bathtub vortex exhibits an interesting two-celled structure with an inner Ekman pumping (EP) and an outer up-drafting motion, termed Taylor upwelling (TU). The two regions of up-drafting motion are separated by a notable finite-thickness structure, identified as a (thin-walled) Taylor column. The thickness ${ \delta }_{T}^{\ast } $ of the Taylor column is found to be well correlated to the height ratio and the Ekman number by ${\delta }_{T} = { \delta }_{T}^{\ast } / R= {(1- h/ H)}^{- 0. 32} {\mathit{Ek}}^{0. 095} $. The Taylor column presents a barrier to the fluid flow such that the fluid from the inlet may only flow into the inner region through the narrow gaps, one above the Taylor column and one beneath it (conveniently called Ekman gaps). As a result, five types of routes along which the fluid may flow to and exit at the drain hole could be identified for the multi-celled vortex structure. Moreover, the flow rates associated with the five routes were calculated and compared to help understand the relative importance of the component flow structures. The weaker influence of the Taylor column effect on the bathtub vortex at $\mathit{Ro}\sim 1$ or even higher $\mathit{Ro}\sim 1{0}^{2} $ is also discussed.

A core–annular flow, the concurrent axial flow of two immiscible fluids in a circular tube or pore with one fluid in the core and the other in the wetting annular region, is frequently used to model technologically important flows, e.g. in liquid–liquid displacements in secondary oil recovery. Most of the existing literature assumes that the pores in which such flows occur are uniform circular cylinders, and examine the interfacial stability of such systems as a function of fluid and interfacial properties. Since real rock pores possess a more complex geometry, the companion paper examined the linear stability of core–annular flows in axisymmetric, corrugated pores in the limit of asymptotically weak corrugation. It found that short-wave disturbances that were stable in straight tubes could couple to the wall's periodicity to excite unstable long waves. In this paper, we follow the evolution of the axisymmetric, linearly unstable waves for fluids of equal densities in a corrugated tube into the weakly nonlinear regime. Here, we ask whether this continual generation of new disturbances by the coupling to the wall's periodicity can overcome the nonlinear saturation mechanism that relies on the nonlinear (kinematic-condition-derived) wave steepening of the Kuramoto–Sivashinsky (KS) equation. If it cannot, and the unstable waves still saturate, then do these additional excited waves make the KS solutions more likely to be chaotic, or does the dispersion introduced into the growth rate correction by capillarity serve to regularize otherwise chaotic motions?

We find that in the usual strong surface tension limit, the saturation mechanism of the KS mechanism remains able to saturate all disturbances. Moreover, an additional capillary-derived nonlinear term seems to favour regular travelling waves over chaos, and corrugation adds a temporal periodicity to the waves associated with their periodical traversing of the wall's crests and troughs. For even larger surface tensions, capillarity dominates over convection and a weakly nonlinear version of Hammond's no-flow equation results; this equation, with or without corrugation, suggests further growth. Finally, for a weaker surface tension, the leading-order base flow interface follows the wall's shape. The corrugation-derived excited waves appear able to push an otherwise regular travelling wave solution to KS to become chaotic, whereas its dispersive properties in this limit seem insufficiently strong to regularize chaotic motions.

Interfacial instability arising from evaporation of a single component liquid is investigated using linear and weakly nonlinear analysis. Evaporative convection is studied taking into account the fluid dynamics of both liquid and vapour phases as well as lateral rigid sidewalls. Both open and closed systems are addressed. The nature of the bifurcation and the change in heat flux in the nonlinear regime are determined. It is shown that depending upon the aspect ratio of the geometry, either supercritical or subcritical branching behaviour is possible.

Mild-slope (MS) type equations are depth-integrated models, which predict under appropriate conditions refraction and diffraction of linear time-harmonic water waves. Among these equations, the complementary mild-slope equation (CMSE) was shown to give better agreement with exact two-dimensional linear theory compared to other MS-type equations. Nevertheless, it has a disadvantage of being a vector equation, i.e. it requires solving a system of two coupled partial differential equations. In addition, for three-dimensional problems, there is a difficulty in constructing the additional boundary condition needed for the solution. In the present work, it is shown how the vector CMSE can be transformed into an equivalent scalar equation using a pseudo-potential formulation. The pseudo-potential mild-slope equation (PMSE) preserves the accuracy of the CMSE while avoiding the need of an additional boundary condition. Furthermore, the PMSE significantly reduces the computational effort relative to the CMSE, since it is a scalar equation. The accuracy of the new model was tested numerically by comparing it to laboratory data and analytical solutions.

A solution for the weak bore is found in which the mean profile is dominated by viscosity, so that the velocity variation is given essentially by a quasi-uniform Poiseuille flow. It is found that such a transition between flows of different depths is possible provided the Froude number is less than 1·58. The possibility of superposing an inviscid perturbation on such a flow is then investigated. Under favourable circumstances the effect of this perturbation is to add to the profile of the free surface a term which decays exponentially in front of the bore, but is oscillatory behind it.

Moderately favourable pressure gradient turbulent boundary layers are investigated within a theoretical framework based on the unintegrated two-dimensional mean momentum equation. The present theory stems from an observed exchange of balance between terms in the mean momentum equation across different regions of the boundary layer. This exchange of balance leads to the identification of distinct physical layers, unambiguously defined by the predominant mean dynamics active in each layer. Scaling domains congruent with the physical layers are obtained from a multi-scale analysis of the mean momentum equation. Scaling behaviours predicted by the present theory are evaluated using direct measurements of all of the terms in the mean momentum balance for the case of a sink-flow pressure gradient generated in a wind tunnel with a long development length. Measurements also captured the evolution of the turbulent boundary layers from a non-equilibrium state near the wind tunnel entrance towards an equilibrium state further downstream. Salient features of the present multi-scale theory were reproduced in all the experimental data. Under equilibrium conditions, a universal function was found to describe the decay of the Reynolds stress profile in the outer region of the boundary layer. Non-equilibrium effects appeared to be manifest primarily in the outer region, whereas differences in the inner region were attributed solely to Reynolds number effects.

Just after formation, optical fibres are wetted stably with acrylate at capillary numbers routinely exceeding 1000. It is hypothesized that this is possible because of dissolution of air into the liquid coating. A lubrication/boundary integral analysis that includes gas diffusion and solubility is developed. It is applied using conservatively estimated solubility and diffusivity coefficients and solutions are found that are consistent with industry practice and with the hypothesis. The results also agree with the claim of Deneka, Kar & Mensah (1988) that the use of high-solubility gases to bathe a wetting line allows significantly greater wetting speeds. The solutions indicate a maximum speed of wetting which increases with gas solubility and with reduction in wetting-channel diameter.

We present an experimental and theoretical study of the dynamics of laterally confined marine ice sheets in the natural limit in which the long, narrow channel into which they flow is wider than the depth of the ice. A marine ice sheet comprises a grounded ice sheet in contact with bedrock that floats away from the bedrock at a ‘grounding line’ to form a floating ice shelf. We model the grounded ice sheet as a viscous gravity current resisted dominantly by vertical shear stresses owing to the no-slip boundary condition applied at the bedrock. We model the ice shelf as a floating viscous current resisted dominantly by horizontal shear stresses owing to no-slip boundary conditions applied at the sidewalls of the channel. The two shear-dominated regions are coupled by jump conditions relating force and fluid flux across a short transition region downstream of the grounding line. We find that the influence of the stresses within the transition region becomes negligible at long times and we model the transition region as a singular interface across which the ice thickness and mass flux can be discontinuous. The confined shelf buttresses the sheet, causing the grounding line to advance more than it would otherwise. In the case that the sheet flows on a base of uniform slope, we find asymptotically that the grounding line advances indefinitely as $t^{1/3}$, where $t$ is time. This contrasts with the two-dimensional counterpart, for which the shelf provides no buttressing and the grounding line reaches a steady state (Robison, J. Fluid Mech., vol. 648, 2010, pp. 363–380).

The contact line of a volatile liquid on a flat substrate is studied theoretically. We show that a remarkable result obtained for a pure-vapour atmosphere (Phys. Rev. E, vol. 87, 2013, 010401) also holds for an isothermal diffusion-limited vapour exchange with air. Namely, for both zero and finite Young’s angles, the motion- and phase-change-related contact-line singularities can in principle be regularised solely by the Kelvin effect (curvature dependence of saturation conditions). The latter prevents the curvature from diverging and rather leads to its versatile self-adjustment. To illustrate the point, the problem is resolved for a distinguished vicinity of the contact line (‘microregion’) in a ‘minimalist’ way, i.e. without any disjoining pressure, precursor film, Navier slip or any other microphysics. This also leads to the determination of the ‘Kelvin-only’ evaporation- and motion-induced apparent contact angles. With the Kelvin-only microscales actually turning out to be quite nanoscopic, other microphysics effects may nonetheless interfere too in reality. The Kelvin-only results will then yield a limiting case within such a more general formulation.

The novel viscous fingering instability recently found in the experiments of Kowal & Worster (J. Fluid Mech., vol. 766, 2015, pp. 626–655), involving two superposed currents of viscous fluid, has been shown to originate at the lubrication front when the fluids are of equal density. However, when the densities are unequal, additional buoyancy forces associated with the underlying layer act to suppress this instability and are largest at the lubrication front, which is where the instability originates. In this paper, we investigate the interaction between the mechanism of the instability and the stabilising influence of these buoyancy forces by performing a global and fully time-dependent analysis, which does not use the frozen-time approximation. We determine a critical condition for instability in terms of the viscosity ratio and the density difference between the two layers. Consistently with the local analysis of the companion paper, instabilities occur when the jump in hydrostatic pressure gradient across the lubrication front is negative, or, equivalently, when the intruding fluid is less viscous than the overlying fluid, provided the two fluids are of equal densities. Once there is a non-zero density difference, these driving buoyancy forces suppress the instability for large wavelengths, giving rise to wavelength selection. As the density difference increases, the instability criterion requires higher viscosity ratios for any instability to occur, and the band of unstable wavenumbers becomes bounded. Large enough density differences suppress the instability completely.

Over two decades ago, some numerical studies and laboratory experiments identified the phenomenon of leapfrogging internal solitary waves located on separated pycnoclines. We revisit this problem to explore the behaviour of the near resonance phenomenon. We have developed a numerical code to follow the long-time inviscid evolution of isolated mode-two disturbances on two separated pycnoclines in a three-layer stratified fluid bounded by rigid horizontal top and bottom walls. We study the dependence of the solution on input system parameters, namely the three fluid densities and the two interface thicknesses, for fixed initial conditions describing isolated mode-two disturbances on each pycnocline. For most parameter values, the initial disturbances separate immediately and evolve into solitary waves, each with a distinct speed. However, in a narrow region of parameter space, the waves pair up and oscillate for some time in leapfrog fashion with a nearly equal average speed. The motion is only quasi-periodic, as each wave loses energy into its respective dispersive tail, which causes the spatial oscillation magnitude and period to increase until the waves eventually separate. We record the separation time, oscillation period and magnitude, and the final amplitudes and celerity of the separated waves as a function of the input parameters, and give evidence that no perfect periodic solutions occur. A simple asymptotic model is developed to aid in interpretation of the numerical results.

Small perturbation evolution in compressible Poiseuille flow is contrasted against the incompressible case using direct simulations and non-modal linear analysis. The onset of compressibility effects leads to a profound change in the behaviour of pressure and its interaction with the velocity field. Linear analysis shows that the most significant compressibility outcome is the harmonic coupling between pressure and wall-normal velocity perturbations. Oscillations in normal perturbations can lead to periods of negative production causing suppression of perturbation growth. The extent of the influence of compressibility can be characterized in terms of an effective gradient Mach number ($M_{g}^{e}$). Analysis shows that $M_{g}^{e}$ diminishes as the angle of the perturbation increases with respect to the shear plane. Direct numerical simulations show that streamwise perturbations, which would lead to Tollmien–Schlichting instability in the incompressible case, are completely suppressed in the compressible case and experience the highest $M_{g}^{e}$. At the other extreme, computations reveal that spanwise perturbations, which experience negligible $M_{g}^{e}$, are entirely unaltered from the incompressible case. Perturbation behaviour at intermediate obliqueness angles is established. Moreover, the underlying pressure–velocity interactions are explicated.

The nonlinear interactions between a free surface and a shed vortex shear layer in the inviscid wake of a vertical surface-piercing plate are studied numerically using a mixed-Eulerian-Lagrangian method. For a plate with initial submergence d starting abruptly from rest to constant horizontal velocity U, the problem is governed by a single parameter, the Froude number Fn = U/(gd)½, where g is the gravitational acceleration. Depending on Fn, three classes of interaction dynamics (subcritical, transcritical and supercritical) are identified. For subcritical Fn ([lsim ] 0.7), the free surface plunges on both the forward and lee sides of the plate before significant interactions with the vortex sheet occur. For transcritical and supercritical Fn, interactions between the free surface and the starting vortex result in a stretching of the vortex sheet which eventually rolls up into double-branched spirals as a result of Kelvin-Helmholtz instability. In the transcritical range (Fn ∼ 0.7–1.0), the effect of the free surface on the double-branched spirals remains weak, while for supercritical Fn ([gsim ] 1.0), strong interactions lead to entrainment of the double-branched spiral into the free surface resulting in prominent surface features.

We present wall-resolved large-eddy simulation (LES) of flow with free-stream velocity $\boldsymbol{U}_{\infty }$ over a cylinder of diameter $D$ rotating at constant angular velocity $\unicode[STIX]{x1D6FA}$, with the focus on the lift crisis, which takes place at relatively high Reynolds number $Re_{D}=U_{\infty }D/\unicode[STIX]{x1D708}$, where $\unicode[STIX]{x1D708}$ is the kinematic viscosity of the fluid. Two sets of LES are performed within the ($Re_{D}$, $\unicode[STIX]{x1D6FC}$)-plane with $\unicode[STIX]{x1D6FC}=\unicode[STIX]{x1D6FA}D/(2U_{\infty })$ the dimensionless cylinder rotation speed. One set, at $Re_{D}=5000$, is used as a reference flow and does not exhibit a lift crisis. Our main LES varies $\unicode[STIX]{x1D6FC}$ in $0\leqslant \unicode[STIX]{x1D6FC}\leqslant 2.0$ at fixed $Re_{D}=6\times 10^{4}$. For $\unicode[STIX]{x1D6FC}$ in the range $\unicode[STIX]{x1D6FC}=0.48{-}0.6$ we find a lift crisis. This range is in agreement with experiment although the LES shows a deeper local minimum in the lift coefficient than the measured value. Diagnostics that include instantaneous surface portraits of the surface skin-friction vector field $\boldsymbol{C}_{\boldsymbol{f}}$, spanwise-averaged flow-streamline plots, and a statistical analysis of local, near-surface flow reversal show that, on the leeward-bottom cylinder surface, the flow experiences large-scale reorganization as $\unicode[STIX]{x1D6FC}$ increases through the lift crisis. At $\unicode[STIX]{x1D6FC}=0.48$ the primary-flow features comprise a shear layer separating from that side of the cylinder that moves with the free stream and a pattern of oscillatory but largely attached flow zones surrounded by scattered patches of local flow separation/reattachment on the lee and underside of the cylinder surface. Large-scale, unsteady vortex shedding is observed. At $\unicode[STIX]{x1D6FC}=0.6$ the flow has transitioned to a more ordered state where the small-scale separation/reattachment cells concentrate into a relatively narrow zone with largely attached flow elsewhere. This induces a low-pressure region which produces a sudden decrease in lift and hence the lift crisis. Through this process, the boundary layer does not show classical turbulence behaviour. As $\unicode[STIX]{x1D6FC}$ is further increased at constant $Re_{D}$, the localized separation zone dissipates with corresponding attached flow on most of the cylinder surface. The lift coefficient then resumes its increasing trend. A logarithmic region is found within the boundary layer at $\unicode[STIX]{x1D6FC}=1.0$.

Surface stimulation of any physical origin (electrohydrodynamic, thermocapillary, etc.) has the goal of generating localized perturbations on the free surface or the velocity field of a capillary jet. Among these perturbations, only the axisymmetric ones are determinant for the jet breakup. Often, the stimulation is weak enough for a linear model to be applicable. Then, the stimulation can be described by means of the Green functions for stresses, both normal and tangential to the interface, the calculations of which are, in addition, uncoupled from the hydrodynamic variables. If a harmonic forcing is applied, these Green functions are combinations of the spatial modes whose associated poles lie inside the appropriate integration contour of the complex wavenumber plane. This is the motivation for a comprehensive enumeration and description of the spatial modes, which has not been done up to now. Modes familiar from a temporal analysis, the dominant and subdominant capillary modes and the hydrodynamic modes, are present, along with modes specific to a spatial analysis. Most of the latter have already been mentioned in the literature for inviscid jets, but not analysed. A mode not previously found is reported. In addition, a description of the velocity field associated with each mode is provided, as a tool to understand their physical origin and behaviour. The relative importance of each mode in both normal- and tangential-stress stimulations is discussed. Finally, the well-known merging of poles below a critical jet velocity, leading to absolute instability, is analysed in the light of the modal description.

For a small sphere suspended in a background fluid flow near an obstacle, we calculate the hydrodynamic force on the sphere in the direction normal to the boundary of the obstacle. Using the Lorentz reciprocal theorem, we obtain analytical expressions for the normal force in the Stokes flow limit, valid for arbitrary separations of the particle from the obstacle, both for solid obstacles and those with free surfaces. The main effect of the boundary is to produce a normal force proportional to extensional flow gradients in the vicinity of the particle. The strength of this force is greatest when the separation between the surfaces of the particle and the obstacle is small relative to the particle size. While the magnitude of the force weakens for large separations between the sphere and the obstacle (decaying quadratically with separation distance), it can significantly modify Faxén’s law even at modestly large separation distances. In addition, we find a second force contribution due to the curvature of the background flow normal to the obstacle, which is also important when the sphere is close to the obstacle. The results of the theory are of importance to the dynamics of particles in confined geometries, whether bounded by a solid obstacle, the wall of a channel or a gas bubble.