Laminar–turbulent transition in shear flow is complicated and follows many possible routes. In this study, we seek to examine a scenario based on three-dimensional (3-D) waves (Jiang et al., 2020, J. Fluid Mech., vol. 890, A11) in compressible mixing layers, and elucidate the role of 3-D waves in generating streamwise vorticity. The Eulerian–Lagrangian coupled method is used to track the evolution of flow structures. Qualitative evidence shows that localised 3-D waves travel coherently with vortex structures at the early transition stage, which is consistent with the behaviours of 3-D waves in boundary layer transitions. To examine the local flow events surrounding 3-D waves and investigate the cause and effect relationships inherent in wave–vortex interaction, the finite-time Lyapunov exponent and components of the strain rate tensor are integrated into evolving Lagrangian material surfaces. The formation of high-shear layers in the flanks of the 3-D waves is observed, driven by fluid ejection and sweep motions induced by the amplification of 3-D waves. The $\Lambda$-shaped vortices are found born in the vicinity of high-shear regions and then stretched into hairpin-shaped vortices farther downstream. Statistical findings reveal that streamwise vorticity develops concurrently with the significant growth of the oblique mode, while the normal motion of wave structures induces a high strain rate layer in the surrounding region. In addition, conditional statistics underscore the significance of high shear in enstrophy generation. Finally, a conceptual model is proposed to depict the evolution of coherent structures based on the relationship among the 3-D waves, high-shear/strain layers, and $\varLambda$-vortices, providing insights into their collective dynamics within transitional mixing layers.

Finite-amplitude spiral vortex flows are obtained numerically for the Taylor–Couette system in the narrow limit of the gap between two concentric rotating cylinders. These spiral vortex flows bifurcate from circular Couette flow before axisymmetric Taylor vortex flow sets in when the ratio $\mu$ of the angular velocities of the outer to the inner cylinder is less than −0.78, consistent with the results of linear stability analysis by Krueger et al. (J. Fluid Mech., vol. 24, 1966, pp. 521–538), while the boundary of existence of spiral vortex flows is determined not by the linear critical point, but by the saddle-node point of the subcritical spiral vortex flow branch for $\mu \lessapprox -0.75$, when the axial wavenumber $\beta =2.0$. It is found that the nonlinear spiral vortex flows exhibit the mean flow in the axial direction as well as in the azimuthal direction, and that the profiles of both mean-flow components are asymmetric about the centre plane between the gap.

In this work the fascinating dynamics of a two-layered channel flow characterised by the dispersion in composite media within its layers is investigated in depth. The top layer comprises of a fluid zone that allows the fluid to travel along its surface easily (with relatively higher velocity), while the bottom layer is packed with porous media. The primary objective of this research is to do an in-depth investigation of the complex two-dimensional concentration distribution of a passive solute discharged from the inflow region. A multi-scale perturbation analysis approach has been implemented to address the system’s inherent complexity. This accurate determination of the dispersion coefficient, mean concentration distribution and two-dimensional concentration distribution is accomplished deftly using Mei’s homogenisation approach up to second-order approximation, which satisfactorily capture the minor variations in the solute dynamics also. The influence of various flow and porous media elements on these basic parameters is thoroughly investigated, expanding our comprehension of the complex interaction between flow dynamics and porous media’s properties. The effect of Darcy number and the ratio of two viscosities ($M$) on the dispersion coefficient depends on the height of the porous layer. As the Péclet number ratio increases, the dispersion coefficient experiences a concurrent increase, resulting in a decline in the concentration peak. The results of the analytical studies have also been compared with those results obtained using a purely computational method to establish the validity of our studies. Both the sets of results show quite good agreement with each other. In this study, alternate flow models have been used for the porous region, and the outcomes are compared to determine which approach yields more suitable results under different conditions.

The paper discusses the stochastic dynamics of the vortex shedding process in the presence of external harmonic excitation and coloured multiplicative noise. The situation is encountered in a turbulent practical combustor experiencing combustion instability. Acoustic feedback and turbulent flow are imitated by the harmonic and stochastic excitations, respectively. The Ornstein–Uhlenbeck process is used to generate the noise. A low-order model for vortex shedding is used. The Fokker–Planck framework is used to obtain the evolution of the probability density function of the shedding time period. Stochastic lock-in and resonance characteristics are studied for various parameters associated with the harmonic (amplitude, frequency) and noise (amplitude, correlation time, multiplicative noise factor) excitations. We observed that: (i) the stochastic lock-in (s-lock-in) boundary strongly depends on the noise correlation time; (ii) the parameter sites for s-lock-in can be approximately identified from the noise-induced shedding statistics; and (iii) stochastic resonance is significant for some intermediate correlation times. The effects of the above-mentioned observations are discussed in the context of combustion instability.

