The phenomenon of bulge evolution under the action of gravity on shallow water is prevalent both in natural occurrences and engineering industries. However, despite its ubiquity, its physical process remains largely unexplored. The evolution of bulge contains two fundamental physical processes: collapse and propagation. The collapse process can be further divided into two sub-processes: squeezing process and diffusion process. Based on the weakly nonlinear shallow water assumption with the classical perturbation method, the governing equations controlling the surface elevations in the diffusion process and the propagation process have been theoretically derived, where a bulge-induced surface pressure is modeled for the propagation process. Moreover, their scaling laws for the decay of wave height are also established, which have been validated by direct numerical simulation results. The derived scaling laws for wave height attenuation of bulge evolution provide profound insights, which hold the potential to applications in the engineering industry.

This paper presents a method to stabilise oscillations occurring in a mixed convective flow in a nearly hemispherical cavity, using actuation based on the receptivity map of the unstable mode. This configuration models the continuous casting of metallic alloys, where hot liquid metal is poured at the top of a hot sump with cold walls pulled in a solid phase at the bottom. The model focuses on the underlying fundamental thermohydrodynamic processes without dealing with the complexity inherent to the real configuration. This flow exhibits three branches of instability. The solution of the adjoint eigenvalue problem for the convective flow equations reveals that the regions of highest receptivity for unstable modes of each branch concentrate near the inflow upper surface. Simulations of the linearised governing equations show that a thermomechanical actuation modelled on the adjoint eigenmode asymptotically suppresses the unstable mode. If the actuation’s amplitude is kept constant in time, which is easier to implement in an industrial environment, the suppression is still effective but only over a finite time, after which it becomes destabilising. Based on this phenomenology, we apply the same actuation during the stabilising phase only in the nonlinear evolution of the unstable mode. It turns out stabilisation persists, even when the unstable mode is left to evolve freely after the actuation period. These results not only demonstrate the effectiveness of receptivity-informed actuation in stabilising convective oscillations but also suggest a simple strategy for their long-term control.

The dynamics of wall-mounted flexible structures, such as aquatic vegetation, is essential for analysing collective behaviours, flow distributions and vortex formation across different scales. To accurately model these structures under various flow conditions, we develop a novel numerical method that couples the immersed boundary method (IBM) with the vector form intrinsic finite element (VFIFE) method, referred to as the IBM–VFIFE method. We simulate both flexible and rigid stems, each with a constant aspect ratio of 10, mounted on an impermeable bottom in uniform flow with the Reynolds number ranging from 200 to 1000. In the rigid case, we identify three distinct flow regimes based on the vortex dynamics and lift spectral characteristics. Due to the influences of downwash flow at the free end and upwash flow near the junction, vortex shedding varies significantly along the vertical direction. For the flexible case, we examine a wide range of stem stiffness values to explore potential dynamic responses. The results reveal that stiffness plays a key role in stem behaviour, leading to three distinct classifications based on amplitude magnitude and displacement spectra respectively. Notably, the vortex dynamics of a flexible stem differs significantly from that of a rigid stem due to shape deformation and stem oscillation. A flexible stem with relatively high stiffness experiences greater hydrodynamic loads compared with its rigid counterpart. This study highlights the unique stem behaviours and vortex dynamics associated with flexible stems. We find that stem oscillation, combined with a near-wake base vortex, contributes to an upwash flow near the stem bottom, which significantly weakens (or, in some cases, eliminates) the downwash flow. Additionally, low-frequency oscillations in the streamwise and vertical directions are observed, while the transverse oscillation exhibits a dominant frequency one order of magnitude higher. Overall, this study provides valuable insights into the response and vortex dynamics of a single stem in uniform flow.

The classical water-wave theory often neglects water compressibility effects, assuming acoustic and gravity waves propagate independently due to their disparate spatial and temporal scales. However, nonlinear interactions can couple these wave modes, enabling energy transfer between them. This study adopts a dynamical systems approach to investigate acoustic–gravity wave triads in compressible water flow, employing phase-plane analysis to reveal complex bifurcation structures and identify steady-state resonant configurations. Through this framework, we identify specific parameter conditions that enable complete energy exchange between surface and acoustic modes, with the triad phase (also known as the dynamical phase) playing a crucial role in modulating energy transfer. Further, incorporating spatial dependencies into the triad system reveals additional dynamical effects that depend on the wave velocity and resonance conditions: we observe that travelling-wave solutions emerge, and their stability is governed by the Hamiltonian structure of the system. The phase-plane analysis shows that, for certain velocity regimes, the resonance dynamics remains similar to the spatially independent case, while in other regimes, bifurcations modify the structure of resonant interactions, influencing the efficiency of energy exchange. Additionally, modulated periodic solutions appear, exhibiting changes in wave amplitudes over time and space, with implications for wave-packet stability and energy localisation. These findings enhance the theoretical understanding of acoustic–gravity wave interactions, offering potential applications in geophysical phenomena such as oceanic microseisms.

