Drag for wall-bounded flows is directly related to the spatial flux of spanwise vorticity outward from the wall. In turbulent flows a key contribution to this wall-normal flux arises from nonlinear advection and stretching of vorticity, interpretable as a cascade. We study this process using numerical simulation data of turbulent channel flow at friction Reynolds number $Re_\tau =1000$. The net transfer from the wall of spanwise vorticity created by downstream pressure drop is due to two large opposing fluxes, one which is ‘down-gradient’ or outward from the wall, where most vorticity concentrates, and the other which is ‘up-gradient’ or toward the wall and acting against strong viscous diffusion in the near-wall region. We present evidence that the up-gradient/down-gradient transport occurs by a mechanism of correlated inflow/outflow and spanwise vortex stretching/contraction that was proposed by Lighthill. This mechanism is essentially Lagrangian, but we explicate its relation to the Eulerian anti-symmetric vorticity flux tensor. As evidence for the mechanism, we study (i) statistical correlations of the wall-normal velocity and of wall-normal flux of spanwise vorticity, (ii) vorticity flux cospectra identifying eddies involved in nonlinear vorticity transport in the two opposing directions and (iii) visualizations of coherent vortex structures which contribute to the transport. The ‘D-type’ vortices contributing to down-gradient transport in the log layer are found to be attached, hairpin-type vortices. However, the ‘U-type’ vortices contributing to up-gradient transport are detached, wall-parallel, pancake-shaped vortices with strong spanwise vorticity, as expected by Lighthill's mechanism. We discuss modifications to the attached eddy model and implications for turbulent drag reduction.

A long-wave asymptotic model is developed for the flow of an axisymmetric viscous film lining the interior of a tube for the case where slip occurs at the tube wall. Both the case of a falling film with a passive air core and that of a film driven up the tube by pressure-driven airflow are considered. The impact of slip on the net liquid volume flux is discussed, and linear stability analysis of the evolution equation is conducted to identify the impact of slip on the phase speed and growth rates of disturbances in each case. The presence of slip enhances the growth rates, though its impact on phase speed depends on the film thickness and the strength of the core airflow. For some parameter combinations, slip can modify the phase speed without altering the base flow. The nonlinear evolution of the free surface is then studied numerically. For falling films, increasing the slip length reduces the critical thickness required for plug formation to occur. Families of travelling wave solutions are found via continuation and are used to derive a simple formula for the dependence of this critical thickness on the slip length; this formula is shown to hold for small slip length. For air-driven films, the topology of streamlines in the film can be altered by slip at the wall; if the slip length is large enough, it can prevent regions of recirculation from forming at the wave crest.

We study the effect of continuum breakdown on the two-dimensional thermoacoustic radiation of a thin plate set in a perfect monatomic gas. The plate is heated harmonically in time and the acoustic field is investigated in the entire range of gas rarefaction rates. Analytical approximations are obtained in the limits of high (free-molecular) and low (continuum limit) gas rarefaction, accompanied by direct simulation Monte Carlo calculations at intermediate flow conditions. While the source acoustic field is of a monopole type in the continuum regime, it turns non-isotropic at non-continuum conditions, exhibiting dipole directivity and exponential decay rate in the collisionless limit. The combined effects of source heating frequency and gas Knudsen number on the far-field acoustic radiation are illustrated and rationalized.

This study investigates the effect of vibration on the flow structure transitions in thermal vibrational convection (TVC) systems, which occur when a fluid layer with a temperature gradient is excited by vibration. Direct numerical simulation (DNS) of TVC in a two-dimensional enclosed square box is performed over a range of dimensionless vibration amplitudes $0.001 \le a \le 0.3$ and angular frequencies $10^{2} \le \omega \le 10^{7}$, with a fixed Prandtl number of 4.38. The flow visualisation shows the transition behaviour of flow structure upon the varying frequency, characterising three distinct regimes, which are the periodic-circulation regime, columnar regime and columnar-broken regime. Different statistical properties are distinguished from the temperature and velocity fluctuations at the boundary layer and mid-height. Upon transition into the columnar regime, columnar thermal coherent structures are formed, in contrast to the periodic oscillating circulation. These columns are contributed by the merging of thermal plumes near the boundary layer, and the resultant thermal updrafts remain at almost fixed lateral position, leading to a decrease in fluctuations. We further find that the critical point of this transition can be described nicely by the vibrational Rayleigh number ${{Ra}}_{vib}$. As the frequency continues to increase, entering the so-called columnar-broken regime, the columnar structures are broken, and eventually the flow state becomes a large-scale circulation (LSC), characterised by a sudden increase in fluctuations. Finally, a phase diagram is constructed to summarise the flow structure transition over a wide range of vibration amplitude and frequency parameters.

The general problem of tracer diffusion in non-equilibrium baths is important in a wide range of systems, from the cellular level to geographical length scales. In this paper, we revisit the archetypical example of such a system: a collection of small passive particles immersed in a dilute suspension of non-interacting dipolar microswimmers, representing bacteria or algae. In particular, we consider the interplay between thermal (Brownian) diffusion and hydrodynamic (active) diffusion due to the persistent advection of tracers by microswimmer flow fields. Previously, it has been argued that even a moderate amount of Brownian diffusion is sufficient to significantly reduce the persistence time of tracer advection, leading to a significantly reduced value of the effective active diffusion coefficient $D_A$ compared to the non-Brownian case. Here, we show by large-scale simulations and kinetic theory that this effect is in fact practically relevant only for microswimmers that effectively remain stationary while still stirring up the surrounding fluid – so-called shakers. In contrast, for moderate and high values of the swimming speed, relevant for biological microswimmer suspensions, the effect of Brownian motion on $D_A$ is negligible, leading to the effects of advection by microswimmers and Brownian motion being additive. This conclusion contrasts with previous results from the literature, and encourages a reinterpretation of recent experimental measurements of $D_A$ for tracer particles of varying size in bacterial suspensions.

