In this work, smoothed particle hydrodynamics (SPH) is employed to investigate the segregation evolution in granular flows. We first provide the Lagrangian description-based governing equations, including the linear momentum conservation and the segregation–diffusion equation. Then the hybrid continuum surface reaction scheme is introduced to formulate the concentration-related inhomogeneous Neumann boundary condition on the free and wall surfaces. We follow a two-stage strategy to advance boundary particle searching and normal direction identification. Moreover, $C^1$ consistency is considered based on the Taylor series to obtain accurate segregation flux gradient along the boundary. Our SPH model is validated with a shear box experiment. The model is then applied to investigate the segregation mechanism in bidisperse-sized granular flows in a rotating drum.

Growth of a fluid-infused patch on a thin porous layer, e.g. on a piece of paper or cloth, is related to the transmission of virus particles through exhaled droplets and aerosols. We present a theoretical model to describe how a wet patch develops gradually through imbibition, once a sessile droplet attaches at a permeable surface and drains gradually into a thin porous layer. Two limiting cases are considered based on different assumptions on the motion of the contact line during the coupled process of drop drainage and patch growth: (i) the apparent contact angle remains unchanged, so the radius of a sessile droplet decreases with time; and (ii) the location of the contact line remains pinned, so the contact angle decreases as time progresses. The model leads to evolution pathways for both the droplet and the fluid film within the porous layer, without introducing arbitrary fitting parameters. Potential implications of the model and its solutions are also discussed briefly in the context of the outspread of COVID-19, employing physical parameters for exhaled droplets, paper and cloth.

For flows in microchannels, a slip on the walls may be efficient in reducing viscous dissipation. A related issue, addressed in this article, is to decrease the effective viscosity of a dilute monodisperse suspension of spheres in Poiseuille flow by using two parallel slip walls. Extending the approach developed for no-slip walls in Feuillebois et al. (J. Fluid Mech., vol. 800, 2016, pp. 111–139), a formal expression is obtained for the suspension intrinsic viscosity $[\mu ]$ solely in terms of a stresslet component and a quadrupole component exerted on a single freely suspended sphere. In the calculation of $[\mu ]$, the hydrodynamic interactions between a sphere and the slip walls are approximated using either the nearest wall model or the wall-superposition model. Both the stresslet and quadrupole are derived and accurately calculated using bipolar coordinates. Results are presented for $[\mu ]$ in terms of $H/(2a)$ and $\tilde{\lambda}=\lambda /a\leq 1$, where $H$ is the gap between walls, $a$ is the sphere radius and $\lambda$ is the wall slip length using the Navier slip boundary condition. As compared with the no-slip case, the intrinsic viscosity strongly depends on $\tilde{\lambda}$ for given $H/(2a)$, especially for small $H/(2a)$. For example, in the very confined case $H/(2a)=2$ (a lower bound found for practical validity of single-wall models) and for $\tilde{\lambda}=1$, the intrinsic viscosity is three times smaller than for a suspension bounded by no-slip walls and five times smaller than for an unbounded suspension (Einstein, Ann. Phys., vol. 19, 1906, pp. 289–306). We also provide a handy formula fitting our results for $[\mu ]$ in the entire range $2\leq H/(2a)\leq 100$ and $\tilde{\lambda}\leq 1$.

The column collapse experiment is a simplified version of natural and industrial granular flows. In this set-up, a column built with grains collapses and spreads over a horizontal plane. Granular flows are often studied with a monodisperse distribution; however, this is not the case in natural granular flows where a variety of grain sizes, known as polydispersity, is a common feature. In this work, we study the effect of polydispersity, and of the inherent changes that polydispersity causes in the initial packing fraction, in dry and immersed columns. We show that dry columns are not significantly affected by polydispersity, reaching similar distances at similar times. In contrast, immersed columns are strongly affected by the polydispersity and packing fraction, and the collapse sequence is linked to changes of the basal pore fluid pressure $P$. At the collapse initiation, negative changes of $P$ beneath the column produce a temporary increase of the column strength. The negative change of $P$ lasts longer in polydisperse columns than in monodisperse columns, delaying the collapse sequence. Conversely, during the column spreading, positive changes of $P$ lead to a decrease of the shear strength. For polydisperse collapses, the excess of $P$ lasts longer, allowing the material to reach farther distances, compared with the collapses of monodisperse materials. Finally, we show that a mobility model that scales the final runout with the collapse kinetic energy remains true for different polydispersity levels in a three-dimensional configuration, capturing the scaling between the micro to macro controlling features.

