High-Reynolds-number data in turbulent boundary layers are analysed to examine statistical moments of streamwise velocity fluctuations ${u}^{\prime } $ . Prior work has shown that the variance of ${u}^{\prime } $ exhibits logarithmic behaviour with distance to the surface, within an inertial sublayer. Here we extend these observations to even-order moments. We show that the $2p$ -order moments, raised to the power $1/ p, $ also follow logarithmic behaviour according to $\langle \mathop{({u}^{\prime + } ){}^{2p} \rangle }\nolimits ^{1/ p} = {B}_{p} - {A}_{p} \ln (z/ \delta )$ , where ${u}^{\prime + } $ is the velocity fluctuation normalized by the friction velocity, $\delta $ is an outer length scale and ${B}_{p} $ are non-universal constants. The slopes ${A}_{p} $ in the logarithmic region appear quite insensitive to Reynolds number, consistent with universal behaviour for wall-bounded flows. The slopes differ from predictions that assume Gaussian statistics, and instead are consistent with sub-Gaussian behaviour.

The stabilization of a swept-wing boundary layer by distributed surface roughness elements is studied by performing direct numerical simulations. The configuration resembles experiments studied by Saric and coworkers at Arizona State University, who employed this control method in order to delay transition. An array of cylindrical roughness elements are placed near the leading edge to excite subcritical cross-flow modes. Subcritical refers to the modes that are not critical with respect to transition. Their amplification to nonlinear amplitudes modifies the base flow such that the most unstable cross-flow mode and secondary instabilities are damped, resulting in downstream shift of the transition location. The experiments by Saric and coworkers were performed at low levels of free stream turbulence, and the boundary layer was therefore dominated by stationary cross-flow disturbances. Here, we consider a more complex disturbance field, which comprises both steady and unsteady instabilities of similar amplitudes. It is demonstrated that the control is robust with respect to complex disturbance fields as transition is shifted from 45 to 65 % chord.

We investigate numerically, in spherical geometry, the interaction of stratification with precession. Both stable stratification and unstable stratification are studied. In the parameter regime we are concerned with, stable stratification suppresses the precessional instability, whereas unstable stratification and precession can either stabilize or destabilize each other at different precession rates.

In this paper we study the self-propulsion of a dumbbell micro-robot submerged in a viscous fluid. The micro-robot consists of two rigid spherical beads connected by a rod or a spring; the rod/spring length changes periodically. The constant density of each sphere differs from the density of the fluid, while the whole micro-robot has neutral buoyancy. An effective oscillating gravity field is created via rigid-body oscillations of the fluid. Our calculations show that the micro-robot undertakes both translational and rotational motion. Using an asymptotic procedure containing a two-time method and a distinguished limit, we obtain analytic expressions for the averaged self-propulsion velocity and averaged angular velocity. The important special case of zero angular velocity represents rectilinear self-propulsion with constant velocity. In particular, we have shown that: (a) no unidirectional oscillations of a fluid result in self-propulsion; and (b) for the oscillations of a fluid in two directions rectilinear motion of a micro-robot can be achieved.

The breakup of droplets due to creeping motion in a confined microchannel geometry is studied using three-dimensional numerical simulations. Analogously to unconfined droplets, there exist two distinct breakup phases: (i) a quasi-steady droplet deformation driven by the externally applied flow; and (ii) a surface-tension-driven three-dimensional rapid pinching that is independent of the externally applied flow. In the first phase, the droplet relaxes back to its original shape if the externally applied flow stops; if the second phase is reached, the droplet will always break. Also analogously to unconfined droplets, there exist two distinct critical conditions: (i) one that determines whether the droplet reaches the second phase and breaks, or it reaches a steady shape and does not break; and (ii) one that determines when the rapid autonomous pinching starts. We analyse the second phase using stop–flow simulations, which reveal that the mechanism responsible for the autonomous breakup is similar to the end-pinching mechanism for unconfined droplets reported in the literature: the rapid pinching starts when, in the channel mid-plane, the curvature at the neck becomes larger than the curvature everywhere else. The same critical condition is observed in simulations in which we do not stop the flow: the breakup dynamics and the neck thickness corresponding to the crossover of curvatures are similar in both cases. This critical neck thickness depends strongly on the aspect ratio, and, unlike unconfined flows, depends only weakly on the capillary number and the viscosity contrast between the fluids inside and outside the droplet.

Liquid rings can be generated in the Leidenfrost state using liquid oxygen of low boiling point ( $- 183~\textdegree \mathrm{C} $ ) and high magnetic susceptibility, allowing one to ‘sculpt’ the liquid into a ring shape using an annular magnet. When the magnetic field is turned off, the ring shrinks back into a puddle with a constant acceleration. A potential flow approach accurately describes the dynamics of closure with an equation reminiscent of the Rayleigh–Plesset equation for the collapse of transient cavities.

