We consider the linearized form of the regularized 13-moment equations (R13) to model the slow, steady gas dynamics surrounding a rigid, heat-conducting sphere when a uniform temperature gradient is imposed far from the sphere and the gas is in a state of rarefaction. Under these conditions, the phenomenon of thermophoresis, characterized by forces on the solid surfaces, occurs. The R13 equations, derived from the Boltzmann equation using the moment method, provide closure to the mass, momentum and energy conservation laws in the form of constitutive, transport equations for the stress and heat flux that extend the Navier–Stokes–Fourier model to include non-equilibrium effects. We obtain analytical solutions for the field variables that characterize the gas dynamics and a closed-form expression for the thermophoretic force on the sphere. We also consider the slow, streaming flow of gas past a sphere using the same model resulting in a drag force on the body. The thermophoretic velocity of the sphere is then determined from the balance between thermophoretic force and drag. The thermophoretic force is compared with predictions from other theories, including Grad’s 13-moment equations (G13), variants of the Boltzmann equation commonly used in kinetic theory, and with recently published experimental data. The new results from R13 agree well with results from kinetic theory up to a Knudsen number (based on the sphere’s radius) of approximately 0.1 for the values of solid-to-gas heat conductivity ratios considered. However, in this range of Knudsen numbers, where for a very high thermal conductivity of the solid the experiments show reversed thermophoretic forces, the R13 solution, which does result in a reversal of the force, as well as the other theories predict significantly smaller forces than the experimental values. For Knudsen numbers between 0.1 and 1 approximately, the R13 model of thermophoretic force qualitatively shows the trend exhibited by the measurements and, among the various models considered, results in the least discrepancy.

In isothermal non-coalescence behaviours of a droplet against a wall, an air film of micrometre thickness plays a crucial role. We experimentally study this phenomenon by letting a droplet levitate over a moving glass wall. The three-dimensional shape of the air film is measured using an interferometric method. The mean curvature distribution of the deformed free surface and the distributions of the lubrication pressure are derived from the experimental measurements. We vary experimental parameters, namely wall velocity, droplet diameter and viscosity of the droplets, over a wide range; for example, the droplet viscosity is varied over two orders of magnitude. For the same wall velocity, the air film of low-viscosity droplets shows little shape oscillation with constant film thickness (defined as the steady state), while that of highly viscous droplets shows a significant shape oscillation with varying film thickness (defined as the unsteady state). The droplet viscosity also affects the surface velocity of a droplet. Under our experimental conditions, where the air film shape can be assumed to be steady, we present experimental evidence showing that the lift force generated inside the air film balances with the droplet’s weight. We also verify that the lubrication pressure locally balances with the surface tension and hydrostatic pressures. This indicates that lubrication pressure and the shape of the free surface are mutually determined. Based on the local pressure balance, we discuss a process of determining the steady shape of an air film that has two areas of minimum thickness in the vicinity of the downstream rim.

The theoretical framework developed by Rayleigh and Plateau in the 19th century has been remarkably accurate in describing macroscale experiments of liquid cylinder instability. Here we re-evaluate and revise the Rayleigh–Plateau instability for the nanoscale, where molecular dynamics experiments demonstrate its inadequacy. A new framework based on the stochastic lubrication equation is developed that captures nanoscale flow features and highlights the critical role of thermal fluctuations at small scales. Remarkably, the model indicates that classically stable (i.e. ‘fat’) liquid cylinders can be broken at the nanoscale, and this is confirmed by molecular dynamics.

The steady behaviour of a rarefied gas around a rotating sphere is studied numerically on the basis of the linearised ellipsoidal statistical model of the Boltzmann equation, also known as the ES model, and the Maxwell diffuse–specular boundary condition. It is demonstrated numerically that the normal derivative of the circumferential component of the flow velocity and that of the heat flux diverge on the boundary with a rate $s^{-1/2}$, where $s$ is the normal distance from the boundary. Further, it is demonstrated that the diverging term is proportional to the magnitude of the jump discontinuity of the velocity distribution function on the boundary, which originates from the mismatch of the incoming and outgoing data on the boundary. The moment of force exerted on the sphere is also obtained for a wide range of the Knudsen number and for various values of the accommodation coefficient.