A relation among invariants of filtered velocity gradients with two different filter sizes is derived. Based on this relation and physical reasoning, it is shown analytically that strain self-amplification contributes more to energy transfer than vortex stretching in homogeneous turbulence, as observed in recent numerical investigations of homogeneous isotropic turbulence. We note that the invariant relation studied and hence the inequality between strain self-amplification and vortex stretching apply to all homogeneous flows, not restricted to isotropic turbulence.

When a liquid jet plunges into a pool, it can generate a bubble-laden jet flow underneath the surface. This common and simple phenomenon is investigated experimentally for circular jets to illustrate and quantify the role played by the net gas/liquid void fraction on the maximum bubble penetration depth. It is first shown that an increase in either the impact diameter or the jet fall height to diameter ratio at constant impact momentum leads to a reduction in the bubble cloud size. By measuring systematically the local void fraction using optical probes in the biphasic jet, it is then demonstrated that this effect is a direct consequence of the increase in air content within the cloud. A simple momentum balance model, including only inertia and the buoyancy force, is shown to predict the bubble cloud depth without any fitting parameters. Finally, a Froude number based on the bubble terminal velocity, the cloud depth and also the net void fraction is introduced to propose a simple criterion for the threshold between the inertia-dominated and buoyancy-dominated regimes.

This study focuses on numerically investigating the impact of fluid viscoelasticity on the flow dynamics around a transversely forced oscillating cylinder operating in the laminar vortex shedding regime at a fixed Reynolds number of $Re = 100$. Specifically, we explore how fluid viscoelasticity affects the boundary between the lock-in and no lock-in regions and the corresponding wake characteristics compared with a simple Newtonian fluid. Our findings reveal that fluid viscoelasticity enables the synchronization of the vortex street with the cylinder motion at lower oscillation frequencies than those required for a Newtonian fluid. Consequently, the lock-in region boundary for a viscoelastic fluid differs from that of a Newtonian fluid and expands in the non-dimensional cylinder oscillation amplitude and frequency parameter space. In the primary synchronization region, the wake of a Newtonian fluid exhibits ‘2S’ (two single vortices) and ‘P+S’ (a pair of vortices and a single vortex) shedding modes. In contrast, a ‘2P’ (two pairs of vortices) vortex mode is observed for a viscoelastic fluid within the same region. To gain a deeper understanding of the differences in the coherent flow structures and their associated frequencies between the two fluids, we employ the data-driven reduced-order modelling technique, known as the dynamic mode decomposition (DMD) technique. Utilizing this technique, we successfully extract and visualize the two competing fundamental frequencies (cylinder oscillation and natural vortex shedding frequencies) and their associated flow structures in the case of the no lock-in state, whereas only the dominant cylinder oscillation frequency and associated flow structure in the case of the lock-in state. Furthermore, we propose that the presence of excess strain resulting from the stretching of polymer molecules in viscoelastic fluids leads to a distinct difference in the wake structure compared with Newtonian fluids. This observation aligns with the findings obtained from the $Q$-criterion and vorticity transport analysis of the wake.

