Roughness of the surface underlying the atmospheric boundary layer causes departures of the near-surface scalar and momentum transport in comparison with aerodynamically smooth surfaces. Here, we investigate the effect of $56\times 56$ homogeneously distributed roughness elements on bulk properties of a turbulent Ekman flow. Direct numerical simulation in combination with an immersed boundary method is performed for fully resolved, three-dimensional roughness elements. The packing density is approximately $10\,\%$ and the roughness elements have a mean height in wall units of $10 \lesssim H^+ \lesssim 40$. According to their roughness Reynolds numbers, the cases are transitionally rough, although the roughest case is on the verge of being fully rough. We derive the friction of velocity and of the passive scalar through vertical integration of the respective balances. Thereby, we quantify the enhancement of turbulent activity with increasing roughness height and find a scaling for the friction Reynolds number that is verified up to $Re_\tau \approx 2700$. The higher level of turbulent activity results in a deeper logarithmic layer for the rough cases and an increase of the near-surface wind veer in spite of higher $Re_\tau$. We estimate the von Kármán constant for the horizontal velocity $\kappa _{m}=0.42$ (offset $A=5.44$) and for the passive scalar $\kappa _{h}=0.35$ (offset $\mathbb {A}=4.2$). We find an accurate collapse of the data under the rough-wall scaling in the logarithmic layer, which also yields a scaling for the roughness parameters $z$-nought for momentum ($z_{0{m}}$) and the passive scalar ($z_{0{h}}$).

Numerous studies have indicated that turbulence typically initiates along the boundary layer of the stationary disk within a rotor–stator cavity. To describe the transition process to turbulence on the stationary side of a closed rotor–stator cavity, a comprehensive approach combining global linear stability analysis with direct numerical simulation was adopted in the present study. The proposed model aligns with that of Yim et al. (J. Fluid Mech., vol. 848, 2018, pp. 631–647), who investigated the stability characteristics of the rotating-disk boundary layer in a rotor–stator cavity. In order to achieve a stable inflow for the stationary-disk boundary layer, we rotate the shroud together with the rotating disk. Through careful global stability analysis, the predominant spiral mode exhibiting the highest instability in the boundary layer of the stationary disk was discerned, corroborating observations from simulations. Initially, the spiral mode undergoes linear amplification, reaches a state of linear saturation and enters the nonlinear regime. Following nonlinear saturation in the flow field, a circular wave mode arises due to the influence of mean flow distortion. As the Reynolds number attained a sufficiently high level, the interplay between the downstream-propagating circular mode and spiral mode amplified disturbances in the boundary layer of the stationary disk, ultimately leading to the development of localised turbulence at the mid-radius of the rotor–stator cavity. Notably, the present study is the first to elucidate the coexistence of laminar–transitional–turbulent flow states in the stationary-disk boundary layer through direct numerical simulations.

A two-dimensional body that moves suddenly in a viscous fluid can instantly generate vortices at its sharp edges. Recent work using inviscid flow theory, based on the Birkhoff–Rott equation and the Kutta condition, predicts that the ‘starting vortices’ generated by the sharp and straight edges of a body – i.e. the vortices formed immediately after motion begins – can be one of three distinct self-similar types. We explore the existence of these starting vortices for a flat plate and two symmetric Joukowski aerofoils immersed in a viscous fluid, using high-fidelity direct numerical simulations (DNS) of the Navier–Stokes equations. A lattice Green's function method is employed and simulations are performed for chord Reynolds numbers ranging from 5040 to 45 255. Vortices generated at the leading and trailing edges of the flat plate show agreement with the derived inviscid theory, for which a detailed assessment is reported. Agreement is also observed for the two symmetric Joukowski aerofoils, demonstrating the utility of the inviscid theory for arbitrary bodies. While this inviscid theory predicts an abrupt transition between the starting-vortex types, DNS shows a smooth transition. This behaviour occurs for all Reynolds numbers and is related to finite-time effects – there is a maximal time for which the (self-similar) starting vortices exist. We confirm the inviscid prediction that the leading-edge starting vortex of a flat plate can be suppressed dynamically. This has implications for the performance of low-speed aircraft such as model aeroplanes, micro air vehicles and unmanned air vehicles.

