This article extends the previous investigation of the spatial evolution of energy-containing motions in wall-bounded turbulent flows (Kannadasan et al., J. Fluid Mech., vol. 955, 2023, R1) by examining their scale-interactions through spectral analysis based on the spanwise scale decomposition of turbulent kinetic energy and the Reynolds stress transport equation. The energy-containing motions located at the inflow of a turbulent channel flow are artificially removed and the interscale transport mechanisms involved in their spatial evolution are studied. This scale interaction analysis reveals the presence of a significant inverse transfer of streamwise Reynolds stress from the near-wall streaks to larger scales in the spatial evolution of energy-containing motions. This transfer is due to the spanwise variation of streamwise velocity fluctuations, represented by $\partial u'/\partial z$, which is the primary mechanism of streak instability. The analysis presented in this study also shows that the inverse cascade of spanwise energy may correspond to the regeneration of streamwise vortices in the process of reactivating the self-sustaining mechanism in the spatial evolution of energy-containing motions.

The topology of turbulent coherent structures is known to regulate the transport of energy, mass and momentum in the atmospheric boundary layer (ABL). While previous research has primarily focused on characterizing the structure of turbulence in stationary ABL flows, real-world scenarios frequently deviate from stationarity, giving rise to nuanced and poorly understood changes in the turbulence geometry and associated transport mechanisms. This study sheds light on this problem by examining topological changes in ABL turbulence induced by non-stationarity and their effects on momentum transport. Results from a large-eddy simulation of pulsatile open channel flow over an array of surface-mounted cuboids are examined. The analysis reveals that the flow pulsation triggers a phase-dependent shear rate, and the ejection-sweep pattern varies with the shear rate during the pulsatile cycle. From a turbulence structure perspective, it is attributed to the changes in the geometry of hairpin vortices. An increase (decrease) in the shear rate intensifies (relaxes) these structures, leading to an increase (decrease) in the frequency of ejections and an amplification (reduction) of their percentage contribution to the total momentum flux. Furthermore, the size of the hairpin packets undergoes variations, which depend on the geometry of the constituting hairpin vortices, yet the packet inclination preserves its orientation throughout the pulsatile cycle. These observations reinforce the important role non-stationarity holds in shaping the structure of ABL turbulence and the momentum transport mechanisms it governs.

The molten sand that is a mixture of calcia, magnesia, alumina and silicate, known as CMAS, is characterized by its high viscosity, density and surface tension. The unique properties of CMAS make it a challenging material to deal with in high-temperature applications, requiring innovative solutions and materials to prevent its buildup and damage to critical equipment. Here, we use multiphase many-body dissipative particle dynamics simulations to study the wetting dynamics of highly viscous molten CMAS droplets. The simulations are performed in three dimensions, with varying initial droplet sizes and equilibrium contact angles. We propose a parametric ordinary differential equation (ODE) that captures the spreading radius behaviour of the CMAS droplets. The ODE parameters are then identified based on the physics-informed neural network (PINN) framework. Subsequently, the closed-form dependency of parameter values found by the PINN on the initial radii and contact angles are given using symbolic regression. Finally, we employ Bayesian PINNs (B-PINNs) to assess and quantify the uncertainty associated with the discovered parameters. In brief, this study provides insight into spreading dynamics of CMAS droplets by fusing simple parametric ODE modelling and state-of-the-art machine-learning techniques.

