We present herein the derivation of a lubrication-mediated (LM) quasi-potential model for droplet rebounds off deep liquid baths, assuming the presence of a persistent dynamic air layer which acts as a lubricating pressure transfer. We then present numerical simulations of the LM model for axisymmetric rebounds of solid spheres and compare quantitatively to current results in the literature, including experimental data in the low-speed impact regime. In this regime the LM model has the advantage of being far more computationally tractable than direct numerical simulation (DNS) and is also able to provide detailed behaviour within the micro-metric thin lubrication region. The LM system has an interesting mathematical structure, with the lubrication layer providing a free-boundary elliptic problem mediating the drop and bath free-boundary evolutionary equations.

Microscopic irregularity (roughness) of bounding surfaces affects macroscopic dynamics of fluid flows. Its effect on bulk flow is usually quantified empirically by means of a roughness coefficient. A new approach, which treats rough surfaces (e.g. parallel plates) as random fields whose statistical properties can be inferred from measurements, is presented. The mapping of a random flow domain onto its deterministic counterpart, and the subsequent stochastic averaging of the transformed Stokes equations, yield expressions for the effective viscosity and roughness coefficient in terms of the statistical characteristics of the irregular geometry of the boundaries. The analytical nature of the solutions allows one to handle surface roughness characterized by short correlation lengths, a challenging feature for numerical stochastic simulations.

Radial unstable stratification is a potential source of turbulence in the cold regions of accretion disks. To investigate this thermal effect, here we focus on two-dimensional Rayleigh–Bénard convection in an annulus subject to radially dependent gravitational acceleration $g \propto 1/r$. Next to the Rayleigh number $Ra$ and Prandtl number $Pr$, the radius ratio $\eta$, defined as the ratio of inner and outer cylinder radii, is a crucial parameter governing the flow dynamics. Using direct numerical simulations for $Pr=1$ and $Ra$ in the range from $10^7$ to $10^{10}$, we explore how variations in $\eta$ influence the asymmetry in the flow field, particularly in the boundary layers. Our results show that in the studied parameter range, the flow is dominated by convective rolls and that the thermal boundary-layer (TBL) thickness ratio between the inner and outer boundaries varies as $\eta ^{1/2}$. This scaling is attributed to the equality of velocity scales in the inner ($u_i$) and outer ($u_o$) regions. We further derive that the temperature drops in the inner and outer TBLs scale as $1/(1+\eta ^{1/2})$ and $\eta ^{1/2}/(1+\eta ^{1/2})$, respectively. The scalings and the temperature drops are in perfect agreement with the numerical data.

The flow of a nematic liquid crystal in a Hele-Shaw cell with an electrically controlled viscous obstruction is investigated using both a theoretical model and physical experiments. The viscous obstruction is created by temporarily electrically altering the viscosity of the nematic in a region of the cell across which an electric field is applied. The theoretical model is validated experimentally for a circular cylindrical obstruction, demonstrating user-controlled flow manipulation of an anisotropic liquid within a heterogeneous single-phase microfluidic device.

The recent work of Siegelman & Young (Proc. Natl Acad. Sci. USA, vol. 120, issue 44, 2023, e2308018120) revealed two extreme states reached by the evolution of unforced and weakly damped two-dimensional turbulence above random rough topography, separated by a critical kinetic energy $E_\#$. The low- and high-energy solutions correspond to topographically locked and roaming vortices, surrounded by non-uniform and homogeneous background potential vorticity (PV), respectively. However, we found that these phenomena are restricted to the particular intermediate length scale where the energy was initially injected into the system. Through simulations initialized by injecting the energy at larger and smaller length scales, we found that the long-term state of topographic turbulence is also dependent on the initial length scale and thus the initial enstrophy. If the initial length scale is comparable to the domain size, the long-term flow field resembles the minimum-enstrophy state proposed by Bretherton & Haidvogel (J. Fluid Mech., vol. 78, issue 1, 1976, pp. 129–154), with very few topographically locked vortices; the long-term enstrophy is quite close to the minimum value, especially when the energy is no larger than $E_\#$. As the initial length scale becomes smaller, more vortices nucleate and become more mobile; the long-term enstrophy increasingly deviates from the minimum value. Simultaneously, the background PV tends to homogenization, even if the energy is below $E_\#$. These results complement the phenomenology of topographic turbulence documented by Siegelman & Young, by showing that the minimum-enstrophy and background PV homogenization states can be adequately approached by large- and small-scale initial fields, respectively, with relatively arbitrary energy.