The motion of several plates in an inviscid and incompressible fluid is studied numerically using a vortex sheet model. Two to four plates are initially placed in line, separated by a specified distance, and actuated in the vertical direction with a prescribed oscillatory heaving motion. The vertical motion induces the plates’ horizontal acceleration due to their self-induced thrust and fluid drag forces. In certain parameter regimes, the plates adopt equilibrium ‘schooling modes’, wherein they translate at a steady horizontal velocity while maintaining a constant separation distance between them. The separation distances are found to be quantised on the flapping wavelength. As either the number of plates increases or the flapping amplitude decreases, the schooling modes destabilise via oscillations that propagate downstream from the leader and cause collisions between the plates, an instability that is similar to that observed in recent experiments on flapping wings in a water tank (Newbolt et al., 2024, Nat. Commun., vol. 15, 3462). A simple control mechanism is implemented, wherein each plate accelerates or decelerates according to its velocity relative to the plate directly ahead by modulating its own flapping amplitude. This mechanism is shown to successfully stabilise the schooling modes, with remarkable impact on the regularity of the vortex pattern in the wake. Several phenomena observed in the simulations are obtained by a reduced model based on linear thin-aerofoil theory.

We present a numerical scheme that solves for the self-similar viscous fingers that emerge from the Saffman–Taylor instability in a divergent wedge. This is based on the formulation by Ben Amar (1991, Phys. Rev. A, vol. 44, pp. 3673–3685). It is demonstrated that there exists a countably infinite set of selected solutions, each with an associated relative finger angle, and furthermore, solutions can be characterised by the number of ripples located at the tip of their finger profiles. Our numerical scheme allows us to observe these ripples and measure them, demonstrating that the amplitudes are exponentially small in terms of the surface tension; the selection mechanism is driven by these exponentially small contributions. A recently published paper derived the selection mechanism for this problem using exponential asymptotic analytical techniques, and obtained bifurcation diagrams that we compare with our numerical results.

Waves propagating over oscillating periodic structures can be reflected and attenuated either by Bragg scattering or by local resonance. In this work, we focus on the interplay between surface gravity waves and submerged resonators, investigating the effect of the local resonance on wave propagation. The study is performed using a state of the art numerical simulation of the Navier–Stokes equation in two-dimensional form with free boundary and moving bodies. A volume of fluid interface technique is employed for tracking the free surface, and an immersed boundary method for the fluid–structure interaction. A wave maker is placed at one end of the flume and an absorbing beach at the other. The evolution in space of a monochromatic wave interacting with up to four resonators coupled only fluid mechanically is presented. We evaluate the efficiency of the system in terms of wave amplitude attenuation and energy transfers between the fluid and the solid phase. The results indicate that, near resonance conditions, both wave reflection and energy dissipation increase significantly. Conversely, far from resonance, waves can propagate through the system with minimal dissipation, even in the presence of numerous resonators. Moreover, when the time scale associated with the resonator’s restoring force is longer than the wave period, the resonators tend to follow the wave motion, oscillating with an amplitude comparable to that of the wave. In contrast, when the two time scales are similar, the resonator motion becomes amplified, resulting in stronger velocity gradients and enhanced viscous dissipation.

This paper presents a theoretical model for the electro-osmotic flow (EOF) of semi-dilute polyelectrolyte (PE) solutions in nanochannels. We use mean-field theories to describe the properties of electric double layer and viscosity of PE solutions that are prerequisites for constructing the EOF model. The EOF model is validated via a good match to the existing experimental results. Based on the validated EOF model, we conduct a comprehensive analysis of EOF of semi-dilute PE solutions in nanochannels. First, we observe considerable EOF of PE solutions in the uncharged nanochannels, which is in stark contrast to EOF of simple electrolyte solutions. The analyses show that the EOF of PE solutions in uncharged nanochannels is triggered by the external electric field acting on the near-wall non-electroneutral regions resulting from the confinement-induced inhomogeneous distribution of PE monomers. Although the solutions are electroneutral as a whole, the presence of local non-electroneutral regions and the mismatch between non-electroneutral regions and high-viscosity regions lead to the net EOF in uncharged nanochannels. Furthermore, we reveal that the EOF mobility $\mu _{{eof}}$ in uncharged nanochannels exhibits a scaling law $\mu _{{eof}} \propto a^{-0.44}$ (wherein $a$ denotes monomer Kuhn length) and is inversely proportional to the PE chain length, while it decreases nonlinearly with the charge fraction of the PE chains. Moreover, the EOF mobility reaches its maximum at specific bulk monomer concentration, and increases with the nanochannel height before converging to that under no confinement. Second, we analyse the EOF of PE solutions in nanochannels with various wall effects, such as surface charge density, slip length and adsorption length. When the surface charge is absent, the adsorption length significantly influences the direction and magnitude of the EOF, whereas the slip length has no effect. When the wall becomes increasingly charged, the influence of adsorption length on EOF gradually diminishes, while the importance of the slip length progressively intensifies and the EOF is highly influenced by the co-action of various wall effects in a complicated manner. When the surface wall is oppositely charged to polymer monomers, the EOF mobility varies nonlinearly with the surface charge density, while a zero net flow of EOF followed by a direction reversal is discovered when the wall is likely charged to polymer monomers.

