Understanding the interplay between buoyancy and fluid motions within stably stratified shear layers is crucial for unravelling the contribution of flow structures to turbulent mixing. In this study, we examine statistically the local relationship between stratification and fluid deformation rate in wave and turbulent regimes, using experimental datasets obtained from a stratified inclined duct (SID) containing fluids of different densities that form an exchange flow. We introduce rotational and shear components of varying strength within the vorticity and a family of coherent gradient Richardson numbers ($Ri_C$), ratios related to the buoyancy frequency and the strength of either the rotational or shearing motion. Conditional statistical analysis reveals that both shear and stratification intensity affect the probability distribution of the $Ri_C$, with extreme events occurring more frequently in areas of weak stratification. In the wave regime, we identify the persistence of fast-spin vortices within the strongly stratified density interface. However, scouring of the density interface is primarily driven by shearing motions, with baroclinic torque making a notable contribution to enstrophy transport. In the turbulent regime, rigid-body rotations occur at significantly shorter time scales than that associated with the local buoyancy frequency, making them more disruptive to stratification than shear. Additionally, correlation analysis reveals that irrotational strain distorts stable stratification similarly to shearing motions, but is weaker than both shearing and rotational motions and less likely to have a time scale longer than that related to the buoyancy frequency. Moreover, we observed that the interplay between rotational and shearing motions intensifies as stratification increases. Finally, a comparison of length scales along the shear layers highlights the $Ri_C$ as a valuable measure of the relative sizes of different motions compared with the Ozmidov scale and shows that stratification can influence sub-Ozmidov scales through baroclinic torque. This study highlights the critical impact of the type, strength and location of fluid deformations on localised mixing, providing new insights into the role of rotational motions in shear-driven stratified flows.

Steady flow at low Reynolds (Re) number through a planar channel with converging or diverging width is investigated in this study. Along the primary direction of flow, the small dimension of the channel cross-section remains constant while the sidewalls bounding the larger dimension are oriented at a constant angle. Due in part to ease of manufacturing, parallel-plate geometries such as this have found widespread use in microfluidic devices for mixing, heat exchange, flow control and flow patterning at small length scales. Previous analytical solutions for flows of this nature have required the converging or diverging aspect of the channel to be gradual. In this work, we derive a matched asymptotic solution, validated against numerical modelling results, that is valid for any sidewall angle, without requiring the channel width to vary gradually. To accomplish this, a cylindrical coordinate system defined by the angle of convergence between the channel sidewalls is considered. From the mathematical form of the composite expansion, a delineation between two secondary flow components emerges naturally. The results of this work show how one of these two components, originating from viscous shear near the channel sidewalls, corresponds to convective mixing, whereas the other component impresses the sidewall geometry on streamlines in the outer flow.

The Lamb–Oseen vortex is a model for practical vortical flows with a finite vortex core. Vortices with a Lamb–Oseen vortex velocity profile are stable according to the Rayleigh criterion in an infinite domain. Practical situations introduce boundary conditions over finite domains. Direct numerical simulations are performed on the evolution of perturbations to a viscous Lamb–Oseen vortex with uniform inlet axial velocity in a pipe of finite length. Linear stability boundaries are determined in the $(\textit{Re},\omega )$ plane. For a given swirl ratio $\omega$, the flow is found to become linearly unstable when the Reynolds number $\textit{Re}$ is above a critical value. The complete evolution history of the flow is followed until it reaches its final state. For small swirl ratios, the axisymmetric mode is linearly unstable and evolves to a final steady axisymmetric but non-columnar accelerated flow state after nonlinear saturation. For large swirl ratios, the spiral mode is linearly unstable. The spiral mode is found to force growth of an axisymmetric component due to nonlinear interaction. The flow evolves to a final unsteady spiral vortex breakdown state after it undergoes nonlinear saturation. The energy transfer between the mean flow and perturbations is studied by the Reynolds–Orr equation. The pressure work at the exit of the finite pipe is a major source of energy production. Finite-domain boundary conditions also modify the perturbation mode shapes, which can render the vortex core from absorbing energy to producing energy, and thus lead to instabilities. As the pipe length increases, the stability behaviour of the flow is found to approach that predicted by the classical Rayleigh criterion.

