We study the coupled settling, deformation and mixing dynamics of a dense blob of fluid falling in an axially (vertically) linearly stratified Taylor–Couette cell (operated in a laminar stable regime). This configuration allows the independent analysis of stretching dynamics, driven by radial (horizontal) velocity variations, and settling dynamics, driven by buoyancy forces associated with vertical density variations. As the blob settles, it is stretched in the horizontal plane and forms an elongated lamella. Through the competing effects of transverse compression of the lamella due to this shear-induced stretching and broadening due to diffusion, the lamella irreversibly mixes with ambient fluid, thus progressively adjusting its own density towards that of the ambient fluid. Eventually, the lamella settling stops at a final equilibrium position that depends on the ambient vertical density gradient and the rate at which it has been deformed by the horizontal shear. We show how this final position is determined by stretching-enhanced diffusion, i.e. mixing. We demonstrate that a theoretical mixing model compares favourably with experiments with various Froude numbers (quantifying the relative strength of the horizontal shear and the vertical stratification) and construct a new criterion for the energetic ‘efficiency’ of this mixing process that explicitly captures its inherently diffusive character.

Ding (J. Fluid Mech., vol. 970, 2023, A27) analysed the dispersion of a multicomponent electrolyte solution flowing through a channel. An inequality in ionic diffusion coefficients induces a spontaneous electric field that leads to nonlinear coupling of species fluxes. Ding presented an effective equation for the long-time evolution of the ionic concentration distributions, which revealed novel features compared with the classic case of an uncharged solute. This work highlights the rich physics of ion diffusion and electro-migration in non-uniform flows.

This paper focuses on the linear evolution of Mack instability modes in a hypersonic boundary layer over a flat plate that is partially coated by a compliant section. The compliant section is a thin, flexible membrane covering on a porous wall consisting of micro holes. The instability pressure could induce a vibration of the membrane, leading to a feedback to the boundary-layer fluids through the transverse velocity fluctuation. Such a process is formulated by an admittance boundary condition for the boundary-layer perturbation, which is dependent on the thickness and tension of the membrane, the properties of the porous wall, and the frequency of the Mack-mode perturbation. Using this admittance condition, the impact of the compliant coating on the Mack growth rate is studied systematically by solving the compressible Orr–Sommerfeld equations. It is found that the compliant coating could suppress the Mack instability with a frequency band in the neighbourhood of the most unstable frequency, and the stabilising frequency band widens as the membrane thickness and tension decrease, indicating a more favourable effect of a softer membrane. For a Mack mode with a specified dimensional frequency – since its dimensionless frequency, normalised by the local boundary-layer thickness and oncoming velocity, increases as it propagates downstream – the second-mode frequency band usually appears in downstream locations, and so does the stabilising effect of the membrane. Thus it is favourable to apply a compliant panel at a downstream region. In this situation, the solid–compliant junction could produce an additional scattering effect on the evolution of the Mack mode due to the sudden change of its boundary condition. The scattering effect is quantified by a transmission coefficient defined by the equivalent amplitude of the compliant-wall perturbation to the solid-wall perturbation, which can be obtained by the harmonic linearised Navier–Stokes (HLNS) approach. If the admittance is weak, then the transmission coefficient can also be predicted by an analytical solution based on the residue theorem. It is found that most of the second modes are suppressed by the scattering effect as long as the argument of the admittance is in the interval $[150^\circ, 210^\circ ]$, agreeing with most of the physical situation. The analytical predictions agree well with the HLNS calculations when the modulus of the admittance is less than $O(0.1)$.

Soft porous materials, such as biological tissues and soils, are exposed to periodic deformations in a variety of natural and industrial contexts. The detailed flow and mechanics of these deformations have not yet been systematically investigated. Here, we fill this gap by identifying and exploring the complete parameter space associated with periodic deformations in the context of a one-dimensional model problem. We use large-deformation poroelasticity to consider a wide range of loading periods and amplitudes. We identify two distinct mechanical regimes, distinguished by whether the loading period is slow or fast relative to the characteristic poroelastic timescale. We develop analytical solutions for slow loading at any amplitude and for infinitesimal amplitude at any period. We use these analytical solutions and a full numerical solution to explore the localisation of the deformation near the permeable boundary as the period decreases and the emergence of nonlinear effects as the amplitude increases. We show that large deformations lead to asymmetry between the loading and unloading phases of each cycle in terms of the distributions of porosity and fluid flux.

