Viscous fingering, a classic hydrodynamic instability, is governed by the the competition between destabilising viscosity ratios and stabilising surface tension or thermal diffusion. We show that the channel confinement can induce ‘diffusion’-like stabilising effects on viscous fingering even in the absence of interfacial tension and thermal diffusion, when a clear oil invades the mixture of the same oil and non-colloidal particles. The key lies in the generation of long-range dipolar disturbance flows by highly confined particles that form a monolayer inside a Hele-Shaw cell. We develop a coarse-grained model whose results correctly predict universal fingering dynamics that is independent of particle concentrations. This new mechanism offers insights into manipulating and harnessing collective motion in non-equilibrium systems.

The impact of several ‘flavours’ of free-stream turbulence (FST) on the structural response of a cantilever cylinder, subjected to a turbulent cross-flow is investigated. At high enough Reynolds numbers, the cylinder generates a spectrally rich turbulent wake that contributes significantly to the experienced loads. The presence of FST introduces additional complexity through two primary mechanisms: directly, by imposing a fluctuating velocity field on the cylinder’s surface, and indirectly, by altering the vortex shedding dynamics, modifying the experienced loads. We employ concurrent temporally resolved particle image velocimetry and distributed strain measurements using Rayleigh backscattering fibre optic sensors to instrument the surrounding velocity field and the structural strain respectively. By using various turbulence-generating grids, and manipulating their distance to the cylinder, we assess a broad FST parameter space allowing us to explore individually the influence of the transverse integral length scale ($\mathcal{L}_{13}/D$) and turbulence intensity of the FST on the developing load dynamics. The FST enhances the magnitude of the loads acting on the cylinder. This results from a decreased vortex formation length, increased coherence of regular vortex shedding, and energy associated with this flow structure in the near wake. The cylinder’s structural response is driven mainly by the vortex shedding dynamics, and its modification induced by the presence of FST, i.e. the indirect effect outweighs the direct effect. From the explored FST parameter space, turbulence intensity was seen to be the main driver of enhanced loading conditions, presenting a positive correlation with the fluctuating loads magnitude at the root.

Complex materials with internal microstructure such as suspensions and emulsions exhibit time-dependent rheology characterised by viscoelasticity and thixotropy. In many large-scale applications such as turbulent pipe flow, the elastic response occurs on a much shorter time scale than the thixotropy, hence these flows are purely thixotropic. The fundamental dynamics of thixotropic turbulence is poorly understood, particularly the interplay between microstructural state, rheology and turbulence structure. To address this gap, we conduct direct numerical simulations (DNS) of fully developed turbulent pipe flow of a model thixotropic (Moore) fluid as a function of the thixoviscous number $\Lambda$, which characterises the thixotropic kinetic rate relative to turbulence eddy turnover time, ranging from slow ($\Lambda \ll 1$) to fast ($\Lambda \gg 1$) kinetics. Analysis of DNS results in the Lagrangian frame shows that, as expected, in the limits of slow and fast kinetics, these time-dependent flows behave as time-independent purely viscous (generalised Newtonian) analogues. For intermediate kinetics ($\Lambda \sim 1$), the rheology is governed by a path integral of the thixotropic fading memory kernel over the distribution of Lagrangian shear history, the latter of which is modelled via a simple stochastic model for the radially non-stationary pipe flow. The DNS computations based on this effective viscosity closure exhibit excellent agreement with the fully thixotropic model for $\Lambda =1$, indicating that the purely viscous (generalised Newtonian) analogue persists for arbitrary values of $\Lambda \in (0,\infty ^+)$ and across nonlinear rheology models. These results significantly simplify our understanding of turbulent thixotropic flow, and provide insights into the structure of these complex time-dependent flows.

