Vortex rings are critical for thrust production underwater. In the ocean, self-propelled mesozooplankton generate vortices while swimming within a weakly stratified fluid. While large-scale biogenic transport has been observed during vertical migration in the wild and lab experiments, little focus has been given to the evolution of induced vortex rings as a function of their propagation direction relative to the density gradient. In this study, the evolution of an isolated vortex ring crossing the interface of a stable two-layer system is examined as a function of its translation direction with respect to gravity. The vortex ring size and position are visualized using planar laser-induced fluorescence (PLIF) and the induced vorticity field derived from particle image velocimetry (PIV) is examined. It is found that the production of baroclinic vorticity significantly affects the propagation of vortex rings crossing the density interface. As a result, any expected symmetry between vortex rings travelling from dense to light fluids and from light to dense fluids breaks down. In turn, the maximum penetration depth of the vortex ring occurs in the case in which the vortex propagates against the density gradient due to the misalignment of the pressure and density gradients. Our results have far-reaching implications for the characterization of local ecosystems in marine environments.

To explore the direction of inter-scale transfer of scalar variance between subgrid scale (SGS) and resolved scalar fields, direct numerical simulation data obtained earlier from two complex-chemistry lean hydrogen–air flames are analysed by applying Helmholtz–Hodge decomposition (HHD) to the simulated velocity fields. Computed results show backscatter of scalar (combustion progress variable $c$) variance, i.e. its transfer from SGS to resolved scales, even in a highly turbulent flame characterized by a unity-order Damköhler number and a ratio of Kolmogorov length scale to thermal laminar flame thickness as low as 0.05. Analysis of scalar fluxes associated with the solenoidal and potential velocity fields yielded by HHD shows that the documented backscatter stems primarily from the potential velocity perturbations generated due to dilatation in instantaneous local flames, with the backscatter being substantially promoted by a close alignment of the spatial gradient of mean scalar progress variable and the potential-velocity contribution to the local SGS scalar flux. The alignment is associated with the fact that combustion-induced thermal expansion increases local velocity in the direction of $\boldsymbol {\nabla } c$. These results call for development of SGS models capable of predicting backscatter of scalar variance in turbulent flames in large eddy simulations.

Coupling surface deformations with active stresses in two-dimensional nematic liquid crystal films leads to a rich area of investigation, particularly in biological fluid mechanics across multiple scales from tissue mechanics to cell membrane mechanics. In Al-Izzi & Morris (J. Fluid Mech., vol. 957, 2023, A4), the authors derive the complete set of governing equations for such systems. Their results provide an extended theoretical framework with which active nematic fluid films with in-plane flow and out-of-plane deformation can be analysed. To illustrate the potential applications of this framework, a few specific biologically inspired examples are discussed.

We conduct three-dimensional direct numerical simulations to investigate the mixing, entrainment and energy budgets of gravity currents emerging from two-layer stratified locks. Depending on the density and layer thickness ratios, we find that either the upper layer or lower layer fluid can propagate faster, and that the density structure of the overall gravity current can range from strongly stratified to near-complete mixing. We furthermore observe that intermediate values of the density ratio can maximise mixing between the gravity current layers. Based on the vorticity budget, we propose a theoretical model for predicting the overall gravity current height, along with the front velocity of the two layers, for situations in which the lower layer moves faster than the upper layer. The model identifies the role of the height and thickness ratios in determining the velocity structure of the current, and it clarifies the dynamics of the ambient counter-current. A detailed analysis of the energy budget quantifies the conversion of potential into kinetic energy as a function of the governing parameters, along with the energy transfer between the different layers of the gravity current and the ambient fluid. Depending on the values of the density and layer thickness ratios, we find that the lower lock layer can gain or lose energy, whereas the upper layer always loses energy.