Quantum turbulence in thermal counterflow of superfluid ${\text{} }^{4} \mathrm{He} $ is studied at length scales comparable to the mean distance $\ell $ between quantized vortices. The Lagrangian dynamics of solid deuterium particles, of radius ${R}_{p} $ about one order of magnitude smaller than $\ell $ , is analysed in a planar section of the experimental volume by using the particle tracking velocimetry technique. We show that the average amplitude of the acceleration of the particles seems to increase as the temperature decreases and applied heat flux increases and this can be explained by exploiting the two-fluid model of superfluid ${\text{} }^{4} \mathrm{He} $ . We also report that, at the probed length scales, the normalized distribution of the acceleration of the particles appears to follow an unexpected classical-like behaviour.

The lift-generating capabilities of a translating and passively pitching rectangular plate are assessed experimentally. The plate pitch dynamics are generated by a rotational spring, and a genetic algorithm isolates a set of spring parameters maximizing the average lift. Our experiments identified a range of parameters that produce kinematic trajectories associated with optimal lift production. The stroke length and the dynamic response of the spring at the driving frequency are revealed to play crucial roles in the generation of such trajectories. Measurements taken with digital particle image velocimetry are used to analyse the results.

The dynamics on the laminar–turbulent separatrix is investigated numerically for boundary-layer flows in the subcritical regime. Constant homogeneous suction is applied at the wall, resulting in a parallel asymptotic suction boundary layer (ASBL). When the numerical domain is sufficiently extended in the spanwise direction, the coherent structures found by edge tracking are invariably localized and their dynamics shows bursts that drive a remarkable regular or irregular spanwise dynamics. Depending on the parameters, the asymptotic dynamics on the edge can be either periodic in time or chaotic. A clear mechanism for the regeneration of streaks and streamwise vortices emerges in all cases and is investigated in detail.

We prove that periodic solutions of the steady Prandtl equations do not exist on a stationary boundary. On a moving boundary, there are no solutions with a monotone velocity profile.

The time evolution of pressure statistics in freely decaying homogeneous isotropic turbulence (HIT) is investigated via eddy-damped quasi-normal Markovian (EDQNM) computations. The present results show that the time decay rate of pressure-based statistical quantities, such as pressure variance and pressure gradient variance, are sensitive to the breakdown of permanence of large eddies. New formulae for the associated time-decay exponents are proposed, which extend previous relations proposed in Lesieur, Ossia & Metais (Phys. Fluids, vol. 11, 1999, p. 1535). Particular attention is paid to finite-Reynolds-number (FRN) effects on the pressure spectrum and pressure statistics. The results also suggest that $R{e}_{\lambda } = O(1{0}^{4} )$ must be considered to observe a one-decade inertial range in the pressure spectrum with Kolmogorov $- 7/ 3$ scaling. This threshold value is larger than almost all existing direct numerical simulation (DNS) and experimental data, justifying the discussion about other possible scaling laws. The $- 5/ 3$ slope reported in some DNS data is also recovered by the EDQNM model, but it is observed to be a low-Reynolds-number effect. Another important result is that FRN effects yield a departure from asymptotic theoretical behaviours which appear similar to some effects attributed to intermittency by most authors. This is exemplified by the ratio between pressure-based and velocity-based Taylor microscales. Therefore, high-Reynolds-number DNS or experiments such that $R{e}_{\lambda } = O(1{0}^{4} )$ would be required in order to remove FRN effects and to analyse pure intermittency effects.

The transient dynamics of the linearized Euler–Boussinesq equations governing parallel stratified shear flows is presented and analysed. Solutions are expressed as integral superpositions of generalized eigenfunctions associated with the continuous-spectrum component of the Taylor–Goldstein linear stability operator, and reveal intrinsic dynamics not captured by its discrete-spectrum counterpart. It is shown how continuous-spectrum perturbations are generally characterized by non-normal energy growth and decay with algebraic asymptotic behaviour in either time or space. This behaviour is captured by explicit long-time/far-field expressions from rigorous asymptotic analysis, and it is illustrated with direct numerical simulations of the whole (non-Boussinesq) stratified Euler system. These results can be helpful in understanding recent numerical observations for parallel and non-parallel perturbed stratified shear flows.

We consider the development of Dean vortices in a curved channel of finite aspect ratio. Solutions to the axisymmetric Navier–Stokes equations are obtained through a finite-element analysis, allowing us to explore the complex and rich bifurcation pattern of the flow as the aspect ratio and Dean number vary. We demonstrate a new class of finite-amplitude vortices and discuss their relationship to similar structures seen in finite-length Taylor–Couette flow.