Previous work suggests that the arrangement of elements in an obstruction may influence the bulk flow velocity through the obstruction, but the physical mechanisms for this influence are not yet clear. This is the motivation for this study, where direct numerical simulation is used to investigate flow through an array of cylinders at a resolution sufficient to observe interactions between wakes of individual elements. The arrangement is altered by varying the gap ratio $G/d$ (1.2 – 18, G is the distance between two adjacent cylinders, d is the cylinder diameter), array-to-element diameter ratio $D/d$ (3.6 – 200, D is the array diameter), and incident flow angle ($0^{\circ} - 30^{\circ}$). Depending on the element arrangement, it is found that the average root-mean-square lift and drag coefficients can vary by an order of magnitude, whilst the average time-mean drag coefficient of individual cylinders ($\overline{C_{d}}$), and the bulk velocity are found to vary by up to $50\,\%$ and a factor of 2, respectively. These arrangement effects are a consequence of the variation in flow and drag characteristics of individual cylinders within the array. The arrangement effects become most critical in the intermediate range of flow blockage parameter $\mathit{\Gamma_{D}^{\prime}} = 0.5-1.5$ ($\mathit{\Gamma_{D}^{\prime}}=\overline{C_{d}}aD/(1-\phi)$, where a is frontal element area per unit volume, and $\phi$ is solid volume fraction), due to the high variability in element-scale flow characteristics. Across the full range of arrangements modelled, it is confirmed that the bulk velocity is governed by flow blockage parameter but only if the drag coefficient incorporates arrangement effects. Using these results, this paper proposes a framework for describing and predicting flow through an array across a variety of arrangements.

The $k^{-23/6}$ wave action spectrum with an inverse cascade is one of the fundamental Kolmogorov–Zakharov solutions for gravity wave turbulence, which is part of the citation for the Dirac Medal in 2003. Instead of confirming this solution, however, several existing simulations and experiments suggest a spectrum of $k^{-3}$ in set-ups corresponding to the inverse cascade. We provide a theoretical explanation for the latter, considering the condensate that naturally forms in finite domains of experiments/simulations. Our new theory hinges on: (1) derivation of a spectral diffusion equation when non-local interactions with the condensate become dominant, for the first time systematically formulated for quartet-interaction systems; and (2) careful analysis of the asymptotics of interaction coefficient with a remarkable cancellation of all leading-order terms.

This work is devoted to a theoretical and numerical study of the dynamics of a two-phase system vapour bubble in equilibrium with its liquid phase under translational vibrations in the absence of gravity. The bubble is initially located in the container centre. The liquid and vapour phases are considered as viscous and incompressible. Analysis focuses on the vibrational conditions used in experiments with the two-phase system SF$_6$ in the MIR space station and with the two-phase system para-Hydrogen (p-H$_2$) under magnetic compensation of Earth's gravity. These conditions correspond to small-amplitude high-frequency vibrations. Under vibrations, additionally to the forced oscillations, an average displacement of the bubble to the wall is observed due to an average vibrational attraction force related to the Bernoulli effect. Vibrational conditions for SF$_6$ correspond to much smaller average vibrational force (weak vibrations) than for p-H$_2$ (strong vibrations). For weak vibrations, the role of the initial vibration phase is crucial. The difference in the behaviour at different initial phases is explained using a simple mechanical model. For strong vibrations, the average displacement to the wall stops when the bubble reaches a quasi-equilibrium position where the resulting average force is zero. At large vibration velocity amplitudes this position is near the wall where the bubble performs only forced oscillations. At moderate vibration velocity amplitudes the bubble average displacement stops at a finite distance from the wall, then large-scale damped oscillations around this position accompanied by forced oscillations are observed. Bubble shape oscillations and the parametric resonance of forced oscillations are also studied.

Thin-film equations are utilised in many different areas of fluid dynamics when there exists a direction in which the aspect ratio can be considered small. We consider thin free films with Marangoni effects in the extensional flow regime, where velocity gradients occur predominantly along the film. In practice, because of the local deposition of surfactants or input of energy, asymmetric distributions of surfactants or surface tension more generally, are possible. Such examples include the surface of bubbles and the rupture of thin films. In this study, we consider the asymmetric thin-film equations for extensional flow with Marangoni effects. Concentrating on the case of small Reynolds number $ Re $, we study the deposition of insoluble surfactants on one side of a liquid sheet otherwise at rest and the resulting thinning and rupture of the sheet. The analogous problem with a uniformly thinning liquid sheet is also considered. In addition, the centreline deformation is discussed. In particular, we show analytically that if the surface tension isotherm $\sigma = \sigma (\varGamma )$ is nonlinear (surface tension $\sigma$ varies with surfactant concentration $\varGamma$), then accounting for top–bottom asymmetry leads to slower (faster) thinning and pinching if $\sigma = \sigma (\varGamma )$ is convex (concave). The analytical progress reported in this paper allows us to discuss the production of satellite drops from rupture via Marangoni effects, which, if relevant to surface bubbles, would be an aerosol production mechanism that is distinct from jet drops and film drops.