Deformation occurs in a thin liquid film when it is subjected to a non-uniform electric field, which is referred to as the electrohydrodynamic patterning. Due to the development of a non-uniform electrical force along the surface, the film would evolve into microstructures/nanostructures. In this work, a linear and a nonlinear model are proposed to thoroughly investigate the steady state (i.e. equilibrium state) of the electrohydrodynamic deformation of thin liquid film. It is found that the deformation is closely dependent on the electric Bond number BoE. Interestingly, when BoE is larger than a critical value, the film would be deformed remarkably and get in contact with the top template. To model the ‘contact’ between the liquid film and the solid template, the disjoining pressure is incorporated into the numerical model. From the nonlinear numerical model, a hysteresis deformation is revealed, i.e. the film may have different equilibrium states depending on whether the voltage is increased or decreased. To analyse the stability of these multiple equilibrium states, the Lyapunov functional is employed to characterise the system’s free energy. According to the Lyapunov functional analysis, at most three equilibrium states can be formed. Among them, one is stable, another is metastable and the third one is unstable. Finally, the model is extended to study the three-dimensional deformation of the electrohydrodynamic patterning.

This study presents a novel approach for constructing turbulence models using the kinetic Fokker–Planck equation. By leveraging the inherent similarities between Brownian motion and turbulent dynamics, we formulate a Fokker–Planck equation tailored for turbulence at the hydrodynamic level. In this model, turbulent energy plays a role analogous to temperature in molecular thermodynamics, and the large-scale structures are characterised by a turbulent relaxation time. This model aligns with the framework of Pope’s generalised Langevin model, with the first moment recovering the Reynolds-averaged Navier–Stokes (RANS) equations, and the second moment yielding a partially modelled Reynolds stress transport equation. Utilising the Chapman–Enskog expansion, we derive asymptotic solutions for this turbulent Fokker–Planck equation. With an appropriate choice of relaxation time, we obtain a linear eddy viscosity model at first order, and a quadratic Reynolds stress constitutive relationship at second order. Comparative analysis of the coefficients of the quadratic expression with typical nonlinear viscosity models reveals qualitative consistency. To further validate this kinetic-based nonlinear viscosity model, we integrate it as a RANS model within computational fluid dynamics codes, and calculate three typical cases. The results demonstrate that this quadratic eddy viscosity model outperforms the linear model and shows comparability to a cubic model for two-dimensional flows, without the introduction of ad hoc parameters in the Reynolds stress constitutive relationship.

Power minimisation in branched fluidic networks has gained significant attention in biology and engineering. The optimal network is defined by channel radii that minimise the sum of viscous dissipation and the volumetric energetic cost of the fluid. For limit cases including laminar flows, high-Reynolds-number turbulence or smooth-channel approximations, optimal solutions are known. However, current methods do not allow optimisation for a large intermediate part of the parameter space which is typically encountered in realistic fluidic networks that exhibit turbulent flow. Here, we present a unifying optimisation approach based on the Darcy friction factor, which has been determined for a wide range of flow regimes and fluid models and is applicable to the entire parameter space: (i) laminar and turbulent flows, including networks that exhibit both flow types, (ii) non-Newtonian fluids (powerlaw, Bingham and Herschel–Bulkley) and (iii) networks with arbitrary wall roughness, including non-uniform relative roughness. The optimal channel radii are presented analytically and graphically. All existing limit cases are recovered, and a concise framework is presented for systematic optimisation of fluidic networks. Finally, the parameter $x$ in the optimisation relationship $Q\propto R^{x}$, with $Q$ the flow rate and $R$ the channel radius, was approximated as a function of the Reynolds number, revealing in which case the entire network can be optimised based on one optimal channel radius, and in which case all radii must be optimised individually. Our approach can be extended to a wide range of fluidic networks for which the friction factor is known, such as different channel curvatures, bubbly flows or specific wall slip conditions.

Convective boundary layers are governed by an interplay of vertical turbulent convection and shear-driven turbulence. Here, we investigate vertical velocity and buoyancy fields in convective boundary layers for varying atmospheric conditions by combining probability density function methods and direct numerical simulations. The evolution equations for the probability density functions of vertical velocity and buoyancy contain unclosed terms in the form of conditional averages. We estimate these terms from our direct numerical simulations data, and discuss their physical interpretation. Furthermore, using the method of characteristics, we investigate how these unclosed terms jointly determine the average evolution of a fluid element in a convective boundary layer, and how it relates to the evolution of the probability density functions of vertical velocity and buoyancy as a function of height. Thereby, our work establishes a connection between the turbulent dynamics of convective boundary layers and the resulting statistics.

