A classic lift decomposition (Von Kármán & Sears, J. Aeronaut. Sci., vol. 5, 1938, pp. 379–390) is conducted on potential flow simulations of a near-ground pitching hydrofoil. It is discovered that previously observed stable and unstable equilibrium altitudes are generated by a balance between positive wake-induced lift and negative quasi-steady lift while the added mass lift does not play a role. Using both simulations and experiments, detailed analyses of each lift component's near-ground behaviour provide further physical insights. When applied to three-dimensional pitching hydrofoils the lift decomposition reveals that the disappearance of equilibrium altitudes for ${A{\kern-4pt}R}\ {\rm (aspect\ ratio)} <1.5$ occurs due to the magnitude of the quasi-steady lift outweighing the magnitude of the wake-induced lift at all ground distances. Scaling laws for the quasi-steady lift, wake-induced lift and the stable equilibrium altitude are discovered. A simple scaling law for the lift of a steady foil in ground effect is derived. This scaling shows that both circulation enhancement and the velocity induced at a foil's leading edge by the bound vortex of its ground image foil are the essential physics to understand steady ground effect. The scaling laws for unsteady pitching foils can predict the equilibrium altitude to within $20\,\%$ of its value when $St\ {\rm (Strouhal\ number)} < 0.45$. For $St \ge 0.45$ there is a wake instability effect, not accounted for in the scaling relations, that significantly alters the wake-induced lift. These results not only provide key physical insights and scaling laws for steady and unsteady ground effects, but also for two schooling hydrofoils in a side-by-side formation with an out-of-phase synchronization.

A computational study is presented on cross-stream migration and focusing of deformable capsules in curved microchannels of square and rectangular sections under inertial and non-inertial regimes. The numerical methodology is based on immersed boundary methods for fluid–structure coupling, a finite-volume-based flow solver and finite-element method for capsule deformation. Different focusing behaviours in the two regimes are predicted that arise due to the interplay of inertia, deformation, altered shear gradient, streamline curvature effect and secondary flow. In the non-inertial regime, a single-point focusing occurs on the central plane, and at a radial location between the interior face (i.e. face with highest curvature) of the channel and the location of zero shear. The focusing position is nearly independent of capsule deformability (represented by the capillary number, $Ca$). A two-step migration is observed that is comprised of a faster radial migration, followed by a slower migration toward the centre plane. The focusing location progressively moves further toward the interior face with increasing curvature and width, but decreasing height. In the inertial regime, single-point focusing is observed near the interior face for channel Reynolds number $Re_{C}\sim {O}(1)$, that is also highly sensitive to $Re_{C}$ and $Ca$, and moves progressively toward the exterior face with increasing $Re_{C}$ but decreasing $Ca$. As $Re_{C}$ increases by an order, secondary flow becomes stronger, and two focusing locations appear close to the centres of the Dean vortices. This location becomes practically independent of $Ca$ at even higher inertia. The inertial focusing positions move progressively toward the exterior face with increasing channel width and decreasing height. For wider channels, the equilibrium location is further toward the exterior face than the vortex centre.

Many key environmental, industrial and energy processes rely on controlling fluid transport within subsurface porous media. These media are typically structurally heterogeneous, often with vertically layered strata of distinct permeabilities – leading to uneven partitioning of flow across strata, which can be undesirable. Here, using direct in situ visualization, we demonstrate that polymer additives can homogenize this flow by inducing a purely elastic flow instability that generates random spatio-temporal fluctuations and excess flow resistance in individual strata. In particular, we find that this instability arises at smaller imposed flow rates in higher-permeability strata, diverting flow towards lower-permeability strata and helping to homogenize the flow. Guided by the experiments, we develop a parallel-resistor model that quantitatively predicts the flow rate at which this homogenization is optimized for a given stratified medium. Thus, our work provides a new approach to homogenizing fluid and passive scalar transport in heterogeneous porous media.

The subtle but crucial effects of large-scale forcing on the small-scale velocity-gradient (VG) dynamics is examined using direct numerical simulation (DNS) data of incompressible turbulence. The interplay among large-scale forcing, inertia, pressure and viscous effects is characterised as a function of local streamline geometry and VG magnitude (Frobenius norm). When conditioned on local topology, forcing: (i) counteracts inertial and viscous action in the strain-dominated nodal topologies; and (ii) balances pressure action in the rotation-dominated unstable focal topologies. Unexpectedly, forcing acts to reduce VG magnitudes in regions of strong dissipation. In these regions, forcing balances the non-local pressure effects whereas viscous action offsets the nonlinear inertial effects. In regions of very low dissipation, forcing combines with inertia and pressure effects to offset viscous action. With regard to the probability distribution of the normalised VG invariants, the primary role of forcing is to nullify certain features (dilatational probability currents) of inertia, pressure and viscous action. This results in the emergence of universal statistical features (solenoidal probability currents) that are determined largely by inertia–pressure–viscous interactions. These findings serve to enhance our understanding of small-scale processes in turbulence and guide VG model development.