Following our previous work on periodic ray paths (He et al., J. Fluid Mech., vol. 939, 2022, A3), we study asymptotically and numerically the structure of internal shear layers for very small Ekman numbers in a three-dimensional spherical shell and in a two-dimensional cylindrical annulus when the rays converge towards an attractor. We first show that the asymptotic solution obtained by propagating the self-similar solution generated at the critical latitude on the librating inner core describes the main features of the numerical solution. The internal shear layer structure and the scaling for its width and velocity amplitude in $E^{1/3}$ and $E^{1/12}$, respectively, are recovered. The amplitude of the asymptotic solution is shown to decrease to $E^{1/6}$ when it reaches the attractor, as is also observed numerically. However, some discrepancies are observed close to the particular attractors along which the phase of the wave beam remains constant. Another asymptotic solution close to those attractors is then constructed using the model of Ogilvie (J. Fluid Mech., vol. 543, 2005, pp. 19–44). The solution obtained for the velocity has an $O(E^{1/6})$ amplitude, but a self-similar structure different from the critical-latitude solution. It also depends on the Ekman pumping at the contact points of the attractor with the boundaries. We demonstrate that it reproduces correctly the numerical solution. Surprisingly, the numerical solution close to an attractor with phase shift (that is, an attractor touching the axis in three or two dimensions with a symmetric forcing) is also found to be $O(E^{1/6})$, but its amplitude is much weaker. However, its asymptotic structure remains a mystery.

Most gravitational currents occur on sloping topographies, often in the presence of particles that can settle during the current propagation. Yet an exhaustive exploration of associated parameters in experimental studies is still lacking. Here, we present an extensive experimental investigation of the slumping regime of turbidity (particle-laden) currents in two lock-release (dam-break) systems with inclined bottoms. We identify three regimes controlled by the ratio between settling and current inertia. (i) For negligible settling, the turbidity current morphodynamics corresponds to that of saline homogeneous gravity currents, in terms of velocity, slumping (constant-velocity) regime duration and current morphology. (ii) For intermediate settling, the slumping regime duration decreases to become fully controlled by a particle settling characteristic time. (iii) When settling overcomes the current initial inertia, the slumping (constant-velocity) regime is no longer detected. In the first two regimes, the current velocity increases with the bottom slope, of approximately $35\,\%$ between $0^\circ$ and $15^\circ$. Finally, our experiments show that the current propagates during the slumping regime with the same shape in the frame of the moving front. Strikingly, the current head is found to be independent of all experimental parameters covered in the present study. We also quantify water entrainment coefficients $E$ and compare them with previous literature, hence finding that $E$ increases rather linearly with the current Reynolds number.

This study experimentally investigates the influence of structural damping on the transverse flow-induced vibration (FIV) of an elastically mounted thin elliptical cylinder. The cylinder tested has an elliptical ratio of $\varepsilon = b/a = 5$, where $a$ and $b$ are the streamwise and cross-flow dimensions, respectively, and a mass ratio (i.e. the total oscillating mass/the displaced fluid mass) of $17.4$. The FIV response was characterised over a reduced velocity range of $2.30 \leq U^* = U/(\,{{f_{{nw}}}} b) \leq 10.00$ (corresponding to a Reynolds number range of $300 \leq \textit {Re} =(U b)/\nu \leq 1300$) and a structural damping ratio range of $3.62\times 10^{-3}\leq \zeta \leq 1.87\times 10^{-1}$. Here, $U$ is the free stream velocity, ${{f_{{nw}}}}$ is the natural frequency of the system in quiescent fluid (water) and $\nu$ is the kinematic viscosity of the fluid. The FIV response was characterised by four wake–body synchronisation regimes (defined as the matching of the dominant fluid forcing and oscillation frequencies, and labelled regime I, regime II, regime III and the hyper branch) and a desynchronisation region, with the hyper branch representing a high amplitude regime not observed for a circular cylinder. Interestingly, the major vortex shedding mode was predominately two single opposite-signed vortices shed per body vibration cycle. Moreover, hydrogen-bubble-based flow visualisations revealed a secondary vortex street forming in the elongated shear layers associated with largest-scale vibration amplitudes ($A^* = A/b$ up to $7.7$) in the hyper branch and regime II. As the structural damping ratio was increased beyond $1.92 \times 10^{-2}$, the hyper branch was found to be suppressed. The results have potential ramifications for the efficient extraction of energy from free-flowing water sources, which has become increasingly topical over the last decade.

