Mixing-induced reactions play an important role in a wide range of porous media processes. Recent advances have shown that fluid flow through porous media leads to chaotic advection at the pore scale. However, how this impacts Darcy-scale reaction rates is unknown. Here, we measure the reaction rates in steady mixing fronts using a chemiluminescence reaction in index-matched three-dimensional porous media. We consider two common mixing scenarios for reacting species, flowing either in parallel in a uniform flow or towards each other in a converging flow. We study the reactive properties of these fronts for a range of Péclet numbers. In both scenarios, we find that the reaction rates significantly depart from the prediction of hydrodynamic dispersion models, which obey different scaling laws. We attribute this departure to incomplete mixing effects at the pore scale, and propose a mechanistic model describing the pore-scale deformations of the front triggered by chaotic advection and their impact on the reaction kinetics. The model shows good agreement with the effective Darcy-scale reaction kinetics observed in both uniform and converging flows, opening new perspectives for upscaling reactive transport in porous media.

We present new unconstrained simulations and constrained experiments of a pair of pitching hydrofoils in a leader–follower in-line arrangement. Free-swimming simulations with matched pitching amplitudes show self-organisation into stable formations at a constant gap distance without any control. Over a wide range of phase synchronisation, amplitude and Lighthill number typical of biology, we discover that the stable gap distance scales with the actual wake wavelength of an isolated foil rather than the nominal wake wavelength. A scaling law for the actual wake wavelength is derived and shown to collapse data across a wide Reynolds number range of $200 \leqslant Re \leqslant 59\,000$. Additionally, vortex analysis uncovers that the leader’s wake wavelength-to-chord ratio, $\lambda /c$, is the key dimensionless variable to maximise the follower’s/collective efficiency. When $\lambda /c \approx 2$ it ensures that the follower’s leading edge suction force and the net force from a nearby vortex pair act in the direction with the foil’s motion thereby reducing the follower’s power. Moreover, in both simulations and experiments mismatched foil amplitudes are discovered to increase the efficiency of hydrofoil schools by 70 % while maintaining a stable formation without closed-loop control. This occurs by (i) increasing the stable gap distance between foils to push them into a high-efficiency zone and (ii) raising the level of efficiency in these zones. This study bridges the gap between constrained and unconstrained studies of in-line schooling by showing that constrained-foil measurements can map out the potential efficiency benefits of schooling. These findings can aid in the design of high-efficiency biorobot schools.

Aerodynamic loads play a central role in many fluid dynamics applications, and we present a method for identifying the structures (or modes) in a flow that make dominant contributions to the time-varying aerodynamic loads in a flow. The method results from the combination of the force partitioning method (Menon & Mittal, 2021, J. Fluid Mech., vol. 907, A37) and modal decomposition techniques such as Reynolds decomposition, triple decomposition and proper orthogonal decomposition, and is applied here to three distinct flows – two-dimensional flows past a circular cylinder and an aerofoil, and the three-dimensional flow over a revolving rectangular wing. We show that the force partitioning method applied to modal decomposition of velocity fields results in complex, and difficult to interpret inter-modal interactions. We therefore propose and apply modal decomposition directly to the $Q$-field associated with these flows. The variable $Q$ is a nonlinear observable that is typically used to identify vortices in a flow, and we find that the direct decomposition of $Q$ leads to results that are more amenable to interpretation. We also demonstrate that this modal force partitioning can be extended to provide insights into the far-field aeroacoustic loading noise of these flows.

In this study, we explore the evolution of instabilities in magneto-quasi-geostrophic (MQG) modons on the $f$-plane using a magnetohydrodynamic rotating shallow water model. The numerical experiments have been conducted using a recently proposed second-order flux-globalisation-based path-conservative central-upwind scheme. Our focus is on the evolution and interactions of three key configurations: singular, regular and hollow MQG modons, which represent cases where the magnetic field is confined within the separatrix, evenly distributed inside and outside the separatrix and localised outside the separatrix, respectively. The singular MQG modon emerges as the most stable configuration, demonstrating the greatest resilience to destabilising forces. A notable observation is its transition from a quadrupolar to a tripolar magnetic field structure before reverting to a quadrupole adjusted magnetic modon, accompanied by a clockwise rotation of the system. In terms of stability, singular modons are the most stable ones, while hollow modons are the least stable. As instabilities develop, southward or northward displacements become significantly more pronounced than eastward or westward movements, primarily due to the Coriolis force. Among the configurations, the hollow (singular) modons experience the biggest (smallest) displacements. Additionally, we investigate modon collisions and highlight three scenarios: interactions between cyclonic and anticyclonic components that form a composite modon with meridional bifurcation; collisions of cyclonic vortices that produce a tripolar structure with counterclockwise rotation; and collisions between anticyclonic components that result in a stable, quasi-stationary tripolar configuration. The resulting magnetic poles exhibit a checkered pattern, with their amplitude decreasing with increasing distance from the central vortex.

