Evaporation of multi-component liquid mixtures in confined geometries, such as capillaries, is crucial in applications such as microfluidics, two-phase cooling devices and inkjet printing. Predicting the behaviour of such systems becomes challenging because evaporation triggers complex spatio-temporal changes in the composition of the mixture. These changes in composition, in turn, affect evaporation. In the present work, we study the evaporation of aqueous glycerol solutions contained as a liquid column in a capillary tube. Experiments and direct numerical simulations show three evaporation regimes characterised by different temporal evolutions of the normalised mass transfer rate (or Sherwood number $Sh$), namely $Sh (\tilde{t} ) = 1$, $Sh \sim 1/\sqrt {\tilde{t} }$ and $Sh \sim \exp (-\tilde{t} )$, where $\tilde {t}$ is a normalised time. We present a simplistic analytical model that shows that the evaporation dynamics can be expressed by the classical relation $Sh = \exp ( \tilde{t} )\,\mathrm {erfc} ( \sqrt {\tilde{t} })$. For small and medium $\tilde{t}$, this expression results in the first and second of the three observed scaling regimes, respectively. This analytical model is formulated in the limit of pure diffusion and when the penetration depth $\delta (t)$ of the diffusion front is much smaller than the length $L(t)$ of the liquid column. When $\delta \approx L$, finite-length effects lead to $Sh \sim \exp (-\tilde{t} )$, i.e. the third regime. Finally, we extend our analytical model to incorporate the effect of advection and determine the conditions under which this effect is important. Our results provide fundamental insights into the physics of selective evaporation from a multi-component liquid column.

It is generally believed that the velocity and passive scalar fields share many similarities and differences in wall-bounded turbulence. In the present study, we conduct a series of direct numerical simulations of compressible channel flows with passive scalars and employ the two-dimensional spectral linear stochastic estimation and the correlation function as diagnostic tools to shed light on these aspects. Particular attention is paid to the relevant multiphysics couplings in the spectral domain, i.e. the velocity–temperature ($u-T$), scalar–temperature ($g-T$) and velocity–scalar ($u-g$) couplings. These couplings are found to be utterly different at a given wall-normal position in the logarithmic and outer regions. Specifically, in the logarithmic region, the $u-T$ and $u-g$ couplings are tight at the scales that correspond to the attached eddies and the very large-scale motions (VLSMs), whereas the $g-T$ coupling is robust in the whole spectral domain. In the outer region, $u-T$ and $u-g$ couplings are only active at the scales corresponding to the VLSMs, whereas the $g-T$ coupling is diminished but still strong at all scales. Further analysis indicates that although the temperature field in the vast majority of zones in a channel can be roughly treated as a passive scalar, its physical properties gradually deviate from those of a pure passive scalar as the wall-normal height increases due to the enhancement of the acoustic mode. Furthermore, the deep involvement of the pressure field in the self-sustaining process of energy-containing motions also drives the streamwise velocity fluctuation away from a passive scalar. The current work is an extension of our previous study (Cheng & Fu, J. Fluid Mech., vol. 964, 2023, A15), and further uncovers the details of the multiphysics couplings in compressible wall turbulence.

Plane Couette flow at Reynolds number $Re=1200$ (based on the channel half-height and half the velocity difference between the top and bottom plates) is investigated with a spatial domain designed to retain only two spanwise integral length scales. In this system, the computation of invariant solutions that are physically representative of the turbulent state has been understood to be challenging. To address this challenge, our approach is to employ an accurate reduced-order model with 600 degrees of freedom (Cavalieri & Nogueira, Phys. Rev. Fluids, vol. 7, 2022, L102601). Using the two-scale energy budget and the temporal cross-correlation of key observables, it is first demonstrated that the model contains most of the multi-scale physical processes identified recently (Doohan et al., J. Fluid Mech., vol. 913, 2021, A8); i.e. the large- and small-scale self-sustaining processes, the energy cascade for turbulent dissipation, and an energy-cascade mediated small-scale production mechanism. Invariant solutions of the reduced-order model are subsequently computed, including 96 equilibria and 43 periodic orbits. It is found that none of the computed equilibrium solutions are able to reproduce an accurate energy balance associated with the multi-scale dynamics of the turbulent state. Incorporation of unsteadiness into invariant solutions is seen to be essential for a sensible description of the multi-scale turbulent dynamics and the related energetics, at least in this type of flow, as periodic orbits with a sufficiently long period are mainly able to describe the complex spatio-temporal dynamics associated with the known multi-scale phenomena.