In this study, we performed a set of experiments of the thermal plume on an open cylinder heated from below and a simple scaling analysis. The flow structure of the plume was visualized using the shadowgraph technique, and the temperature at specific points was measured using a thermistor. Transient plumes in a developing stage, an equilibrating stage and a fully developed stage were described. The new scaling laws of the stem radius and the velocity of the rising plume were presented, which are different from those on a two-dimensional heated plate due to the cylindrical effect and agree with experimental results. In the fully developed stage, there is a transition route of the plume from a steady to chaotic state with an increase in the Rayleigh number, which involves a series of bifurcations. A reverse bifurcation from a periodic (Ra = 1.17 × 106) back to a steady (Ra = 1.30 × 106) state has been observed in the experiment, which is referred to as the period bubbling bifurcation. Thus, the present experimental results validate the previous numerical results presented by Zhang et al. (Phys. Fluids, vol. 33 (6), 2021, 064110). The bifurcation diagram, spectral analysis and attractor are used to characterize the transitional plume. In addition, the heat and mass transfer have also been quantified.

From biological tissues to microactuators and absorption of solvents into layers of paint, macroscopically non-porous materials with the capacity to swell when in contact with a solvent are ubiquitous. In these systems, owing to strong solid–fluid interactions, chemically driven flows can yield large geometric changes. We study experimentally and theoretically the canonical problem of the swelling of a thin hydrogel layer by a single water drop. Using a bespoke experimental set-up, we observe fast absorption leading to a radially spreading axisymmetric blister. We use a fully three-dimensional linear poroelastic framework with nonlinear kinematic equations to obtain governing equations, which we then reduce with thin-layer scalings to a one-dimensional nonlinear diffusion equation for the evolution of the blister geometry. In the limits of large and small deformations, the evolution of the blister characteristic height and radius are self-similar, following power laws in time. Our experimental measurements show that the evolution of the blister is broadly within the theoretical predictions in the large and small deformation regimes. In the general intermediate deformation regime, the measurements are well captured by our reduced one-dimensional diffusion model, which does not require the sophisticated and computationally expensive numerical approaches necessary for the original two-dimensional nonlinear coupled transport problem. By adapting modelling techniques from the fluid dynamics of thin porous elastic layers to a polymer swelling problem, our modelling framework extends the range of polymer swelling problems that can be treated with semi-analytical methods. Moreover, our detailed experimental data can serve as a test case for future nonlinear poroelastic frameworks of swelling polymer materials.

Accelerated glacier melt and the loss of perennial snowfields have been associated with increased warming in polar regions, at rates up to four times faster than the rest of the world, thereby reinforcing the critical need for improved models (and predictions) of glacier melt. An essential requirement for such models is an improved understanding of the sensible heat fluxes over glaciers. Since their complexity makes them difficult to model, and direct measurements of sensible turbulent heat fluxes over real glaciers are both rare and impractical, the present work involves simultaneous hot-wire anemometry and cold-wire thermometry measurements of two components of velocity and temperature above a melting glacier model in a series of wind-tunnel experiments. Both single- and multi-variable statistics were used to compare the turbulent velocity field measured over melting ice with that of a similar flow in the absence of ice. The results demonstrate that the ice's presence reduces the magnitude of the Reynolds stresses and vertical velocity variance, but also increases the streamwise velocity variance. The transient evolution of temperature statistics throughout the melt process was also investigated and found to be similar when suitably non-dimensionalized. The velocity and temperature fields were furthermore evaluated at an equivalent non-dimensional time during the melt process, in which statistics of the temperature field, and joint statistics of the vertical velocity and temperature, were studied. The present work lays the foundation for future laboratory-scale replications of the flow above melting glaciers, and provides additional insight into turbulent heat transfer over melting ice.

