Plastic pollution is recognised as one of this century’s most significant environmental challenges and has the characteristics of a super wicked problem. Though researchers and governments around the globe are coming up with promising technological interventions, awareness among citizens and stringent policies are the need of the hour to tackle this issue. A few countries have issued postage stamps and postal materials showcasing the various dimensions of plastic pollution. Historically, stamps depicted every progress, problem and various milestones of humanity spanning multiple fields. We contend that the plastic pollution problems and impact should be depicted through postage stamps from all countries. Through this feat, the message of the need for sustainable usage of plastics for the common good of all species can be spread by showcasing various dimensions of the sustainability of plastic usage in postage materials. This article discusses the rise of plastic pollution, its emerging impacts, and contemporary issues and mitigation strategies through postage stamps and materials. Philately can be a medium for providing environmental awareness, considering the case of plastic pollution. It can be a strong driver to promote consciousness regarding various environmental problems among students undergoing multiple levels of education and the general public.

We investigate the effects of thermal boundary conditions and Mach number on turbulence close to walls. In particular, we study the near-wall asymptotic behaviour for adiabatic and pseudo-adiabatic walls, and compare to the asymptotic behaviour recently found near isothermal cold walls (Baranwal et al. 2022. J. Fluid Mech. 933, A28). This is done by analysing a new large database of highly-resolved direct numerical simulations of turbulent channels with different wall thermal conditions and centreline Mach numbers. We observe that the asymptotic power-law behaviour of Reynolds stresses as well as heat fluxes does change with both centreline Mach number and thermal condition at the wall. Power-law exponents transition from their analytical expansion for solenoidal fields to those for non-solenoidal field as the Mach number is increased, though this transition is found to be dependent on the thermal boundary conditions. The correlation coefficients between velocity and temperature are also found to be affected by these factors. Consistent with recent proposals on universal behaviour of compressible turbulence, we find that dilatation at the wall is the key scaling parameter for these power-law exponents, providing a universal functional law that can provide a basis for general models of near-wall behaviour.

With coastal populations rising at three times the global average, sustainable ways of safeguarding human needs around access and use of the coast alongside lasting ecosystem health of coastal environments must be developed. At the same time, human populations are facing the challenge of managing coastal access on the back of a legacy of human interventions that have already altered – and have often had unintended or unforeseen impacts on – the coastal system and its functioning.

We chart the history of the evolution of North Bull Island in Dublin Bay as an example of major unforeseen sedimentation in a coastal estuarine bay following the construction of river mouth training walls. We investigate the impact of a constructed causeway on the evolved ‘naturescape’ by comparing accretion and elevation change on the mid-marsh either side of the access road over a 32-month period (autumn 2021 to summer 2024) and measuring water levels either side of the causeway on six spring tides on consecutive days characterised by varying meteorological conditions in early September 2023. The results allow us to consider the potential implications a lack of physical connectivity may cause for the future of the two artificially separated back-barrier lagoon environments.

Debris flows are a growing natural hazard as a result of climate change and population density. To effectively assess this hazard, simulating field-scale debris flows at a reasonable computational cost is crucial. We enhance existing debris flow models by rigorously deriving a series of depth-averaged shallow models with varying complexities describing the behaviour of grain–fluid flows, considering granular mass dilatancy and pore fluid pressure feedback. The most complete model includes a mixture layer with an upper fluid layer, and solves for solid and fluid velocity in the mixture and for the upper fluid velocity. Simpler models are obtained by assuming velocity equality in the mixture or single-layer descriptions with a virtual thickness. Simulations in a uniform configuration mimicking submarine landslides and debris flows reveal that these models are extremely sensitive to the rheology, the permeability (grain diameter) and initial volume fraction, parameters that are hard to measure in the field. Notably, velocity equality assumptions in the mixture hold true only for low permeability (corresponding to grain diameter $d=10^{-3}$ m). The one-layer models’ results can strongly differ from those of the complete model, for example, the mass can stop much earlier. One-layer models, however, provide a rough estimate of two-layer models when permeability is low, initial volume fraction is far from critical and the upper fluid layer is very thin. Our work with uniform settings highlights the need of developing two-layer models accounting for dilatancy and for an upper layer made either of fluid or grains.

An analytical theory is developed that describes acoustic microstreaming produced by the interaction of an oscillating gas bubble with a viscoelastic particle. The bubble is assumed to undergo axisymmetric oscillation modes, which can include radial oscillation, translation and shape modes. The oscillations of the particle are excited by the oscillations of the bubble. No restrictions are imposed on the ratio of the bubble and the particle radii to the viscous penetration depth and the separation distance, as well as on the ratio of the viscous penetration depth to the separation distance. Capabilities of the developed theory are illustrated by computational examples. The shear stress produced by the acoustic microstreaming on the particle’s surface is calculated. It is shown that this stress is much higher than the stress predicted by Nyborg’s formula (1958 J. Acoust. Soc. Am. 30, 329–339), which is commonly used to evaluate the time-averaged shear stress produced by a bubble on a rigid wall.