We studied the flow organization and heat transfer properties in two-dimensional and three-dimensional Rayleigh–Bénard cells that are imposed with different types of wall shear. The external wall shear is added with the motivation of manipulating flow mode to control heat transfer efficiency. We imposed three types of wall shear that may facilitate the single-roll, the horizontally stacked double-roll, and the vertically stacked double-roll flow modes, respectively. Direct numerical simulations are performed for fixed Rayleigh number $Ra = 10^{8}$ and fixed Prandtl number $Pr = 5.3$, while the wall-shear Reynolds number ($Re_{w}$) is in the range $60 \leqslant Re_{w} \leqslant 6000$. Generally, we found enhanced heat transfer efficiency and global flow strength with the increase of $Re_{w}$. However, even with the same magnitude of global flow strength, the heat transfer efficiency varies significantly when the cells are under different types of wall shear. An interesting finding is that by increasing the wall-shear strength, the thermal turbulence is relaminarized, and more surprisingly, the heat transfer efficiency in the laminar state is higher than that in the turbulent state. We found that the enhanced heat transfer efficiency at the laminar regime is due to the formation of more stable and stronger convection channels. We propose that the origin of thermal turbulence laminarization is the reduced amount of thermal plumes. Because plumes are mainly responsible for turbulent kinetic energy production, when the detached plumes are swept away by the wall shear, the reduced number of plumes leads to weaker turbulent kinetic energy production. We also quantify the efficiency of facilitating heat transport via external shearing, and find that for larger $Re_{w}$, the enhanced heat transfer efficiency comes at a price of a larger expenditure of mechanical energy.

Large-eddy simulation is used to simulate the flow around a 6 : 1 prolate spheroid at $10^\circ$ and $20^\circ$ angles of attack, and Reynolds number $4.2 \times 10^6$. Flows with and without trip are compared to understand the relative effect of the trip on the state of the boundary layer and separation. For the tripped case, the geometry of the trip is resolved to better predict its effect on the downstream flow. The simulations employ overset grids that allow adequate resolution of the trips without significant increase in the overall computational cost. Results suggest that while the trip accelerates transition to turbulence at $10^\circ$, it does not induce a fully developed turbulent boundary layer as intended at $20^\circ$. Rather, the influence of the trip is localized, and the near-wall flow converges towards a solution similar to that of the non-tripped case upstream of separation. This is due to two distinct phenomena: directly downstream of the trip, favourable pressure gradient and streamline curvature effects suppress the disturbance on the windward side. Further along the spheroid, the boundary layer receives a small fraction of the initial perturbation due to spanwise and wall-normal streamline curvatures inducing a secondary flow that advects the low-momentum trip wake to the leeward side. The locations of transition and separation are insensitive to the presence of the trip. The simulation results are used to construct a regime map that identifies different regions characterized by distinct boundary layer properties and flow features. The present results underscore the difficulty associated with tripping smooth bodies at angle of attack, and the importance of accounting for transition in simulations of such flows, even on tripped geometries.

In this paper, we propose artificial-neural-network-based (ANN-based) nonlinear algebraic models for the large-eddy simulation (LES) of compressible wall-bounded turbulence. An innovative modification is applied to the invariants and the tensor bases of the nonlinear algebraic models through using the local grid widths along each direction to normalise the corresponding gradients of the flow variables. Furthermore, the dimensionless model coefficients are determined by the ANN method. The modified ANN-based nonlinear algebraic model (MANA model) has much higher correlation coefficients and much lower relative errors than the dynamic Smagorinsky model (DSM), Vreman model and wall-adapting local eddy-viscosity model in the a priori test. The significantly more accurate estimations of the mean subgrid-scale (SGS) fluxes of the kinetic energy and temperature variance are also obtained by the MANA models in the a priori test. Furthermore, in the a posteriori test, the MANA model can give much more accurate predictions of the flow statistics and the mean SGS fluxes of the kinetic energy and the temperature variance than other traditional eddy-viscosity models in compressible turbulent channel flows with untrained Reynolds numbers, Mach numbers and grid resolutions. The MANA model has a better performance in predicting the flow statistics in supersonic turbulent boundary layer. The MANA model can well predict both direct and inverse transfer of the kinetic energy and temperature variance, which overcomes the inherent shortcoming that the traditional eddy-viscosity models cannot predict the inverse energy transfer. Moreover, the MANA model is computationally more efficient than the DSM.