The effect of variations in the integral length scale of incoming free-stream turbulence on a NACA0012 wing is investigated with the use of force, moment and particle image velocimetry measurements. At a chord-based Reynolds number ($Re = U_\infty c/\nu$ where c is the chord length, $U_\infty$ is the free-stream velocity and $\nu$ is the kinematic viscosity) of $2\times 10^5$, an active grid generates turbulence intensities of 15 % at normalised integral length scales ranging from 0.5$c$ to 1$c$. The introduction of turbulence improves the time-averaged performance characteristics of the wing by delaying stall and increasing the peak lift coefficient. It is found that for half-chord integral length scales, the magnitude of the fluctuations in forces and moments is larger than that of full-chord integral length scales, as the former amplifies the naturally occurring unsteadiness in the flow (when there is no free-stream turbulence). The increase in magnitude is ascribed to a larger density of smaller-scale vortices within the separated flow and wake region of the wing.

High-speed vehicles experience a highly challenging environment in which the freestream Mach number and surface temperature greatly influence aerodynamic drag and heat transfer. The interplay of these two parameters strongly affects the near-wall dynamics of high-speed turbulent boundary layers (TBLs) in a non-trivial way, breaking similarity arguments on velocity and temperature fields, typically derived for adiabatic cases. We present direct numerical simulations of flat-plate zero-pressure-gradient TBLs spanning three freestream Mach numbers $[2,4,6]$ and four wall temperature conditions (from adiabatic to very cold walls), emphasising the choice of the wall-cooling parameter to recover a similar flow organisation at different Mach numbers. We link qualitative observations on flow patterns to first- and second-order statistics to explain the decoupling of temperature–velocity fluctuations that occurs at reduced wall temperatures and high Mach numbers. For these cases, we discuss the formation of a secondary peak of thermal production in the viscous sublayer, which is in contrast with the monotonic behaviour of adiabatic profiles. We propose different physical mechanisms induced by wall-cooling and compressibility that result in apparently similar flow features, such as a higher peak in the streamwise velocity turbulence intensity, and distinct features, such as the separation of turbulent scales.

The irreversible behaviour of a highly confined non-Brownian suspension of spherical particles at low Reynolds number in a Newtonian fluid is studied experimentally and numerically. In the experiment, the suspension is confined in a thin rectangular channel that prevents complete particle overlap in the narrow dimension and is subjected to an oscillatory pressure-driven flow. In the small cross-sectional dimension, particles rapidly separate to the walls, whereas in the large dimension, features reminiscent of shear-induced migration in bulk suspensions are recovered. Furthermore, as a consequence of the channel geometry and the development and application of a single-camera particle tracking method, three-dimensional particle trajectories are obtained that allow us to directly associate relative particle proximity with the observed migration. Companion simulations of a steadily flowing suspension highly confined between parallel plates are conducted using the force coupling method, which also show rapid migration to the walls as well as other salient features observed in the experiment. While we consider relatively low volume fractions compared to most prior work in the area, we nevertheless observe significant and rapid migration, which we attribute to the high degree of confinement.