We consider the vortex–wedge interaction problem, taking as a departure point Howe’s model of a point vortex interacting with a semi-infinite half-plane, where the vortex path is influenced by its image and a closed-form analytical solution is obtained for the sound field. We generalise Howe’s model to consider wedges of arbitrary angles and explore the influence of vortex circulation, distance from the edge and the wedge half-angle. The effect of wedge angle on sound emission involves a reduced amplitude of the latter as the former is increased. An extension of the model is proposed to account for convection effects by a non-zero ambient flow. We identify a non-dimensional parameter that characterises the vortex kinematics close to the edge and the associated acoustic effect: high and low values of the parameter correspond, respectively, to high- and low-amplitude sound emission of high and low frequency.

We report our finding from direct numerical simulations that polygonal cell structures are formed by inertial particles in turbulent Rayleigh–Bénard convection in a large aspect ratio channel at Rayleigh numbers of $10^6, 10^7$ and $10^8$, and Prandtl number of 0.7. The settling of small particles modified the flow structures only through momentum interactions. From the simulations performed for various sizes and mass loadings of particles, we discovered that for small- and intermediate-sized particles, cell structures such as square, pentagonal or hexagonal cells were observed, whereas a roll structure was formed by large particles. As the mass loading increased, the sizes of the cells or rolls decreased for all particle sizes. The Nusselt number increased with the mass loading of intermediate and large particles, whereas it decreased with the mass loading of small particles compared with the value for particle-free convection. A detailed investigation of the effective feedback forces of the settling particles inside the hot and cold plumes near the walls revealed that the feedback forces break the up–down symmetry between the hot and cold plumes near the surfaces. This enhances the hot plume ascent while not affecting the cold plume, which leads to the preferred formation of cellular structures. The energy budget analysis provides a detailed interaction between particles and fluid, revealing that the net energy is transferred from the fluid to particles when the particles are small, while settling intermediate and large particles drag the fluid so strongly that energy is transferred from particles to fluid.

This study suggests that partial changes in adverse pressure gradient (APG) turbulent boundary layers (TBLs) relative to zero pressure gradient (ZPG) conditions can be obtained quantitatively by the wall-normal integral, while clarifying the partial influence of non-equilibrium effects. Specifically, the term $u_{\tau }^{2}/ ( {U_{e}V_{e}} )$, which is found to describe the degree of scale separation under non-equilibrium conditions, is decomposed into three terms. Here, $u_{\tau }$ is the frictional velocity, $U_{e}$ is the streamwise velocity at the boundary layer edge, and $V_{e}$ is the normal velocity at the boundary layer edge. This equation includes a ZPG term, a pressure gradient term and a streamwise variation term, indicating that the pressure gradient promotes scale separation. The equation can be applied to ZPG TBLs and equilibrium APG TBLs by separately ignoring the pressure gradient term and the streamwise variation term. By using this equation to simplify the integral of the inertia term of the mean momentum equation, an expression for the Reynolds shear stress in the outer region can be obtained, which indicates how APG affects the Reynolds shear stress through the mean velocity. The above quantitative results support further study of non-equilibrium APG TBLs.

The interaction between the flow in a channel with multiple obstructions on the bottom and an elastic ice sheet covering the liquid is studied for the case of steady flow. The mathematical model employs velocity potential theory and fully accounts for the nonlinear boundary conditions at the ice/liquid interface and on the channel bottom. The integral hodograph method is used to derive the complex velocity potential of the flow, explicitly containing the velocity magnitude at the interface. This allows the boundary-value problem to be reduced to a system of nonlinear equations for the unknown velocity magnitude at the ice/liquid interface, which is solved using the collocation method. Case studies are carried out for a widened rectangular obstruction, whose width exceeds the wavelength of the interface, and for arrays of triangular ripples forming the undulating bottom shape. The influence of the bottom shape on the interface is investigated for three flow regimes: the subcritical regime, $F \lt F_{{cr}}$, for which the depth-based Froude number is less than the critical Froude number, and the interface perturbation decays upstream and downstream of the obstruction; the ice-supercritical and channel-subcritical regime, $F_{cr} \lt F \lt 1$, for which two waves of different wavelengths extend upstream and downstream to infinity; and the channel-supercritical regime, $F \gt 1$, for which the hydroelastic wave extends downstream to infinity. The results revealed a trapped interface wave above the rectangular obstruction and the ripple patch. The resonance behaviour of the interface over the undulating bottom occurs when the period of ripples approaches the wavelength of the ice/liquid interface.