We discuss the modal, linear stability analysis of generalized Couette–Poiseuille (GCP) flow between two parallel plates moving with relative speed in the presence of an applied pressure gradient vector inclined at an angle $0\leqslant \phi \leqslant 90^\circ$ to the plate relative velocity vector. All possible GCP flows can be described by a global Reynolds number $Re$, $\phi$ and an angle $0\leqslant \theta \leqslant 90^\circ$, where $\cos \theta$ is a measure of the relative weighting of Couette flow to the composite GCP flow. This provides a novel and uncommon group of generally three-dimensional base velocity fields with wall-normal twist, for which Squire's theorem does not generally apply, requiring study of oblique perturbations with wavenumbers $(\alpha,\beta )$. With $(\theta,\phi )$ fixed, the neutral surface $f(\theta,\phi ;Re,\alpha,\beta )=0$ in $(Re, \alpha,\beta )$ space is discussed. A mapping from GCP to plane Couette–Poiseuille flow stability is found that suggests a scaling relation $Re^*\alpha /k = H(\theta ^*)$ that collapses all critical parameters, where ${Re}^*= Re\,({\alpha _1}/{\alpha })\,({\sin \theta }/{\sin {\theta }^*})$ and $\tan \theta ^*=({\alpha _{1}}/\alpha )\tan \theta$, with $\alpha _1=\alpha \cos \phi +\beta \sin \phi$. This analysis does not, however, directly reveal global critical properties for GCP flow. The global $Re_{cr}(\theta,\phi )$ shows continuous variation, while $\alpha _{cr}(\theta,\phi )$ and $\beta _{cr}(\theta,\phi )$ show complex behaviour, including discontinuities owing to jumping of critical states across neighbouring local valleys (in $Re$) or lobes of the neutral surface. The discontinuity behaviour exists for all low $\phi$. For $\phi \gtrsim 21^\circ$, variations of $\alpha _{cr}(\theta )$ and $\beta _{cr}(\theta )$ are generally smooth and monotonic.

The first experimental results on pattern transitions in the co-rotation regime (i.e. the rotation ratio $\varOmega = \omega _o/\omega _i > 0$, where $\omega _i$ and $\omega _o$ are the angular speeds of the inner and outer cylinders, respectively) of the Taylor–Couette flow (TCF) are reported for a neutrally buoyant suspension of non-colloidal particles, up to a particle volume fraction of $\phi = 0.3$. While the stationary Taylor vortex flow (TVF) is the primary bifurcating state in dilute suspensions ($\phi \leq ~0.05$), the non-axisymmetric oscillatory states, such as the spiral vortex flow (SVF) and the ribbon (RIB), appear as primary bifurcations with increasing particle loading, with an overall de-stabilization of the primary bifurcating states (TVF/SVF/RIB) being found with increasing $\phi$ for all $\varOmega \geq ~0$. At small co-rotations ($\varOmega \sim 0$), the particles play the dual role of stabilization ($\phi < 0.1$) and destabilization ($\phi \geq ~0.1$) on the secondary/tertiary oscillatory states. The distinctive features of the ‘particle-induced’ spiral vortices are identified and contrasted with those of the ‘fluid-induced’ spirals that operate in the counter-rotation regime.

Linear unsteady aerofoil theory, while successfully used for the prediction of unsteady aerofoil lift for many decades, has yet to be proven adequate for predicting the propulsive performance of oscillating aerofoils. In this paper we test the hypothesis that the central shortcoming of linear small-amplitude models, such as the Garrick function, is the failure to account for the flow acceleration caused by aerofoil thrust. A new analytical model is developed by coupling the Garrick function to a cycle-averaged actuator disc model, in a manner analogous to the blade-element momentum theory for wind turbines and propellers. This amounts to assuming the Garrick function to be locally valid and, in combination with a global control volume analysis, enables the prediction of flow acceleration at the aerofoil. The new model is demonstrated to substantially improve the agreement with large-eddy simulations of an aerofoil in combined heave and pitch motion.