We study the capillary rise of viscous liquids into sharp corners formed by two surfaces whose geometry is described by power laws $h_i(x) = c_i x^n$, $i = 1,2$, where $c_2 > c_1$ for $n \geq 1$. Prior investigations of capillary rise in sharp corners have shown that the meniscus altitude increases with time as $t^{1/3}$, a result that is universal, i.e. applies to all corner geometries. The universality of the phenomenon of capillary rise in sharp corners is revisited in this work through the analysis of a partial differential equation for the evolution of a liquid column rising into power-law-shaped corners, which is derived using lubrication theory. Despite the lack of geometric similarity of the liquid column cross-section for $n>1$, there exist a scaling and a similarity transformation that are independent of $c_i$ and $n$, which gives rise to the universal $t^{1/3}$ power law for capillary rise. However, the prefactor, which corresponds to the tip altitude of the self-similar solution, is a function of $n$, and it is shown to be bounded and monotonically decreasing as $n\to \infty$. Accordingly, the profile of the interface radius as a function of altitude is also independent of $c_i$ and exhibits slight variations with $n$. Theoretical results are compared against experimental measurements of the time evolution of the tip altitude and of profiles of the interface radius as a function of altitude.

The present research will provide an analytical explanation to experiments destabilising turbulence in pipe flow reported in Kuehnen et al. (Nat. Phys., vol. 14, 2018, 386–390). Those experiments show four methods by which turbulence vanishes from steady-state pipe flow, without decreasing its bulk velocity, until it becomes completely laminar. The explanation is based on our theory of underlying laminar flow (TULF), which has already been successfully applied to account for other uncommon experiments reported in the literature. The TULF is founded on the Reynolds-averaged Navier–Stokes equations and thus is a theory of ensemble-averaged flows. The zero theorem for steady-state flow is introduced as a universal result that will help explain the laminarisation process described in experiments. After presenting the most comprehensive solution for the mean pipe flow governing equation that, to our knowledge, has ever been reported, we uncover a general sequence for laminarisation, called the laminarisation pattern, and we introduce a mathematical model for it. We show that a drastic decrease in a flow's mean-pressure gradient, while maintaining constant its Reynolds number, is necessary and sufficient to erase turbulence. Equations derived from our model are used to calculate the minimum pressure gradient necessary to cause complete laminarisation in each experiment. Results are then contrasted with reported experimental data, with noticeable agreement. We also propose a figure of merit to assess the efficiency of each laminarisation method. Having disclosed the intrinsic mechanism leading to complete laminarisation, we expect researchers will propose other ingenious methods to achieve it.

Recent progress in the direct measurement of turbulent dissipation in the Arctic Ocean has highlighted the need for an improved parametrization of the turbulent diapycnal diffusivities of heat and salt that is suitable for application in the turbulent environment characteristic of this polar region. In support of this goal we describe herein a series of direct numerical simulations of the turbulence generated in the process of growth and breaking of Kelvin–Helmholtz billows. These simulations provide the data sets needed to serve as basis for a study of the stratified turbulent mixing processes that are expected to exist in the Arctic Ocean environment. The mixing properties of the turbulence are studied using a previously formulated procedure in which the temperature and salinity fields are sorted separately in order to enable the separation of irreversible Arctic mixing from reversible stirring processes and thus the definition of turbulent diffusivities for both heat and salt that depend solely upon irreversible mixing. These analyses allow us to demonstrate that the irreversible diapycnal diffusivities for heat and salt are both solely dependent on the buoyancy Reynolds number in the Arctic Ocean environment. These are found to be in close agreement with the functional forms inferred for these turbulent diffusivities in the previous work of Bouffard & Boegman (Dyn. Atmos. Oceans, vol. 61, 2013, pp. 14–34). Based on a detailed comparison of our simulation data with this previous empirical work, we propose an algorithm that can be used for inferring the diapycnal diffusivities from turbulent dissipation measurements in the Arctic Ocean.