A deep-learning-based closure model to address energy loss in low-dimensional surrogate models based on proper-orthogonal-decomposition (POD) modes is introduced. Using a transformer-encoder block with an easy-attention mechanism, the model predicts the spatial probability density function of fluctuations not captured by the truncated POD modes. The methodology is demonstrated on the wake of the Windsor body at yaw angles of $\delta = [2.5^\circ ,5^\circ ,7.5^\circ ,10^\circ ,12.5^\circ ]$, with $\delta = 7.5^\circ$ as a test case, and in a realistic urban environment at wind directions of $\delta = [-45^\circ ,-22.5^\circ ,0^\circ ,22.5^\circ ,45^\circ ]$, with $\delta = 0^\circ$ as a test case. Key coherent modes are identified by clustering them based on dominant frequency dynamics using Hotelling’s $T^2$ on the spectral properties of temporal coefficients. These coherent modes account for nearly $60 \,\%$ and $75 \,\%$ of the total energy for the Windsor body and the urban environment, respectively. For each case, a common POD basis is created by concatenating coherent modes from training angles and orthonormalising the set without losing information. Transformers with different size on the attention layer, (64, 128 and 256), are trained to model the missing fluctuations in the Windsor body case. Larger attention sizes always improve predictions for the training set, but the transformer with an attention layer of size 256 slightly overshoots the fluctuation predictions in the Windsor body test set because they have lower intensity than in the training cases. A single transformer with an attention size of 256 is trained for the urban flow. In both cases, adding the predicted fluctuations close the energy gap between the reconstruction and the original flow field, improving predictions for energy, root-mean-square velocity fluctuations and instantaneous flow fields. For instance, in the Windsor body case, the deepest architecture reduces the mean energy error from $37 \,\%$ to $12 \,\%$ and decreases the Kullback–Leibler divergence of velocity distributions from ${\mathcal{D}}_{\mathcal{KL}}=0.2$ to below ${\mathcal{D}}_{\mathcal{KL}}=0.026$.

This paper explores decaying turbulence beneath surface waves that is initially isotropic and shear free. We start by presenting phenomenology revealed by wave-averaged numerical simulations: an accumulation of angular momentum in coherent vortices perpendicular to the direction of wave propagation, suppression of kinetic energy dissipation and the development of depth-alternating jets. We interpret these features through an analogy with rotating turbulence (Holm 1996 Physica D. 98, 415–441), wherein the curl of the Stokes drift, ${\boldsymbol{\nabla}} \times {\boldsymbol{u^{S}}}$, takes on the role of the background vorticity (for example, $(f_0 + \beta y) {\boldsymbol{\hat{z}}}$ on the beta plane). We pursue this thread further by showing that a two-equation model proposed by Bardina et al. (1985 J. Fluid Mech. 154, 321–336) for rotating turbulence reproduces the simulated evolution of volume-integrated kinetic energy. This success of the two-equation model – which explicitly parametrises wave-driven suppression of kinetic energy dissipation – carries implications for modelling turbulent mixing in the ocean surface boundary layer. We conclude with a discussion about a wave-averaged analogue of the Rossby number appearing in the two-equation model, which we term the ‘pseudovorticity number’ after the pseudovorticity ${\boldsymbol{\nabla }} \times {\boldsymbol{u}}^S$. The pseudovorticity number is related to the Langmuir number in an integral sense.

Compliant walls made from homogeneous viscoelastic materials may attenuate the amplification of Tollmien–Schlichting waves (TSWs) in a two-dimensional boundary-layer flow, but they also amplify travelling-wave flutter (TWF) instabilities at the interface between the fluid and the solid, which may lead to a premature laminar-to-turbulent transition. To mitigate the detrimental amplification of TWF, we propose to design compliant surfaces using phononic structures that aim at avoiding the propagation of elastic waves in the solid in the frequency range corresponding to the TWF. Thus, stiff inserts are periodically incorporated into the viscoelastic wall in order to create a band gap in the frequency spectrum of the purely solid modes. Fluid–structural resolvent analysis shows that a significant reduction in the amplification peak related to TWF is achieved while only marginal deterioration in the control of TSWs is observed. This observation suggests that the control of TSWs is still achieved by the overall compliance of the wall, while the periodic inserts inhibit the amplification of TWF. Bloch analysis is employed to discuss the propagation of elastic waves in the phononic surface to deduce design principles, accounting for the interaction with the flow.

