The acoustic field radiated by a system of contra-rotating propellers in wetted conditions (with no cavitation) is reconstructed by exploiting the Ffowcs Williams–Hawkings acoustic analogy and a database of instantaneous realizations of the flow. They were generated by high-fidelity computations using a large eddy simulation approach on a cylindrical grid of 4.6 billion points. Results are also compared against the cases of the front and rear propellers working alone. The analysis shows that the importance of the quadrupole component of sound, originating from wake turbulence and instability of the tip vortices, is reinforced, relative to the linear component radiated from the surface of the propeller blades. The sound from the contra-rotating propellers decays at a slower rate for increasing radial distances, compared with the cases of the isolated front and rear propellers, again due to the quadrupole component. The quadrupole sound is often neglected in the analysis of the acoustic signature of marine propellers, by considering the only linear component. In contrast, the results of this study point out that the quadrupole component becomes the leading one in the case of contra-rotating propulsion systems, due to the increased complexity of their wake. This is especially the result of the mutual inductance phenomena between the tip vortices shed by the front and rear propellers of the contra-rotating system.

The hydrodynamic behaviours of finite-size microorganisms in turbulent channel flows are investigated using a direct-forcing fictitious domain method. The classical ‘squirmer’ model, characterized by self-propulsion through tangential surface waves at its boundaries, is employed to mimic the swimming microorganisms. We adopt various simulation parameters, including a friction Reynolds number Reτ = 180, two squirmer volume fractions 𝜑0 = 12.7 % and 2.54 % and a blocking ratio (squirmer radius/half-channel width) κ = 0.125. Results show that pushers (propelled from the rear) induce a more pronounced decrease in the velocity profile than neutral squirmers and pullers (propelled from the front). This hindrance and the induced particle inner stress τpI positively correlate with the quantity of squirmers accumulated in the near-wall region. Notably, the increase in τpI primarily occurs at the expense of diminishing the fluid Reynolds stress τfR. Compared with passive spheres, a low volume fraction (𝜑0 = 2.54 %) of pullers results in a slightly enhanced velocity profile across the channel. In the near-wall region, the swimming direction of the squirmers shows no significant tendency with respect to the flow direction. In the bulk-flow region, pushers and neutral squirmers tend to align their axes more along the flow direction, whereas pullers exhibit a slight preference for alignment with the normal direction.

Within the frameworks of the amplitude method and the linear stability theory, a statistical model of the initial stage of laminar–turbulent transition caused by atmospheric particulates (aerosols) penetrating into the boundary layer is developed. The model accounts for the process of boundary layer receptivity to particulates, asymptotic behaviour of unstable wave packets propagating downstream from particle–wall collisions and the amplitude criterion for the transition onset. The resulting analytical relationships can be used for quick predictions of the transition onset on bodies of relatively simple shape, where the undisturbed boundary layer is quasi-two-dimensional. The model allows us to explore the transition onset at realistic distributions of the particle concentration selected based on an analysis of known empirical data. As an example, a 14° half-angle sharp wedge flying in atmosphere at 20 km altitude and Mach number 4 is considered. It is shown that the transition onset corresponds to an N-factor of 15.3 for a flight under normal atmospheric conditions and 12.2 for a flight in a cloud after volcanic eruption. In accordance with physical restrictions, these values are below the upper limit $N\approx 16.8$ predicted for transition due to thermal fluctuations (perfectly ‘clean’ case). Nevertheless, they are significantly greater than $N=10$ which is commonly recommended for estimates of the transition onset in flight.