Turbulent mixing driven by the reshocked Richtmyer–Meshkov (RM) instability plays a critical role in numerous natural phenomena and engineering applications. As the most fundamental physical quantity characterizing the mixing process, the mixing width transitions from linear to power-law growth following the initial shock. However, there is a notable absence of quantitative models for predicting the pronounced compression of initial interface perturbations or mixing regions at the moment of shock impact. This gap has restricted the development of integrated algebraic models to only the pre- and post-shock evolution stages. To address this limitation, the present study develops a predictive model for the compression of the mixing width induced by shocks. Based on the general principle of growth rate decomposition proposed by Li et al. (Phy. Rev. E, vol. 103, issue 5, 2021, 053109), two distinct types of shock-induced compression processes are identified, differentiated by the dominant mechanism governing their evolution: light–heavy and heavy–light shock-induced compression. For light–heavy interactions, both stretching (compression) and penetration mechanisms are influential, whereas heavy–light interactions are governed predominantly by the stretching (compression) mechanism. To characterize these mechanisms, the average velocity difference between the extremities of the mixing zone is quantified, and a physical model of RM mixing is utilized. A quantitative theoretical model is subsequently formulated through the independent algebraic modelling of these two mechanisms. The proposed model demonstrates excellent agreement with numerical simulations of reshocked RM mixing, offering valuable insights for the development of integrated algebraic models for mixing width evolution.

Predicting the temperature distribution in laminar two-phase flows is essential in a wide range of engineering applications, like heat dissipation of electronic equipment and thermal design of biological reactors. Motivated by this, we extend the classical Graetz problem, studying the heat transfer between two flowing phases in a core-annular flow configuration. Using a rigorous two-scale asymptotic analysis, we derived two coupled one-dimensional advection–diffusion heat-transfer equations (one for each phase) embedding the effects of advection, diffusion (both axial and transverse) and viscous dissipation. Specifically, the heat-transfer mechanisms are described through effective velocity and effective diffusion coefficients, while the interaction between the phases is accounted for via ad hoc coupling and source terms, respectively. The dynamics of the problem is controlled by seven dimensionless groups: the Péclet and Brinkman numbers, the heat flux, the viscosity, thermal diffusivity and thermal conductivity ratios, and the volume fraction. Our analysis reveals the existence of two main regimes, depending on the disparity in thermal conductivity between the phases. When the conductivity ratio is of order one, the problem is strongly coupled; otherwise, the phases are thermally decoupled. Interestingly, we investigate the evolution of the heat-transfer coefficient in the thin-film limit, shedding light on the most common assumptions underlying extensively used models in the context of film flows. Finally, we derived closed-form scaling laws for the Nusselt number clarifying the impact of the phases topology on heat-transfer dynamics. Since our model has been derived by first principles, we hope that it will improve the understanding of two-phase forced convection.

The dispersion behaviour of solutes in flow is crucial to the design of chemical separation systems and microfluidics devices. These systems often rely on coupled electroosmotic and pressure-driven flows to transport and separate chemical species, making the transient dispersive behaviour of solutes highly relevant. However, previous studies of Taylor dispersion in coupled electroosmotic and pressure-driven flows focused on the long-term dispersive behaviour and the associated analyses cannot capture the transient behaviour of solute. Further, the radial distribution of solute has not been analysed. In the current study, we analyse the Taylor dispersion for coupled electroosmotic and pressure-driven flows across all time regimes, assuming a low zeta potential (electric potential at the shear plane), the Debye–Hückel approximation and a finite electric double layer thickness. We first derive analytical expressions for the effective dispersion coefficient in the long-time regime. We also derive an unsteady, two-dimensional (radial and axial) solute concentration field applicable in the latter regime. We next apply Aris’ method of moments to characterise the unsteady propagation of the mean axial position and the unsteady growth of the variance of the solute zone in all time regimes. We benchmark our predictions with Brownian dynamics simulations across a wide and relevant dynamical regime, including various time scales. Lastly, we derive expressions for the optimal relative magnitudes of electroosmotic versus pressure-driven flow and the optimum Péclet number to minimise dispersion across all time scales. These findings offer valuable insights for the design of chemical separation systems, including the optimisation of capillary electrophoresis devices and electrokinetic microchannels and nanochannels.