Although poroelastic clusters in nature, such as bristled wings and plumed seeds, exhibit remarkable flight performances by virtue of their porous structure, the effects of another key feature, elasticity, on aerodynamic loading remain elusive. For a poroelastic cluster, we investigate the aerodynamic effects of elastic deformation that occurs through the collective rearrangement of many elastic components and the fluid-dynamic interactions between them. As a simple two-dimensional model, an array of multiple cylinders which are individually and elastically mounted is employed with diverse values of porosity and elasticity. Under a uniform free stream, the poroelastic cluster enlarges its frontal area and augments the total drag force in the quasi-steady state; this is in contrast to the general reconfiguration of fixed elastic structures, which tends to reduce the frontal area and drag. The rearrangement of the poroelastic cluster is dominated by the virtual fluid barrier that develops in a gap between the elastic components, interrupting the flow penetrating between them. The effects of this hydrodynamic blockage on changes in the frontal area and drag force are analysed in terms of porosity and elasticity, revealing the fluid-dynamic mechanism underlying the appearance of peak drag at an intermediate porosity. Moreover, to represent the coupled effects of porosity and elasticity on the rearrangement, a scaled elastic energy is derived through a consideration of the energy balance.

The cubic interactions in a discrete system of four weakly nonlinear waves propagating in a conservative dispersive medium are studied. By reducing the problem to a single ordinary differential equation governing the motion of a classical particle in a quartic potential, the complete explicit branches of solutions are presented, either steady, periodic, breather or pump, thus recovering or generalizing some already published results in hydrodynamics, nonlinear optics and plasma physics, and presenting some new ones. Various stability criteria are also formulated for steady equilibria. Theory is applied to deep-water gravity waves for which models of isolated quartets are described, including bidirectional standing waves and quadri-directional travelling waves, steady or not, resonant or not.

Equations describing mushy systems, in which solid and liquid are described as a single continuum, have been extensively studied. Most studies of mushy layers have assumed them to be ‘ideal’, such that the liquid and solid were in perfect thermodynamic equilibrium. It has become possible to simulate flows of passive porous media at the pore scale, where liquid and solid are treated as separate continua. In this contribution, we study the simplest possible mushy layers at the pore scale, modelling a single straight cylindrical pore surrounded by a cylindrical annulus representing the solid matrix. Heat and solute can be exchanged at the solid–liquid boundary. We consider harmonic temperature and concentration perturbations and examine their transport rates due to advection and diffusion and the melting and solidification driven by this transport. We compare the results of this numerical model with those of a one-dimensional ideal mushy layer and with analytical solutions valid for ideal mushy layers for small temperature variations. We demonstrate that for small values of an appropriately defined Péclet number, the results of the pore-scale model agree well with ideal mushy layer theory for both transport rates and melting. As this Péclet number increases, the temperature and concentration profiles with radius within the pore differ significantly from constant, and the behaviour of the pore-scale model differs significantly from that of an ideal mushy layer. Some effects of mechanical dispersion arise naturally in our pore-scale model and are shown to be important at high Péclet number.

This paper examines the rigid body motion of a spheroid sedimenting in a Newtonian fluid with a spatially varying viscosity field. The fluid is at zero Reynolds number, and the viscosity varies linearly in space in an arbitrary direction with respect to the external force. First, we obtain the correction to the spheroid's rigid body motion in the limit of small viscosity gradients, using a perturbation expansion combined with the reciprocal theorem. Next, we determine the general form of the particle's mobility tensor relating its rigid body motion to an external force and torque. The viscosity gradient does not alter the force/translation and torque/rotation relationships, but introduces new force/rotation and torque/translation couplings that are determined for a wide range of particle aspect ratios. Finally, we discuss results for the spheroid's rotation and centre-of-mass trajectory during sedimentation. A steady orientation arises at long time whose value depends on the viscosity gradient direction and particle shape. These results are significantly different than when no viscosity gradient is present, where the particle stays at its initial orientation for all times. We summarize the observations for prolate and oblate spheroids for different viscosity gradient directions and provide plots for the orientation and centre-of-mass trajectory versus time. We also provide guidelines to extend the analysis when the viscosity gradient exhibits a more complicated spatial behaviour.