We study the dynamic deflation of a hydraulic fracture subject to fluid withdrawal through a narrow conduit located at the centre of the fracture. Recent work revealed a self-similar dipole-flow regime, when the influence of material toughness is negligibly small. The focus of the current work is on the influence of material toughness, which leads to an additional self-similar regime of fracture deflation with fixed frontal locations in the toughness-dominated regime. The two limiting regimes can be distinguished by a dimensionless material toughness $\Pi _k$, defined based on a comparison with the influence of the viscous thin film flow within the fracture: $\Pi _k \to 0$ indicates the dipole-flow regime, while $\Pi _k \to \infty$ indicates the fixed-length regime. For intermediate $\Pi _k$, the fracture’s front continues to propagate during an initial period of deflation before it remains pinned at a fixed location thereafter. A regime diagram is then derived, with key scaling behaviours for the frontal dynamics, pressure and volume evolution summarised in a table for the self-similar stage. A comparison is also attempted between theoretical predictions and available experimental observations of viscous backflows from transparent solid gelatins.

We use direct numerical simulations to investigate the energy pathways between the velocity and the magnetic fields in a rotating plane layer dynamo driven by Rayleigh–Bénard convection. The kinetic and magnetic energies are divided into mean and turbulent components to study the production, transport and dissipation in large- and small-scale dynamos. This energy balance-based characterisation reveals distinct mechanisms for large- and small-scale magnetic field generation in dynamos, depending on the nature of the velocity field and the conditions imposed at the boundaries. The efficiency of a dynamo in converting the kinetic energy to magnetic energy, apart from the energy redistribution inside the domain, is found to depend on the kinematic and magnetic boundary conditions. In a small-scale dynamo with a turbulent velocity field, the turbulent kinetic energy converts to turbulent magnetic energy via small-scale magnetic field stretching. This term also represents the amplification of the turbulent magnetic energy due to work done by stretching the small-scale magnetic field lines owing to fluctuating velocity gradients. The stretching of the large-scale magnetic field plays a significant role in this energy conversion in a large-scale turbulent dynamo, leading to a broad range of energetic scales in the magnetic field compared with a small-scale dynamo. This large-scale magnetic field stretching becomes the dominant mechanism of magnetic energy generation in a weakly nonlinear dynamo. We also find that, in the weakly nonlinear dynamo, an upscale energy transfer from the small-scale magnetic field to the large-scale magnetic field occurs owing to the presence of a gradient of the mean magnetic field.

The dependence of the Richtmyer–Meshkov instability (RMI) on post-shock Atwood number ($A_1$) is experimentally investigated for a heavy–light single-mode interface. We create initial interfaces with density ratios of heavy to light gases ranging from 1.73 to 34.07, and achieve the highest $|A_1|$ value reported to date for gaseous-interface experiments (0.95). For the first time, spike acceleration is observed in experiments with a heavy–light configuration. The models for the start-up, linear and weakly nonlinear evolution stages are evaluated over a wide range of $A_1$ conditions. Specifically, the models proposed by Li et al. (Phys. Fluids, vol. 36, 2024, 056104) and Wouchuk & Nishihara (Phys. Plasmas, vol. 4, 1997, 1028–1038) effectively describe the start-up and linear stages, respectively, across all cases. None of the considered nonlinear models is valid under all $A_1$ conditions. Based on the dependence of spike and bubble evolutions on $A_1$ provided by the present work and previous study (Chen et al., J. Fluid Mech., vol. 975, 2023, A29), the SEA model (Sadot et al., Phys. Rev. Lett., vol. 80, 1998, pp. 1654–1657), whose expression has clear physical meanings, is modified by revising the coefficient that governs its prediction for early-time evolution. The modified model applies to prediction of the weakly nonlinear evolution of RMI with $A_1$ ranging from −0.95 to −0.35 and from 0.30 to 0.86. Based on this model, an approximation of the critical $A_1$ for the occurrence of spike acceleration is obtained.