We investigate systematically the free settling of a single Platonic polyhedron in an unbounded domain filled with an otherwise quiescent Newtonian fluid. We consider a particle–fluid density mimicking a rock in water. Five Platonic polyhedrons of increasing sphericity are studied for a range of Galileo numbers $10 \leqslant \mathcal {G}a \leqslant 300$. We construct a regime map in the parameter space of Galileo number and particle volume fraction ($\mathcal {G}a, \phi$), highlighting how the angularity of the Platonic polyhedron impacts its settling path and the onset of instabilities. We find that the initial angular position solely affects the transient settling process. All the Platonic polyhedrons maintain a stable settling angular position at low $\mathcal {G}a$. Higher angularity leads to path unsteadiness at lower $\mathcal {G}a$. Path instability progresses from steady vertical to unsteady vertical, followed by oblique settling observed for highly spherical particles, but helical settling (HS) for more angular particles. The particle autorotation is found to be the pivotal factor influencing path instability and the regime transition of angular particles. Beginning in the unsteady oblique and helical regimes, particle autorotation becomes more prevalent, escalating further in the chaotic regime as $\mathcal {G}a$ increases. The particle angular velocity vector is shown to be predominantly situated in the horizontal plane. A thorough force balance in the horizontal plane reveals that the Magnus force is the primary driving force of the HS regime. Additionally, we establish two new empirical correlations to predict the particle settling velocity and the disturbed wake length that solely require the physical properties of the system ($\mathcal {G}a$ and $\phi$). Our numerical results suggest that an increase of the density ratio from $2$ to $3$ exerts only a marginal impact on the path instability of the most angular particle, the settling tetrahedron.

The aeroacoustics of a boundary layer ingesting (BLI) ducted fan is investigated experimentally. The study examines a ducted fan immersed in an adverse streamwise pressure gradient turbulent boundary layer developed over a curved wall. Aeroacoustics measurements indicate that the noise from the BLI ducted fan results from a complex interaction among the fan, duct and the incoming boundary layer. The fundamental mechanisms of noise generation are explained using a general source separation strategy. A detailed noise comparison is made at varying fan rotational speeds and across a wide range of axial inflow velocities. In a low thrust regime, the noise is found to be driven by the fan loading, coupled with duct acoustics and the haystacking phenomenon. In a high thrust regime, the contribution from duct acoustics diminishes, and the noise is predominantly driven by the fan loading coupled with the haystacking phenomenon.

Using low-dimensional numerical simulations, we investigate the characteristics of complex and three-dimensional surface waves in a liquid film flowing over a rotating disk, focusing on large flow rates from a nozzle. Existing integral boundary layer (IBL) models, which are based on spatially averaged variables along the direction normal to the disk surface, have primarily focused on the formation of axisymmetric waves under relatively small flow rates. In this study, an extended IBL model that accounts for both laminar and turbulent regimes is developed by considering the non-uniformity of the local flow rate in the spreading film flow and incorporating closure models dependent on the local Reynolds number. Our numerical results successfully capture the generation of concentric waves by an impinging circular liquid jet and their transition into three-dimensional solitary waves. These findings are in good agreement with visualization images and time-series data of free-surface fluctuations from a displacement sensor. The backscattering of small-scale three-dimensional turbulence into large-scale horizontal turbulence inside the film plays a critical role in determining the transition of wave modes and the nonlinear dynamics of the waves in the turbulent regime. Furthermore, the behaviour of three-dimensional waves in the downstream region, including frequent wave coalescence in the transition region and the breakup of small-scale solitons, is distinct from that of gravity-driven falling film flows. The amplitude of the three-dimensional waves is inversely related to the generalized Reynolds number defined for rotating films.

A revolution in chemical biology occurred with the introduction of click chemistry. Click chemistry plays an important role in protein chemistry modifications, providing specific, sensitive, rapid, and easy-to-handle methods. Under physiological conditions, click chemistry often overlaps with bioorthogonal chemistry, defined as reactions that occur rapidly and selectively without interfering with biological processes. Click chemistry is used for the posttranslational modification of proteins based on covalent bond formations. With the contribution of click reactions, selective modification of proteins would be developed, representing an alternative to other technologies in preparing new proteins or enzymes for studying specific protein functions in different biological processes. Click-modified proteins have potential in diverse applications such as imaging, labeling, sensing, drug design, and enzyme technology. Due to the promising role of proteins in disease diagnosis and therapy, this review aims to highlight the growing applications of click strategies in protein chemistry over the last two decades, with a special emphasis on medicinal applications.

