We compute the full set of second-order inertial corrections to the instantaneous force and torque acting on a small spherical rigid particle moving unsteadily in a general steady linear flow. This is achieved by using matched asymptotic expansions and formulating the problem in a coordinate system co-moving with the background flow. Effects of unsteadiness and fluid-velocity gradients are assumed to be small, but to dominate in the far field over those of the velocity difference between the body and fluid, making the results essentially relevant to weakly positively or negatively buoyant particles. The outer solution (which at first order is responsible for the Basset–Boussinesq history force at short time and for shear-induced forces such as the Saffman lift force at long time) is expressed via a flow-dependent tensorial kernel. The second-order inner solution brings a number of different contributions to the force and torque. Some are proportional to the relative translational or angular acceleration between the particle and fluid, while others take the form of products of the rotation/strain rate of the background flow and the relative translational or angular velocity between the particle and fluid. Adding the outer and inner contributions, the known added-mass force or the spin-induced lift force are recovered, and new effects involving the velocity gradients of the background flow are revealed. The resulting force and torque equations provide a rational extension of the classical Basset–Boussinesq–Oseen equation incorporating all first- and second-order fluid inertia effects resulting from both unsteadiness and velocity gradients of the carrying flow.

Understanding the statistics of bedload particle motions is of great importance. To model the hop events which are defined as trajectories of particles moving successively from the start to the end of their motions, recently, Wu et al. (Water Resour. Res., vol. 56, 2020, p. e2019WR025116) have successfully performed individual-based simulations according to the Fokker–Planck equation for particle velocities. However, analytical solutions are still not available due to (i) difficulties in treating the velocity-dependent diffusivity, and (ii) a knowledge gap in incorporating the termination of particle motions for the equation. To tackle the above-mentioned challenges, we first specify a Robin boundary condition representing the deposition of particles. Second, for analytical solutions of hop statistics, a variable transformation is devised to deal with the velocity-dependent diffusivity. The original bedload transport problem is thus found to be governed by the classic equation for the solute transport in tube flows with a constant diffusivity after the transformation. Finally, through solving the spatial and temporal moments of the governing equation, we investigate the influence of the deposition rate on three key characteristics of particle hops. Importantly, we have related the deposition rate to the mean travel times and hop distances, enabling a direct determination of this physical parameter based on measured particle motion statistics. The analytical solutions are validated by experimental observations with different bedload particle diameters and transport conditions. Based on the limited experimental datasets, the deposition frequency is shown to decrease as the shear stress increases when the flow rate is not small.

Using laboratory experiments and numerical simulations, we examine the transfer of soluble material from small, spherical particles sinking in homogeneous turbulence at large Péclet number. A theoretical analysis predicts two distinct mechanisms of convective mass transfer: strain due to turbulence and slip due to gravitational settling. Their relative strength is parametrised by the sinking ratio, ${Sr} = w_0 \tau _\eta /a$, where $w_0$ is the quiescent settling velocity, $a$ is the particle radius and $\tau _\eta$ is the Kolmogorov time scale. This analysis predicts that the topology of the concentration wake changes from a symmetric topology at ${Sr} \ll 1$ to an asymmetric topology at ${Sr} \gg 1$ as the dominant mechanism of mass transfer changes. Particle tracking flow visualisations of small spheres releasing dye in turbulence confirm the existence of this change in mechanism at ${Sr} = O(1)$. We complement these experiments with numerical simulations of the mass transfer from sinking particles. The transfer rate predicted by the simulation is found to be in good agreement with literature data for mass transfer to turbulent suspensions of solid particles and is consistent with asymptotic expressions for mass transfer in uniform flow when ${Sr} \gg 1$. A decomposition of the convective fluxes confirms the transition in the transfer mechanism. At ${Sr} = O(1)$, both mechanisms provide comparable contributions to the transfer rate. Cross-correlation analysis reveals that particle-scale knowledge of both the recent strain and velocity history is required to predict the instantaneous transfer rate. Turbulence-induced particle rotation has a modest suppression effect upon convective transfer by sinking.