The system composed of a circular cylinder free to move along a transverse rectilinear path within a cross-current has often served as a canonical problem to study the vortex-induced vibrations (VIV) developing in the absence of structural restoring force, thus without structural natural frequency. The object of the present work is to extend the exploration of the behaviour of this system when the path is set to an arbitrary orientation, varying from the transverse to the streamwise direction, and the cylinder is forced to rotate about its axis. The investigation is conducted numerically at a Reynolds number equal to $100$, based on the body diameter and oncoming flow velocity, for structure to displaced fluid mass ratios down to $0.01$ and values of the rotation rate (ratio between body surface and oncoming flow velocities) ranging from $0$ to $1$. When the transverse symmetry is broken by the orientation of the trajectory or the forced rotation, the cylinder drifts along the rectilinear path, at a velocity that can be predicted by a quasi-steady approach. Three distinct regimes are encountered: a pure drift regime, where the body translates at a constant velocity, and two oscillatory regimes, characterised by contrasted forms of displacement fluctuation about the drifting motion, but both closely connected to flow unsteadiness. VIV, nearly sinusoidal, persist over a wide range of path orientations, for all rotation rates. On the other hand, irregular jumps of the body, triggered by the rotation and named saccades, emerge when the trajectory is aligned, or almost aligned, with the current. The two forms of response differ by their regularity, but also by their amplitudes and frequencies, which deviate by one or more orders of magnitude. The rotation attenuates both VIV and saccades. Yet, an increase of the rotation rate enhances the erratic nature of the saccade regime.

Spontaneous flow reversals in buoyancy-driven flows are ubiquitous in many fields of science and engineering, often characterized by violent, intermittent occurrences. In this study, we present a complex-network-based reduced-order model to analyse intermittent events in turbulent flows, using temporal and spatial snapshot data. This framework combines elements of dynamical system theory with network science. We demonstrate its utility by applying it to data sequences from intermittent flow reversal events in two-dimensional thermal convection. This approach has proven robust in detecting and quantifying structures and predicting reversals. Additionally, it provides a perspective on the physical mechanisms underlying flow reversals through cluster evolution. This purely data-driven methodology shows the potential to enhance our understanding, prediction and control of turbulent flows and complex systems.

The interaction of a swimmer with unsteady vortices in complex flows remains a topic of interest and open discussion. The present study, employing the immersed boundary method with a flexible fin model, explores swimming behaviours behind a circular cylinder with vortex-induced vibration (VIV). Five distinct swimming modes are identified on the $U_r$–$G_0$ plane, where $U_r$ denotes the reduced velocity, and $G_0$ represents the fin’s initial position. These modes include drifting upstream I/II (DU-I/II), Karman gait I/II (KG-I/II), and large oscillation (LO), with the DU-II, KG-II and LO modes being newly reported. The fin can either move around or cross through the vortex cores in the KG-I and KG-II modes, respectively, for energy saving and maintaining a stable position. When the upstream cylinder vibrates with its maximum amplitude, a double-row vortex shedding forms in the wake, allowing the DU-II mode to occur with the fin to achieve high-speed locomotion. This is attributed to a significant reduction in the streamwise velocity caused by vortex-induced velocity. Furthermore, a symmetry breaking is observed in the fin’s wake in the DU-II mode, potentially also contributing to high-speed locomotion. Overall, compared to the case without an upstream cylinder, we demonstrate that a self-propelled fin gains hydrodynamic advantages with various swimming modes in different VIV wakes. Interestingly, increased power transferred from flows by the oscillating cylinder leads to a more favourable environment for the downstream fin’s propulsion, indicating that a fin in VIV wakes obtains more advantages compared to the vortex street generated by a stationary cylinder.