Recent experiments and simulations have sparked growing interest in the study of Rayleigh–Bénard convection in very slender cells. One pivotal inquiry arising from this interest is the elucidation of the flow structure within these very slender cells. Here we employ tomographic particle image velocimetry, for the first time, to capture experimentally the full-field three-dimensional and three-component velocity field in a very slender cylindrical cell with aspect ratio $\Gamma =1/10$. The experiments cover a Rayleigh number range $5.0 \times 10^8 \leqslant Ra \leqslant 5.0 \times 10^9$ and Prandtl number 5.7. Our experiments reveal that the flow structure in the $\Gamma =1/10$ cell is neither in the multiple-roll form nor in the simple helical form; instead, the ascending and descending flows can intersect and cross each other, resulting in the crossing events. These crossing events separate the flow into segments; within each segment, the ascending and descending flows ascend or descend side by side vertically or in the twisting manner, and the twisting is not unidirectional, while the segments near the boundary can also be in the form of a donut like structure. By applying the mode decomposition analyses to the measured three-dimensional velocity fields, we identified the crossing events as well as the twisting events for each instantaneous flow field. Statistical analysis of the modes reveals that as $Ra$ increases, the average length of the segments becomes smaller, and the average number of segments increases from 2.5 to 3.9 in the $Ra$ range of our experiments.

The interaction between a turbulent flow and a porous boundary is analysed with focus on the sensitivity of the roughness function, $\Delta U^+$, to the upscaled coefficients characterizing the wall. The study is aimed at (i) demonstrating that imposing effective velocity boundary conditions at a virtual plane boundary, next to the physical one, can efficiently simplify the direct numerical simulations (DNS); and (ii) pursuing correlations to estimate $\Delta U^+$ a priori, once the upscaled coefficients are calculated. The homogenization approach employed incorporates near-interface advection via an Oseen-like linearization, and the macroscopic coefficients thus depend on both the microstructural details of the wall and a slip-velocity-based Reynolds number, $Re_{slip}$. A set of homogenization-simplified DNS is run to study the channel flow over transversely isotropic porous beds, testing values of the grains’ pitch within $0\lt \ell ^+\lt 40$. Reduction of the skin-friction drag is attainable exclusively over streamwise-aligned inclusions for $\ell ^+$ values up to $20{-}30$. The drag increase over spanwise-aligned inclusions (or streamwise-aligned ones at large $\ell ^+$) is accompanied by enhanced turbulence levels, including intensified sweep and ejection events. The root-mean-square of the transpiration velocity fluctuations at the virtual plane, $\tilde V_{rms}$, is the key control parameter of $\Delta U^+$; our analysis shows that, provided $\tilde V_{rms} \lesssim 0.25$, then $\tilde V_{rms}$ is strongly correlated to a single macroscopic quantity, $\Psi$, which comprises the Navier-slip and interface/intrinsic permeability coefficients. Fitting relationships for $\Delta U^+$ are proposed, and their applicability is confirmed against reference results for the turbulent flow over impermeable walls roughened with three-dimensional protrusions or different geometries of riblets.

A new temporal vortex tracking algorithm allows the first long-term temporal observation of the lives of the intense vorticity structures (IVS). The algorithm is applied to direct numerical simulations of statistically stationary isotropic turbulence, with Taylor-based Reynolds numbers in the range $54 \leqslant Re_{\lambda } \leqslant 239$. In the highest Reynolds number case, the continuous time tracking of millions of ‘worms’ is achieved for more than seven integral time scales and close to 200 Kolmogorov time scales. Within an integral scale volume, more than 66 structures exist, and approximately 20 new structures are created per Kolmogorov time. More than $80\, \%$ of the structures live a solitary ‘life’ without any visible interaction with the other structures, while approximately $15\, \%$ break into new structures. Less than $2\, \%$ of the structures merge with others to form new vortices. A ‘population model’ is developed to estimate the numbers of existing vortices for very long simulated times, and it is observed that the birth rate density of these structures slowly increases with the Reynolds number. The survival functions of the IVS lives exhibit an exponential distribution, with some structures living for more than $35$ Kolmogorov time scales (more than four integral time scales). The mean lifetime of the IVS scales with the mean turnover time scale of the worms, defined by their radii and tangential velocity, attaining $\approx 6.5$ turnover time scales at high Reynolds numbers.

