We consider an internally heated fluid between parallel plates with fixed thermal fluxes. For a large class of heat sources that vary in the direction of gravity, we prove that $\smash { \smash {{\langle {\delta T} \rangle _h}} } \geq \sigma R^{-1/3} - \mu$, where $\smash { \smash {{\langle {\delta T} \rangle _h}} }$ is the average temperature difference between the bottom and top plates, $R$ is a ‘flux’ Rayleigh number and the constants $\sigma,\mu >0$ depend on the geometric properties of the internal heating. This result implies that mean downward conduction (for which $\smash { \smash {{\langle {\delta T} \rangle _h}} }< 0$) is impossible for a range of Rayleigh numbers smaller than a critical value $R_0:=(\sigma /\mu )^{3}$. The bound demonstrates that $R_0$ depends on the heating distribution and can be made arbitrarily large by concentrating the heating near the bottom plate. However, for any given fixed heating profile of the class we consider, the corresponding value of $R_0$ is always finite. This points to a fundamental difference between internally heated convection and its limiting case of Rayleigh–Bénard convection with fixed-flux boundary conditions, for which $\smash {{\langle {\delta T} \rangle _h}}$ is known to be positive for all $R$.

Two-dimensional free-surface flow over localised topography is examined, with the emphasis on the stability of hydraulic-fall solutions. A Gaussian topography profile is assumed with a positive or negative amplitude modelling a bump or a dip, respectively. Steady hydraulic-fall solutions to the full incompressible, irrotational Euler equations are computed, and their linear and nonlinear stability is analysed by computing eigenspectra of the pertinent linearised operator and by solving an initial value problem. The computations are carried out numerically using a specially developed computational framework based on the finite-element method. The Hamiltonian structure of the problem is demonstrated, and stability is determined by computing eigenspectra of the pertinent linearised operator. It is found that a hydraulic-fall flow over a bump is spectrally stable. The corresponding flow over a dip is found to be linearly unstable. In the latter case, time-dependent simulations show that ultimately, the flow settles into a time-periodic motion that corresponds to an invariant solution in an appropriately defined phase space. Physically, the solution consists of a localised large-amplitude wave that pulsates above the dip while simultaneously emitting nonlinear cnoidal waves in the upstream direction and multi-harmonic linear waves in the downstream direction.

Consider the motion of a thin layer of electrically conducting fluid, between two closely spaced parallel plates, in a classical Hele-Shaw geometry. Furthermore, let the system be immersed in a uniform external magnetic field (normal to the plates) and let electrical current be driven between conducting probes immersed in the fluid layer. In the present paper, we analyse the ensuing fluid flow at low Hartmann numbers. Physically, the system is particularly interesting because it allows for circulation in the flow, which is not possible in the standard pressure-driven Hele-Shaw cell. We first elucidate the mechanism of flow generation both physically and mathematically. After formulating the problem using complex variables, we present mathematical solutions for a class of canonical multiply connected geometries in terms of the prime function framework developed by Crowdy (Solving Problems in Multiply Connected Domains, SIAM, 2020). We then demonstrate how recently developed fast numerical methods may be applied to accurately determine the flow field in arbitrary geometries.

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.

Diffuse interface models are an important class of models used to describe multi-phase flows. In the case of incompressible viscous fluids there are a number of different diffuse interface models which have been known for many years. Nevertheless, a model based on mixture theory with a full set of conservation laws for the conservation of linear momentum of each constituent was not yet available. This gap was filled by ten Eikelder et al. (J. Fluid Mech., in press) recently and a first comparison with known models of Navier–Stokes/Cahn–Hilliard type is given. A detailed understanding of the relations between these models remains an important question.

Liquid metal buoyant flow around two differentially heated horizontal cylinders in the presence of a uniform vertical magnetic field is investigated experimentally. While magneto-convection in pipes or ducts has been studied theoretically and experimentally in recent years, data for heat transfer at immersed obstacles are rare and, to our knowledge, detailed experimental investigations on this fundamental magnetohydrodynamic problem do not exist. In the present work, two horizontal cylinders inserted into an adiabatic rectangular cavity filled with gallium–indium–tin are kept at constant temperatures to establish a driving temperature gradient in the surrounding liquid metal. The buoyancy-driven flow, quantified by the Grashof number $Gr$, is varied in the range ${10^{6} \leq Gr \leq ~5\times 10^{7}}$. With increasing magnetic field, expressed via the Hartmann number $Ha$, different flow regimes are identified from measurements for $0 \leq Ha \leq ~3000$. The effect of the electromagnetic force primarily consists in suppressing turbulence and damping the convective flow. The heat transfer is quantified in terms of the non-dimensional Nusselt number $Nu$, and its dependence on $Gr/{Ha}^{2}$, which is identified as the important group governing the flow, is discussed.