Riparian vegetation along riverbanks and seagrass along coastlines interact with water currents, significantly altering their flow. To characterise the turbulent fluid motion along the streamwise-edge of a region covered by submerged vegetation (canopy), we perform direct numerical simulations of a half-channel partially obstructed by flexible stems, vertically clamped to the bottom wall. An intense streamwise vortex forms along the canopy edge, drawing high-momentum fluid into the side of the canopy and ejecting low-momentum fluid from the canopy tip, in an upwelling close to the canopy edge. This mechanism has a profound impact on the mean flow and on the exchange of momentum between the fluid and the structure, which we thoroughly characterise. The signature of the canopy-edge vortex is also found in the dynamical response of the stems, assessed for two different values of their flexibility. Varying the flexibility of the stems, we observe how different turbulent structures over the canopy are affected, while the canopy-edge vortex does not exhibit major modifications. Our results provide a better understanding of the flow in fluvial and coastal environments, informing engineering solutions aimed at containing the water flow and protecting banks and coasts from erosion.

This paper provides direct experimental evidence for the coexistence of both a laminar separation bubble and a secondary vortex on the advancing side of a rotating sphere when subjected to the inverse Magnus effect. Detailed experiments were conducted in a wind tunnel on two spheres of varying surface roughness to investigate both ordinary and inverse Magnus effects. Experiments took place for $0.5\times 10^{5}\leqslant {\textit{Re}}\leqslant 3\times 10^{5}$ and rotation rates $0\leqslant \alpha \leqslant 0.45$, where the spheres were rotated via a shaft that was oriented perpendicularly to the free stream flow. Static pressure measurements were made on the non-shaft hemisphere using a spline of taps spanning from the equator to the pole. The ordinary Magnus effect was generally observed at the lowest ${\textit{Re}}$ tested, with a transition to the inverse Magnus effect occurring as ${\textit{Re}}$ increased. Time-averaged pressure coefficient distributions across the equatorial plane were obtained for the smooth and rough spheres. Cross-flow particle image velocimetry was used to visualise the downstream wake velocity field. A pair of counter-rotating wing-tip-like vortices were detected when the sphere experienced the ordinary Magnus effect, generated by flow leakage from the advancing to the retreating side. When the sphere experienced the inverse Magnus effect, the polarity of the counter-rotating vortex pair reversed. This is the first experimental observation of the vortex polarity reversal associated with the inverse Magnus effect in the wake of a rotating sphere. The results provide qualitative visualisation of the complex fluid dynamics and inform future applications of the Magnus effect.

Within the context of machine learning-based closure mappings for Reynolds-averaged Navier Stokes turbulence modelling, physical realisability is often enforced using ad hoc postprocessing of the predicted anisotropy tensor. In this study, we address the realisability issue via a new physics-based loss function that penalises non-realisable results during training, thereby embedding a preference for realisable predictions into the model. Additionally, we propose a new framework for data-driven turbulence modelling which retains the stability and conditioning of optimal eddy viscosity-based approaches while embedding equivariance. Several modifications to the tensor basis neural network to enhance training and testing stability are proposed. We demonstrate the conditioning, stability and generalisation of the new framework and model architecture on three flows: flow over a flat plate, flow over periodic hills and flow through a square duct. The realisability-informed loss function is demonstrated to significantly increase the number of realisable predictions made by the model when generalising to a new flow configuration. Altogether, the proposed framework enables the training of stable and equivariant anisotropy mappings, with more physically realisable predictions on new data. We make our code available for use and modification by others. Moreover, as part of this study, we explore the applicability of Kolmogorov–Arnold networks to turbulence modelling, assessing its potential to address nonlinear mappings in the anisotropy tensor predictions and demonstrating promising results for the flat plate case.

Thixotropic fluids with a non-monotonic flow curve display viscosity bifurcations at certain stresses. It has been proposed that these transitions can introduce interfaces (or shear bands) into thin films that can destabilize inertialess flows over inclined planes. This proposition is confirmed in the present paper by formulating a thin-film model, then using this model to construct sheet-like base flows and test their linear stability. It is also found that viscosity bifurcations, and the associated interfaces, are not necessary for instability, but that the time-dependent relaxation of the microstructure responsible for thixotropy within the bulk of the film can promote instability instead. Computations with the thin-film model demonstrate that instabilities saturate supercritically into steadily propagating nonlinear waves that travel faster than the mean flow.