Liquid metal flows are important for many industrial processes, including liquid metal batteries (LMBs), whose efficiency and lifetime can be affected by fluid mixing. We experimentally investigate flows driven by electrical currents in an LMB model. In our cylindrical apparatus, we observe a poloidal flow that descends near the centreline for strong currents, and a poloidal flow that rises near the centreline for weak currents. The first case is consistent with electrovortex flow, which is an interaction between current and its own magnetic field, whereas the second case is consistent with an interaction between current and the external field, which drives Ekman pumping. Notably, we also observe an intermediate case where the two behaviours appear to compete. Comparing results with Frick et al. (2022 J. Fluid Mech. 949, A20), we test prior estimates of the scaling of flow speed with current to predict the observed reversal. Based on these data, we propose two different ways to apply the Davidson et al. (1999 J. Fluid Mech. 245, 669–699) poloidal suppression theory that explain both experimental results simultaneously: either taking the wire radius into account to scale the Lorentz force, or taking viscous dissipation into account to scale the swirl velocity, following Herreman et al. (2021 J. Fluid Mech. 915, A17).

Self-similar fractal tree models are numerically investigated to elucidate the drag coefficient, non-equilibrium dissipation behaviour and various turbulence statistics of fractal trees. For the simulation, a technique based on the lattice Boltzmann method with a cumulant collision term is used. Self-similar fractal tree models for aerodynamic computations are constructed using parametric L-system rules. Computations across a range of tree-height-based Reynolds numbers $Re_H$, from 2500 to 120 000, are performed using multiple tree models. As per the findings, the drag coefficients ($C_D$) of these models match closely with those of the previous literature at high Reynolds numbers ($Re_H \geqslant 60\,000$). A normalization process that collapses the turbulence intensity across various tree models is formulated. For a single tree model, a consistent centreline turbulence intensity trend is maintained in the wake region beyond a Reynolds number of 60 000. The global and local isotropy analysis of the turbulence generated by fractal trees indicates that, at high Reynolds numbers ($Re_H \geqslant 60\,000$), the distant wake can be considered nearly locally isotropic. The numerical results confirm the non-equilibrium dissipation behaviour demonstrated in previous studies involving space-filling fractal square grids. The non-dimensional dissipation rate $C_\epsilon$ does not remain constant; instead, it becomes approximately inversely proportional to the local Taylor-microscale-based Reynolds number, $C_\epsilon \propto 1/Re_\lambda$. We find significant one-point inhomogeneity, production and transverse transport of turbulent kinetic energy within the non-equilibrium dissipation near wake region.

We present a linear analysis of a minimal model of moist convection under a variety of atmospheric conditions. The stationary solutions that we analyse include both fully saturated and partially unsaturated atmospheres in both unconditionally and conditionally unstable cases. We find that all of the solutions we consider are linearly unstable via exchange of stability when sufficiently driven. The critical Rayleigh numbers vary by over an order of magnitude between unconditionally unstable and conditionally unstable atmospheres. The unsaturated atmospheres are notable for the presence of linear gravity wave-like oscillations even in unstable conditions. We study their eigenfunction structure and find that the buoyancy and moisture perturbations are anticorrelated in $z$, such that regions of negative buoyancy have positive moisture content. We suggest that these features in unsaturated atmospheres may explain the phenomenon of gravity wave shedding by moist convective plumes.

The bi-stable dynamics of a one-degree-of-freedom disk pendulum swept by a flow and allowed to rotate in the cross-flow direction is investigated experimentally. For increasing flow velocity, a subcritical bifurcation is observed from a Pendulum state, characterised by an increasing time-averaged pendulum angle with large amplitude fluctuations, to a rotating state with a non-zero mean rotation velocity at a critical free stream velocity $U_{P2W}$. The rotating state, referred to as Windmill state, presents a strong hysteresis: once initiated, it is sustained down to velocities $U_{W2P}\lt U_{P2W}$ before bifurcating towards the Pendulum state. A thorough experimental characterisation of the dynamical features of each state is reported, with a particular focus on the influence of the static yaw angle of the disk $\beta _0$ and the free stream velocity. In the Pendulum state, the system behaves differently depending on whether $\beta _0$ lies below or above the stall angle of the disk, with more regular dynamics below. We demonstrate that the bifurcation between the Pendulum state and the Windmill state is triggered by aerodynamic fluctuations, while the bifurcation between the Windmill state and the Pendulum state is deterministic. A stochastic model faithfully reproduces the dynamical features of both states, as well as the characteristics of the bifurcations.