This paper presents detailed analyses of the Reynolds stresses and their budgets in temporally evolving stratified wakes using direct numerical simulation. Ensemble averaging is employed to mitigate statistical errors in the data, and the results are presented as functions of both the transverse and vertical coordinates – at time instants across the near-wake, non-equilibrium, and quasi-two-dimensional regimes for wakes in weakly and strongly stratified environments. Key findings include the identification of dominant terms in the Reynolds stress transport equations and their spatial structures, the generation and destruction processes of the Reynolds stresses, and the energy transfer between the Reynolds stress and the mean flow. The study also clarifies the effects of the Reynolds number and the Froude number. Additionally, we assess the validity of the eddy-viscosity type models and some existing closures for the Reynolds stress model, highlighting the limitations of isotropy and return-to-isotropy hypotheses in stratified flows.

Turbulent mixing in a supercritical CO$_2$ shear layer is examined using both experimental and numerical methods. Boundary conditions are selected to focus on the rarely studied near-critical regime, where thermophysical properties vary nonlinearly with respect to temperature and pressure. Experimental results are obtained via Raman spectroscopy and shadowgraphy, while numerical results are obtained via direct numerical simulation. The shear layer growth rate is found to be 0.2. Additionally, density profiles indicate a relaxation of density gradients between the mixed fluid and heavy fluid as the flow evolves downstream, which runs counter to existing supercritical shear layer data in the literature. The computational results identify significant anisotropy in the turbulence in the shear layer, which is discussed in terms of the development of regions of high density gradient magnitude. The Reynolds-averaged enstrophy budget at various streamwise locations indicates no significant dilatational or baroclinic contribution within the shear layer.

An important parameter characterising the synchronisation of turbulent flows is the threshold coupling wavenumber. This study investigates the relationship between the threshold coupling wavenumber and the leading Lyapunov vector using large eddy simulations and the SABRA model. Various subgrid-scale stress models, Reynolds numbers and different coupling methods are examined. A new scaling relation is identified for the leading Lyapunov exponents in large eddy simulations, showing that they approximate those of filtered direct numerical simulations. This interpretation provides a physical basis for results related to the Lyapunov exponents of large eddy simulations, including those related to synchronisation. Synchronisation experiments show that the peak wavenumber of the energy spectrum of the leading Lyapunov vector coincides with the threshold coupling wavenumber, in large eddy simulations of box turbulence with standard Smagorinsky or dynamic mixed models as well as in the SABRA model, replicating results from direct numerical simulations of box turbulence. Although the dynamic Smagorinsky model exhibits different behaviour, the totality of the results suggests that the relationship is an intrinsic property of a certain class of chaotic systems. We also confirm that conditional Lyapunov exponents characterise the synchronisation process in indirectly coupled systems as they do in directly coupled ones, with their values insensitive to the nature of the master flow. These findings advance the understanding of the role of the Lyapunov vector in the synchronisation of turbulence.

Estimation of near-wall turbulence in channel flow from outer observations is investigated using adjoint-variational data assimilation. We first consider fully resolved velocity data, starting at a distance from the wall. By enforcing the estimated flow to exactly satisfy the Navier–Stokes equations, we seek a statistically stationary turbulent state that reproduces the instantaneous outer measurements. Such an estimated state provides full access to the unknown near-wall turbulence, including the wall shear stresses and pressure. When the first observation is within 50 viscous units from the wall, the correlation coefficient between the true and estimated state exceeds 95 %. As the observations are further separated from the wall, at 90 viscous units, the accuracy of the assimilated wall stresses decreases to 40 % at the wall. This trend is nearly independent of the Reynolds number. The Fourier spectrum of the estimation error is qualitatively consistent with the coherence spectrum between the outer and the inner state variables: observed long wavelength structures in the outer flow have deeper coherence into the unobserved near-wall region and, therefore, the error is lowest at large scales. Nevertheless, the adjoint-variational approach provides a more rigorous quantification of the capacity to accurately predict the instantaneous near-wall turbulence from outer measurements. Lastly, we demonstrate the robustness of the estimation accuracy using filtered and sub-sampled outer observations.