We study the dynamics of capillary waves at the interface of a two-layer stratified turbulent channel flow. We use a combined pseudo-spectral/phase field method to solve for the turbulent flow in the two liquid layers and to track the dynamics of the liquid–liquid interface. The two liquid layers have same thickness and same density, but different viscosity. We vary the viscosity of the upper layer (two different values) to mimic a stratified oil–water flow. This allows us to study the interplay between inertial, viscous and surface tension forces in the absence of gravity. In the present set-up, waves are naturally forced by turbulence over a broad range of scales, from the larger scales, whose size is of order of the system scale, down to the smaller dissipative scales. After an initial transient, we observe the emergence of a stationary capillary wave regime, which we study by means of temporal and spatial spectra. The computed frequency and wavenumber power spectra of wave elevation are in line with previous experimental findings and can be explained in the frame of the weak wave turbulence theory. Finally, we show that the dispersion relation, which gives the frequency ($\omega$) as a function of the wavenumber ($k$), is in good agreement with the well-established theoretical prediction, $\omega (k) \sim k^{3/2}$.

The dissolution of a horizontal salt body is theoretically and numerically investigated by considering the coupled dynamics of the dissolution reaction and the diffusive and convective mass transfer at the salt surface. To take into account the nature of the recessing salt surface and the dissolution-driven convection, a dynamic moving boundary condition coupled with a dissolution reaction are employed By employing a newly derived moving interface condition and parameters, a theoretical analysis predicts the onset of gravitational instability, suggesting scaling relationships between parameters to describe the onset of instability (in transport-dominant (Damkholer number $Da \to \infty $ and Rayleigh number $Ra > {10^5}$), onset time ${\tau _d}\sim R{a^{ - 2/3}}$ and ${\tau _d}\sim (RaDa)^{1/4}$ for reaction-dominant regimes $(Da \to 0)$). The scaling relationships on dissolution (the average height change, $\Delta {h_{avg}}\sim \sqrt \tau $ for the diffusion-dominant regime and $\Delta {h_{avg}}\sim R{a^{0.79/3}}{\tau ^{0.86}}$ for the convection-dominant regime) are also suggested to explain the pattern formation. Under a convective mass transfer condition, the concentration gradient of the solute developed during the recession of the salt surface results in an adverse density gradient, which causes gravitational instability motion and finally determines the formation of a dissolution pattern. In addition, two-dimensional and three-dimensional numerical simulations visualize dissolution-driven convective flow fields, the recession of the liquid-solid interface and pattern formation on the dissolving interface. Theoretical and numerical analyses of the dissolution process suggest that the control of gravitational instability is important to prevent or enhance dissolution pattern formation. This systematic parameter study provides a deep understanding of the effect of gravitational instability on convection-driven dissolution in a horizontal geometry.

Exact budget equations are derived for the coherent and stochastic contributions to the second-order structure function tensor. They extend the anisotropic generalised Kolmogorov equations (AGKE) by considering the coherent and stochastic parts of the Reynolds stress tensor, and are useful for the statistical description of turbulent flows with periodic or quasi-periodic features, like, for example, the alternate shedding after a bluff body. While the original AGKE describe production, transport, inter-component redistribution and dissipation of the Reynolds stresses in the combined space of scales and positions, the new equations, called $\varphi$AGKE, contain the phase $\varphi$ as an additional independent variable, and describe the interplay among the mean, coherent and stochastic fields at the various phases. The newly derived $\varphi$AGKE are then applied to a case where an exactly periodic external forcing drives the flow: a turbulent plane channel flow modified by harmonic spanwise oscillations of the wall to reduce drag. The phase-by-phase action of the oscillating transversal Stokes layer generated by the forcing on the near-wall turbulent structures is observed, and a detailed description of the scale-space interaction among mean, coherent and stochastic fields is provided thanks to the $\varphi$AGKE.

The instability mechanism behind a geometrically simple cylindrical roughness element continues to be a challenging topic in fluid mechanics. Considerable progress has been made in understanding the phenomena in recent years, but more research is needed to predict the temporal nature and spatial structure of the dominant instability in a given flow configuration. This is of particular interest, as these instabilities dictate the transition to turbulence and thus are significant for large-scale effects such as skin friction drag. A smoke-flow visualization study with a large variation of parameters, featuring a cylindrical roughness element connected to a linear traverse, has been performed. Results show good agreement with previous investigations and provide further insights into the stability properties, revealing several unexpected effects. For a low roughness aspect ratio $\eta$, no global instability is detected even at the highest roughness Reynolds number $Re_{kk}$, whereas a high aspect ratio indicates a delay in the onset of instability. From the acquired visualizations, we constructed the, so far, richest instability diagram of the wake behind an isolated roughness element in the $Re_{kk}\unicode{x2013}\eta$ space, sampled in the same measurement campaign. Furthermore, information regarding the dominant frequency in the wake can be extracted from the visualization images. Our results suggest a new scaling of the frequency as the velocity is increased. Finally, it is shown that the dominant frequency in a certain flow regime can be well predicted using a Strouhal number based on the cylinder diameter and the roughness velocity.