In this work, we analyse wall-bounded flows in the continuum to transition regime with the help of higher-order transport equations (super-set of the Navier–Stokes equation). Towards this, we incorporate second-order in Knudsen number accurate terms in the single-particle distribution function, and with this complete representation, we first derive the second-order accurate extended-OBurnett (EOBurnett) and third-order accurate super-OBurnett (SOBurnett) equations. We then demonstrate that these newly derived equations exhibit unconditional linear stability. We finally validate the equations by solving for plane Poiseuille flow and derive closed-form analytical solutions for the pressure and velocity fields. The pressure and velocity results thus obtained have been compared with direct simulation Monte Carlo (DSMC) data in the transition regime. The results from both the EOBurnett and SOBurnett equations are found to yield better agreement with DSMC data than that obtained from the Navier–Stokes equations. This improved agreement is attributed to the presence of additional terms in the proposed equations, which effectively capture the effect of the Knudsen layer near the wall. The obtained higher-order transport equations and the closed-form solution presented in this work are novel. The ability of the equations to describe the flow in the transition regime should form the basis for conducting further realistic analytical studies of wall-bounded flows in the future.

We present evidence for layered anisotropic stratified turbulence (LAST) and mixing produced in a freely evolving uniformly stratified shear layer where the direction of shear is orthogonal to gravity. As originally reported by Basak & Sarkar (J. Fluid. Mech., vol. 568, 2006, pp. 19–54), such a flow develops a rich three-dimensional structure in the form of interlocking columnar vortices formed by horizontal shear instability that remain coherent at large scales due to the stabilising vertical stratification. Here, we modify the initial velocity field by introducing additional small-amplitude vertical perturbations designed to be representative of pre-existing horizontal layers often observed in strongly stratified ocean environments. This reveals a novel finite amplitude, non-normal mode growth mechanism through which the vertical shear between layers may be rapidly amplified by its interaction with the horizontal shear layer prior to the growth of shear instability, leading to a rapid turbulent transition instigated by the subsequent interaction of the layers with the emerging columnar vortices. Through a consideration of relevant flow statistics and associated dimensionless parameters, we demonstrate that turbulence can enter the LAST regime, thereby indicating a generic mechanism leading to the transient development of regions of strongly stratified turbulence in the ocean. We discuss the properties of mixing and the parameterisation of mixing efficiency in terms of the relationship between turbulent length scales in the flow, in particular highlighting links to models based on the classical vertical shear instability paradigm typically associated with more weakly stratified flows that produce isolated turbulent ‘patches’.

Direct numerical simulations are utilised to investigate mass-transfer processes at gas-evolving electrodes that experience successive formation and detachment of bubbles. The gas–liquid interface is modelled employing an immersed boundary method. We simulate the growth phase of the bubbles followed by their departure from the electrode surface in order to study the mixing induced by these processes. We find that the growth of the bubbles switches from a diffusion-limited mode at low to moderate fractional bubble coverages of the electrode to a reaction-limited growth dynamics at high coverages. Furthermore, our results indicate that the net transport within the system is governed by the effective buoyancy driving induced by the rising bubbles and that mechanisms commonly subsumed under the term ‘microconvection’ do not significantly affect the mass transport. Consequently, the resulting gas transport for different bubble sizes, current densities and electrode coverages can be collapsed onto one single curve and only depends on an effective Grashof number. The same holds for the mixing of the electrolyte when additionally taking the effect of surface blockage by attached bubbles into account. For the gas transport to the bubble, we find that the relevant Sherwood numbers also collapse onto a single curve when accounting for the driving force of bubble growth, incorporated in an effective Jakob number. Finally, linking the hydrogen transfer rates at the electrode and the bubble interface, an approximate correlation for the gas-evolution efficiency has been established. Taken together, these findings enable us to deduce parametrisations for all response parameters of the systems.

The impact of a liquid droplet with another droplet or onto a solid surface are important basic processes that occur in many applications such as agricultural sprays and inkjet printing, and in nature such as pathogens transport by raindrops. We investigated the head-on collision of unequal-size droplets of the same liquid on wetting surfaces using the direct numerical simulations technique at different size ratios. The unsteady Navier–Stokes equations are solved and the liquid–gas interface is tracked using the geometric volume-of-fluid method. The numerical model is validated by comparing simulation results of two extreme cases of droplets bouncing with the experimental data from previous studies and the agreement is quite accurate. The validated model is employed to simulate droplets bouncing at several size ratios at different Weber numbers and Ohnesorge number. Two distinct regimes are identified, namely, the inertial regime, where the restitution coefficient is a constant value close to 0.3, the viscous regime, where the restitution coefficient declines. To understand the bouncing behaviour, the velocity field is analysed and an energy budget calculation is performed. The distribution of the sessile droplet energy is found to be important and the sessile droplet surface energy is calculated by its deformation characteristics such as crater depth. Finally, a scaling analysis is performed to rationalize the insensitivity of the coefficient of restitution in the inertial regime, and its decline in the viscous regime, at large size ratios.

