Numerical simulation is used to document the statistical structure and better understand energy transfers in a low-Reynolds-number turbulent flow generated by negative axial buoyancy in a long circular tilted pipe under the Boussinesq approximation. The flow is found to exhibit specific features which strikingly contrast with the familiar characteristics of pressure-driven pipe and channel flows. The mean flow, dominated by an axial component exhibiting a uniform shear in the core, also comprises a weak secondary component made of four counter-rotating cells filling the entire cross-section. Within the cross-section, variations of the axial and transverse velocity fluctuations are markedly different, the former reaching its maximum at the edge of the core while the latter two decrease monotonically from the axis to the wall. The negative axial buoyancy component generates long plumes travelling along the pipe, yielding unusually large longitudinal integral length scales. The axial and crosswise mean density variations are shown to be respectively responsible for a quadratic variation of the crosswise shear stress and density flux which both decrease from a maximum on the pipe axis to near-zero values throughout the near-wall region. Although the crosswise buoyancy component is stabilizing everywhere, the crosswise density flux is negative in some peripheral regions, which corresponds to apparent counter-gradient diffusion. Budgets of velocity and density fluctuations variances and of crosswise shear stress and density flux are analysed to explain the above features. A novel two-time algebraic model of the turbulent fluxes is introduced to determine all components of the diffusivity tensor, revealing that they are significantly influenced by axial and crosswise buoyancy effects. The eddy viscosity and eddy diffusivity concepts and the Reynolds analogy are found to work reasonably well within the central part of the section whereas non-local effects cannot be ignored elsewhere.

The three-dimensional flow around two spherical bubbles moving in a viscous fluid is studied numerically by solving the full Navier–Stokes equations. The study considers the interaction between two bubbles for moderate Reynolds numbers (50 ≤ Re ≤ 500, Re being based on the bubble diameter) and for positions described by the separation S (2.5 ≤ S ≤ 10, S being the distance between the bubble centres normalised by the bubble radius) and the angle θ (0° ≤ θ ≤ 90°) formed between the centreline and the direction perpendicular to the direction of the motion. We provide a general description of the interaction extending the results obtained for two bubbles moving side by side (θ = 0°) by Legendre, Magnaudet & Mougin (J. Fluid Mech., vol. 497, 2003, p. 133) and for two bubbles moving in line (θ = 90°) by Yuan & Prosperetti J. Fluid Mech., vol. 278, 1994, p. 325). Simple models based on physical arguments are given for the drag and lift forces experienced by each bubble. The interaction is the combination of three effects: a potential effect, a viscous correction (Moore's correction) and a significant wake effect observed on both the drag and the transverse forces of the second bubble when located in the wake of the first one.

We study numerically the viscous phase of horizontal gravity currents of immiscible fluids in the lock-exchange configuration. A numerical technique capable of dealing with stiff density gradients is used, allowing us to mimic high-Schmidt-number situations similar to those encountered in most laboratory experiments. Plane two-dimensional computations with no-slip boundary conditions are run so as to compare numerical predictions with the ‘short reservoir’ solution of Huppert (J. Fluid Mech., vol. 121, 1982, pp. 43–58), which predicts the front position lf to evolve as t1/5, and the ‘long reservoir’ solution of Gratton & Minotti (J. Fluid Mech., vol. 210, 1990, pp. 155–182) which predicts a diffusive evolution of the distance travelled by the front xf ~ t1/2. In line with dimensional arguments, computations indicate that the self-similar power law governing the front position is selected by the flow Reynolds number and the initial volume of the released heavy fluid. We derive and validate a criterion predicting which type of viscous regime immediately succeeds the slumping phase. The computations also reveal that, under certain conditions, two different viscous regimes may appear successively during the life of a given current. Effects of sidewalls are examined through three-dimensional computations and are found to affect the transition time between the slumping phase and the viscous regime. In the various situations we consider, we make use of a force balance to estimate the transition time at which the viscous regime sets in and show that the corresponding prediction compares well with the computational results.