Turbulent pipe flow is of substantial importance in practical applications, and it remains challenging to depict the characteristic complex multiscale dynamics by a unified theoretical framework, hindered by its inherent intermittency. Inspired by a recent study of velocity circulation in turbulent channel flows from Duan, Chen & Sreenivasan (2025 J. Fluid Mech., vol. 1009, p. R4), in this study, we investigate the statistical characteristics of velocity circulation (or equally the area integral of wall-normal vorticity) over rectangular loops in concentric cylindrical shells, parallel to the pipe wall. The statistics are implemented using direct numerical simulation data at friction Reynolds numbers of $ \textit{Re}_\tau =1057$ and $2000$. Close to the pipe wall, the circulation in the inertial range resides on space-filling unifractal sets, with the Hölder exponent smaller than Kolmogorov’s $4/3$. Away from the pipe wall, the circulation displays bifractal characteristics and the Hölder exponents for high moment orders are very close to those reported in channel flows and homogeneous isotropic turbulence. The circulation statistics are only dependent on the area enclosed by the loops, and are invariant to the loop aspect ratio, once both edge lengths of the loops are in the inertial range.

Traditional Reynolds-averaged Navier–Stokes (RANS) closures, based on the Boussinesq eddy-viscosity hypothesis and calibrated on canonical flows, often yield inaccurate predictions of both mean flow and turbulence statistics. Here, we consider flow past a circular cylinder over a range of Reynolds numbers ($3900$–$100\,000$) and Mach numbers ($0$–$0.3$), encompassing incompressible and weakly compressible regimes, with the goal of improving predictions of mean velocity and Reynolds forces. To this end, we assemble a cross-validated dataset comprising hydrodynamic particle image velocimetry (PIV) in a towing tank, aerodynamic PIV in a wind tunnel and high-fidelity spectral element direct numerical simulation and large eddy simulation. Analysis of these data reveals a universal distribution of Reynolds stresses across the parameter space, which provides the foundation for a data-driven closure. We employ physics-informed neural networks (PINNs), trained with the unclosed RANS equations, to infer the velocity field and Reynolds-stress forcing from boundary information alone. The resulting closure, embedded in a forward PINN solver and the numerical solver OpenFOAM, significantly improves RANS predictions of both mean flow and turbulence statistics relative to conventional models.

Electrohydrodynamic (EHD) instabilities at polymer–porous interfaces play a pivotal role in determining interfacial morphology, wettability and pattern formation, with implications for energy storage, diagnostics and flexible electronics. This study presents a comprehensive general linear stability analysis to examine electric-field-induced instabilities at a confined interface between a viscoelastic polymer gel and a saturated porous medium. By coupling Maxwell stresses with a modified Darcy–Brinkman–Kelvin–Voigt framework, the model captures how porous medium-moderated EHD instabilities influence both the onset and dominant instability modes. Key parameters – including the electric Rayleigh number, Darcy number, dielectric contrast and geometric filling ratio – govern the spatio-temporal features of emerging patterns. The analysis reveals a sigmoidal dependence of characteristic length and time scales on permeability, i.e. Darcy number, establishing three regimes: impermeable, transitional and highly permeable, with a shift toward shorter wavelengths. The length and time scale transitions, triggered by the solid-saturated porous medium, are further moderated by the dielectric contrast – instabilities are suppressed when the contrast is low and amplified when it is high, enabling sub-micron patterning. Geometric confinement, i.e. increasing filling ratio, further intensifies pattern length scales, suggesting the feasibility of fabricating complex ultra-fine nanoscale encapsulated porous patterns. The elasticity of the viscoelastic layer imposes a threshold for instability onset and is critical for identifying wettability transitions at the interface. This framework offers predictive insight into tuning instability modes through permeability–viscoelasticity–electrostatics interplay, laying the foundation for wettability-controlled interfaces and self-organised interfacial patterns in next-generation EHD-driven systems.

This paper investigates asymmetric shock reflection in a supersonic overexpanded jet. A theoretical model that can predict the size and shape of flow field structures such as shocks, expansion waves and sliplines is established. For symmetric Mach reflection, the current model exhibits better agreement with numerical simulation results. For asymmetric Mach reflection, the current model also shows good agreement with numerical simulation results in predicting Mach stem height. The research indicates that the Mach stem height decreases with increasing nozzle pressure ratio and nozzle exit length difference, and decreases as the nozzle exit Mach number decreases. In addition, the critical geometric condition for complete misalignment of the upper and lower slipline interference segments (i.e. when the interference between one side of the expansion wave and the slipline ends, the other side of the expansion wave has not yet begun to interfere with the slipline) is given, which increases approximately linearly with nozzle pressure ratio and decreases as the nozzle exit Mach number increases. This study provides important theoretical support for the engineering application of asymmetric shock reflection in supersonic overexpanded jet.