The squirmer is a popular model to analyse the fluid mechanics of a self-propelled object, such as a micro-organism. We demonstrate that some fore–aft symmetric squirmers can spontaneously self-propel above a critical Reynolds number. Specifically, we numerically study the effects of inertia on spherical squirmers characterised by an axially and fore–aft symmetric ‘quadrupolar’ distribution of surface-slip velocity; under creeping-flow conditions, such squirmers generate a pure stresslet flow, the stresslet sign classifying the squirmer as either a ‘pusher’ or ‘puller’. Assuming axial symmetry, and over the examined range of the Reynolds number $Re$ (defined based upon the magnitude of the quadrupolar squirming), we find that spontaneous symmetry breaking occurs in the puller case above $Re \approx 14.3$, with steady swimming emerging from that threshold consistently with a supercritical pitchfork bifurcation and with the swimming speed growing monotonically with $Re$.

Fine-scale flow motions are measured in a Newtonian and polymer drag-reduced turbulent boundary layer (TBL) at a common momentum thickness Reynolds number $Re_{\theta }$ of 2300. Relative to the Newtonian TBL, the polymer-laden flow has a 33 % lower skin-friction coefficient. Three-dimensional (3-D) particle tracking velocimetry is used to measure the components of the velocity gradient tensor (VGT), rate of deformation tensor (RDT) and rate of rotation tensor (RRT). The invariants in these tensors are then used to distinguish the different types of fine-scale flow motions – a method called the $\varDelta$-criterion. Joint probability density functions (j.p.d.f.s) of the VGT invariants, $Q$ and $R$, for the Newtonian TBL produce the familiar tear-drop pattern, commonly seen in direct numerical simulations of Newtonian turbulence. Relative to the Newtonian TBL, the polymer-laden flow has significantly attenuated values of $R$, implying an overall reduction in fluid stretching. The invariants in the RDT, $Q_D$ and $R_D$, imply that straining motions of the polymeric flow are more two dimensional compared with the Newtonian flow. Moreover, j.p.d.f.s of $Q_D$ and the invariant in the RRT $Q_W$, suggest that the flow consists of fewer biaxial extensional events and more shear-dominated flow. Few, if any, experimental investigations have measured the 3-D structure of fine-scale motions in a Newtonian and polymer drag-reduced TBL using the $\varDelta$-criterion. We provide the first experimental evidence that supports the notion that an attenuation of fluid stretching, particularly biaxial straining motions, is central to the mechanism of polymer drag reduction.

Flow physics vary in different regimes across the full Mach number range, with our knowledge being particularly poor about the hypersonic regime. An Eulerian realization of the particles on demand method, a kinetic model formulated in the comoving reference frame, is proposed to simulate hypersonic compressible flows. The present model allows for flux evaluation in different reference frames, in this case rescaled and shifted by local macroscopic quantities, i.e. fluid speed and temperature. The resulting system of coupled hyperbolic equations is discretized in physical space with a finite volume scheme ensuring exact conservation properties. Regularization via Grad expansion is introduced to implement distribution function and flux transformation between different reference frames. It is shown that the proposed method possesses Galilean invariance at a Mach number up to $100$. Different benchmarks including both inviscid and viscous flows are reproduced with the Mach number up to $198$ and pressure ratio up to $10^5$. Finally, the new model is demonstrated to be capable of simulating hypersonic reactive flows, including one-dimensional and two-dimensional detonations. The developed methodology opens up possibilities for the simulation of the full range of compressible flows, without or with chemical reactions, from the subsonic to hypersonic regimes, leading to enhanced understanding of flow behaviours across the full Mach number range.

The familiar process of bubbles generated via breaking waves in the ocean is foundational to many natural and industrial applications. In this process, large pockets of entrained gas are successively fragmented by the ambient turbulence into smaller and smaller bubbles. The key question is how long it takes for the bubbles to reach terminal sizes for a given system. Despite decades of effort, the reported breakup time from multiple experiments differs significantly. Here, to reconcile those results, rather than focusing on one scale, we measure multiple time scales associated with the process through a unique experiment that resolves bubbles’ local deformation and curvature. The results emphasize that the scale separation among various time scales is controlled by the Weber number, similar to how the Reynolds number determines the scale separation in single-phase turbulence, but shows a distinct transition at a critical Weber number.