Lagrangian data assimilation is a complex problem in oceanic and atmospheric modelling. Tracking drifters in large-scale geophysical flows can involve uncertainty in drifter location, complex inertial effects and other factors which make comparing them to simulated Lagrangian trajectories from numerical models extremely challenging. Temporal and spatial discretisation, factors necessary in modelling large scale flows, also contribute to separation between real and simulated drifter trajectories. The chaotic advection inherent in these turbulent flows tends to separate even closely spaced tracer particles, making error metrics based solely on drifter displacements unsuitable for estimating model parameters. We propose to instead use error in the coherent structure colouring (CSC) field to assess model skill. The CSC field provides a spatial representation of the underlying coherent patterns in the flow, and we show that it is a more robust metric for assessing model accuracy. Through the use of two test cases, one considering spatial uncertainty in particle initialisation, and one examining the influence of stochastic error along a trajectory and temporal discretisation, we show that error in the coherent structure colouring field can be used to accurately determine single or multiple simultaneously unknown model parameters, whereas a conventional error metric based on error in drifter displacement fails. Because the CSC field enhances the difference in error between correct and incorrect model parameters, error minima in model parameter sweeps become more distinct. The effectiveness and robustness of this method for single and multi-parameter estimation in analytical flows suggest that Lagrangian data assimilation for real oceanic and atmospheric models would benefit from a similar approach.

Freak or rogue waves are so called because of their unexpectedly large size relative to the population of smaller waves in which they occur. The 25.6 m high Draupner wave, observed in a sea state with a significant wave height of 12 m, was one of the first confirmed field measurements of a freak wave. The physical mechanisms that give rise to freak waves such as the Draupner wave are still contentious. Through physical experiments carried out in a circular wave tank, we attempt to recreate the freak wave measured at the Draupner platform and gain an understanding of the directional conditions capable of supporting such a large and steep wave. Herein, we recreate the full scaled crest amplitude and profile of the Draupner wave, including bound set-up. We find that the onset and type of wave breaking play a significant role and differ significantly for crossing and non-crossing waves. Crucially, breaking becomes less crest-amplitude limiting for sufficiently large crossing angles and involves the formation of near-vertical jets. In our experiments, we were only able to reproduce the scaled crest and total wave height of the wave measured at the Draupner platform for conditions where two wave systems cross at a large angle.

This paper presents a comprehensive investigation into flow past a circular cylinder where compressibility and rarefaction effects play an important role. The study focuses on steady subsonic flow in the Reynolds-number range 0.1–45. Rarefaction, or non-equilibrium, effects in the slip and early transition regime are accounted for using the method of moments and results are compared to data from kinetic theory obtained from the direct simulation Monte Carlo method. Solutions obtained for incompressible continuum flow serve as a baseline to examine non-equilibrium effects on the flow features. For creeping flow, where the Reynolds number is less than unity, the drag coefficient predicted by the moment equations is in good agreement with kinetic theory for Knudsen numbers less than one. When flow separation occurs, we show that the effects of rarefaction and velocity slip delay flow separation and will reduce the size of the vortices downstream of the cylinder. When the Knudsen number is above 0.028, the vortex length shows an initial increase with the Reynolds number, as observed in the standard no-slip continuum regime. However, once the Reynolds number exceeds a critical value, the size of the downstream vortices decreases with increasing Reynolds number until they disappear. An existence criterion, which identifies the limits for the presence of the vortices, is proposed. The flow physics around the cylinder is further analysed in terms of velocity slip, pressure and skin friction coefficients, which highlights that viscous, rarefaction and compressibility effects all play a complex role. We also show that the local Knudsen number, which indicates the state of the gas around the cylinder, can differ significantly from its free-stream value and it is essential that computational studies of subsonic gas flows in the slip and early transition regime are able to account for these strong non-equilibrium effects.

The fate of a contracting liquid filament depends on the Ohnesorge number ($Oh$), the initial aspect ratio ($\unicode[STIX]{x1D6E4}$) and surface perturbation. Generally, it is believed that there exists a critical aspect ratio $\unicode[STIX]{x1D6E4}_{c}(Oh)$ such that longer filaments break up and shorter ones recoil into a single drop. Through computational and experimental studies, we report a transitional regime for filaments with a broad range of intermediate aspect ratios, where there exist multiple $\unicode[STIX]{x1D6E4}_{c}$ thresholds at which a novel breakup mode alternates with no-break mode. We develop a simple model considering the superposition of capillary waves, which can predict the complicated new phase diagram. In this model, the breakup results from constructive interference between the capillary waves that originate from the ends of the filament.