Modons, or dipolar vortices, are common and long-lived features of the upper ocean, consisting of a pair of counter-rotating monopolar vortices moving through self-advection. Such structures remain stable over long times and may be important for fluid transport over large distances. Here, we present a semi-analytical method for finding fully nonlinear modon solutions in a multi-layer quasi-geostrophic model with arbitrarily many layers. Our approach is to reduce the problem to a multi-parameter linear eigenvalue problem which can be solved using numerical techniques from linear algebra. The method is shown to replicate previous results for one- and two-layer models and is applied to a three-layer model to find a solution describing a mid-depth propagating, topographic vortex.

The stimulation of instability and transport in the bottom boundary layer by internal solitary waves has been documented for over twenty years. However, the challenge of shallow slopes and a disparity of scales between the large-scale wave and the small-scale boundary layer has proven challenging for simulations. We present laboratory scale simulations that resolve the three-dimensionalisation in the boundary layer during the entire shoaling process. We find that the late stage, in which the incoming wave fissions into boluses, provides the most consistent source of three-dimensionalisation. In the early stage of shoaling, three-dimensionalisation occurs not so much due to separation bubble instability, but due to the interaction of vortices shed from the separation bubble with the overlying pycnocline. This interaction overturns the pycnocline, and creates bursts in kinetic energy and viscous dissipation, suggesting that the shed vortices induce turbulent motion and sediment resuspension in the water column above and behind the separation bubble.

We present a systematic study on the effects of small aspect ratios $\varGamma$ on heat transport in liquid metal convection with a Prandtl number of $Pr=0.029$. The study covers $1/20\le \varGamma \le 1$ experimentally and $1/50\le \varGamma \le 1$ numerically, and a Rayleigh number $Ra$ range of $4\times 10^3 \le Ra \le 7\times 10^{9}$. It is found experimentally that the local effective heat transport scaling exponent $\gamma$ changes with both $Ra$ and $\varGamma$, attaining a $\varGamma$-dependent maximum value before transition-to-turbulence and approaches $\gamma =0.25$ in the turbulence state as $Ra$ increases. Just above the onset of convection, Shishkina (Phys. Rev. Fluids, vol 6, 2021, 090502) derived a length scale $\ell =H/(1+1.49\varGamma ^{-2})^{1/3}$. Our numerical study shows $Ra_{\ell }$, i.e. $Ra$ based on $\ell$, serves as a proper control parameter for heat transport above the onset with $Nu-1=0.018(1+0.34/\varGamma ^2)(Ra/Ra_{c,\varGamma }-1)$. Here $Ra_{c,\varGamma }$ represents the $\varGamma$-dependent critical $Ra$ for the onset of convection and $Nu$ is the Nusselt number. In the turbulent state, for a general scaling law of $Nu-1\sim Ra^\alpha$, we propose a length scale $\ell = H/(1+1.49\varGamma ^{-2})^{1/[3(1-\alpha )]}$. In the case of turbulent liquid metal convection with $\alpha =1/4$, our measurement shows that the heat transport will become weakly dependent on $\varGamma$ with $Ra_{\ell }\equiv Ra/(1+1.49\varGamma ^{-2})^{4/3} \ge 7\times 10^5$. Finally, once the flow becomes time-dependent, the growth rate of $Nu$ with $Ra$ declines compared with the linear growth rate in the convection state. A hysteresis is observed in a $\varGamma =1/3$ cell when the flow becomes time-dependent. Measurements of the large-scale circulation suggest the hysteresis is caused by the system switching from a single-roll-mode to a double-roll-mode in an oscillation state.

Bluff-body wakes generally become three-dimensional (3-D) and then turbulent when the Reynolds number exceeds a few hundred. Other than an alternate shedding of the spanwise vortices behind the body and a gradual decay and annihilation of the vortices with distance downstream, whether a secondary vortex street would develop in the relatively far wake has been a long-standing argument in the literature. This argument is addressed in the present study. Specifically, direct numerical simulations and transient growth analysis are performed to examine the two-dimensional and 3-D wakes of different bluff bodies, including circular cylinder, square cylinder, diamond cylinder and rectangular cylinders with different cross-sectional aspect ratios. We found that a secondary vortex street is absent for most 3-D and turbulent wakes. The root cause is the weakening of spanwise vortices by 3-D wake instability modes and streamwise circulation/vorticity. The weakened spanwise vortices induce reduced mean shear in the intermediate wake, which then induces much smaller perturbation energy growth that is below the threshold for the emergence of a secondary vortex street. This finding suggests that the 3-D and turbulence characteristics, and the momentum, mass and heat transport in the relatively far wake of bluff bodies, would not be influenced by extra anisotropy or inhomogeneity caused by a secondary vortex street.