We report on near-turbulent thermal convection of a nematic liquid crystal heated from below in a cylindrical cell with an aspect ratio (diameter/height) equal to 0.50 for Rayleigh numbers $2\times 1{0}^{7} \lesssim \mathit{Ra}\lesssim 3\times 1{0}^{8} $ and a Prandtl number of about 355. The Nusselt number $\mathit{Nu}$ as a function of $\mathit{Ra}$ did not differ significantly from that of an isotropic fluid. In a vertical magnetic field $\mathbi{H}$ , we found $\mathit{Nu}(H)/ \mathit{Nu}(0)= 1+ a(\mathit{Ra}){H}^{2} $ , with $a(\mathit{Ra})= 0. 24{\mathit{Ra}}^{0. 75} ~{\mathrm{G} }^{- 2} $ . We present a model that describes the $H$ dependence in terms of a change of the thermal conductivity in the thermal boundary layers due to a field-induced director alignment.

We investigate the shape of bubbles in liquid oxygen under magnetic levitation conditions: a magnetic field is applied that polarizes bulk oxygen, and its spatial variation induces a body force opposed to its weight. In these conditions, bubbles appear to have a smooth ellipsoidal shape, which may be prolate (elongated in the vertical direction), oblate (elongated in the horizontal plane) or perfectly spherical. The dependence of the elongation ratio $\eta $ on the volume and levitation position is explored. It is found that the bubble shape is prescribed by the minimization of the sum of surface tension, demagnetization and magnetic–gravitational potential energies.

We analyse axisymmetric gravity currents of power-law fluids theoretically and experimentally. We use aqueous suspensions of Xanthan gum in laboratory experiments of constant-volume and constant-flux release to resolve the rheological parameters of the fluid, which we then compare with measurements made using a strain-controlled rheometer. We find that the constant-volume release of highly shear-thinning fluids involves an early-time evolution dominated by inertia, and non-convex free surfaces that make the application of similarity solutions of the late-time viscously dominated evolution inefficient at resolving material properties. In contrast, constant-flux release of the same fluids can be viscously dominated and consistent with the self-similar solution from early in the evolution, which makes it a more useful method for measuring rheological parameters.

We examine pulse interaction and bound-state formation in interfacial turbulence using the problem of a falling liquid film as a model system. We perform direct numerical simulations (DNSs) of the full Navier–Stokes equations and associated wall and free-surface boundary conditions and we examine both analytically and numerically a low-dimensional (LD) model for the film. For a two-pulse system, DNSs reveal the existence of very rich and complex pulse interactions, characterized by either overdamped, underdamped or self-sustained oscillating behaviours, all of them found to be in excellent agreement with LD results. Having demonstrated the reliability of the LD model for two-pulse systems/smaller domains, we use it to investigate larger domains with many interacting pulses, where DNSs are computationally expensive. We demonstrate that such systems are likely to be dominated by a self-sustained oscillatory dynamics.

Using high-speed photography coupled with optical interference, we experimentally study the air entrapment during a liquid drop impacting a solid substrate. We observe the formation of a compressed air film before the liquid touches the substrate, with internal pressure considerably higher than the atmospheric value. The degree of compression highly depends on the impact velocity, as explained by balancing the liquid deceleration with the large pressure of the compressed air. After contact, the air film expands vertically at the edge, reducing its pressure within a few tens of microseconds and producing a thick rim on the perimeter. This thick-rimmed air film subsequently contracts into an air bubble, governed by the complex interaction between surface tension, inertia and viscous drag. Such a process is universally observed for impacts above a few centimetres high.

We report results of a combined numerical and experimental study on axisymmetric and non-axisymmetric flow states in a finite-length, co-rotating Taylor–Couette system in the Taylor vortex regime but also in the Rayleigh stable regime for moderate Reynolds numbers ( ${\leq }1000$ ). We found the dominant boundary-driven axisymmetric circulation to play a crucial role in the mode selection and the bifurcation behaviour in this flow. A sequence of partially hysteretic transitions to other axisymmetric multi-cell flow states is observed. Furthermore, we observed spiral states bifurcating via a supercritical Hopf bifurcation out of these multi-cell states which strongly determine the shape of the spiral. Finally, an excellent agreement between experimental and numerical results of the full Navier–Stokes equations is found.

Three-dimensional direct numerical simulation of horizontal convection is reported for a large Rayleigh number, $\mathit{Ra}\sim O(1{0}^{12} )$ , and boundary conditions that allow comparison with previous laboratory experiments. The convection is forced by heating over half of the horizontal base of a long channel and cooling over the other half of the base. The solutions are consistent with the experiments, including small-scale streamwise roll instability developing into a convectively mixed layer within the bottom thermal boundary layer, and a turbulent endwall plume. The mechanical energy budget is shown to be dominated by conversions of available potential energy to kinetic energy by buoyancy flux in the plume and the reverse in the interior of the circulation. These local conversions are three orders of magnitude greater than the total rate of viscous dissipation. The total irreversible mixing is exactly equal to the generation of available potential energy by buoyancy forcing, and one order of magnitude larger than the viscous dissipation. This confirms that dissipation rate is not an indicator of the strength of the circulation and explains why horizontal convection is more energetic than might be expected.