We investigate the deformation of a linear viscoelastic compliant coating in a turbulent flow for a wide range of coating parameters. A one-way coupling model is proposed in which the turbulent surface stresses are expressed as a sum of streamwise-travelling waves with amplitudes determined from the stress spectra of the corresponding flow over a rigid wall. The analytically calculated coating deformation is analysed in terms of the root-mean-square (r.m.s.) surface displacement and the corresponding point frequency spectra. The present study systematically investigates the influence of five coating properties namely density, stiffness, thickness, viscoelasticity and compressibility. The surface displacements increase linearly with the fluid/solid density ratio. They are linearly proportional to the coating thickness for thin coatings, while they become independent of the thickness for thick coatings. Very soft coatings show resonant behaviour, but the displacement for stiffer coatings is proportional to the inverse of the shear modulus. The viscoelastic loss angle has only a significant influence when resonances occur in the coating response, while Poisson’s ratio has a minor effect for most cases. The modelled surface displacement is qualitatively compared with recent measurements on the deformation of three different coatings in a turbulent boundary-layer flow. The model predicts the order of magnitude of the surface displacement, and it captures the increase of the coating displacement with the Reynolds number and the coating softness. Finally, we propose a scaling that collapses all the experimental data for the r.m.s. of the vertical surface displacement onto a single curve.

The objective of this work is to investigate linear modal and algebraic instability in Poiseuille flows with fluids close to their vapour–liquid critical point. Close to this critical point, the ideal gas assumption does not hold and large non-ideal fluid behaviours occur. As a representative non-ideal fluid, we consider supercritical carbon dioxide ($\text{CO}_{2}$) at a pressure of 80 bar, which is above its critical pressure of 73.9 bar. The Poiseuille flow is characterized by the Reynolds number ($Re=\unicode[STIX]{x1D70C}_{w}^{\ast }u_{r}^{\ast }h^{\ast }/\unicode[STIX]{x1D707}_{w}^{\ast }$), the product of the Prandtl ($Pr=\unicode[STIX]{x1D707}_{w}^{\ast }C_{pw}^{\ast }/\unicode[STIX]{x1D705}_{w}^{\ast }$) and Eckert numbers ($Ec=u_{r}^{\ast 2}/C_{pw}^{\ast }T_{w}^{\ast }$) and the wall temperature that in addition to pressure determine the thermodynamic reference condition. For low Eckert numbers, the flow is essentially isothermal and no difference with the well-known stability behaviour of incompressible flows is observed. However, if the Eckert number increases, the viscous heating causes gradients of thermodynamic and transport properties, and non-ideal gas effects become significant. Three regimes of the laminar base flow can be considered: the subcritical (temperature in the channel is entirely below its pseudo-critical value), transcritical and supercritical temperature regimes. If compared to the linear stability of an ideal gas Poiseuille flow, we show that the base flow is modally more unstable in the subcritical regime, inviscid unstable in the transcritical regime and significantly more stable in the supercritical regime. Following the principle of corresponding states, we expect that qualitatively similar results will be obtained for other fluids at equivalent thermodynamic states.

We study air entrainment by a solid plate plunging into a viscous liquid, theoretically and numerically. At dimensionless speeds $Ca=U\unicode[STIX]{x1D702}/\unicode[STIX]{x1D6FE}$ of order unity, a near-cusp forms due to the presence of a moving contact line. The radius of curvature of the cusp’s tip scales with the slip length multiplied by an exponential of $-Ca$. The pressure from the air flow drawn inside the cusp leads to a bifurcation, at which air is entrained, i.e. there is ‘wetting failure’. We develop an analytical theory of the threshold to air entrainment, which predicts the critical capillary number to depend logarithmically on the viscosity ratio, with corrections coming from the slip in the gas phase.

Direct numerical simulations are used to characterize wind-shear effects on entrainment in a barotropic convective boundary layer (CBL) that grows into a linearly stratified atmosphere. We consider weakly to strongly unstable conditions $-z_{enc}/L_{Ob}\gtrsim 4$, where $z_{enc}$ is the encroachment CBL depth and $L_{Ob}$ is the Obukhov length. Dimensional analysis allows us to characterize such a sheared CBL by a normalized CBL depth, a Froude number and a Reynolds number. The first two non-dimensional quantities embed the dependence of the system on time, on the surface buoyancy flux, and on the buoyancy stratification and wind velocity in the free atmosphere. We show that the dependence of entrainment-zone properties on these two non-dimensional quantities can be expressed in terms of just one independent variable, the ratio between a shear scale $(\unicode[STIX]{x0394}z_{i})_{s}\equiv \sqrt{1/3}\unicode[STIX]{x0394}u/N_{0}$ and a convective scale $(\unicode[STIX]{x0394}z_{i})_{c}\equiv 0.25z_{enc}$, where $\unicode[STIX]{x0394}u$ is the velocity increment across the entrainment zone, and $N_{0}$ is the buoyancy frequency of the free atmosphere. Here $(\unicode[STIX]{x0394}z_{i})_{s}$ and $(\unicode[STIX]{x0394}z_{i})_{c}$ represent the entrainment-zone thickness in the limits of weak convective instability (strong wind) and strong convective instability (weak wind), respectively. We derive scaling laws for the CBL depth, the entrainment-zone thickness, the mean entrainment velocity and the entrainment-flux ratio as functions of $(\unicode[STIX]{x0394}z_{i})_{s}/(\unicode[STIX]{x0394}z_{i})_{c}$. These scaling laws can also be expressed as functions of only a Richardson number $(N_{0}z_{enc}/\unicode[STIX]{x0394}u)^{2}$, but not in terms of only the stability parameter $-z_{enc}/L_{Ob}$.