The effect of single perforations and their location on the drag and reconfiguration of flexible plates was explored through laboratory experiments and direct numerical simulations. The plates were subjected to uniform flows with negligible turbulence, and the perforations had a square cross-section resulting in a low porosity ratio of $\gamma \approx 0.028$. Rigid plates with and without perforations and flexible plates without perforations served as the baseline cases. The perforated plates exhibited distinct jets through the openings in the wake, significantly impacting the aerodynamic force and plate deformation. The velocity and position of the centre jet velocity in relation to downwind distance were influenced by both the incoming flow and the location of the perforations. The centre jet velocity profiles were normalized using an effective velocity and corrected perforation half-width, revealing their dependence on these factors. A simple first-order formulation was developed to predict the change in drag for various perforated plates under a wide range of incoming velocities. This formulation was supported by numerical simulations across a wider range of Cauchy number to confirm the proposed model and separate the effect of the Cauchy and Reynolds numbers. The results of this study may inform the design of flexible structures, define effective porosity and serve as an initial step towards modelling the complex interaction between flow and structures with low porosity.

A microscale lubrication flow of a gas between eccentric circular cylinders is studied on the basis of kinetic theory. The dimensionless curvature, defined by the mean clearance divided by the radius of the inner cylinder, is small, and the rotation speed of the inner cylinder is also small. The Knudsen number, defined by the mean free path divided by the mean clearance, is arbitrary. The Boltzmann equation is studied analytically using the slowly varying approximation following the method proposed in the author's previous study (Doi, Phys. Rev. Fluids, vol. 7, 2022, 034201). A macroscopic lubrication equation, which is a microscale generalization of the Reynolds lubrication equation, is derived. To assess this, a direct numerical analysis of the Boltzmann equation in a bipolar coordinate system is conducted using the Bhatnagar–Gross–Krook–Welander kinetic equation. It is demonstrated that the solution of the derived lubrication equation approximates that of the Boltzmann equation over a wide range of the eccentricity and the whole range of the Knudsen number. It is also demonstrated that another lubrication equation derived by a formal application of the slowly varying approximation produces a non-negligible error of the order of the square root of the dimensionless curvature for large Knudsen numbers.

Deep-water surface wave breaking affects the transfer of mass, momentum, energy and heat between the air and sea. Understanding when and how the onset of wave breaking will occur remains a challenge. The mechanisms that form unforced steep waves, i.e. nonlinearity or dispersion, are thought to have a strong influence on the onset of wave breaking. In two dimensions and in deep water, spectral bandwidth is the main factor that affects the roles these mechanism play. Existing studies, in which the relationship between spectral bandwidth and wave breaking onset is investigated, present varied and sometimes conflicting results. We perform potential-flow simulations of two-dimensional focused wave groups on deep water to better understand this relationship, with the aim of reconciling existing studies. We show that the way in which steepness is defined may be the main source of confusion in the literature. Locally defined steepness at breaking onset reduces as a function of bandwidth, and globally defined (spectral) steepness increases. The relationship between global breaking onset steepness and spectral shape (using the parameters bandwidth and spectral skewness) is too complex to parameterise in a general way. However, we find that the local surface slope of maximally steep non-breaking waves, of all spectral bandwidths and shapes that we simulate, approaches a limit of $1/\tan ({\rm \pi} /3)\approx 0.5774$. This slope-based threshold is simple to measure and may be used as an alternative to existing kinematic breaking onset thresholds. There is a potential link between slope-based and kinematic breaking onset thresholds, which future work should seek to better understand.

Modal and non-modal linear stability analyses are employed to investigate the effect of internal and external heating on disturbance temporal growth for the Darcy–Bénard convection with throughflow. A matrix-forming approach is employed for both purposes, where the generalised eigenvalue problem is built using the generalised integral transform technique. Although the disturbance equations are not self-adjoint, the non-modal analysis indicates that there is no transient growth. Hence, any disturbance growth in time must be induced by modal mechanisms. An absolute instability analysis reveals that viscous dissipation has a destabilising effect and introduces new modes that are eventually destabilised by increasing the Péclet number. Beyond critical values of the Péclet number, where codimension-two absolutely unstable points exist, these modes become more unstable than the classical mode found in the absence of viscous dissipation, which is stabilised by an increasing Péclet number. This internal heating mechanism generated by viscous dissipation is so strong at high enough Péclet numbers that instability becomes possible through heating from above.