The settling of highly elastic non-Brownian closed fibres (called loops) under gravity in a viscous fluid is investigated numerically. The loops are represented using a bead–spring model with harmonic bending potential and finitely extensible nonlinear elastic stretching potential. Numerical solutions to the Stokes equations are obtained with the use of HYDROMULTIPOLE numerical codes, which are based on the multipole method corrected for lubrication to calculate hydrodynamic interactions between spherical particles with high precision. Depending on the elasto-gravitation number $B$, a ratio of gravitation to bending forces, the loop approaches different attracting dynamical modes, as described by Gruziel-Słomka et al. (2019 Soft Matt. 15, 7262–7274) with the use of the Rotne–Prager mobility of the elastic loop made of beads. Here, using a more precise method, we find and characterise a new mode, analyse typical time scales, velocities and orientations of all the modes, compare them and investigate their coexistence. We analyse numerically the transitions (bifurcations) to a different mode at certain critical values of the elasto-gravitation number.

The settling dynamics of finite size, slightly heavier-than-fluid Kolmogorov-scale particles in homogeneous, isotropic turbulence at moderate volume loadings is investigated. A thoroughly validated two-way-coupled, point-particle model based on the complete Maxey–Riley–Gatignol equation of motion is used with closure models for all forces, including the history force, together with corrections for the self-disturbance field created by the particle using a novel zonal-advection-diffusion-reaction method. Settling dynamics is investigated by varying turbulence intensities relative to the particle settling speed in quiescent flow for multiple Stokes numbers. The length scales associated with the turbulence structures that strongly interact with and influence the settling dynamics are investigated using multiscale statistical analysis of the fluid velocity and second invariant of the velocity gradient tensors sampled by the particles. The time scales are investigated using trajectory curvature angle statistics of inertial and fluid particles. Low-to-moderate Stokes number particles tend to sample strain-rate dominated regions of the flow, tend to follow the curvature of the flow paths and show enhanced settling at higher turbulence intensities due to fast tracking and preferential sampling. Higher Stokes number particles, on the other hand, have a tendency to travel in straight lines relative to the flow and result in reduced settling speeds due to loitering. For the low mass loadings considered in this work, there is minimal global effect on the turbulent flow characteristics; however, it is found that the Kolmogorov-scale particles interact with and locally modify flow structures approximately twice their size, whereas they sample flow velocities from scales up to ten times the particle size, influencing preferential sampling and settling characteristics.

Capsules are widely used in bioengineering, chemical engineering and industry. The development of drug delivery systems using deformable capsules is progressing, yet the regulation of drug release within a capsule remains a challenge. Meanwhile, a microswimmer enclosed in a capsule can generate a large lubrication force on the capsule membrane, which could result in deformation and mechanical damage to the membrane. In this study, we numerically investigate how a capsule can be damaged by an enclosed microswimmer. The capsule membrane is modelled as a two-dimensional neo-Hookean material, with its deformability parametrised by capillary number. An isotropic brittle damage model is applied to express membrane rupture, with the Lighthill–Blake squirmer serving as the microswimmer model. In a sufficiently small capillary number regime, pusher-type squirmers exhibit stable swimming along the capsule membrane, while neutral-type and puller-type squirmers exhibit swimming towards the membrane and remain stationary. As capillary number increases, the damage to the membrane increases and rupture occurs in all swimming modes. For pusher-type squirmers, the critical capillary number leading to rupture is dependent on the initial incidence angle, whereas neutral-type and puller-type squirmers are independent of the initial value. Furthermore, we present methods for controlling membrane damage by magnetically orienting the microswimmer. The findings reveal that a static magnetic field can orient the microswimmer, leading to membrane damage and rupture even for a capsule that cannot be damaged by free swimming, while controlling the swimming path with a rotating magnetic field enables soft membranes to maintain deformation without rupture.

The flow-induced oscillations of a clamped flexible ring in a uniform flow were explored using the penalty immersed boundary method. Both inverted and conventional ring configurations were examined, with systematic analysis focused on the effects of bending rigidity and eccentricity. Four distinct oscillation modes were identified across parameter variations: flapping (F), deflected oscillation (DO), transverse oscillation (TO) and equilibrium (E) modes. Each mode exhibited a 2S wake pattern. The inverted ring sustained the DO mode under low bending rigidity with a deflected shape, transitioning to the TO mode at higher bending rigidity. In the TO mode, a lock-in phenomenon emerged, enabling the inverted ring to achieve a high power coefficient due to a simultaneous rise in both oscillation amplitude and frequency. By contrast, the conventional ring exhibited the F mode at low bending rigidity and transitioned to the E mode as rigidity increased, although its power coefficient remained lower because of reduced critical bending rigidity. For the inverted ring, low eccentricity enhanced oscillation intensity but limited the operational range of the TO mode. In contrast, for the conventional ring, reducing eccentricity led to an increase in oscillation amplitude. Among the investigated configurations, the inverted-clamped ring achieved the highest energy-harvesting efficiency, surpassing those of the conventional clamped ring and a buckled filament.