Both in rheometry and in fundamental fluid mechanics studies, the Taylor–Couette geometry is used frequently to investigate viscoelastic fluids. In order to ensure a constant shear rate in the gap between the inner and outer cylinders, such studies are usually restricted to the small-gap limit where the assumption of a linear velocity distribution is well justified. In conjunction with a sufficiently large aspect ratio $\unicode[STIX]{x1D6EC}$ (i.e. ratio of cylinder height to gap), the flow is then assumed to be viscometric. Here we demonstrate, using a perturbation technique with the curvature ratio (i.e. ratio of the half-gap to the mid-radius of the cylinders) as the perturbation parameter, full nonlinear simulations using a finite-volume technique, and supporting experiments, that, even in the creeping-flow (inertialess) narrow-gap limit, for viscoelastic fluids end effects due to finite aspect ratio always give rise to a secondary motion. Using the constant-viscosity Oldroyd-B model we are able to show that this secondary motion, as has been observed in related pressure-driven flows with curvature, such as the viscoelastic Dean flow, is solely a consequence of the combination of gradients of the first normal-stress difference and curvature. Our results show that end effects can significantly change the flow characteristics, especially for small aspect ratios, and this may have important consequences in some situations such as the onset criteria for purely elastic instabilities.

The presence and the microscopic origin of normal stress differences in dense suspensions under simple shear flows are investigated by means of inertialess particle dynamics simulations, taking into account hydrodynamic lubrication and frictional contact forces. The synergic action of hydrodynamic and contact forces between the suspended particles is found to be the origin of negative contributions to the first normal stress difference $N_{1}$, whereas positive values of $N_{1}$ observed at higher volume fractions near jamming are due to effects that cannot be accounted for in the hard-sphere limit. Furthermore, we found that the stress anisotropy induced by the planarity of the simple shear flow vanishes as the volume fraction approaches the jamming point for frictionless particles, while it remains finite for the case of frictional particles.

Internal gravity wave energy contributes significantly to the energy budget of the oceans, affecting mixing and the thermohaline circulation. Hence it is important to determine the internal wave energy flux $\boldsymbol{J}=p\,\boldsymbol{v}$, where $p$ is the pressure perturbation field and $\boldsymbol{v}$ is the velocity perturbation field. However, the pressure perturbation field is not directly accessible in laboratory or field observations. Previously, a Green’s function based method was developed to calculate the instantaneous energy flux field from a measured density perturbation field $\unicode[STIX]{x1D70C}(x,z,t)$, given a constant buoyancy frequency $N$. Here we present methods for computing the instantaneous energy flux $\boldsymbol{J}(x,z,t)$ for an internal wave field with vertically varying background $N(z)$, as in the oceans where $N(z)$ typically decreases by two orders of magnitude from the pycnocline to the deep ocean. Analytic methods are presented for computing $\boldsymbol{J}(x,z,t)$ from a density perturbation field for $N(z)$ varying linearly with $z$ and for $N^{2}(z)$ varying as $\tanh (z)$. To generalize this approach to arbitrary $N(z)$, we present a computational method for obtaining $\boldsymbol{J}(x,z,t)$. The results for $\boldsymbol{J}(x,z,t)$ for the different cases agree well with results from direct numerical simulations of the Navier–Stokes equations. Our computational method can be applied to any density perturbation data using the MATLAB graphical user interface ‘EnergyFlux’.

Oceanic internal tides and other inertia–gravity waves propagate in an energetic turbulent flow whose length scales are similar to the wavelengths. Advection and refraction by this flow cause the scattering of the waves, redistributing their energy in wavevector space. As a result, initially plane waves radiated from a source such as a topographic ridge become spatially incoherent away from the source. To examine this process, we derive a kinetic equation which describes the statistics of the scattering under the assumptions that the flow is quasigeostrophic, barotropic and well represented by a stationary homogeneous random field. Energy transfers are quantified by computing a scattering cross-section and shown to be restricted to waves with the same frequency and identical vertical structure, hence the same horizontal wavelength. For isotropic flows, scattering leads to an isotropic wave field. We estimate the characteristic time and length scales of this isotropisation, and study their dependence on parameters including the energy spectrum of the flow. Simulations of internal tides generated by a planar wavemaker carried out for the linearised shallow-water model confirm the pertinence of these scales. A comparison with the numerical solution of the kinetic equation demonstrates the validity of the latter and illustrates how the interplay between wave scattering and transport shapes the wave statistics.