Large numbers of relative periodic orbits (RPOs) have been found recently in doubly periodic, two-dimensional Kolmogorov flow at moderate Reynolds numbers ${\textit{Re}} \in \{40, 100\}$. While these solutions lead to robust statistical reconstructions at the ${\textit{Re}}$ values where they were obtained, it is unclear how their dynamical importance changes with ${\textit{Re}}$. Arclength continuation on this library of solutions reveals that large numbers of RPOs quickly become dynamically irrelevant, reaching dissipation values either much larger or smaller than the values typical of the turbulent attractor at high ${\textit{Re}}$. The scaling of the high-dissipation RPOs is shown to be consistent with a direct connection to solutions of the unforced Euler equation, and is observed for a wide variety of states beyond the ‘unimodal’ solutions considered in previous work (Kim & Okamoto, Nonlinearity vol. 28, 2015, p. 3219). However, the weakly dissipative states have properties indicating a connection to exact solutions of a forced Euler equation. The dynamical irrelevance of many solutions leads to poor statistical reconstruction at higher ${\textit{Re}}$, raising serious questions for the future use of RPOs for estimating probability densities. Motivated by the Euler connection of some of our RPOs, we also show that many of these states can be well described by exact relative periodic solutions in a system of point vortices. The point vortex RPOs are converged via gradient-based optimisation of a scalar loss function which (i) matches the dynamics of the point vortices to the turbulent vortex cores and (ii) insists the point vortex evolution is itself time-periodic.

In the present study, we observe interesting profiles and fluctuations in a quasi-two-dimensional thermal convection system filled with low-Prandtl-number liquid metal. A high-precision thermistor, which can be precisely controlled to move up and down, is used to measure the temperature distribution along the centreline of a convection cell. As the thermistor probes move away from the heated wall surface, the measured temperatures initially decrease to values below the central temperature of the cell, then recover to the central temperature, indicating an inverse temperature gradient. Furthermore, by analysing the root-mean-square temperature ($\sigma _T (z)$) along the centreline, we find a second peak away from the wall location, which has never been reported before, in addition to the first peak associated with the thermal boundary thickness. This phenomenon is also confirmed by the results of third- and fourth-order moments of temperature. Experimental results, together with insights from previous studies, suggest that in liquid metal, the distinct flow organisation arising from the large thermal diffusivity plays an important role in shaping the observed temperature distribution.

Pendant drops appear in many engineering applications, such as inkjet printing and optical tensiometry, and they have also been the subject of studies of droplet–particle interaction. While the hydrostatics of pendant drops has been studied extensively, the influence of external flow disturbances has received limited attention. This research aims to incorporate aerodynamic factors into the understanding of pendant drop behaviour. Employing a simplified model, an irrotational flow aligned with the drop’s axis is derived from a distribution of singularity elements within the drop. The drop’s equilibrium shape is then determined using a numerical model that couples the flow field with the Young–Laplace equation. The model’s predictions are compared to droplet images captured via high-speed shadowgraph in a vertical wind tunnel, showing good agreement with the experimentally observed shapes. Additionally, under certain flow conditions, the drop exhibits instability in the form of periodic pendulum-like motion. This instability was linked to two distinct critical drop heights, and the corresponding stability criterion was mathematically derived from the numerical model. Our theoretical and experimental findings provide the first quantitative description of the equilibrium shape and stability criterion of pendant drops under the influence of external flow.

Experiments have shown that ultrasound-stimulated microbubbles can translate through gel phantoms and tissues, leaving behind tunnel-like degraded regions. A computational model is used to examine the tunnelling mechanisms in a model material with well-defined properties. The high strain rates motivate the neglect of weak elasticity in favour of viscosity, which is taken to degrade above a strain threshold. The reference parameters are motivated by a 1 $\unicode{x03BC}$m diameter bubble in a polysaccharide gel tissue phantom. This is a reduced model and data are scarce, so close quantitative agreement is not expected, but tunnels matching observations do form at realistic rates, which provides validation sufficient to analyse potential mechanisms. Simulations of up to 100 acoustic cycles are used to track tunnelling over 10 bubble diameters, including a steady tunnelling phase during which tunnels extend each forcing cycle in two steps: strain degrades the tunnel front during the bubble expansion, and then the bubble is drawn further along the tunnel during its subsequent inertial collapse. Bubble collapse jetting is damaging, though it is only observed during a transient for some initial conditions. There is a threshold behaviour when the viscosity of the undamaged material changes the character of the inertial bubble oscillation. Apart from that, the tunnel growth rate is relatively insensitive to the high viscosity of the material. Higher excitation amplitudes and lower frequencies accelerate tunnelling. That acoustic radiation force, elasticity and bubble jetting are not required is a principal conclusion.