Dense suspensions of solid particles in viscous liquid are ubiquitous in both industry and nature, and there is a clear need for efficient numerical routines to simulate their rheology and microstructure. Particles of micron size present a particular challenge: at low shear rates, colloidal interactions control their dynamics while at high rates, granular-like contacts dominate. While there are established particle-based simulation schemes for large-scale non-Brownian suspensions using only pairwise lubrication and contact forces, common schemes for colloidal suspensions generally are more computationally costly and thus restricted to relatively small system sizes. Here, we present a minimal particle-based numerical model for dense colloidal suspensions that incorporates Brownian forces in pairwise form alongside contact and lubrication forces. We show that this scheme reproduces key features of dense suspension rheology near the colloidal-to-granular transition, including both shear thinning due to entropic forces at low rates and shear thickening at high rates due to contact formation. This scheme is implemented in LAMMPS, a widely used open source code for parallelised particle-based simulations, with a runtime that scales linearly with the number of particles, making it amenable for large-scale simulations.

Detailed observations of the velocities of Jovian vortices exist at only one height in the atmosphere, so their vertical structures are poorly understood. This motivates this study that computes stable three-dimensional, long-lived planetary vortices that satisfy the equations of motion. We solve the anelastic equations with a high-resolution pseudo-spectral method using the observed Jovian atmospheric temperatures and zonal flow. We examine several families of vortices and find that constant-vorticity vortices, which have nearly uniform vorticity as a function of height, and horizontal areas that go to zero at their tops and bottoms, converge to stable vortices that look like the Great Red Spot (GRS) and other Jovian anticyclones. In contrast, the constant-area vortices proposed in previous studies, which have nearly uniform areas as a function of height, and vertical vorticities that go to zero at their tops and bottoms, are far from equilibrium, break apart, and converge to constant-vorticity vortices. Our late-time vortices show unexpected properties. Vortices that are initially non-hollow become hollow (i.e. have local minima of vertical vorticity at their centres), which is a feature of the GRS that cannot be explained with two-dimensional simulations. The central axes of the final vortices align with the planetary spin axis even if they align initially with the local direction of gravity. We present scaling laws for how vortex properties change with the Rossby number and other non-dimensional parameters. We prove analytically that the horizontal mid-plane of a stable vortex must lie at a height above the top of the convective zone.

A dense particle suspension under shear may lose its uniform state to large local density and stress fluctuations, which challenge the mean-field description of the system. Here, we explore the novel dynamics of a non-Brownian suspension under orbital oscillations, where localized density waves along the flow direction appear beyond an excitation frequency threshold and self-organize into a hexagonal pattern across the system. The spontaneous occurrence of the inhomogeneity pattern arises from a coupling between particle advection and the shear-thickening nature of the suspension. Through linear stability analysis, we show that they overcome the stabilizing effects of particle pressure at sufficient particle volume fraction and oscillation frequency. In addition, the long-standing density waves degenerate into random fluctuations when replacing the free surface with rigid confinement. It indicates that the shear-thickened state is intrinsically heterogeneous, and the boundary conditions are crucial for developing local disturbance.

The generation of water waves by wind is a fundamental and much-studied problem of scientific and operational concern. One mechanism that has obtained a wide degree of acceptance and use is the Miles air shear flow instability theory in which water waves grow due to energy transfer from the air at a critical level where the wave phase speed matches the wind speed. In this paper we revisit and extend this theory by examining how the predicted wave growth rate depends on the water depth and the wind shear parameters. At the same time we include frictional effects due to water bottom stress and at the water surface, and add an additional driving term due to wave stress in the air near the sea surface, using a turbulent parametrisation. The theory is developed using a frictional modification of the usual potential flow theory for water waves, the modified Euler equations, coupled to linearised air-flow equations. To assist our analysis we develop a long-wave approximation for the key Miles growth parameter, which explicitly shows how this depends on the water depth and the wind shear parameters. Since our focus is on the development of wave groups, we present a forced nonlinear Schrödinger equation in which the forcing term describes wave growth or decay. For very shallow water the modified Euler equations are reduced to their shallow-water counterpart, which are briefly discussed using Riemann invariants.