We study how channel width variations influence the dynamics of free-surface granular flows. For this purpose, we extend a continuum model framework to granular flows passing through channels that narrow or widen. Our theory uses a linearized approximation to an established dense granular flow rheology and a Coulomb friction law to model interaction between flow and sidewalls. We test the theoretical predictions using two novel 40 cm-diameter drums (convex and concave) filled halfway with 2 mm diameter particles rotated at rates in which the shear layer remains shallow and dense. We apply particle tracking velocimetry to enable quantitative comparisons between experimental data and theoretical predictions. We find that our experimental kinematics and energy profiles largely agree with the theoretical predictions. In general, flows through narrowing channels are faster and deeper than flows through widening channels. The influence of width variations grows with increasing flow speed, and the form of the rate dependence changes fundamentally as the regime changes from one in which kinetic energy is dissipated locally to one in which it is advected downstream. For both regimes, theoretical scaling analysis leads us to experimentally validated power laws, in which the exponent depends on the flow regime, and the multiplicative coefficient depends on channel geometry alone. Finally, we discuss how the differences between theoretical predictions and experimental data may be useful for improving our understanding of flows through non-uniform channels.

The highly nonlinear evolution of the single-mode stratified compressible Rayleigh–Taylor instability (RTI) is investigated via direct numerical simulation over a range of Atwood numbers ($A_T=0.1$–$0.9$) and Mach numbers ($Ma=0.1$–$0.7$) for characterising the isothermal background stratification. After the potential stage, it is found that the bubble is accelerated to a velocity which is well above the saturation value predicted in the potential flow model. Unlike the bubble re-acceleration behaviour in quasi-incompressible RTI with uniform background density, the characteristics in the stratified compressible RTI are driven by not only vorticity accumulation inside the bubble but also flow compressibility resulting from the stratification. Specifically, in the case of strong stratification and high $A_T$, the flow compressibility dominates the bubble re-acceleration characters. To model the effect of flow compressibility, we propose a novel model to reliably describe the bubble re-acceleration behaviours in the stratified compressible RTI, via introducing the dilatation into the classical model that takes into account only vorticity accumulation.

Classical theories for suspensions have been formulated by starting from the Navier–Stokes equations describing pure liquid flow and then introducing additional dependencies to account for the presence of suspended particles. These models are often accurate for low particle concentrations but have lacked a convincing description of the frictional interactions of particles, which are important at larger solid volume fractions. The $\mu (J), \varPhi (J)$ rheology, which draws a direct analogy between suspension flow and dry granular flow, is a recent theory that addresses this issue, but is shown here to be dynamically ill-posed for large solid volume fractions. An alternative well-posed theory is introduced that includes additional dependence on the particle-phase dilation and compression. The new theory, denoted vCIDR, is tested numerically to show grid convergence for problems in which the $\mu (J), \varPhi (J)$ rheology instead suffers from catastrophic blow-up. A further well-posed extension provides a framework for handling the transition between viscous and inertial flows.

Experiments are conducted in water films falling along the bottom wall of a weakly inclined rectangular channel of height $5$ mm in the presence of a laminar counter-current air flow. Boundary conditions have been specifically designed to avoid flooding at the liquid outlet, thus allowing us to focus on the wave dynamics in the core of the channel. Surface waves are excited via coherent inlet forcing before they come into contact with the air flow. The effect of the air flow on the height, shape and speed of two-dimensional travelling nonlinear waves is investigated and contrasted with experiments of Kofman, Mergui & Ruyer-Quil (Intl J. Multiphase Flow, vol. 95, 2017, pp. 22–34), which were performed in a weakly confined channel. We observe a striking difference between these two cases. In our strongly confined configuration, a monotonic stabilizing effect or a non-monotonic trend (i.e. the wave height first increases and then diminishes upon increasing the gas flow rate) is observed, in contrast to the weakly confined configuration where the gas flow is always destabilizing. This stabilizing effect implies the possibility of attenuating waves via the gas flow and it confirms recent numerical results obtained by Lavalle et al. (J. Fluid Mech., vol. 919, 2021, R2) in a superconfined channel.