Viscous flow through high-permeability channels occurs in many environmental and industrial applications, including carbon sequestration, groundwater flow and enhanced oil recovery. In this work, we study the displacement of a less-viscous fluid by a more-viscous fluid in a layered porous medium in a rectilinear configuration, where two low-permeability layers sandwich a higher-permeability layer. We derive a theoretical model that is validated using corroborative laboratory experiments, when the influence of the density difference is negligible. We find that the location of the propagating front increases with time according to a power-law form $x_f \propto t^{1/2}$, while the fluid–fluid interface exhibits a self-similar shape, when the motion of the displaced fluid is negligible in an unconfined porous medium. In the experimental set-up, distinct permeability layers were constructed using various sizes of spherical glass beads. The working fluids comprised fresh water as the less-viscous ambient fluid, and a glycerine–water mixture as the more-viscous injecting fluid. Our experimental measurement show a better match with the theory for the experiments performed at low Reynolds numbers and with permeable boundaries in the far field.

We study the homogeneous isotropic turbulence of a shear-thinning fluid modelled by the Carreau model, and show how the variable viscosity affects the multiscale behaviour of the turbulent flow. We show that Kolmogorov theory can be extended to such non-Newtonian fluids, provided that the correct choice of average is taken when defining the mean Kolmogorov scale and dissipation rate, to properly capture the effect of the variable viscosity. Thus the classical phenomenology à la Kolmogorov can be observed in the inertial range of scale, with the energy spectra decaying as $k^{-5/3}$, with $k$ being the wavenumber, and the third-order structure function obeying the $4/5$ law. The changing viscosity instead strongly alters the small scale of turbulence, leading to an enhanced intermittent behaviour of the velocity field.

A nonlinear Schrödinger equation for pure capillary waves propagating at the free surface of a vertically sheared current has been used to study the stability and bifurcation of capillary Stokes waves on arbitrary depth. A linear stability analysis of weakly nonlinear capillary Stokes waves on arbitrary depth has shown that (i) the growth rate of modulational instability increases as the vorticity decreases whatever the dispersive parameter $kh$, where $k$ is the carrier wavenumber and $h$ the depth; (ii) the growth rate is significantly amplified for shallow water depths; and (iii) the instability bandwidth widens as the vorticity decreases. Particular attention has been paid to damping due to viscosity and forcing effects on modulational instability. In addition, a linear stability analysis to transverse perturbations in deep water has been carried out, demonstrating that the dominant modulational instability is two-dimensional whatever the vorticity. Near the minimum of linear phase velocity in deep water, we have shown that generalised capillary solitary waves bifurcate from linear capillary Stokes waves when the vorticity is positive. Moreover, we have shown that the envelope of pure capillary waves in deep water is unstable to transverse perturbations. Consequently, deep-water generalised capillary solitary waves are expected to be unstable to transverse perturbations.

Görtler vortices induced by concave curvature in supersonic turbulent flows are investigated using resolvent analysis and large-eddy simulations at Mach 2.95 and Reynolds number $ Re_{\delta }=63\,500$ based on the boundary-layer thickness $ \delta$. Resolvent analysis reveals that the most amplified coherent structures manifest as streamwise counter-rotating vortices with optimal spanwise wavelength $ 2.4\delta$ at cut-off frequency $f\delta /{u}_{\infty } =0.036$, where $ {u}_{\infty }$ is the freestream velocity. The leading spectral proper orthogonal decomposition modes with spanwise wavelength approximately $ 2\delta$ align well with the predicted coherent structures from resolvent analysis at $f\delta /{u}_{\infty } =0.036$. These predicted and extracted coherent structures are identified as Görtler vortices, driven by the Görtler instability. The preferential spanwise scale of the Görtler vortices is further examined under varying geometric and freestream parameters. The optimal spanwise wavelength is insensitive to the total turning angle beyond a critical value, but sensitive to the concave curvature $ K$ at the same turning angle. A limit spanwise wavelength $ 1.96\delta$, corresponding to an infinite concave curvature as $ K\rightarrow \infty$, is identified and validated. Increasing the freestream Mach number or decreasing the ratio of wall temperature to freestream temperature reduces the optimal wavelength normalised by $ \delta$, while variations in freestream Reynolds number have negligible impact. Additionally, a modified definition of the turbulent Görtler number $ G_{T}$ based on the peak eddy viscosity in boundary layers is proposed and employed to assess the occurrence of Görtler instability.