We perform a comprehensive linear non-modal stability analysis of the Rayleigh–Bénard convection with and without a Poiseuille/Couette flow in Oldroyd-B fluids. In the absence of shear flow, unlike the Newtonian case in which the perturbation energy decays monotonically with time, the interaction between temperature gradient and polymeric stresses can surprisingly cause a transient growth up to 104. This transient growth is maximized at the Hopf bifurcation when the stationary instability dominant in the weakly elastic regime transitions to the oscillatory instability dominant in the strongly elastic regime. In the presence of a Poiseuille/Couette flow, the streamwise-uniform disturbances may achieve the greatest energy amplification, and similar to the pure bounded shear flows, Gmax ∝ Re2 and tmax ∝ Re, where Gmax is the maximum energy growth, tmax the time to attain Gmax, Re the Reynolds number. It is noteworthy that there exist two peaks during the transient energy growth at high-Re cases. Different from the first one which is less affected by the temperature gradient and elasticity, the second peak, at which the disturbance energy is the largest, is simultaneously determined by the temperature gradient, elasticity and shear intensity. Specifically, the polymeric stresses field absorbs energy from the temperature field and base flow, which is partially transferred into the perturbed hydrodynamic field eventually, driving the transient amplification of the perturbed wall-normal vorticity.

The influence of parametric forcing on a viscoelastic fluid layer, in both gravitationally stable and unstable configurations, is investigated via linear stability analysis. When such a layer is vertically oscillated beyond a threshold amplitude, large interface deflections are caused by Faraday instability. Viscosity and elasticity affect the damping rate of momentary disturbances with arbitrary wavelength, thereby altering the threshold and temporal response of this instability. In gravitationally stable configurations, calculations show that increased elasticity can either stabilize or destabilize the viscoelastic system. In weakly elastic liquids, higher elasticity increases damping, raising the threshold for Faraday instability, whereas the opposite is observed in strongly elastic liquids. While oscillatory instability occurs in Newtonian fluids for all gravity levels, we find that parametric forcing below a critical frequency will cause a monotonic instability for viscoelastic systems at microgravity. Importantly, in gravitationally unstable configurations, parametric forcing above this frequency stabilizes viscoelastic fluids, until the occurrence of a second critical frequency. This result contrasts with the case of Newtonian liquids, where under the same conditions, forcing stabilizes a system for all frequencies below a single critical frequency. Analytical expressions are obtained under the assumption of long wavelength disturbances predicting the damping rate of momentary disturbances as well as the range of parameters that lead to a monotonic response under parametric forcing.

The attached-eddy model (AEM) predicts that the mean streamwise velocity and streamwise velocity variance profiles follow a logarithmic shape, while the vertical velocity variance remains invariant with height in the overlap region of high Reynolds number wall-bounded turbulent flows. Moreover, the AEM coefficients are presumed to attain asymptotically constant values at very high Reynolds numbers. Here, the AEM predictions are examined using sonic anemometer measurements in the near-neutral atmospheric surface layer, with a focus on the logarithmic behaviour of the streamwise velocity variance. Utilizing an extensive 210-day dataset collected from a 62 m meteorological tower located in the Eastern Snake River Plain, Idaho, USA, the inertial sublayer is first identified by analysing the measured momentum flux and mean velocity profiles. The logarithmic behaviour of the streamwise velocity variance and the associated ‘$-1$’ scaling of the streamwise velocity energy spectra are then investigated. The findings indicate that the Townsend–Perry coefficient ($A_1$) is influenced by mild non-stationarity that manifests itself as a Reynolds number dependence. After excluding non-stationary runs, and requiring the bulk Reynolds number defined using the atmospheric boundary layer height to be larger than $4 \times 10^{7}$, the inferred $A_1$ converges to values ranging between 1 and 1.25, consistent with laboratory experiments. Furthermore, nine benchmark cases selected through a restrictive quality control reveal a close relation between the ‘$-1$’ scaling in the streamwise velocity energy spectrum and the logarithmic behaviour of streamwise velocity variance. However, additional data are required to determine whether the plateau value of the pre-multiplied streamwise velocity energy spectrum is identical to $A_1$.