In this study, we focus on Langmuir turbulence in the deep ocean with the presence of a large macroalgal farm using a large eddy simulation method. The wave–current interactions are modelled by solving the wave-averaged equations. The hydrodynamic process over the farm is found to drive a persistent flow pattern similar to Langmuir circulations but is locked in space across the farm. These secondary circulations are generated because the cross-stream shear produced by the rows of canopy elements leads to a steady vertical vorticity field, which is then rotated to the downstream direction under the effect of vortex force. Since the driving mechanism is similar to the Craik–Leibovich type 2 instability theory, these secondary circulations are also termed as attached Langmuir circulations. We then apply a triple decomposition on the flow field to unveil the underlying kinematics and energy transfer between the mean flow, the secondary flow resulting from the farm drag and the transient eddies. Flow visualizations and statistics suggest that the attached Langmuir circulations result from the adjustment of the upper ocean mixed layer to the macroalgal farm, and they will weaken (if not disappear) when the flow reaches an equilibrium state within the farm. The triple-decomposed energy budgets reveal that the energy of the secondary flow is transferred from the mean flow under the action of canopy drag, while the transient eddies feed on wave energy transferred by the Stokes drift and energy conversion from the secondary flow.

Turbulent mixing layers in nature are often characterised by the presence of a mean shear and an unstable buoyancy gradient between two streams of different velocities. Depending on the relative strength of shear versus buoyancy, either the former or the latter may dominate the turbulence and mixing between the two streams. In this paper, we present a phenomenological theory that leads to the identification of two distinct turbulent regimes: an early regime, dominated by mean shear, and a later regime dominated by buoyancy. The main theoretical result consists of the identification of a cross-over timescale that distinguishes between the shear- and the buoyancy-dominated turbulence. This cross-over time depends on three large-scale constants of the flow, namely, the buoyancy difference, the velocity difference between the two streams and the gravitational acceleration. We validate our theory against direct numerical simulations of a temporal turbulent mixing layer compounded with an unstable stratification. We observe that the cross-over time correctly predicts the transition from shear- to buoyancy-driven turbulence, in terms of turbulent kinetic energy production, energy spectra scaling and mixing layer thickness.

Particle–wall interactions have broad biological and technological applications. In particular, some artificial microswimmers capitalize on their translation–rotation coupling near a wall to generate directed propulsion. Emerging biomedical applications of these microswimmers in complex biological fluids prompt questions on the impact of non-Newtonian rheology on their propulsion. In this work, we report some intriguing effects of shear-thinning rheology, a ubiquitous non-Newtonian behaviour of biological fluids, on the translation–rotation coupling of a particle near a wall. One particularly interesting feature revealed here is that the wall-induced translation by rotation can occur in a direction opposite to what might be intuitively expected for an object rolling on a solid substrate. We elucidate the underlying physical mechanism and discuss its implications on the design of micromachines and bacterial motion near walls in complex fluids.

The propagation of a normal shock wave along a coupled convex–concave surface of equal radii has been analysed experimentally and numerically in this study. The experimental and numerical studies were conducted using a similar geometry as of that used by Ram et al. (J. Fluid Mech., vol. 768, 2015, pp. 219–239) for studying the shock wave transition from regular reflection to Mach reflection. Many interesting flow features such as shock wave transitions over the ramp, characteristics of the induced flow behind the shock wave and the development of a stationary separation shock wave have been observed in the study. The numerical results are validated with experimental data. While the shock wave transitions over the ramp are found to depend mainly on the ramp geometry, the characteristics of the stationary shock wave and the flow separation in the concave region of the ramp surface have been found to vary with the shock wave Mach numbers.

Direct numerical simulation of spanwise-rotation-driven flow transitions in viscoelastic plane Couette flow from a drag-reduced inertial to a drag-enhanced elasto-inertial turbulent flow state followed by full relaminarization is reported for the first time. Specifically, this novel flow transition begins with a drag-reduced inertial turbulent flow state at a low rotation number $0\leqslant Ro \leqslant 0.1$, and then transitions to a rotation/polymer-additive-driven drag-enhanced inertial turbulent regime, $0.1\leqslant Ro \leqslant 0.3$. In turn, the flow transitions to a drag-enhanced elasto-inertial turbulent state, $0.3\leqslant Ro \leqslant 0.9$, and eventually relaminarizes at $Ro=1$. In addition, two novel rotation-dependent drag enhancement mechanisms are proposed and substantiated. (1) The formation of large-scale roll cells results in enhanced convective momentum transport along with significant polymer elongation and stress generated in the extensionally dominated flow between adjacent roll cells at $Ro\leqslant 0.2$. (2) Coriolis-force-generated turbulent vortices cause strong incoherent transport and homogenization of significant polymer stress in the bulk via their vortical circulations at $Ro=0.5 - 0.9$.

The problem of interaction of a uniform current with a submerged horizontal circular cylinder in an ice-covered channel is considered. The fluid flow is described by linearized velocity potential theory and the ice sheet is treated as a thin elastic plate. The potential due to a source or the Green function satisfying all boundary conditions apart from that on the body surface is first derived. This can be used to derive the boundary integral equation for a body of arbitrary shape. It can also be used to obtain the solution due to multipoles by differentiating the Green function with its position directly. For a transverse circular cylinder, through distributing multipoles along its centre line, the velocity potential can be written in an infinite series with unknown coefficients, which can be determined from the impermeable condition on a body surface. A major feature here is that different from the free surface problem, or a channel without the ice sheet cover, this problem is fully three-dimensional because of the constraints along the intersection of the ice sheet with the channel wall. It has been also confirmed that there is an infinite number of critical speeds. Whenever the current speed passes a critical value, the force on the body and wave pattern change rapidly, and two more wave components are generated at the far-field. Extensive results are provided for hydroelastic waves and hydrodynamic forces when the ice sheet is under different edge conditions, and the insight of their physical features is discussed.