Turbulence can have a strong effect on the fall speed of snowflakes and ice crystals. In this experimental study, the behaviour of thin disks falling in homogeneous turbulence is investigated, in a range of parameters relevant to plate crystals. Disks ranging in diameter from 0.3 to 3 mm, and in Reynolds number $Re = 10\unicode{x2013}435$, are dispersed in two air turbulence levels, with velocity fluctuations comparable to the terminal velocity. For each case, thousands of trajectories are captured and reconstructed by high-speed laser imaging, allowing for statistical analysis of the translational and rotational dynamics. Air turbulence reduces the disk terminal velocities by up to 35 %, with the largest diameters influenced most significantly, which is primarily attributed to drag nonlinearity. This is evidenced by large lateral excursions of the trajectories, which correlate with cross-flow-induced drag enhancement as reported previously for falling spheres and rising bubbles. As the turbulence intensity is increased, flat-falling behaviour is progressively eliminated and tumbling becomes prevalent. The rotation rates of the tumbling disks, however, remain similar to those displayed in still air. This is due to their large moment of inertia compared to the surrounding fluid, in stark contrast with studies conducted in water. In fact, the observed reduction of settling velocity is opposite to previous findings on disks falling in turbulent water. This emphasizes the importance of the solid-to-fluid density ratio in analogous experiments that aim to mimic the behaviour of frozen hydrometeors.

Mixing describes the process by which solutes evolve from an initial heterogeneous state to uniformity under the stirring action of a fluid flow. Fluid stretching forms thin scalar lamellae that coalesce due to molecular diffusion. Owing to the linearity of the advection–diffusion equation, coalescence can be envisioned as an aggregation process. Here, we demonstrate that in smooth two-dimensional chaotic flows, mixing obeys a correlated aggregation process, where the spatial distribution of the number of lamellae in aggregates is highly correlated with their elongation, and is set by the fractal properties of the advected material lines. We show that the presence of correlations makes mixing less efficient than a completely random aggregation process because lamellae with similar elongations and scalar levels tend to remain isolated from each other. We show that correlated aggregation is uniquely determined by a single exponent that quantifies the effective number of random aggregation events. These findings expand aggregation theories to a larger class of systems, which have relevance to various fundamental and applied mixing problems.

Recently, direct numerical simulations (DNS) of stably stratified turbulence have shown that as the Prandtl number ($Pr$) is increased from 1 to 7, the mean turbulent potential energy dissipation rate (TPE-DR) drops dramatically, while the mean turbulent kinetic energy dissipation rate (TKE-DR) increases significantly. Through an analysis of the equations governing the fluctuating velocity and density gradients we provide a mechanistic explanation for this surprising behaviour and test the predictions using DNS. We show that the mean density gradient gives rise to a mechanism that opposes the production of fluctuating density gradients, and this is connected to the emergence of ramp cliffs. The same term appears in the velocity gradient equation but with the opposite sign, and is the contribution from buoyancy. This term is ultimately the reason why the TPE-DR reduces while the TKE-DR increases with increasing $Pr$. Our analysis also predicts that the effects of buoyancy on the smallest scales of the flow become stronger as $Pr$ is increased, and this is confirmed by our DNS data. A consequence of this is that the standard buoyancy Reynolds number does not correctly estimate the impact of buoyancy at the smallest scales when $Pr$ deviates from 1, and we derive a suitable alternative parameter. Finally, an analysis of the filtered gradient equations reveals that the mean density gradient term changes sign at sufficiently large scales, such that buoyancy acts as a source for velocity gradients at small scales, but as a sink at large scales.

We systematically study the dissipative anomaly in compressible magnetohydrodynamic (MHD) turbulence using direct numerical simulations, and show that the total dissipation remains finite as viscosity diminishes. The dimensionless dissipation rate $\mathcal {C}_{\varepsilon }$ fits well with the model $\mathcal {C}_{\varepsilon } = \mathcal {C}_{\varepsilon,\infty } + \mathcal {D}/R_L^-$ for all levels of flow compressibility considered here, where $R_L^-$ is the generalized large-scale Reynolds number. The asymptotic value $\mathcal {C}_{\varepsilon,\infty }$ describes the total energy transfer flux, and decreases with increase of the flow compressibility, indicating non-universality of the dimensionless dissipation rate in compressible MHD turbulence. After introducing an empirically modified dissipation rate, the data from compressible cases collapse to a form similar to the incompressible MHD case depending only on the modified Reynolds number.

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’.

Chiral fluids – such as fluids under rotation or a magnetic field as well as synthetic and biological active fluids – flow in a different way than ordinary ones. Due to symmetries broken at the microscopic level, chiral fluids may have asymmetric stress and viscosity tensors, for example giving rise to a hydrostatic torque or non-dissipative (odd) and parity-violating viscosities. In this article, we investigate the motion of rigid bodies in such an anisotropic fluid in the incompressible Stokes regime through the mobility matrix, which encodes the response of a solid body to forces and torques. We demonstrate how the form of the mobility matrix, which is usually determined by particle geometry, can be analogously controlled by the symmetries of the fluid. By computing the mobility matrix for simple shapes in a three-dimensional (3-D) anisotropic chiral fluid, we predict counterintuitive phenomena such as motion at an angle to the direction of applied forces and spinning under the force of gravity.