We explore the instability and oscillation dynamics of barrel-shaped droplets on cylindrical fibres, contributing to a deeper understanding of fibre–droplet interactions critical to both natural systems and industrial applications. Unlike sessile droplets on flat surfaces, droplets on fibres exhibit unique behaviours due to the curvature of the fibre, such as transitions from axisymmetric (barrel) to non-axisymmetric (clamshell) shapes governed by droplet volume, contact angle and fibre radius. Using a linear inviscid theory, we compute the frequency spectrum of barrel-shaped droplets and identify stability thresholds for the barrel-to-clamshell transition by examining the first rocking mode, with a focus on the role of contact line conditions. This analysis resolves experimental anomalies concerning the stability of half-barrel-shaped droplets on hydrophobic fibres. Our findings also reveals diverse frequency spectra: droplets on thin fibres exhibit Rayleigh–Lamb-like spectral features, while those on thicker fibres show reduced sensitivity to azimuthal wavenumber. Interestingly, the instability of sectoral modes on thick fibres resembles the Rayleigh–Plateau instability of static rivulets, with fibre curvature slightly reducing growth rates at small axial wavenumbers but increasing them at larger ones.

The recently proposed near-wall turbulence predictive model quantifies the degree of the superposition and the amplitude modulation exerted by large-scale coherent structures on small scales in the linear and nonlinear terms of the formula, respectively, and achieves the prediction of streamwise velocity in the inner region. However, the multiscale effect and the time shift confirmed in the amplitude modulation have not yet been simultaneously taken into account in the model, which could limit the prediction accuracy especially at high Reynolds numbers. In this study, the role of the nonlinear term in the model is clarified based on high-quality flow data obtained in atmospheric surface layers: it redistributes the energy of the universal signal in the time domain and determines the accuracy of the predictive odd moments. An analysis of the multiscale effect and the time shifts in the nonlinear term is subsequently conducted, followed by a demonstration of the refinement in the quality of the universal signal after separately incorporating them into the model. The amplitude modulation is revealed when the two factors are simultaneously considered, and profiles of the scales that dominate the modulation and time shifts with height is provided. Thus, the nonlinear term of the existing model is modified, proposing an polished scheme that can quantify the nonlinear modulation terms more accurately.

The instability characteristics and laminar–turbulent transition of a series of laminar separation bubbles (LSBs) formed due to a single sinusoidal surface waviness are investigated in the absence of external disturbances or forcing. A scaling based on the geometrical parameters of the waviness and flow Reynolds number is found that enables the prediction of flow separation on the wall leeward side. The analysis of three-dimensional instabilities of two-dimensional base flows reveals a relation between the number of changes in the curvature sign of the recirculating streamlines and the number of unstable centrifugal modes that coexist for the same flow. When multiple curvature changes occur, in addition to the usual steady mode reported for two-dimensional recirculation bubbles, a new self-excited mode with a higher growth rate emerges, localised near the highest streamline curvature, close to the reattachment point. A detailed analysis of the mode growth and saturation using DNS reveals that the localised mode only disturbs the LSB locally, while the usual one leads to a global distortion of the bubble in the spanwise direction; this has a distinctive impact on the self-excited secondary instabilities. Then, the complete transition scenario is studied for two selected LSB cases. The first one only presents an unstable eigenmode, namely the usual centrifugal mode in recirculating flows. The second case presents three unstable eigenmodes: two centrifugal eigenmodes (the usual and the localised ones) and a two-dimensional eigenmode associated with the self-sustained Kelvin–Helmholtz waves. These results show how completely different transition scenarios can emerge from subtle changes in the LSB characteristics.

This paper investigates the linear and nonlinear dynamics of two-dimensional penetrative convection subjected to radiative volumetric thermal forcing, focusing on ice-covered freshwater systems. Linear stability analysis reveals how critical wavenumbers $k_c$ and Rayleigh numbers $Ra_c$ are influenced by the attenuation lengths and incoming heat flux. In this configuration, the system easily becomes unstable with a small $Ra_c$, which is two decades smaller than that of the classical Rayleigh–Bénard convection problem, with typically $O(10)$. Weakly nonlinear analysis figures out that this configuration is supercritical, contrasting with the subcritical case by Veronis (Astrophys. J., vol. 137, 1963, 641–663). Numerical bifurcation solutions are performed from the critical points and over several decades, up to $Ra \sim O(10^6)$. This paper found that the system exhibits multiple steady solutions, and under certain specific conditions, a staircase temperature profile emerges. Meanwhile, we further discuss the influence of incoming heat flux and the Prandtl number $Pr$ on the primary bifurcation. Direct numerical simulations are also carried out, showing that heat is transported more efficiently via unsteady convection.