When turbulent convection interacts with a turbulent shear flow, the cores of convective cells become aligned with the mean current, and these cells (which span the height of the domain) may interact with motions closer to the solid boundary. In this work, we use coupled Eulerian–Lagrangian direct numerical simulations of a turbulent channel flow to demonstrate that, under conditions of turbulent mixed convection, interactions between motions associated with ejections and low-speed streaks near the solid boundary and coherent superstructures in the interior of the flow interact and lead to significant vertical transport of strongly settling Lagrangian particles. We show that the primary suspension mechanism is associated with strong ejection events (canonical low-speed streaks and hairpin vortices characterised by $u'\lt 0$ and $w'\gt 0$, where $u'$ and $w'$ are the streamwise and vertical turbulent velocity fluctuations), whereas secondary suspension is strongly associated with large-scale plume structures aligned with the mean shear (characterised by $w'\gt 0$ and $\theta '\gt 0$, where $\theta$ represents temperature fluctuations). This coupling, which is absent in the limiting cases (pure channel flow or free convection) is shown to lead to a sudden increase in the interior concentration profiles as ${Ri}_\tau$, the friction Richardson number, increases, resulting in concentrations that are larger by roughly an order of magnitude at the channel midplane.

It is generally accepted that the evolution of the deep-water surface gravity wave spectrum is governed by quartet resonant and quasi-resonant interactions. However, it has also been reported in both experimental and computational studies that non-resonant triad interactions can play a role, e.g. generation of bound waves. In this study, we investigate the effects of triad and quartet interactions on the spectral evolution, by numerically tracking the contributions from quadratic and cubic terms in the dynamical equation. In a finite time interval, we find that the contribution from triad interactions follows the trend of that from quartet resonances (with comparable magnitude) for most wavenumbers, except that it peaks at low wavenumbers with very low initial energy. This result reveals two effects of triad interactions. (1) The non-resonant triad interactions can be connected to form quartet resonant interactions (hence exhibiting the comparable trend), which is a reflection of the normal form transformation applied in wave turbulence theory of surface gravity waves. (2) The triad interactions can fill energy into the low-energy portion of the spectrum (low wavenumber part in this case) on a very fast time scale, with energy distributed in both bound and free modes at the same wavenumber. We further analyse the latter mechanism using a simple model with two initially active modes in the wavenumber domain. Analytical formulae describing the distribution of energy in free and bound modes are provided, along with numerical validations.

Thermal Marangoni effects play important roles in bubble dynamics such as bubbles generated by water electrolysis or boiling. As macroscopic bubbles often originate from nucleated nanobubbles, it is crucial to understand how thermocapillarity operates at the nanoscale. In this study, the motion of transient bulk gas nanobubbles in water driven by a vertical temperature gradient between two solid plates is investigated using molecular dynamics simulations and analytical theory. The simulation results show that due to the thermal Marangoni force, nanobubbles move towards the hot plate at a constant velocity, similar to the behaviour of macroscale gas bubbles. However, unlike macroscale gas bubbles whose thermal conductivity and viscosity are negligible compared to those of water, the thermal conductivity and viscosity of nanoscale gas bubbles are significantly increased due to their large gas density. The thermal resistance and the slip length are also found to matter at the liquid–gas interface, though they decrease with increasing gas densities. The previous thermocapillary theory for macroscale bubbles is extended by considering all these nanoscopic effects. Expressions of the Marangoni force and the drag force are derived. By balancing the Marangoni force and the drag force, the theoretical velocity of the nanobubble migration in a thermal gradient is obtained. When using the measured transport properties of liquid, gas, and their interfaces, the theoretically obtained velocity is consistent with the result of the molecular simulations. We find that the slip length is too small to have considerable effects on nanobubble motions in the current liquid–gas system.

An important feature of the dynamics of double-diffusive fluids is the spontaneous formation of thermohaline staircases, where wide regions of well-mixed fluid are separated by sharp density interfaces. Recent developments have produced a number of one-dimensional reduced models to describe the evolution of such staircases in the salt fingering regime relevant to mid-latitude oceans; however, there has been significantly less work done on layer formation in the diffusive convection regime. We aim to fill this gap by presenting a new model for staircases in diffusive convection based on a regularisation of the $\gamma$-instability (Radko 2003 J. Fluid Mech. vol. 805, 147–170), with a range of parameter values relevant to both polar oceans and astrophysical contexts. We use the results of numerical simulations to inform turbulence-closure parametrisations as a function of the horizontally averaged kinetic energy $e$, and ratio of the haline to thermal gradients $R_0^*$. These parametrisations result in a one-dimensional model that reproduces the critical value of $R_0^*$ for the layering instability, and the spatial scale of layers, for a wide range of parameter values, although there is a mismatch between the range of $R_0^*$ for layer formation in the model and observational values from polar oceans. Staircases form in the one-dimensional model, evolving gradually through layer merger events that closely resemble simulations.

