Direct numerical simulations with two-way coupled Lagrangian tracking are carried out to study the bubble preferential concentration and the flow field modification. Simulations are conducted in an upward vertical turbulent channel driven by a constant pressure gradient, corresponding to a friction Reynolds number $Re_{\tau 0}=180$. Micro-sized bubbles with diameters ranging from 0.72 to 1.43 wall units are considered. Competition between lift force and wall-lift force in the wall-normal direction leads to significant near-wall bubble accumulation and directly results in distinct preferential concentration patterns across the channel. Below (above) the peak concentration height, the wall-lift (lift) force dominates, driving bubbles to accumulate in regions of high-speed sweep (low-speed ejection) events. In the vicinity of the wall, the wall-normal lift force exhibits a strong correlation with the local streamwise flow velocity, further reinforcing the preferential concentration of bubbles in high-speed regions. Additionally, bubbles show a strong preference for the low-enstrophy and high-dissipation nodal topologies. Furthermore, small bubbles primarily accumulate in the vicinity of the wall, reducing the work done on the flow and leading to a decrease in bulk velocity and turbulence statistics. In contrast, the turbulence statistics of large bubbles are nearly identical to those of the unladen flow. The impact of large bubbles on the flow field primarily manifests as an effective increase in the mean pressure gradient. These findings demonstrate that bubbles in the upward vertical channel flow exhibit strong preferential concentration behaviours, whereas their ability to modulate turbulence remains limited.

Exact mathematical expressions are derived to predict the exponent $p$ observed in non-equilibrium turbulence, where the classical dissipation law is replaced by a new dissipation scaling law $C_{\varepsilon } \sim \textit{Re}_{\lambda }^p$. Here, $ \textit{Re}_{\lambda }$ is the Taylor-based Reynolds number and $C_{\varepsilon } = \varepsilon L_{11} / u^{\prime 3}$ is the non-dimensional dissipation rate, defined by the viscous dissipation rate, $\varepsilon$, longitudinal integral scale, $L_{11}$, and root-mean-square of the velocity fluctuations $u^{\prime} = \sqrt {\overline {u^{\prime 2}}}$ (Vassilicos, Annu. Rev. Fluid Mech., vol. 47, 2015, pp. 95–114). Assuming homogeneous and isotropic turbulence, it is shown that the exact value of $p$ involves only first-order derivatives of these variables; however, at very high Reynolds numbers, and under particularly strong changes in the power input of the external forcing (without changing the shape of the forcing spectrum), the exact expression simplifies to $p = 3\pi / 4\alpha L_{110} - 5 / 2$, where $L_{110}$ is the initial value of the longitudinal integral scale and $\alpha$ represents an effective forcing wavenumber. Thus, the main finding is that only large-scale effects are involved in the imposition of the non-equilibrium dissipation scaling law. The results are compared with direct numerical simulation (DNS) results of isotropic turbulence under abruptly changing forcing conditions and with experimental data of non-equilibrium decaying isotropic turbulence, showing consistent results.

An oscillating body floating at the water surface produces a field of self-generated waves. When the oscillation induces a difference in fore–aft wave amplitude squared, these self-generated waves can be used as a mechanism to propel the body horizontally across the surface (Longuet-Higgins 1977 Proc. R. Soc. Lond. A, vol. 352, no. 1671, pp. 463–480). The optimisation of this wave-driven propulsion is the interest of this work. To study the conditions necessary to produce optimal thrust we will consider a shallow water set-up where a periodically oscillating pressure source acts as the body. In this framework, an expression for the thrust is derived by relation to the difference in fore–aft amplitude squared. The conditions on the source for maximal thrust are explored both analytically and numerically in two optimal control problems. The first case is where a bound is imposed on the norm of the control function to regularise it. Secondly, a more physically motivated case is studied where the power injected by the source is bounded. The body is permitted to have a drift velocity $U$. When scaled with the wave speed $c$, the dimensionless velocity $v=U/c$ divides the study into subcritical, critical and supercritical regimes and the optimal conditions are presented for each. The result in the bounded power case is then used to demonstrate how the modulation of power injected can slowly change the cruising velocity from rest to supercritical velocities.