In this study, the length scaling for the boundary layer separation induced by two incident shock waves is experimentally and analytically investigated. The experiments are performed in a Mach 2.73 flow. A double-wedge shock generator with two deflection angles ($\alpha _1$ and $\alpha _2$) is employed to generate two incident shock waves. Two deflection angle combinations with an identical total deflection angle are adopted: ($\alpha _1 = 7^\circ$, $\alpha _2 = 5^\circ$) and ($\alpha _1 = 5^\circ$, $\alpha _2 = 7^\circ$). For each deflection angle combination, the flow features of the dual-incident shock wave–turbulent boundary layer interactions (dual-ISWTBLIs) under five shock wave distance conditions are examined via schlieren photography, wall-pressure measurements and surface oil-flow visualisation. The experimental results show that the separation point moves downstream with increasing shock wave distance ($d$). For the dual-ISWTBLIs exhibiting a coupling separation state, the upstream interaction length ($L_{int}$) of the separation region approximately linearly decreases with increasing $d$, and the decrease rate of $L_{int}$ with $d$ increases with the second deflection angle under the condition of an identical total deflection angle. Based on control volume analysis of mass and momentum conservations, the relation between $L_{int}$ and $d$ is analytically determined to be approximately linear for the dual-ISWTBLIs with a coupling separation region, and the slope of the linear relation obtained analytically agrees well with that obtained experimentally. Furthermore, a prediction method for $L_{int}$ of the dual-ISWTBLIs with a coupling separation region is proposed, and the relative error of the predicted $L_{int}$ in comparison with the experimental result is $\sim$10 %.

A better understanding of the fluid dynamics of disease transmission by disintegrated respiratory droplets has been the focus of great attention since the recent outbreak of COVID-19. In particular, human respiratory activities such as coughing, sneezing and even talking and eating expel a large amount of pathogen-laden droplets. Particularly, during eating or drinking, the physical properties of saliva can be changed. In this study, we investigate the atomization morphology of expelled artificial saliva mixtures from the perspective of varying fluid physical properties, specifically surface tension and dynamic viscosity. Using high-speed shadowgraph experiments on artificial saliva, we visualize and analyse the disintegration of saliva liquid sheets into ligaments and droplets. We find that the viscosity and surface tension affect the droplet size formed from expelled saliva and follow scaling laws that have been previously observed and predicted for constant shear viscosity. We conclude that the changes in physical properties of saliva induced by eating and drinking tend to favour the formation of smaller droplets during sneezing or coughing, which could drive the airborne transmission pathway of pathogens. Furthermore, we derive a theoretical model based on scaling arguments that shows the breakup time of ligaments produced from the artificial saliva mixtures is dependent on the capillary number.

This work presents the first experimental characterization of the flow field in the vicinity of periodically spaced discrete roughness elements (DRE) in a swept wing boundary layer. The time-averaged velocity fields are acquired in a volumetric domain by high-resolution dual-pulse tomographic particle tracking velocimetry. Investigation of the stationary flow topology indicates that the near-element flow region is dominated by high- and low-speed streaks. The boundary layer spectral content is inferred by spatial fast Fourier transform (FFT) analysis of the spanwise velocity signal, characterizing the chordwise behaviour of individual disturbance modes. The two signature features of transient growth, namely algebraic growth and exponential decay, are identified in the chordwise evolution of the disturbance energy associated with higher harmonics of the primary stationary mode. A transient decay process is instead identified in the near-wake region just aft of each DRE, similar to the wake relaxation effect previously observed in two-dimensional boundary layer flows. The transient decay regime is found to condition the onset and initial amplitude of modal crossflow instabilities. Within the critical DRE amplitude range (i.e. affecting boundary layer transition without causing flow tripping) the transient disturbances are strongly receptive to the spanwise spacing and diameter of the elements, which drive the modal energy distribution within the spatial spectra. In the super-critical amplitude forcing (i.e. causing flow tripping) the near-element stationary flow topology is dominated by the development of a high-speed and strongly fluctuating region closely aligned with the DRE wake. Therefore, elevated shears and unsteady disturbances affect the near-element flow development. Combined with the harmonic modes transient growth these instabilities initiate a laminar streak structure breakdown and a bypass transition process.