We report on an experimental study in which Lagrangian tracking is applied to millions of microscopic particles floating on the free surface of turbulent water. We leverage a large jet-stirred zero-mean-flow apparatus, where the Reynolds number is sufficiently high for an inertial range to emerge while the surface deformation remains minimal. Two-point statistics reveal specific features of the flow, deviating from the classic description derived for incompressible turbulence. The magnitude of the relative velocity is strongly intermittent, especially at small separations, leading to anomalous scaling of the second-order structure functions in the dissipative range. This is driven by the divergent component of the flow, leading to fast approaching/separation rates of nearby particles. The Lagrangian relative velocity shows strong persistence of the initial state, such that the ballistic pair separation extends to the inertial range of time delays. Based on these observations, we propose a classification of particle pairs based on their initial separation rate. When this is much smaller than the relative velocity prescribed by inertial scaling (which is the case for the majority of the observed particle pairs), the relative velocity transitions to a diffusive growth and the Richardson–Obukhov super-diffusive dispersion is recovered.

At high incidence, low-aspect-ratio wings present a unique set of aerodynamic characteristics, including flow separation, vortex shedding and unsteady force production. Furthermore, low-aspect-ratio wings exhibit a highly impactful tip vortex, which introduces strong spanwise gradients into an already complex flow. In this work, we explore the interaction between leading-edge flow separation and a strong, persistent tip vortex over a Reynolds number range of $600 \leq Re \leq 10{\,}000$. In performing this study, we aim to bridge the insight gained from existing low-Reynolds-number studies of separated flow on finite wings ($Re \approx 10^2$) and turbulent flows at higher Reynolds numbers ($Re \approx 10^4$). Our study suggests two primary effects of the Reynolds number. First, we observe a break from periodicity, along with a dramatic increase in the intensity and concentration of small-scale eddies, as we shift from $Re = 600$ to $Re = 2500$. Second, we observe that many of our flow diagnostics, including the time-averaged aerodynamic force, exhibit reduced sensitivity to Reynolds number beyond $Re = 2500$, an observation attributed to the stabilising impact of the wing tip vortex. This latter point illustrates the manner by which the tip vortex drives flow over low-aspect-ratio wings, and provides insight into how our existing understanding of this flow field may be adjusted for higher-Reynolds-number applications.

The effects of thermal convection on turbulence in accretion discs, and particularly its interplay with the magnetorotational instability (MRI), are of significant astrophysical interest. Despite extensive theoretical and numerical studies, such an interplay has not been explored experimentally. We conduct linear analysis of the azimuthal version of MRI (AMRI) in the presence of thermal convection and compare the results with our experimental data published before. We show that the critical Hartmann number ($Ha$) for the onset of AMRI is reduced by convection. Importantly, convection breaks symmetry between $m = \pm 1$ instability modes ($m$ is the azimuthal wavenumber). This preference for one mode over the other makes the AMRI wave appear as a ‘one-winged butterfly’.

The nonlinear dynamics of a helical vortex disturbed by a long-wave-instability mode is studied by direct numerical simulation. Vortex reconnection or self-reconnection of the helical vortex is shown to play a crucial role depending on the pitch of the helical vortex. For the larger pitch, a vortex ring is created after the vortex reconnection; it detaches from the remaining helical vortex, whose pitch is doubled. A vortex ring is also created for the smaller pitch; however, it forms a linked system with the remaining vortex. The topological constraint due to this linkage forces strong interaction between the different parts of the helical vortex, leading to turbulent transition.