A solid body in a viscous fluid undergoing oscillatory motion naturally produces a steady secondary flow due to convective inertia. This phenomenon is embodied in the streaming flow generated by a sphere in an unbounded fluid executing rectilinear oscillations. We review the considerable literature on this canonical problem and summarise exact and asymptotic formulas in the small-amplitude limit. These analytical formulas are used to explore the characteristic flow structure of this problem and clarify previously unreported features. A single, toroidal-shaped vortex exists in each hemisphere regardless of the oscillation frequency, which can drive a counter-flow away from the sphere. The vortex centre moves monotonically away from the sphere with decreasing oscillation frequency, and engulfs the entire flow domain for $\beta \equiv \omega R^2/\nu < 16.317$, where $\omega$ is the angular oscillation frequency, $R$ the sphere radius, and $\nu$ the fluid kinematic viscosity. This seemingly abrupt change in flow structure at the critical frequency $\beta _{critical} =16.317$, and its quantification, appear to have not been reported previously. We perform a direct numerical simulation of the Navier–Stokes equations, to (1) confirm existence of this critical frequency at finite amplitude, and (2) examine its variation with amplitude. This reveals a universal relationship between the critical frequency and oscillation amplitude, clarifying previous reports on the structure of this streaming flow. The critical frequency is shown to be identical for the streaming flow and the cycle-averaged particle paths, establishing that the critical frequency is accessible directly using standard measurements.

Understanding the mechanism of thermal instability in nanofluids is of fundamental importance to explore the reasons behind the enhancement of heat transfer efficiency. Since Buongiorno (ASME J. Heat Transfer, vol. 128, 2006, pp. 240–250) proposed his theoretical model of nanofluids, most studies focusing on the thermal instability analysis exclusively considered Brownian motion and thermophoresis as the main diffusion mechanisms of nanoparticles. All the analyses concluded that a nanofluid layer is much more unstable than its pure counterpart as it is heated from below. However, a recent experimental observation on Rayleigh–Bénard convection appears to contradict the theoretical prediction, implying that some mechanisms neglected in the previous model may have a significant impact on the onset of thermal convection. In the present study, we revise the convective transport model of nanofluids proposed by Buongiorno and find that the gravitational settling of nanoparticles is a crucial factor influencing the thermal instability behaviour of nanofluids. By performing a linear stability analysis based on the novel model, the effect of gravity settling exhibits a stabilizing mechanism to resist the destabilizing effect of thermophoresis. Furthermore, the onset of instability can be delayed once the nanoparticle diameter exceeds a certain threshold, which explains the phenomenon observed in experiments. Particularly, the oscillatory mode is found to emerge and dominate the flow instability when the gravity settling effect is competitive with the effect of thermophoresis.

For turbulent bubbly flows, multi-phase simulations resolving both the liquid and bubbles are prohibitively expensive in the context of different natural phenomena. One example is breaking waves, where bubbles strongly influence wave impact loads, acoustic emissions and atmospheric-ocean transfer, but detailed simulations in all but the simplest settings are infeasible. An alternative approach is to resolve only large scales, and model small-scale bubbles adopting sub-resolution closures. Here, we introduce a large eddy simulation smoothed particle hydrodynamics (SPH) scheme for simulations of bubbly flows. The continuous liquid phase is resolved with a semi-implicit isothermally compressible SPH framework. This is coupled with a discrete Lagrangian bubble model. Bubbles and liquid interact via exchanges of volume and momentum, through turbulent closures, bubble breakup and entrainment, and free-surface interaction models. By representing bubbles as individual particles, they can be tracked over their lifetimes, allowing closure models for sub-resolution fluctuations, bubble deformation, breakup and free-surface interaction in integral form, accounting for the finite time scales over which these events occur. We investigate two flows: bubble plumes and breaking waves, and find close quantitative agreement with published experimental and numerical data. In particular, for plunging breaking waves, our framework accurately predicts the Hinze scale, bubble size distribution, and growth rate of the entrained bubble population. This is the first coupling of an SPH framework with a discrete bubble model, with potential for cost-effective simulations of wave–structure interactions and more accurate predictions of wave impact loads.