Local shearing motions in turbulence form small-scale shear layers, which are unstable to perturbations at approximately 30 times the Kolmogorov scale. This study conducts direct numerical simulations of passive-scalar mixing layers in a shear-free turbulent front to investigate mixing enhancements induced by such perturbations. The initial turbulent Reynolds number based on the Taylor microscale is $ Re_\lambda = 72$ or 202. The turbulent front develops by entraining outer fluid. Weak sinusoidal velocity perturbations are introduced locally, either inside or outside the turbulent front, or globally throughout the flow. Perturbations at this critical wavelength promote small-scale shear instability, complicating the boundary geometry of the scalar mixing layer at small scales. This increases the fractal dimension and enhances scalar diffusion outward from the scalar mixing layer. Additionally, the promoted instability increases the scalar dissipation rate and turbulent scalar flux at small scales, facilitating faster scalar mixing. The effects manifest locally; external perturbations intensify mixing near the boundary, while internal perturbations affect the entire turbulent region. The impact of perturbations is consistent across different Reynolds numbers when the amplitudes normalised by the Kolmogorov velocity are the same, indicating that even weaker perturbations can enhance scalar mixing at higher Reynolds numbers. The mean scalar dissipation rate increases by up to 50 %, even when the perturbation energy is only 2.5 % of the turbulent kinetic energy. These results underscore the potential to leverage small-scale shear instability for efficient mixing enhancement in applications such as chemically reacting flows.

Using clean numerical simulation (CNS) in which artificial numerical noise is negligible over a finite, sufficiently long interval of time, we provide evidence, for the first time, that artificial numerical noise in direct numerical simulation (DNS) of turbulence is approximately equivalent to thermal fluctuation and/or stochastic environmental noise. This confers physical significance on the artificial numerical noise of DNS of the Navier–Stokes equations. As a result, DNS on a fine mesh should correspond to turbulence under small internal/external physical disturbance, whereas DNS on a sparse mesh corresponds to turbulent flow under large physical disturbance. The key point is that all of them have physical meanings and so are correct in terms of their deterministic physics, even if their statistics are quite different. This is illustrated herein. Our paper provides a positive viewpoint regarding the presence of artificial numerical noise in DNS.

We study the onset of spontaneous dynamics in the follower force model of an active filament, wherein a slender elastic filament in a viscous liquid is clamped normal to a wall at one end and subjected to a tangential compressive force at the other. Clarke et al. (Phys. Rev. Fluids, vol. 9, 2024, 073101) recently conducted a thorough investigation of this model using methods of computational dynamical systems; inter alia, they showed that the filament first loses stability via a supercritical double-Hopf bifurcation, with periodic ‘planar-beating’ states (unstable) and ‘whirling’ states (stable) simultaneously emerging at the critical follower-force value. We complement their numerical study by carrying out a weakly nonlinear analysis close to this unconventional bifurcation, using the method of multiple scales. The main outcome is an ‘amplitude equation’ governing the slow modulation of small-magnitude oscillations of the filament in that regime. Analysis of this reduced-order model provides insights into the onset of spontaneous dynamics, including the creation of the nonlinear whirling states from particular superpositions of linear planar-beating modes as well as the selection of whirling over planar beating in three-dimensional scenarios.

A point force acting on a Brinkman fluid in confinement is always counterbalanced by the force on the porous medium, the force on the walls and the stress at open boundaries. We discuss the distribution of those forces in different geometries: a long pipe, a medium with a single no-slip planar boundary, a porous sphere with an open boundary and a porous sphere with a no-slip wall. We determine the forces using the Lorentz reciprocal theorem and additionally validate the results with explicit analytical flow solutions. We discuss the relevance of our findings for cellular processes such as cytoplasmic streaming and centrosome positioning.