Microorganisms, such as spermatozoa, exhibit rich behaviours when in close proximity to each other. However, their locomotion is not fully understood when coupled mechanically and hydrodynamically. In this study, we develop hydrodynamic models to investigate the locomotion of paired spermatozoa, predicting the fine structure of their swimming. Experimentally, sperm pairs are observed to transition between different modes of flagellar synchronisation: in-phase, anti-phase and lagged synchronisation. Using our models, we assess their swimming performances in these synchronisation modes in terms of average swimming speed, average power consumption, and swimming efficiency. The swimming performances of paired spermatozoa are shown to depend on their flagellar phase lag, flagellar waveforms, and the mechanical coupling between their heads.

In this study, we conducted interface-capturing high-resolution simulations of a bubbly upflow in a vertical channel to investigate the bubble distribution and its interaction with surrounding turbulence, focusing on the effects of the density ratio. A bulk Reynolds number $Re_b=2300$ was used for all simulations. The influence of density ratio on vortex structures and turbulence statistics differed between the near-wall and core regions of the channel. Adding 5.43 $\%$ gas caused an increase in wall friction. By applying a generalised FIK identity to analyse wall friction, it was determined that the drag rise in the bubbly channel was mostly due to the near-wall region. Visualisation of the bubble and vortex structures showed that small bubbles near the wall induced larger magnitude of Reynolds shear stress and increased wall friction. Bubble behaviour near the wall region was similar for density ratios above 30, leading to wall friction saturation. In the core region, large deformable bubbles created wake vortices due to slip velocity between liquid and gas phases. Wake vortices help large bubbles absorb smaller bubbles and maintain their sizes. As the density ratio increased, the slip velocity increased owing to greater difference in the gravitational acceleration between liquid and gas phases, resulting in corresponding increase in wake intensity and velocity fluctuations. However, quadrant analysis showed that Q1 and Q3 events increased together with Q2 and Q4 events in the core region, cancelling out any net effect of wake vortices on Reynolds shear stress or wall friction.

Highly resolved simulations reveal the fundamental influence of a carrier fluid’s flow dynamics on triboelectric powder charging. We found that particles transported through a square-shaped duct charge faster than in a channel flow caused by secondary flows that led to more severe particle–wall collisions. Specifically, particles with a Stokes number of 4.69 achieve 85 % of their equilibrium charge approximately 1.5 times faster in duct flow than in channel flow. Also, charge distribution is more uniform in a duct cross-section compared with a channel cross-section. In channel flow, particles are trapped near the walls and collide frequently due to limited movement in the wall-normal direction, causing localized charge buildup. In contrast, duct flow promotes better mixing through secondary flows, reducing repeating collisions and providing uniform charge distribution across the cross-section. Upon charging, electrostatic forces significantly reshape particle behaviour and distribution. Once the powder achieves half of its equilibrium charge, particles increasingly accumulate at the wall, leading to a reduced concentration in the central region. These changes in particle distribution have a noticeable impact on the surrounding fluid phase and alter the overall flow dynamics. These findings open the possibility for a new measure to control powder charging by imposing a specific pattern.