The essence of the maximum drag reduction (MDR) state of viscoelastic drag-reducing turbulence (DRT) is still under debate, which mainly holds two different types of views: the marginal state of inertial turbulence (IT) and elasto-inertial turbulence (EIT). To further promote its understanding, this paper conducts a large number of direct numerical simulations of DRT at a modest Reynolds number Re with $Re = 6000$ for the FENE-P model that covers a wide range of flow states and focuses on the problem of how nonlinear extension affects the nature of MDR by varying the maximum extension length $L$ of polymers. It demonstrates that the essence of the MDR state can be both IT and EIT, where $L$ is somehow an important parameter in determining the dominant dynamics. Moreover, there exists a critical $L_c$ under which the minimum flow drag can be achieved in the MDR state even exceeding the suggested MDR limit. Systematic analyses on the statistical properties, energy spectrum, characteristic structures and underly dynamics show that the dominant dynamics of the MDR state gradually shift from IT-related to EIT-related dynamics with an increase of $L$. The above effects can be explained by the effective elasticity introduced by different $L$ at a fixed Weissenberg number (Wi) as well as the excitation of pure EIT. It indicates that larger $L$ introduces more effective elasticity and is favourable to EIT excitation. Therefore, we argue that the MDR state is still dominated by IT-related dynamics for the case of small $L$, but replaced by EIT-related dynamics at high $L$. The obtained results can harmonize the seemingly controversial viewpoints on the dominant dynamics of the MDR state and also provide some ideas for breaking through the MDR limit, such as searching for a polymer solution with a proper molecular length and concentration.

This study analyses the stability characteristics of the shear layer vortices (SLV) in a reacting jet in crossflow, analysing effects of flame position, momentum flux ratio ($J$) and density ratio ($S$). It utilizes 40 kHz particle image velocimetry to characterize the dominant SLV frequencies, streamwise evolution and convective/global stability characteristics for three different canonical configurations, one non-reacting and two reacting (‘R1’ and ‘R2’). In the non-reacting case, both convective and global instability is observed, depending upon $S$ and $J$. Qualitatively similar $S$ dependencies occur for the R1 reacting case where the radial flame position lies outside the jet shear layer, albeit with slower SLV growth rates. When the flame lies inside the jet shear layer, the R2 reacting case, a qualitatively different behaviour is observed, as vorticity concentration in the shear layers is suppressed almost completely. Finally, we show that frequency and stability characteristics of the non-reacting and R1 cases can be scaled in a unified manner using a counter-current shear layer model. This model relates these SLV behaviours to a vorticity layer thickness, a velocity scale and an effective density ratio (noting that there are three distinct densities associated with the jet, the crossflow and the burned gases). These parameters were extracted from the data and used to collapse the frequency scaling, and to explain the transition to self-excited oscillatory behaviour.

The vortex shedding topology of a heavy pendulum oscillating in a dense fluid is investigated using time-resolved three-dimensional particle tracking velocimetry (tr-3-D-PTV). A series of experiments with eight different solid to fluid mass ratios $m^*$ in the range $[1.14, 14.95]$ and corresponding Reynolds numbers of up to $Re \sim O(10^4)$ was conducted. The period of oscillation depends heavily on $m^*$. The relation between amplitude decay and oscillation frequency is non-monotonic, having a damping optimum at $m^* \approx 2.50$. Moreover, a novel digital object tracking (DOT) method using vorticity-magnitude iso-surfaces is implemented to analyse vortical structures. A similar vortex shedding topology is observed for various mass ratios $m^*$. Our observations show that first, a vortex ring in the pendulum's wake is formed. Soon after, the initial ring breaks down to two clearly distinguishable structures of similar size. One of the two vortices remains on the circular path of the pendulum, while the other detaches, propagates downwards, and eventually dissipates. The time when the first vortex is shed, and its initial propagation velocity, depend on $m^*$ and the momentum imparted by the spherical bob. The findings further show good agreement between the experimentally determined vortex shedding frequency and the theoretical vortex shedding time scale based on the Strouhal number.