Surface gravity waves induce a drift on objects floating on the water’s surface. This study presents laboratory experiments investigating the drift of large two-dimensional floating objects on deep-water, unidirectional, regular waves, with wave steepness ranging from 0.04 to 0.31 (0.04 $\lt k{a_w}\lt$ 0.31, where $k$ is the wavenumber and $a_w$ the wave amplitude). The objects were carefully designed to have a rectangular cross-section with a constant aspect ratio; their size varied from 2.6 $\%$ to 27 $\%$ of the incident wavelength. We observed Lagrangian behaviour for small objects. Small and large objects exhibited fundamentally different drift behaviour at high compared with low wave steepness, with a regime shift observed at a certain size and wave steepness. The scaling of object drift with steepness depends on the relative size of the object. For small objects, drift scales with steepness squared, whereas drift becomes a linear function of steepness as the object size increases. For objects that are relatively large but smaller than 13–16% of a wavelength (low to high steepness), we provide experimental evidence supporting the mechanisms of drift enhancement recently identified by Xiao et al. (J. Fluid Mech., vol. 980, 2024, p. A27) and termed the ‘diffraction-modified Stokes drift’. This enhanced drift behaviour, compared with the theoretical Stokes drift for infinitely small fluid parcels, is attributed to changes in the objects’ oscillatory motion and local wave amplitude distribution (standing wave pattern) due to the presence of the object. In the case of larger objects, similar to Harms (J. Waterw. Port Coast. Ocean Eng., vol. 113(6), 1987, pp. 606–622), we relate the critical size at which drift is maximised to their vertical bobbing motion. We determine the domain of validity for both Stokes drift and the diffraction-modified Stokes drift model of Xiao et al. (J. Fluid Mech., vol. 980, 2024, A27) in terms of relative size and wave steepness and propose an empirical parametrisation based on our experimental data.

Granular surface flows are frequently encountered in nature as well as during handling of powders in different industries. An experimental study of granular surface flow on a heap is carried out. The heap is formed by pouring nearly monodisperse spherical particles from the rectangular slit orifice of a hopper on a rough rectangular plate. A flow of particles is developed on the heap surface, which is planar in the central region, with particles flowing over the edge of the plate into a collection chamber. The geometry considered in this study is an example of a fully three-dimensional heap without side walls. The surface velocities of the particles are measured using high-speed videography and particle tracking velocimetry for different mass flow rates with steel balls and glass beads of two different sizes, for heaps of different aspect ratios. The flow is uni-directional and fully developed in a central zone on the heap surface. The flowing layer thickness is measured in this zone by immersing a soot-coated blade into the flow. The angle of inclination of the free surface of the heap is found to be nearly constant for a ten-fold increase of the mass flow rate. The scaled flowing layer thickness is found to vary linearly with the scaled flow rate and the data for all the particles collapse to a single line over a ten-fold increase in the scaled flow rate and an increase in the aspect ratio of the heap by a factor of 1.75. The predicted scaled surface velocity and scaled shear rate using this correlation match the measured values.

The laminar flow past rectangular prisms is studied in the space of length-to-height ratio ($1 \leqslant L/H \leqslant 5$), width-to-height ratio ($1.2 \leqslant W/H \leqslant 5$) and Reynolds number ($Re \lessapprox 700$); $L$ and $W$ are the streamwise and cross-flow dimensions of the prisms. The primary bifurcation is investigated with linear stability analysis. For large $W/L$, an oscillating mode breaks the top/bottom planar symmetry. For smaller $W/L$, the flow becomes unstable to stationary perturbations and the wake experiences a static deflection, vertical for intermediate $W/L$ and horizontal for small $W/L$. Weakly nonlinear analysis and nonlinear direct numerical simulations are used for $L/H = 5$ and larger $Re$. For $W/H = 1.2$ and 2.25, the flow recovers the top/bottom planar symmetry but loses the left/right one, via supercritical and subcritical pitchfork bifurcations, respectively. For even larger $Re$, the flow becomes unsteady and oscillates around either the deflected (small $W/H$) or the non-deflected (intermediate $W/H$) wake. For intermediate $W/H$ and $Re$, a fully symmetric periodic regime is detected, with hairpin vortices shed from the top and bottom leading-edge (LE) shear layers; its triggering mechanism is discussed. At large $Re$ and for all $W/H$, the flow approaches a chaotic state characterised by the superposition of different modes: shedding of hairpin vortices from the LE shear layers, and wake oscillations in the horizontal and vertical directions. In some portions of the parameter space the different modes synchronise, giving rise to periodic regimes also at relatively large $Re$.