Many environmental, energy and industrial processes involve the flow of viscoelastic polymer solutions in three-dimensional (3-D) porous media where fluid is confined to navigate through complex pore space geometries. As polymers are transported through the tortuous pore space, elastic stresses accumulate, leading to the onset of unsteady, time-dependent flow fluctuations above a threshold flow rate. How does pore space geometry influence the development and features of this elastic instability? Here, we address this question by directly imaging polymer solution flow in microfabricated 3-D ordered porous media with precisely controlled geometries consisting of simple-cubic (SC) or body-centred cuboid (BC) arrays of spherical grains. In both cases, we find that the flow instability is generated at stagnation points arising at the contacts between grains rather than at the polar upstream/downstream grain surfaces, as is the case for flow around a single grain. The characteristics of the flow instability are strongly dependent on the unit cell geometry: in SC packings, the instability manifests through the formation of time-dependent, fluctuating 3-D eddies; whereas in BC packings, it manifests as continual fluctuating ‘wobbles’ and crossing in the flow pathlines. Despite this difference, we find that characteristics of the transition from steady to unsteady flow with increasing flow rate have commonalities across geometries. Moreover, for both packing geometries, our data indicate that extensional flow-induced polymeric stresses generated by contact-associated stagnation points are the primary contributor to the macroscopic resistance to flow across the entire medium. Altogether, our work highlights the pivotal role of inter-grain contacts – which are typically idealised as discrete points and therefore overlooked, but are inherent in most natural and engineered media – in shaping elastic instabilities in porous media.

We investigate mixing dynamics in porous media at finite times, using pore-scale lattice-Boltzmann simulations combined with Lagrangian particle tracking. We compute fluid deformation in randomly packed beds based on the moving Protean frame approach introduced by Lester et al. (2018 J. Fluid Mech. 855, 770–803). From the extracted Lagrangian kinematics, we construct a mixing model based on lamellar aggregation that well predicts the Eulerian scalar fields obtained from simulations. Our results reveal an early-time mixing regime dominated by shear-driven fluid deformation, where solute mixing arises from the random overlap of diffusive concentration elements. In this regime, mixing proceeds slowly and follows a temporal decay of concentration variance, $\sigma _c^2 \propto \textit{Pe}^{-\alpha /(2\alpha +1)} t^{-1/2}$, where $ \textit{Pe}$ is the Péclet number and $\alpha$ the exponent characterising shear deformation. This dynamic arises when the Péclet number is small relative to the ratio between the exponential-mixing and shear-deformation time scales. This analysis also demonstrates that shear-induced mixing governs the homogenisation of early-stage reactions at the fluid–solid interface in finite-size random packed beds, typically operating at moderate Péclet numbers $ \textit{Pe} =O(10^2)$.

We experimentally study a scallop-like swimmer with reciprocally flapping wings in a nearly frictionless, cohesive granular medium consisting of hydrogel spheres. Significant locomotion is found when the swimmer’s flapping frequency matches the inverse relaxation time of the material. Remarkably, the swimmer moves in the opposite direction compared with its motion in a cohesion-free granular material of hard plastic spheres. At higher or lower frequencies, we observe no motion of the swimmer, apart from a short initial transient phase. X-ray radiograms reveal that the wing motions create low-density zones, which in turn give rise to a hysteresis in drag and propulsion forces. This time-dependent effect, combined with the swimmer’s inertia, accounts for locomotion at intermediate frequencies.

Coarse-grained continuous descriptions for lipid bilayers are typically based on minimising the Helfrich energy. Such models consider the fluid properties of these structures only implicitly and have been shown to nicely reproduce equilibrium properties. Model extensions that also address the dynamics of these structures are surface (Navier–)Stokes–Helfrich models. They explicitly account for membrane viscosity. However, these models also usually treat the lipid bilayer as a homogeneous continuum, neglecting the molecular degrees of freedom of the lipids. Here, we derive refined models that consider in addition a scalar order parameter representing the molecular alignment of the lipids along the surface normal. Starting from hydrodynamic surface liquid crystal models, we obtain a hydrodynamic surface Landau–Helfrich model for asymmetric lipid bilayers and a surface Beris–Edwards model for symmetric lipid bilayers. The fully ordered case for both models leads to the known surface (Navier–)Stokes–Helfrich models. Besides more detailed continuous models for lipid bilayers, we therefore also provide an alternative derivation of surface (Navier–)Stokes–Helfrich models. The impact on the dynamics is demonstrated by numerical simulations.

The inner–outer interaction model (IOIM), first proposed by Marusic et al. (Science, 2010, vol. 329, pp. 193–196), has proven to be an effective turbulence model for canonical and non-canonical wall-bounded flows, where a reference velocity signal from the logarithmic region acts as the input for predicting near-wall velocity fluctuations. Its most recent iteration by Baars et al. (Phys. Rev. Fluids, 2016, vol. 1, p. 054406) further proposes a user-independent scale separation point, refining model parameters. In this study, we compared the long-perceived universal IOIM’s parameters, including the linear transfer kernel, amplitude modulation coefficients and the universal signal for a range of Reynolds and Mach numbers, where mathematical relationships between the parameters are proposed. We observed that while the universal signals exhibit a high degree of similarity, particularly near the wall, the amplitude modulation coefficients and linear transfer kernels display Reynolds and Mach number effects, where varying the reference location also causes them to exhibit significant changes. We have found transformations to collapse amplitude modulation coefficients for incompressible flows and differing reference locations, improving modelling via the IOIM across flow parameters. Despite this, compressibility effects cannot be suitably accounted for currently and remain a future challenge for the IOIM framework.