Large eddy simulation of flow past a cricket ball with its seam at $30^\circ$ to the free stream is carried out for $5 \times 10^4 \le Re \le 4.5 \times 10^5$. Three regimes of flow are identified on the basis of the time-averaged swing force coefficient ($\bar {C}_Z$) – no swing (NS), conventional swing (CS, $\bar {C}_Z>0$) and reverse swing (RS, $\bar {C}_Z<0$). The effect of seam on the boundary layer is investigated. Contrary to the popular belief, the boundary layer does not transition to a turbulent state in the initial stages of CS. The seam energizes the laminar boundary layer and delays its separation. The delay is significantly larger in a region near the poles, whose extent increases with an increase in $Re$ causing $\bar {C}_Z$ to increase. Here $\bar {C}_Z$ assumes a near constant value in the later stage of CS. The boundary layer transitions to a turbulent state via formation of a laminar separation bubble (LSB) in the equatorial region and directly, without a LSB, in the polar region. The extent of the LSB shrinks while the region of direct transition near the poles increases with an increase in $Re$. A LSB forms on the non-seam side of the ball in the RS regime. A secondary vortex is observed in the wake bubble. While it exists on the non-seam side for the entire range of $Re$ considered, the mixing in the flow introduced by the seam causes it to disappear beyond a certain $Re$ on the seam side. The pressure difference between the seam and non-seam sides sets up wing-tip-like vortices. Their polarity reverses with the switch from the CS to RS regime.

This study presents a noise-robust closed-loop control strategy for wake flows employing model predictive control. The proposed control framework involves the autonomous offline selection of hyperparameters, eliminating the need for user interaction. To this purpose, Bayesian optimization maximizes the control performance, adapting to external disturbances, plant model inaccuracies and actuation constraints. The noise robustness of the control is achieved through sensor data smoothing based on local polynomial regression. The plant model can be identified through either theoretical formulation or using existing data-driven techniques. In this work we leverage the latter approach, which requires minimal user intervention. The self-tuned control strategy is applied to the control of the wake of the fluidic pinball, with the plant model based solely on aerodynamic force measurements. The closed-loop actuation results in two distinct control mechanisms: boat tailing for drag reduction and stagnation point control for lift stabilization. The control strategy proves to be highly effective even in realistic noise scenarios, despite relying on a plant model based on a reduced number of sensors.

There is currently considerable interest in the guided-jet mode, as a result of recent works demonstrating it being the upstream component of various resonant systems in high-speed flows. For given jet operating conditions, the mode is known to exist over only a finite-frequency range that, for a twin-jet system, has been observed to vary with both jet separation and solution symmetry. Vortex-sheet and finite-thickness linear stability models are here employed to consider the behaviour of the guided-jet mode as the two jets are brought together, for both a planar and round twin-jet system. It is demonstrated that in both cases as the twin-jet system merges it forms a higher-order mode of an equivalent single-jet geometry. This then imposes a constraint on the guided-jet mode as the finite-frequency range must change to meet that of the equivalent geometry the system merges to, explaining the previously observed dependence on jet separation.

Amplification of velocity gradients, a key feature of turbulent flows, is affected by the non-local character of the incompressible fluid equations expressed by the second derivative (Hessian) of the pressure field. By analysing the structure of the flow in regions where the vorticity is the highest, we propose an approximate expression for the pressure Hessian in terms of the local vorticity, consistent with the existence of intense vortex tubes. Contrary to the often used simplification of an isotropic form for the pressure Hessian, which in effect inhibits vortex stretching, the proposed approximate form of the pressure Hessian enables much stronger vortex stretching. The prediction of the approximation proposed here is validated with results of direct numerical simulations of turbulent flows.

This study investigates the scenario of flow past a circular cylinder in proximity to a moving wall (or equally a body translating in still fluid parallel to a stationary wall). Fifty high-fidelity three-dimensional direct numerical simulations are performed over a parameter space of turbulent Reynolds numbers (Re) of 300 to 1000 combined with gap-to-diameter ratios (G/D) of 0.2 to 3. The flow, hydrodynamic and turbulence characteristics over the (Re, G/D) parameter space are examined in detail. Small-scale rib-like mode B structures and alternate vortex shedding are observed for all cases. The streamwise location for the vortex shedding (quantified by the wake recirculation length Lr) varies strongly with both Re and G/D. The variation of Lr with Re and G/D can be explained by the spanwise circulation $\varGamma_z$ fed into the wake, where the variation trends of Lr and $\varGamma_z$ are inversely correlated. The variations of the mean drag and fluctuating lift coefficients are also governed by the variations in Lr and $\varGamma_z$. The total kinetic energy in the wake region reduces drastically as G/D reduces below 0.8, which is contributed collectively by (i) reduction in the strength of the shed vortices, (ii) downstream movement of the location of vortex shedding and (iii) associated delayed generation of streamwise vortices. The present results on a moving wall also help to explain several flow and hydrodynamic characteristics reported in the literature for a stationary wall, because the moving wall eliminates the complex wall boundary layer and retains a ‘clean’ near-wall effect.