We consider the propagation of inviscid bubbles in a Hele-Shaw cell under a uniform background flow. We focus on the distinguished limit in which the hydrodynamic pressure gradient due to the external flow balances viscous drag effects due to thin liquid films between the bubbles and the cell walls (Bretherton, J. Fluid Mech., vol. 10, issue 2, 1961, pp. 166–188), with the ratio between these two effects measured by a single dimensionless parameter that we label $\delta$. In this regime, we find that each bubble remains approximately circular, and its propagation velocity is determined by a net force balance. The analytical solution for the problem of an isolated bubble in an infinite Hele-Shaw cell is found to agree well with experimental data in the literature. In particular, we find that the bubble may move faster or slower than the background fluid speed, depending on whether $\delta >1$ or $\delta <1$, or precisely with the background flow if $\delta =1$. When the model is generalised to include the effects of multiple bubbles and boundaries in the Hele-Shaw cell, we still find that the sign of $\delta -1$ causes striking changes in the qualitative behaviour. For a train of three or more bubbles moving along a Hele-Shaw channel, we observe longitudinal waves that propagate forwards or backwards along the bubble train, depending on whether $\delta >1$ or $\delta <1$, resembling a Hele-Shaw Newton's cradle.

A two-dimensional depth-averaged model is derived for open-channel flows in the smooth turbulent case. The derivation is consistently obtained with a method of matched asymptotic expansions in the outer and inner layers using a mixing length model of turbulence including the free-surface reduction of the eddy viscosity. The shear effects are taken into account by an extra tensor variable of the model called the enstrophy tensor. The friction coefficient is an explicit expression of the water depth. The three-dimensional velocity field and the friction velocity can be reconstructed from the values of the depth-averaged quantities. Numerical simulations show that the enstrophy can be used to evaluate the development of the turbulent boundary layer. In the case of subcritical unsteady flows, the reconstructed velocity can be described with a logarithmic law modified by Coles’ wake function with apparent von Kármán constant, integration constant and wake-strength parameter, which differ from their values in steady flows. In the viscous sublayer the steady-state relation between the velocity and the vertical coordinate, in the inner scaling, is not valid for unsteady flows. Large errors on the calculation of the von Kármán constant can be made if the validity of the steady-state relation is assumed for unsteady flows. The comparisons of the reconstructed velocity profiles in the case of one-dimensional unsteady open-channel flows and two-dimensional wide trapezoidal channels show a good agreement with experiments.

Convection occurs ubiquitously on and in rotating geophysical and astrophysical bodies. Prior spherical shell studies have shown that the convection dynamics in polar regions can differ significantly from the lower latitude, equatorial dynamics. Yet most spherical shell convective scaling laws use globally-averaged quantities that erase latitudinal differences in the physics. Here we quantify those latitudinal differences by analysing spherical shell simulations in terms of their regionalized convective heat-transfer properties. This is done by measuring local Nusselt numbers in two specific, latitudinally separate, portions of the shell, the polar and the equatorial regions, $Nu_p$ and $Nu_e$, respectively. In rotating spherical shells, convection first sets in outside the tangent cylinder such that equatorial heat transfer dominates at small and moderate supercriticalities. We show that the buoyancy forcing, parameterized by the Rayleigh number $Ra$, must exceed the critical equatorial forcing by a factor of ${\approx }20$ to trigger polar convection within the tangent cylinder. Once triggered, $Nu_p$ increases with $Ra$ much faster than does $Nu_e$. The equatorial and polar heat fluxes then tend to become comparable at sufficiently high $Ra$. Comparisons between the polar convection data and Cartesian numerical simulations reveal quantitative agreement between the two geometries in terms of heat transfer and averaged bulk temperature gradient. This agreement indicates that rotating spherical shell convection dynamics is accessible both through spherical simulations and via reduced investigatory pathways, be they theoretical, numerical or experimental.