Turbulent flows above a solid surface are characterised by a hydrodynamic roughness that represents, for the far velocity field, the typical length scale at which momentum mixing occurs close to the surface. Here, we are theoretically interested in the hydrodynamic roughness induced by a two-dimensional modulated surface, the elevation profile of which is decomposed in Fourier modes. We describe the flow for a sinusoidal mode of given wavelength and amplitude with Reynolds-averaged Navier–Stokes equations closed by means of a mixing-length approach that takes into account a possible surface geometrical roughness as well as the presence of a viscous sublayer. It also incorporates spatial transient effects at the laminar–turbulent transition. Performing a weekly nonlinear expansion in the bedform aspect ratio, we predict the effective hydrodynamic roughness when the surface wavelength is varied and we show that it presents a non-monotonic behaviour at the laminar–turbulent transition when the surface is hydrodynamically smooth. Further, with a self-consistent looped calculation, we are able to recover the smooth–rough transition of a flat surface, for which the hydrodynamic roughness changes from a regime where it is dominated by the viscous length to another one where it scales with the surface corrugation. We finally apply the results to natural patterns resulting from hydrodynamic instabilities such as those associated with dissolution or sediment transport. We discuss in particular the aspect ratio selection of dissolution bedforms and roughness hierarchy in superimposed ripples and dunes.

Direct numerical simulations up to Reλ = 1445 show that the scaling exponents for the enstrophy and the dissipation rate extrema are different and depend on the Reynolds number. A similar Reynolds number dependence of the scaling exponents is observed for the moments of the dissipation rate, but not for the moments of the enstrophy. Significant changes in the exponents occur at approximately Reλ ≈ 250, where Reλ is the Taylor based Reynolds number, which coincides with structural changes in the flow, in particular the development of large-scale shear layers. A model for the probability density functions (PDFs) of the enstrophy and dissipate rate is presented, which is an extension of our existing model (Proc. R. Soc. A, vol. 476, 2020, p. 20200591) and is based on the mentioned development of large-scale layer regions within the flow. This model is able to capture the observed Reynolds number dependencies of the scaling exponents, in contrast to the existing theories which yield constant exponents. Moreover, the model reconciles the scaling at finite Reynolds number with the theoretical limit, where the enstrophy and dissipation rate scale identically at infinite Reynolds number. It suggests that the large-scale shear layers are vital for understanding the scaling of the extrema. Furthermore, to reach the theoretical limit, the scaling exponents must remain Reynolds number dependent beyond the present Reλ range.

A theoretical and experimental investigation of two-dimensional (2-D) liquid curtains (gravitationally thinning liquid sheets) is provided under conditions where the curtain issues from a thin slot whose centreline is inclined with respect to the vertical. This analysis is motivated in part by recent works where it has been proposed that oblique liquid curtains (those exiting a non-vertical slot) may bend upwards against gravity when the relevant Weber number at the slot is less than unity ($We <1$). By contrast, Weinstein et al. (J. Fluid Mech., vol. 876, 2019, R3) have proposed that such $We<1$ curtains must be vertical and downward falling regardless of the inclination of the slot. Under low-Reynolds-number ($Re$) conditions typical of liquid film coating operations, our experiments show that the curtain shape follows the classic ballistic (parabolic) trajectory in the supercritical regime ($We>1$). In subcritical conditions ($We<1$), experiments show that the downward-falling curtain is vertical except in a relatively small region near the slot, where the combined effects of viscosity and surface tension induce the so-called teapot effect. These experimental results are confirmed by 2-D numerical simulations, which predict the curtain behaviour ranging from highly viscous ($Re = O(1)$) to nearly inviscid conditions. The one-dimensional (1-D) inviscid model of Weinstein et al. is recast in a different form to facilitate comparisons with the 2-D model, and 1-D and 2-D results agree favourably for supercritical and subcritical conditions. Despite the large parameter range explored, we have found no evidence that upward-bending curtains exist in an oblique configuration.