The accurate simulation of complex dynamics in fluid flows demands a substantial number of degrees of freedom, i.e. a high-dimensional state space. Nevertheless, the swift attenuation of small-scale perturbations due to viscous diffusion permits in principle the representation of these flows using a significantly reduced dimensionality. Over time, the dynamics of such flows evolves towards a finite-dimensional invariant manifold. Using only data from direct numerical simulations, in the present work we identify the manifold and determine evolution equations for the dynamics on it. We use an advanced autoencoder framework to automatically estimate the intrinsic dimension of the manifold and provide an orthogonal coordinate system. Then, we learn the dynamics by determining an equation on the manifold by using both a function-space approach (approximating the Koopman operator) and a state-space approach (approximating the vector field on the manifold). We apply this method to exact coherent states for Kolmogorov flow and minimal flow unit pipe flow. Fully resolved simulations for these cases require $O(10^3)$ and $O(10^5)$ degrees of freedom, respectively, and we build models with two or three degrees of freedom that faithfully capture the dynamics of these flows. For these examples, both the state-space and function-space time evaluations provide highly accurate predictions of the long-time dynamics in manifold coordinates.

When two fluid drops touch, they coalesce due to surface tension. At early times, there is only a relatively small fluid bridge joining the drops. An asymptotic solution is presented for an inertial regime of early-time coalescence, in which inertial forces balance surface tension at leading order. It is demonstrated that viscosity nevertheless has a leading-order effect. Radial momentum is created at the tightly curved edge of the fluid bridge by the net force $2\gamma$ (per unit length) due to surface tension. This momentum is left behind the radially expanding bridge edge in a thin viscous wake. The divergent volume flux in the wake entrains fluid from above and below the bridge, and drives an inviscid irrotational flow in the drops on the scale of the bridge radius. This flow widens the gap between the drops ahead of the bridge, and the larger gap width results in a lower rate of coalescence. Including viscosity in this way improves the agreement between theory and the available experimental and numerical data.

We consider a thin horizontal layer of a non-magnetic electrolyte containing a bulk solution of salt and carrying an electric current. The layer is bounded by two deformable free surfaces loaded with an insoluble surfactant and is placed in a vertical magnetic field. The arising Lorentz force drives the electrolyte in the plane of the layer. We employ the long-wave approximation to derive general two-dimensional hydrodynamic equations describing symmetric pinching-type deformations of the free surfaces. These equations are used to study the azimuthal flow in an annular film spanning the gap between two coaxial cylindrical electrodes. In weakly deformed films, the base azimuthal flow and its linear stability with respect to azimuthally invariant perturbations are studied analytically. For relatively thick layers and weak magnetic fields, the leading mode with the smallest decay rate is found to correspond to a monotonic azimuthal velocity perturbation. The Marangoni effect leads to further stabilisation of the flow while perturbations of the solute concentration in the bulk of the fluid have no influence on the flow stability. In strongly deformed films in the diffusion-dominated regime, the azimuthal flow becomes linearly unstable with respect to an oscillatory mixed mode characterised by the combination of radial and azimuthal velocity perturbations when the voltage applied between electrodes exceeds the critical value.