The impact of freestream turbulence (FST) on the aerodynamic performance of a flexible finite wing and the produced wingtip vortex was investigated. The wing had a NACA 4412 airfoil profile and the chord-based Reynolds number was $1.4\times 10^{5}$. The experiments were conducted in a closed-loop wind tunnel with four different inflow turbulence intensities ($0.2\,\%$, $3\,\%$, $8\,\%$ and $13\,\%$) generated using an active turbulence grid. Force balance measurements revealed that increasing the scale of the FST increased the maximum lift and delayed stall. Digital image correlation (DIC) measured deflections of the wing’s structure. Spanwise bending was found to be the dominant deformation. While the wing vibrated at its natural frequency in all conditions, FST increased the amplitude of the vibrations. A similar spectral signature was observed in the lift force fluctuations as well. Stereoscopic particle image velocimetry measurements were obtained two chord lengths downstream of the trailing edge simultaneously with DIC. FST decreased the vortex strength, and marginally increased vortex diffusion and size. It also increased the vortex meandering amplitude, while reducing the meandering frequency band. For the cases with a turbulence intensity of $8\,\%$ and $13\,\%$, the frequency of meandering and the wing’s vibration were similar and a phase relation between the two motions was observed. Proper orthogonal decomposition of the vortex (after removing meandering) and the subsequent velocity field reconstruction revealed temporal fluctuations in the vortex strength at the same frequency as the wing’s vibration. This was linked to the lift force fluctuations induced by the wing’s deformations.

The effect of a horizontal magnetic field on heat transport and flow structures in vertical liquid metal convection (Prandtl number $Pr \approx 0.03$) is investigated experimentally. The experiments are carried out for Rayleigh numbers in the range of $1.48 \times 10^6 \leqslant Ra \leqslant 3.54 \times 10^{7}$ and Chandrasekhar numbers in the range of $2 \times 10^2 \leqslant Q \leqslant 1.86 \times 10^6$, as well as for the non-magnetic case ($Q=0$). Measurements of the heat transport show a rise in the Nusselt number at low and moderate magnetic field strengths up to an optimum value of $Q$, before a further increase in the magnetic field leads to a decrease in the transport properties. By applying simultaneous velocity and temperature measurements, we are able to identify three different oscillatory flow regimes for $10^{-5}\lt Q/Ra \lt 0.5$ and assign them to the respective heat transfer characteristics. In the range $10^{-5}\gt Q/Ra\gt 10^{-3}$, first evidence of a transition to anisotropic flow structures caused by the magnetic field is visible. Two strongly oscillatory regimes are identified, where the energy is either distributed around a dominant frequency ($10^{-3}\gt Q/Ra\gt 10^{-2}$), or strongly concentrated on a single frequency ($10^{-2}\gt Q/Ra\gt 0.5$). The dominating frequency increases with the Rayleigh number according to $Ra^{0.71\pm 0.02}$. This flow structure based regime separation correspond to changes of both the heat transfer through the Nusselt number and mass transfer through the Reynolds number.

Magnetohydrodynamic turbulence with Hall effects is ubiquitous in heliophysics and plasma physics. Direct numerical simulations reveal that, when the forcing scale is comparable to the ion inertial scale, the Hall effects induce remarkable cross-helicity. It then suppresses the cascade efficiency, leading to the accumulation of large-scale magnetic energy and helicity. The process is accompanied by the disruption of current sheets through the entrainment by vortex tubes or the excitation of whistler waves. Using the solar wind data from the Parker Solar Probe, the numerical findings are separately confirmed. These findings provide new insights into the emergence of large-scale solar wind turbulence driven by helical fields and Hall effects.

Turbulent separating and reattaching flows are known to exhibit low-frequency fluctuations manifested in a large-scale contraction and expansion of the reverse-flow region. Previous experimental investigations have been restricted to planar measurements, while the computational cost to resolve the low-frequency spectrum with high-fidelity simulations currently appears to be unaffordable. In this article, we make use of volumetric measurements to reveal the low-frequency dynamics of a turbulent separation bubble (TSB) formed in the fully turbulent flow past a smooth backward-facing ramp. The volumetric velocity field measurements cover the entire separated flow region over a domain with a spanwise extent of $S=0.6\, {\textrm{m}}$. Spectral proper orthogonal decomposition (SPOD) of the velocity fluctuations reveals low-rank low-frequency behaviour at Strouhal numbers ${\textit{St}}\lt 0.05$, which was also observed in previous planar measurements. However, in contrast with the interpretation of a two-dimensional contraction/expansion motion, the low-frequency dynamics is shown to be inherently three-dimensional, and governed by large elongated structures with a spanwise wavelength of approximately $S/2$. A low-order model constructed with the leading SPOD mode confirms substantial changes of the TSB extent in the centre plane, linking it to the modal pattern that is strongly non-uniform in the spanwise direction. The findings presented in this study promote a more complete understanding of the low-frequency dynamics in turbulent separated flows, thereby enabling novel modelling and control approaches.