We perform linear stability analysis and direct numerical simulations to study the effect of the radius ratio on the instability and flow characteristics of the sheared annular centrifugal Rayleigh–Bénard convection, where the cold inner cylinder and the hot outer cylinder rotate with a small angular velocity difference. With the shear enhancement, the thermal convection is suppressed and finally becomes stable for different radius ratios $\{\eta \in \mathbb {R}|0.2\leqslant \eta \le 0.95\}$. Considering the inhomogeneous distribution of shear stresses in the base flow, a new global Richardson number $Ri_g$ is defined and the marginal-state curves for different radius ratios are successfully unified in the parameter domain of $Ri_g$ and the Rayleigh number $Ra$. The results are consistent with the marginal-state curve of the wall-sheared classical Rayleigh–Bénard convection in the streamwise direction, demonstrating that the basic stabilization mechanisms are identical. Moreover, systems with small radius ratios exhibit greater geometric asymmetry. On the one hand, this results in a smaller equivalent aspect ratio for the system, accommodating fewer convection roll pairs; fewer roll pairs are more likely to cause a transition in the flow structure during shear enhancement. On the other hand, the shear distribution is more inhomogeneous, allowing for an outward shift of the convection region and the elevation of bulk temperature under strong shear.

Single-flagellated bacteria are ubiquitous in nature. They exhibit various swimming modes using their flagella to explore complex surroundings such as soil and porous polymer networks. Some single-flagellated bacteria swim with two distinct modes, one with the flagellum extended away from its body and another with the flagellum wrapped around it. The wrapped mode has been observed when bacteria swim under tight confinements or in highly viscous polymeric melts. In this study we investigate the hydrodynamics of these two modes inside a circular pipe. We find that the wrapped mode is slower than the extended mode in bulk but more efficient under strong confinement due to a hydrodynamic increase of its flagellum translation–rotation coupling and an Archimedes’ screw-like configuration that helps to move the fluid along the pipe.

The practical implementation of machine learning in flow control is limited due to its significant training expenses. In the present study the convolutional neural network (CNN) trained with the data of the restricted nonlinear (RNL) model is used to predict the normal velocity on a detection plane at $y^+=10$ in a turbulent channel flow, and the predicted velocity is used as wall blowing and suction for drag reduction. An active control test is carried out by using the well-trained CNN in direct numerical simulation (DNS). Substantial drag reduction rates up to 19 % and 16 % are obtained based on the spanwise and streamwise wall shear stresses, respectively. Furthermore, we explore the online control of wall turbulence by combining the RNL model with reinforcement learning (RL). The RL is constructed to determine the optimal wall blowing and suction based on its observation of the wall shear stresses without using the label data on the detection plane for training. The controlling and training processes are conducted synchronously in a RNL flow field. The control strategy discovered by RL has similar drag reduction rates with those obtained previously by the established method. Also, the training cost decreases by over thirty times at $Re_{\tau }=950$ compared with the DNS-RL model. The present results provide a perspective that combining the RNL model with machine learning control for drag reduction in wall turbulence can be effective and computationally economical. Also, this approach can be easily extended to flows at higher Reynolds numbers.

The nature and behaviour of the drag coefficient $C_D$ of irregularly shaped grains within a wide range of Reynolds numbers $Re$ is discussed. The morphology of the grains is controlled by their fractal description, and they differ in shape. Using computational fluid dynamics tools, the characteristics of the boundary layer at high $Re$ has been determined by applying the Reynolds-averaged Navier–Stokes turbulence model. Both grid resolution and mesh size dependence are validated with well-reported previous experimental results applied in flow around isolated smooth spheres. The drag coefficient for irregularly shaped grains is shown to be higher than that for spherical shapes, also showing a strong drop in its value at high $Re$. This drag crisis is reported at lower $Re$ compared to the smooth sphere, but higher critical $C_D$, demonstrating that the morphology of the particle accelerates this crisis. Furthermore, the dependence of $C_D$ on $Re$ in this type of geometry can be represented qualitatively by four defined zones: subcritical, critical, supercritical and transcritical. The orientational dependence for both particles with respect to the fluid flow is analysed, where our findings show an interesting oscillatory behaviour of $C_D$ as a function of the angle of incidence, fitting the results to a sine-squared interpolation, predicted for particles within the Stokes laminar regime ($Re\ll 1$) and for elongated/flattened spheroids up to $Re=2000$. A statistical analysis shows that this system satisfies a Weibullian behaviour of the drag coefficient when random azimuthal and polar rotation angles are considered.