Real-time wave forecasting (RTWF) consists in predicting ocean wave motion or forces, from seconds to minutes in advance, using real-time measurements. For the successful development of RTWF, understanding wave predictability is essential. Usually, a deterministic ‘predictable zone’ (DPZ) is geometrically constructed from the wave group velocities and directions present in the spectrum. DPZs have little experimental evidence, and suffer ambiguities regarding the choice of cutoff frequencies and directions – since actual ocean waves are not band-limited. The present study addresses those shortcomings, by defining probabilistic predictable zones (PPZs) with respect to chosen uncertainty thresholds, using a rigorous statistical framework restricted to near-Gaussian sea states (precisely those where RTWF would be employed). PPZs are examined in idealised spectra and in a stereo dataset of a real wave field. It is shown that the PPZ geometry is quantitatively related to the sea state characteristics, through three physical parameters: two limiting group velocities (similar to the deterministic theory), and a directional spreading effect, which also limits the PPZ extent. While the lower group velocity depends on the chosen uncertainty threshold, the upper group velocity is better approximated by that of the spectrum peak frequency, which is a novel finding. The empirical data support the validity of the present PPZ theory. In contrast, both theoretical and empirical results contradict the fan-shaped predictable zones, constructed in the three-dimensional deterministic theory, thus highlighting the importance of considering stochasticity to understand the predictability of actual ocean waves.

We report laboratory experiments of long-crested irregular water surface waves propagating over a shoal, with attention to the region over the down-slope behind the shoal. We measure the surface elevation field, the horizontal velocity field in the water, and the resulting forces on a horizontal submerged cylinder placed over the down-slope of the shoal. In addition, we calculate the horizontal acceleration field. From this, we find that the presence of the shoal can modify the wave field such that the resulting forces on the submerged cylinder can be enhanced with thicker extreme tails and increased values of skewness and kurtosis depending on the location of the cylinder. The spatial dependence of the statistics of forces is different from the spatial dependence of the statistics of horizontal velocity, horizontal acceleration and surface elevation.

A pressure-gradient-induced laminar separation bubble (LSB) was examined using wind-tunnel experiments, direct numerical simulations (DNS) and linear local/global stability analysis. The LSB was experimentally generated on a flat plate using the favourable-to-adverse pressure gradient imposed by an inverted modified NACA $64_3-618$ airfoil. Direct numerical simulation was performed using boundary conditions extracted from a steady precursor simulation of the entire flow field. Despite good agreement in the upstream boundary layer with the experiment, DNS exhibited an approximately 25 % longer mean separation bubble, attributed to an earlier onset of transition due to the free-stream turbulence (FST) in the experiment. Introducing a very low level of isotropic FST in the DNS, similar to that measured in the experiment, caused earlier transition, decreased the mean bubble length and led to a remarkably good agreement between the DNS and experiments. Differences were observed for the dominant frequencies in the experiment and DNS, but both were within the band of most amplified frequencies predicted by LST. Proper orthogonal decomposition confirmed that dominant coherent structures from DNS and experiments are related to the inviscid Kelvin–Helmholtz instability and have similar characteristics despite slight differences in frequency. Local and global stability and dynamic mode decomposition analysis corroborated the convective nature of the dominant two-dimensional (2-D) instability and showed that the LSB is globally unstable to a range of 3-D wavenumbers, in agreement with 3-D structures observed in experiments. Results confirmed the strong impact of very low FST levels on the LSB and indicate a close agreement of the time-averaged and instability characteristics between the experiments and DNS.

The incompressible Navier–Stokes equations in spherical coordinates are solved using a pseudo-spectral method to simulate the problem of spherical Couette flow. The flow is investigated for a narrow-gap ratio with only the inner sphere rotating. We find that the flow is sensitive to the initial conditions and have used various initial conditions to obtain different branches of the bifurcation curve of the flow. We have identified three different branches dominated respectively by axisymmetric flow, travelling wave instability and equatorial instability. The axisymmetric branch shows unsteadiness at large Reynolds numbers. The travelling wave instability branch shows spiral instability and is prominent the near poles. The travelling wave instability branch further exhibits a reversal in the propagation direction of the spiral instability as the Reynolds number is increased. This branch also exhibits a multi-mode equatorial instability at larger Reynolds numbers. The equatorial instability branch exhibits twin jet streams on either side of the equator, which become unstable at larger Reynolds numbers. The flow topology on the three branches is also investigated in their phase space and found to exhibit chaotic behaviour at large Reynolds numbers on the travelling wave instability branch.