Two-way momentum-coupled direct numerical simulations of a particle-laden turbulent channel flow are addressed to investigate the effect of the particle Stokes number and of the particle-to-fluid density ratio on the turbulence modification. The exact regularised point-particle method is used to model the interphase momentum exchange in presence of solid boundaries, allowing the exploration of an extensive region of the parameter space. Results show that the particles increase the friction drag in the parameter space region considered, namely the Stokes number $St_+ \in [2,80]$, and the particle-to-fluid density ratio $\rho _p/\rho _f \in [90,5760]$ at a fixed mass loading $\phi =0.4$. It is noteworthy that the highest drag occurs for small Stokes number particles. A measurable drag increase occurs for all particle-to-fluid density ratios, the effect being reduced significantly only at the highest value of $\rho _p/\rho _f$. The modified stress budget and turbulent kinetic energy equation provide the rationale behind the observed behaviour. The particles’ extra stress causes an additional momentum flux towards the wall that modifies the structure of the buffer and of the viscous sublayer where the streamwise and wall-normal velocity fluctuations are increased. Indeed, in the viscous sublayer, additional turbulent kinetic energy is produced by the particles’ back-reaction, resulting in a strong augmentation of the spatial energy flux towards the wall where the energy is ultimately dissipated. This behaviour explains the increase of friction drag in particle-laden wall-bounded flows.

Wave turbulence is the study of the long-time statistical behaviour of equations describing a set of weakly nonlinear interacting waves. Such a theory, which has a natural asymptotic closure, allows us to probe the nature of turbulence more deeply than the exact Kolmogorov laws by rigorously proving the direction of the cascade and the existence of an inertial range, predicting stationary spectra for conserved quantities, or evaluating the Kolmogorov constant. An emblematic example is given by fast rotating fluids for which a wave turbulence theory has been derived by Galtier (Phys. Rev. E, vol. 68, issue 1, 2003, p. 015301). This work involves non-trivial analytical developments for a problem that is anisotropic by nature. We propose here a new path for the derivation of the kinetic equation by using the anisotropy at the beginning of the analysis. We show that the helicity basis is not necessary to obtain the wave amplitude equation for the canonical variables that involve a combination of poloidal and toroidal fields. The multiple time scale method adapted to this anisotropic problem is then used to derive the kinetic equation that is the same as the original work when anisotropy is eventually taken into account. This result proves the commutativity between asymptotic closure and anisotropy. In addition, the multiple time scale method informs us that the kinetic equation can be derived without imposing restrictions on the probability distribution of the wave amplitude such as quasi-Gaussianity, or on the phase such as random phase approximation that naturally occurs dynamically.

Understanding the skin friction in an axisymmetric turbulent boundary layer (ATBL) flow is a key to designing and optimising the flow past axisymmetric bodies, for example, a rocket engine nozzle and a submarine hull. In this study, we propose a universal law of the skin-friction coefficient in an ATBL flow. The flow is steady and fully developed with a zero pressure gradient. The governing equation for the ATBL flow is derived by methodically applying the boundary layer approximation. Subsequently, the scaling law of the Reynolds shear stress, caused by turbulent eddies, at the surface tangential to the wall roughness summits is derived by incorporating the role of transverse curvature. The skin-friction coefficient in a smooth ATBL flow is found to depend on two parameters, namely, the Reynolds number based on the cylinder radius, Rea, and the ratio of boundary layer thickness to cylinder radius, δ/a. The analysis predicts a simple form of the skin-friction coefficient as ${C_f} = 4.56 \times {10^{ - 2}}{[R{e_a}\ln (1 + \delta /a)]^{ - 1/4}}$, which agrees satisfactorily with the available experimental data and the numerical simulations in all the axisymmetric flow regimes. The proposed law, in the limit of infinite radius, is consistent with the classical law of the skin-friction coefficient in a plane turbulent boundary layer flow as ${C_{f0}}\sim Re_\delta ^{ - 1/4}$, where Reδ is the Reynolds number based on the boundary layer thickness. This study reveals that, for δ/a < 1, the relative skin-friction coefficient, $({C_f} - {C_{f0}})/{C_{f0}}$, follows a linear law with δ$/$a.

Modal and non-modal linear stability analyses are employed to investigate the effect of internal and external heating on disturbance temporal growth for the Darcy–Bénard convection with throughflow. A matrix-forming approach is employed for both purposes, where the generalised eigenvalue problem is built using the generalised integral transform technique. Although the disturbance equations are not self-adjoint, the non-modal analysis indicates that there is no transient growth. Hence, any disturbance growth in time must be induced by modal mechanisms. An absolute instability analysis reveals that viscous dissipation has a destabilising effect and introduces new modes that are eventually destabilised by increasing the Péclet number. Beyond critical values of the Péclet number, where codimension-two absolutely unstable points exist, these modes become more unstable than the classical mode found in the absence of viscous dissipation, which is stabilised by an increasing Péclet number. This internal heating mechanism generated by viscous dissipation is so strong at high enough Péclet numbers that instability becomes possible through heating from above.