Three-dimensional steady-state Arnold–Beltrami–Childress (ABC) flow has a chaotic Lagrangian structure, and also satisfies the Navier–Stokes (NS) equations with an external force per unit mass. It is well known that, although trajectories of a chaotic system have sensitive dependence on initial conditions, i.e. the famous ‘butterfly effect’, their statistical properties are often insensitive to small disturbances. This kind of chaos (such as governed by the Lorenz equations) is called normal-chaos. However, a new concept, i.e. ultra-chaos, has been reported recently, whose statistics are unstable to tiny disturbances. Thus, ultra-chaos represents higher disorder than normal chaos. In this paper, we illustrate that ultra-chaos widely exists in Lagrangian trajectories of fluid particles in steady-state ABC flow. Moreover, solving the NS equation when $Re=50$ with the ABC flow plus a very small disturbance as the initial condition, it is found that trajectories of nearly all fluid particles become ultra-chaotic when the transition from laminar to turbulence occurs. These numerical experiments and facts highly suggest that ultra-chaos should have a relationship with turbulence. This paper identifies differences between ultra-chaos and sensitivity of statistics to parameters. Possible relationships between ultra-chaos and the Poincaré section, ultra-chaos and ergodicity/non-ergodicity, etc., are discussed. The concept of ultra-chaos opens a new perspective of chaos, the Poincaré section, ergodicity/non-ergodicity, turbulence and their inter-relationships.

Numerical simulations in two space dimensions are used to examine the dynamics, transport and equilibrium behaviours of a neutrally buoyant circular object immersed in an active suspension within a larger circular container. The continuum model of Gao et al. (Phys. Rev. Fluids, vol. 2, issue 9, 2017, 093302) represents the suspension of non-interacting, immotile, extensor-type microscopic agents that have a direction and strength, and align in response to strain rate. Such a suspension is well known to be unstable above an activity strength threshold, which depends upon the length scale of the confinement. Introducing the object leads to additional phenomenology. It can confine fluid between it and the container wall, which suppresses local suspension activity. However, its motion also correlates strain rates near its surface, which induce a correspondingly correlated active-stress response. Depending on the suspension activity strength, these mechanisms lead to either an attraction toward or a repulsion away from the container wall. In addition, a persistent propagating behaviour is found for modest activity strength, which provides a mechanism for long-range transport. When activity is so weak that the mobility of the object is essential to support suspension instability and sustain flow, all flow terminates when its mobility is diminished as it nears the container wall. If activity strength is scheduled in time, then these mechanisms could be used to perform relatively complex tasks with simple active agents.

Wall-pressure and velocity statistics in the turbulent boundary layer (TBL) on a cambered controlled-diffusion aerofoil at $8^{\circ }$ incidence, a Mach number of 0.25 and a chord-based Reynolds number ${Re}_c=1.5\times 10^{5}$ are analysed at four locations on the suction side with zero and adverse pressure gradients (ZPG and APG), characterised by increasing Reynolds numbers based on momentum thickness, ${Re}_{\theta }=319$, 390, 877 and $1036$. The strong APG yields a highly non-equilibrium TBL at the trailing edge that significantly affects the turbulent flow statistics. Different normalisations of the full wall-pressure statistics involved in trailing-edge noise are analysed for the first time in such strong APG with convex curvature, and compared with available experimental and numerical data. Good overall agreement is found in the ZPG region, and most results obtained in previous APG TBL can be extended to the present highly non-equilibrium case. The presence of strong APG augments the intensity of wall-pressure fluctuations noticeably at low frequencies, shortens the streamwise and broadens the spanwise coherence of wall-pressure fluctuations in both time and space, and significantly reduces the convection velocity. The wall-pressure power spectral density are found to scale with the displacement thickness, the Zaragola–Smits velocity and the root-mean-squared pressure, the latter possibly being replaced by the local maximum Reynolds shear stress. The other two key parameters to trailing-edge noise modelling, the spanwise coherence length and the convection velocity, rather scale with displacement thickness and friction velocity, respectively.