Shock waves are of great interest in many fields of science and engineering, but the mechanisms of their formation, maintenance and dissipation are still not well understood. While all transport processes existing in a shock wave contribute to its compression and irreversibility, they are not of equal importance. To figure out the roles of viscosity and heat conduction in shock transition, the existence of smooth shock solutions and the counter-intuitive entropy overshoot phenomenon (the specific entropy is not monotonically increasing and exhibits a peak inside the shock front) are theoretically and numerically investigated, with emphasis on the effects of viscosity and heat conduction. Instead of higher-order hydrodynamics, the Navier–Stokes formalism is employed for its stability and simplicity. Supplemented with nonlinear thermodynamically consistent constitutive relations, the Navier–Stokes equations are adequate to demonstrate the general nature of shock profiles. It is found that heat conduction cannot sustain strong shocks without the presence of viscosity, while viscosity can maintain smooth shock transition at all strengths, regardless of heat conduction. Hence, the critical role in shock compression is played by viscosity rather than heat conduction. Nevertheless, the dispensability of heat conduction would not compromise its essential role in the emergence of an entropy peak. It is the entropy flux resulting from heat conduction that neutralises the positive entropy production and thus prevents the decreasing entropy from violating the second law of thermodynamics. This mechanism of entropy overshoot has not been addressed previously in the literature and is revealed using the entropy balance equation.

Transitional separating flow induced by a rectangular plate subjected to uniform incoming flow at Reynolds number (based on the incoming velocity and half plate height) 2000 is investigated using direct numerical simulation. The objective is to unveil the long-lasting mystery of low-frequency flapping motion (FM) in flow separation. At a fixed streamwise-vertical plane or from the perspective of previous experimental studies using pointwise or planar measurements, FM manifests as a low-frequency periodic switching between low and high velocities covering the entire separation bubble. The results indicate that in three-dimensional space, FM reflects an intricate evolution of streamwise elongated streaky structures under the influence of separated shear layer and mean flow reversal. The FM is an absolute instability, and is initiated through a lift-up mechanism boosted by mean flow deceleration near the crest of the separating streamline. At this particular location, the shear bends the vortex filament abruptly, so that one end is vertically struck into the first half of the separation bubble, whereas the other end is extended in the streamwise direction in the second half of the separation bubble. These two ends of vortex filament are mutually sustained and also stretched by the vertical acceleration and streamwise acceleration in the first and second halves of the separation bubble, respectively. This process periodically switches the low-velocity (or high-velocity) streaky structure to a high-velocity (or low-velocity) streaky structure encompassing the entire separation bubble, and thus flaps the separated shear layer up and down in the vertical direction. A ‘large vortex’ shedding manifests when the streaky structure switches signs. This large vortex is fundamentally different from the spanwise vortex shedding residing in the shear layer originated from the Kelvin–Helmholtz instability and successive vortex amalgamation. It is also believed that the three-dimensional evolution of streaky structures in the form of FM is applicable for both geometry- and pressure-induced separating flows.

We performed numerical simulations of a homogeneous swarm of bubbles rising at large Reynolds number, $Re=760$, with volume fractions ranging from 1 % to 10 %. We consider a simplified model in which the interfaces are not resolved, but which allows us to simulate flows with a large number of bubbles and to emphasize the interactions between bubble wakes. The liquid phase is described by solving, on an Eulerian grid, the Navier–Stokes equations, including sources of momentum which model the effect of the bubbles. The dynamics of each bubble is determined within the Lagrangian framework by solving an equation of motion involving the hydrodynamic forces exerted by the fluid accounting for the correction of the fictitious self-interaction of a bubble with its own wake. The comparison with experiments shows that this coarse-grained simulations approach can reliably describe the dynamics of the resolved flow scales. We use conditional averaging to characterize the mean bubble wakes and obtain in particular the typical shear imposed by the rising bubbles. On the basis of the spectral decomposition of the energy budget, we observe that the flow is dominated by production at large scales and by dissipation at small scales and we rule out the presence of an intermediate range in which the production and dissipation are locally in balance. We propose that the $k^{-3}$ subrange of the energy spectra results from the mean shear rate imposed by the bubbles, which controls the rate of return to isotropy.