Particle-laden horizontal turbulent pipe flow is studied experimentally in the two-way coupling regime with a focus on delineating the effects of particle-to-fluid density ratio $\rho _{p}/\rho _{f}=1$ and 1.05 on the fluid and particle statistics. Particle volume fraction $\phi _{v}$ up to $1\,\%$ and viscous Stokes numbers ranging from $St^+ \approx 1.2$ to $St^+ \approx 3.8$ are investigated at friction Reynolds number $Re_\tau \approx 195$ using time-resolved two-dimensional particle image and tracking velocimetry. Substantial differences are observed between the statistics of neutrally buoyant (i.e. $\rho _{p}/\rho _{f}=1$) and denser (i.e. $\rho _{p}/\rho _{f}=1.05$) settling particles (with settling velocities 0.12–0.32 times the friction velocity), which, at most instances, show opposing trends compared to unladen pipe flow statistics. Neutrally buoyant particles show a slightly increased overall drag and suppressed turbulent stresses, but elevated particle–fluid interaction drag and results in elongated turbulent structures compared to the unladen flow, whereas $\rho _{p}/\rho _{f}=1.05$ particles exhibit a slight overall drag reduction even with increased radial turbulent stresses, and shorter streamwise structures compared to the unladen flow. These differences are enhanced with increasing $St^+$ and $\phi _v$, and can be attributed to the small but non-negligible settling velocity of denser particles, which also leads to differing statistics in the upper and lower pipe halves.

We perform simulations of a two-fluid–structure interaction problem involving liquid–gas flow past a fully submerged stationary circular cylinder. Interactions between the liquid–gas interface with finite surface tension and flow disturbances arising from the cylinder induce a variety of interfacial phenomena and wake structures. We map different interface regimes in a parameter space defined by the Bond number $Bo \in [100, 5000]$ and the submergence depth $h/D \in [1, 2.5]$ of the cylinder while keeping the Reynolds (Re) and Weber (We) numbers fixed at 150 and 1000, respectively. The emerging interface features are classified into three distinct regimes: interfacial waves generated by Strouhal vortices, the entrainment of multi-scale gas bubbles and the reduced deformation state. In the interfacial wave regime, we demonstrate that the frequency of transverse interface fluctuations at a specific streamwise location is identical to the vortex shedding frequency. Additionally, the wavelength of interfacial waves is determined by the size of vortex pairs consisting of alternating Strouhal vortices. In the gas entrainment regime at $ Bo = 1000$, our bubble-size distributions reveal that the entrained bubbles have sizes ranging from one to two orders of magnitude smaller than the cylinder. These multi-scale bubbles are formed primarily through plunging and surfing breakers at $h/D = 2.5$. In contrast, at $h/D = 1$, smaller bubbles initially emerge from the breakup of a gas finger. Over time, some of these bubbles grow in size through coalescence cascades. The influence of $ Re \in [50, 150]$ and $ We \in [700, 1100]$ on gas entrainment is quantified in terms of mean bubble size and count. Lastly, we demonstrate how the deformability of the liquid–gas interface drives the hydrodynamic lift force acting on the cylinder. The net downward lift materializes only in the gas entrainment and reduced deformation regimes due to the broken symmetry of the front stagnation point. While our study focuses on two-dimensional simulations, we also provide insights into the three-dimensional gas entrainment mechanism for one of the extreme cases at $h/D = 1$.

Mangrove restoration efforts have been ongoing, but with varying levels of success, requiring spatial and temporal monitoring to better understand the stocks and drivers of success. Here, we used multi-spectral remote sensing and spatial regression techniques to examine mangrove distribution and restoration potential in the Vietnamese Southern Coastal (VSC) region from 1988 to 2023, an area where multiple episodes of mangrove restoration have been attempted over the past decades. Our results show that 51.5% of the mangrove area has recovered from previous losses, while 48.5% has been lost during the 1988–2023 period. Significant gains were observed between 2018 and 2023, accounting for 77.8% of the total restoration. However, over 40,000 ha of mangroves were lost during each decade between 1988 and 2018, primarily due to land-use changes. Regression analyses estimated a sustainable mangrove cover increase of 9.9% (23,407 ha) and persistence of 22.5% (52,936 ha), mainly in protected areas and low-impact zones. Conversely, 9.8% (23,056 ha) of mangroves in erosion-prone and human-disturbed regions face continued decline. Our study demonstrated the effectiveness of integrating long-term Normalised Difference Vegetation Index time-series analysis with spatial regression to monitor mangrove ecosystems. These techniques offered a scalable framework for global mangrove monitoring and restoration planning, supporting evidence-based conservation policies.