We derive a limit on energy savings in controlled channel flow. For flow in a channel driven by pressure, shear or any combination of the two, and controlled via wall transpiration or spanwise wall motion, the uncontrolled laminar state requires the least net energy (accounting for the energetic cost of control). Thus, the optimal control solution is to laminarize the flow. Additionally, we raise the possibility of beating this limit. By simultaneously applying wall transpiration and spanwise wall motion, we show that it may be possible to attain sustained sub-laminar energy expenditure in a controlled flow. We provide a necessary design criterion for net energy savings.

When life arose from prebiotic molecules 3.5 billion years ago, what came first? Informational molecules (RNA, DNA), functional ones (proteins), or something else? We argue here for a different logic: rather than seeking a molecule type, we seek a dynamical process. Biology required an ability to evolve before it could choose and optimise materials. We hypothesise that the evolution process was rooted in the peptide folding process. Modelling shows how short random peptides can collapse in water and catalyse the elongation of others, powering both increased folding stability and emergent autocatalysis through a disorder-to-order process.

The convergence of free-energy calculations based on importance sampling depends heavily on the choice of collective variables (CVs), which in principle, should include the slow degrees of freedom of the biological processes to be investigated. Autoencoders (AEs), as emerging data-driven dimension reduction tools, have been utilised for discovering CVs. AEs, however, are often treated as black boxes, and what AEs actually encode during training, and whether the latent variables from encoders are suitable as CVs for further free-energy calculations remains unknown. In this contribution, we review AEs and their time-series-based variants, including time-lagged AEs (TAEs) and modified TAEs, as well as the closely related model variational approach for Markov processes networks (VAMPnets). We then show through numerical examples that AEs learn the high-variance modes instead of the slow modes. In stark contrast, time series-based models are able to capture the slow modes. Moreover, both modified TAEs with extensions from slow feature analysis and the state-free reversible VAMPnets (SRVs) can yield orthogonal multidimensional CVs. As an illustration, we employ SRVs to discover the CVs of the isomerizations of N-acetyl-N′-methylalanylamide and trialanine by iterative learning with trajectories from biased simulations. Last, through numerical experiments with anisotropic diffusion, we investigate the potential relationship of time-series-based models and committor probabilities.

Respiratory disease and breathing abnormalities worsen with dehydration of the upper airways. We find that humidification of inhaled air occurs by evaporation of water over mucus lining the upper airways in such a way as to deliver an osmotic force on mucus, displacing it towards the epithelium. This displacement thins the periciliary layer of water beneath mucus while thickening topical water that is partially condensed from humid air on exhalation. With the rapid mouth breathing of dry air, this condensation layer, not previously reported while common to transpiring hydrogels in nature, can deliver an osmotic compressive force of up to around 100 cm H2O on underlying cilia, promoting adenosine triphosphate secretion and activating neural pathways. We derive expressions for the evolution of the thickness of the condensation layer, and its impact on cough frequency, inflammatory marker secretion, cilia beat frequency and respiratory droplet generation. We compare our predictions with human clinical data from multiple published sources and highlight the damaging impact of mouth breathing, dry, dirty air and high minute volume on upper airway function. We predict the hypertonic (or hypotonic) saline mass required to reduce (or amplify) dysfunction by restoration (or deterioration) of the structure of ciliated and condensation water layers in the upper airways and compare these predictions with published human clinical data. Preserving water balance in the upper airways appears critical in light of contemporary respiratory health challenges posed by the breathing of dirty and dry air.