Recent theoretical progress using multiscale asymptotic analysis has revealed various possible regimes of stratified turbulence. Notably, buoyancy transport can either be dominated by advection or diffusion, depending on the effective Péclet number of the flow. Two types of asymptotic models have been proposed, which yield measurably different predictions for the characteristic vertical velocity and length scale of the turbulent eddies in both diffusive and non-diffusive regimes. The first, termed a ‘single-scale model’, is designed to describe flow structures having large horizontal and small vertical scales, while the second, termed a ‘multiscale model’, additionally incorporates flow features with small horizontal scales, and reduces to the single-scale model in their absence. By comparing predicted vertical velocity scaling laws with direct numerical simulation data, we show that the multiscale model correctly captures the properties of strongly stratified turbulence within regions dominated by small-scale isotropic motions, whose volume fraction decreases as the stratification increases. Meanwhile its single-scale reduction accurately describes the more orderly, layer-like, quiescent flow outside those regions.

Direct numerical simulations (DNS) were carried out to investigate flow control for transition delay using steady blowing/suction strips at the wall of a flared cone at Mach 6 and zero angle of attack. For the numerical investigations of the transition control strategy, the flared cone geometry and the flow conditions of the experiments in the Boeing/Air Force Office of Scientific Research (AFOSR) Mach 6 Quiet Tunnel (BAM6QT) at Purdue University were chosen. For the DNS, transition was initiated by introducing random disturbances at the inflow of the computational domain, emulating ‘natural’ transition in wind-tunnel experiments caused by free-stream noise. In both wind-tunnel experiments and numerical simulations, streamwise ‘hot’ streaks were found on the surface of the flared cone, which are caused by a nonlinear interaction of an axisymmetric second-mode wave and a pair of oblique waves of the same frequency (‘fundamental resonance’). The objective of the flow control strategy proposed here is to delay the transition onset, and thus mitigate the negative consequences associated with the nonlinear transition stages, i.e. the development of hot streaks and large wall-pressure amplitudes that were observed in experiments and DNS. Our previous so-called ‘controlled’ transition simulations have shown that flow control using steady blowing and suction strips can lead to a significant delay of the hot streak development on the surface of the flared cone. The simulation results presented in this paper show that this flow control strategy remains effective, even in a natural transition scenario characterized by broadband disturbances.

This study proposes a novel super-resolution (or SR) framework for generating high-resolution turbulent boundary layer (TBL) flow from low-resolution inputs. The framework combines a super-resolution generative adversarial neural network (SRGAN) with down-sampling modules (DMs), integrating the residual of the continuity equation into the loss function. The DMs selectively filter out components with excessive energy dissipation in low-resolution fields prior to the super-resolution process. The framework iteratively applies the SRGAN and DM procedure to fully capture the energy cascade of multi-scale flow structures, collectively termed the SRGAN-based energy cascade reconstruction framework (EC-SRGAN). Despite being trained solely on turbulent channel flow data (via ‘zero-shot transfer’), EC-SRGAN exhibits remarkable generalization in predicting TBL small-scale velocity fields, accurately reproducing wavenumber spectra compared to direct numerical simulation (DNS) results. Furthermore, a super-resolution core is trained at a specific super-resolution ratio. By leveraging this pretrained super-resolution core, EC-SRGAN efficiently reconstructs TBL fields at multiple super-resolution ratios from various levels of low-resolution inputs, showcasing strong flexibility. By learning turbulent scale invariance, EC-SRGAN demonstrates robustness across different TBL datasets. These results underscore the potential of EC-SRGAN for generating and predicting wall turbulence with high flexibility, offering promising applications in addressing diverse TBL-related challenges.

We consider the long-time propagation of a Boussinesq inertia–buoyancy (large-Reynolds- number) gravity current released from a lock over a downslope of angle $\gamma$, affected by entrainment and drag. We show that the shallow-water (depth-averaged) equations with a Benjamin-type front-jump condition admit a similarity solution $x_N(t) = K t^{2/3}$ while $h, \phi, u$ change like $t$ to the power of $2/3, -4/3, -1/3$, respectively; here $x_N, h, \phi, u$ and $t$ are the position of the nose (distance from backwall), thickness, concentration of dense fluid, velocity and time, respectively, and K is a constant. Assuming that $\gamma$ and the coefficients of entrainment and drag are constant, we derive an analytical exact solution for the similarity profiles and show that $K \propto (\tan \gamma )^{1/3}$; the driving of the slope is balanced by entrainment and/or drag. The predicted $t^{2/3}$ propagation is in agreement with previously published experimental data but a conclusive quantitative assessment of the present theory cannot be performed due to various uncertainties (discussed in the paper) that must be resolved by future work.