This study proposes a new mechanism that can lead to layering or convection from the finite amplitude perturbation acting on the double diffusive convection with uniform background shear. We focus on the double diffusive convection in the diffusive regime with the cold fresh water laying above the warm salty water. We demonstrate that, although the unperturbed system is linearly stable, the finite amplitude perturbation can trigger the initial flow motions which subsequently obtain energy from the gravitational potential energy and from the uniform background shear, and evolve to layering or convection. By using the linear stability analysis for the initial growth stage and the energy analysis for the following transitional stage, the critical Richardson number can be predicted theoretically. Here the Richardson number measures the relative strength of stratification to the background shear. The dominant wavenumbers and the growth rates of the corresponding modes given by linear theory agree well with the two-dimensional direct numerical simulations, and so does the critical Richardson number predicted by the theoretical model. The layering state is dominated by the double diffusion process, while the convection state at smaller Richardson number exhibits stronger influences from shear and generates smaller heat and salinity fluxes. The theoretical model is further applied to the parameter range which is relevant to the real oceanic environments and reveals that for the typical density ratio observed in the staircase regions in the Arctic Ocean, the current mechanism can lead to layering for relatively weak shear.

The scaling relations mapping the turbulence statistics in compressible turbulent boundary layers (TBLs) onto their incompressible counterparts are of fundamental significance for turbulence modelling, such as the Morkovin scaling for velocity fields, while for pressure fluctuation fields, a corresponding scaling relation is currently absent. In this work, the underlying scaling relations of pressure fluctuations about Mach number ($M$) contained in their generation mechanisms are explored by analysing a series of direct numerical simulation data of compressible TBLs over a wide Mach number range $(0.5\leq M \leq 8.0)$. Based on the governing equation of pressure fluctuations, they are decomposed into components according to the properties of source terms. It is notable that the intensity of the compressible component, predominantly originating from the acoustic mode, obeys a monotonic distribution about the Mach number and wall distance; further, the intensity of the rest of the pressure components, which are mainly generated by the vorticity mode, demonstrates a uniform distribution consistent with its incompressible counterpart. Moreover, the coupling between the two components is negligibly weak. Based on the scaling relations, semiempirical models for the fluctuation intensity of both pressure and its components are constructed. Hence, a mapping relation is obtained that the profiles of pressure fluctuation intensities in compressible TBLs can be mapped onto their incompressible counterparts by removing the contribution from the acoustic mode, which can be provided by the model. The intrinsic scaling relation can provide some basic insight for pressure fluctuation modelling.

We introduce a low-order dynamical system to describe thermal convection in an annular domain. The model derives systematically from a Fourier–Laurent truncation of the governing Navier–Stokes Boussinesq equations and accounts for spatial dependence of the flow and temperature fields. Comparison with fully resolved direct numerical simulations (DNS) shows that the model captures parameter bifurcations and reversals of the large-scale circulation (LSC), including states of (i) steady circulating flow, (ii) chaotic LSC reversals and (iii) periodic LSC reversals. Casting the system in terms of the fluid's angular momentum and centre of mass (CoM) reveals equivalence to a damped pendulum with forcing that raises the CoM above the fulcrum. This formulation offers a transparent mechanism for LSC reversals, namely the inertial overshoot of a forced pendulum, and it yields an explicit formula for the frequency $f^*$ of regular LSC reversals in the high-Rayleigh-number (Ra) limit. This formula is shown to be in excellent agreement with DNS and produces the scaling law $f^* \sim {Ra}^{0.5}$.

This study conducts particle-resolved direct numerical simulations to analyse how finite-size spherical particles affect the decay rate of turbulent kinetic energy in non-sustained homogeneous isotropic turbulence. The decaying particle-laden homogeneous isotropic turbulence is generated with two set-ups, i.e. (1) releasing particles into a single-phase decaying homogeneous isotropic turbulence and (2) switching off the driving force of a sustained particle-laden homogeneous isotropic turbulence. With both set-ups, the decay of turbulent kinetic energy follows a power-law when the flow is fully relaxed, similar to their single-phase counterparts. The dependence of the power-law decay exponent $n$ on the particle-to-fluid density ratio, particle size and volume fraction is also investigated, and a predictive model is developed. We find that the presence of heavier particles slows down the long-time power-law decay exponent.