We study the shape and motion of gas bubbles in a liquid flowing through a horizontal or slightly inclined thin annulus. Experimental data show that in the horizontal annulus, bubbles develop a unique ‘tadpole-like’ shape with a semi-circular cap and a highly stretched tail. As the annulus is inclined, the bubble tail tends to vanish, resulting in a significant decrease of bubble length. To model the bubble evolution, the thin annulus is conceptualised as a ‘Hele-Shaw’ cell in a curvilinear space. The three-dimensional flow within the cell is represented by a gap-averaged, two-dimensional model, which achieved a close match to the experimental data. The numerical model is further used to investigate the effects of gap thickness and pipe diameter on the bubble behaviour. The mechanism for the semi-circular cap formation is interpreted based on an analogous irrotational flow field around a circular cylinder, based on which a theoretical solution to the bubble velocity is derived. The bubble motion and cap geometry is mainly controlled by the gravitational component perpendicular to the flow direction. The bubble elongation in the horizontal annulus is caused by the buoyancy that moves the bubble to the top of the annulus. However, as the annulus is inclined, the gravitational component parallel to the flow direction becomes important, causing bubble separation at the tail and reduction in bubble length.

Very low Reynolds number propulsion is a topic of enduring interest due to its importance in biological systems such as sperm migration in the female reproductive tract. Motivated by the fibrous nature of cervical mucus, several recent studies have considered the effect of anisotropic rheology; these studies have generally employed the classical swimming sheet model of G. I. Taylor. The models of Cupples et al. (J. Fluid Mech. vol. 812, 2017, pp. 501–524) and Shi & Powers (Phys. Rev. Fluids vol. 2, 2017, 123102) consider related problems which in a common limit (passive, slightly anisotropic) make different predictions regarding how swimming speed depends on alignment angle. In the present paper we find that this discrepancy is due to missing terms in the analysis of Cupples et al., and that when these terms are correctly included, the models agree in their common limit. We further explore the predictions of the corrected model for both passive and active cases; it is found that for certain combinations of alignment angle and activity parameter, propulsion is halted; in other cases the small amplitude asymptotic expansion is no longer valid, motivating future numerical study.

Instability evolution in a transitional hypersonic boundary layer and its effects on aerodynamic heating are investigated over a 260 mm long flared cone. Experiments are conducted in a Mach 6 wind tunnel using Rayleigh-scattering flow visualization, fast-response pressure sensors, fluorescent temperature-sensitive paint (TSP) and particle image velocimetry (PIV). Calculations are also performed based on both the parabolized stability equations (PSE) and direct numerical simulations (DNS). Four unit Reynolds numbers are studied, 5.4, 7.6, 9.7 and $11.7\times 10^{6}~\text{m}^{-1}$. It is found that there exist two peaks of surface-temperature rise along the streamwise direction of the model. The first one (denoted as HS) is at the region where the second-mode instability reaches its maximum value. The second one (denoted as HT) is at the region where the transition is completed. Increasing the unit Reynolds number promotes the second-mode dissipation but increases the strength of local aerodynamic heating at HS. Furthermore, the heat generation rates induced by the dilatation and shear processes (respectively denoted as $w_{\unicode[STIX]{x1D703}}$ and $w_{\unicode[STIX]{x1D714}}$) were investigated. The former item includes both the pressure work $w_{\unicode[STIX]{x1D703}1}$ and dilatational viscous dissipation $w_{\unicode[STIX]{x1D703}2}$. The aerodynamic heating in HS mainly arose from the high-frequency compression and expansion of fluid accompanying the second mode. The dilatation heating, especially $w_{\unicode[STIX]{x1D703}1}$, was more than five times its shear counterpart. In a limited region, the underestimated $w_{\unicode[STIX]{x1D703}2}$ was also larger than $w_{\unicode[STIX]{x1D714}}$. As the second-mode waves decay downstream, the low-frequency waves continue to grow, with the consequent shear-induced heating increasing. The latter brings about a second, weaker growth of surface-temperature HT. A theoretical analysis is provided to interpret the temperature distribution resulting from the aerodynamic heating.