An asymptotic model for the flow of a highly viscous film coating the interior of a slippery, flexible tube is developed and studied. The model is valid for the axisymmetric flow of moderately thick films, and accounts for tube flexibility, wall damping, longitudinal tension, slip length and strength of base flow due either to gravity or airflow. In the absence of base flow, linear stability analysis shows the existence of one unstable mode; the presence of base flow allows for multiple unstable modes arising due to the Plateau–Rayleigh instability and elastic instability, with stronger base flow reducing the maximum growth rate. Numerical solutions in the absence of base flow show that slip decreases the amplitude of wall deformations and can significantly decrease the time to plug formation in weakly flexible or strongly damped tubes. For falling films, the impact of model parameters on the critical thickness required for plug formation was analysed by studying turning points in families of travelling-wave solutions; this thickness decreases with slip, flexibility and tension, while damping had a non-monotonic impact on critical thickness. In contrast to model solutions in rigid tubes, for flexible tubes the critical thickness cannot be made arbitrarily large through simply increasing the strength of the base flow. For air-driven films, both slip and flexibility increase the rate of film transport along the tube.

We investigate theoretically the breakup dynamics of an elasto-visco-plastic filament surrounded by an inert gas. The filament is initially placed between two coaxial disks, and the upper disk is suddenly pulled away, inducing deformation due to both constant stretching and capillary forces. We model the rheological response of the material with the Saramito–Herschel–Bulkley (SHB) model. Assuming axial symmetry, the mass and momentum balance equations, along with the constitutive equation, are solved using the finite element framework PEGAFEM-V, enhanced with adaptive mesh refinement with an underlying elliptic mesh generation algorithm. As the minimum radius decreases, the breakup dynamics accelerates significantly. We demonstrate that the evolution of the minimum radius, velocity and axial stress follow a power-law scaling, with the corresponding exponent depending on the SHB shear-thinning parameter, $n$. The scaling exponents obtained from our axisymmetric simulations under creeping flow are verified through asymptotic analysis of the slender filament equations. Our findings reveal three distinct breakup regimes: (a) elasto-plastic, (b) elasto-plasto-capillary, both with finite-time breakup for $n\lt 1$, and (c) elasto-plasto-capillary with no finite-time breakup for $n=1$. We show that self-similar solutions close to filament breakup can be achieved by appropriate rescaling of length, velocity and stress. Notably, the effect of the yield stress becomes negligible in the late stages of breakup due to the local dominance of high elastic stresses. Moreover, the scaling exponents are independent of elasticity, resembling the breakup behaviour of finite extensible viscoelastic materials.

In air-entraining flows, there is often strong turbulence beneath the free surface. We consider the entrainment of bubbles at the free surface by this strong free-surface turbulence (FST). Our interest is the entrainment size distribution (per unit free surface area) $I(a)/A_{\textit{FS}}$, for bubbles with radius $a$ greater than the capillary scale ($\approx 1.3\ \mathrm{mm}$ for air–water on Earth), where gravity dominates surface tension. We develop a mechanistic model based on entrained bubble size being proportional to the minimum radius of curvature of the initial surface deformation. Using direct numerical simulation of a flow that isolates entrainment by FST, we show that, consistent with our mechanism, $I(a)/A_{\textit{FS}} = C_I \, g^{-3} \varepsilon ^{7/3} (2 a)^{-14/3}$, where $g$ is gravity, and $\varepsilon$ is the turbulence dissipation rate. In the limit of negligible surface tension, $C_I\approx 3.62$, and we describe how $C_I$ decreases with increasing surface tension. This scaling holds for sufficiently strong FST such that near-surface turbulence is nearly isotropic, which we show is true for turbulent Froude number ${\textit{Fr}}^2_T = \varepsilon /u_{\textit{rms}} g \gt 0.1$. While we study FST entrainment in isolation, our model corroborates previous numerical results from shear-driven flow, and experimental results from open-channel flow, showing that the FST entrainment mechanism that we elucidate can be important in broad classes of air-entraining flows.

The propagation of linear waves in non-ideal compressible fluids plays a crucial role in numerous physical and engineering applications, particularly in the study of instabilities, aeroacoustics and turbulence modelling. This work investigates linear waves in viscous and heat-conducting non-ideal compressible fluids, modelled by the Navier–Stokes–Fourier equations and a fully arbitrary equation of state (EOS). The linearised governing equations are derived to analyse the dispersion relations when the EOS differs from that of an ideal gas. Special attention is given to the influence of non-ideal effects and various dimensionless numbers on wave propagation speed and attenuation. By extending classical results from Kovásznay (1953 J. Aeronaut. Sci. vol. 20, no. 10, pp. 657–674) and Chu (1965 Acta Mech. vol. 1, no. 3, pp. 215–234) obtained under the ideal gas assumption, this study highlights the modifications introduced by arbitrary EOSs to the linear wave dynamics in non-ideal compressible flows. This work paves the path for an improved understanding and modelling of wave propagation, turbulence and linear stability in arbitrary viscous and heat-conducting fluids.