Direct numerical simulations of the turbulence of a Herschel–Bulkley (HB) fluid in a rough channel are performed at a shear Reynolds number $Re_{\tau } \approx 300$ and a Bingham number ${Bn} \approx 0.9$. For the type of rough surface used in this study, the results indicate that Townsend's wall similarity hypothesis also holds for HB fluids. However, there are notable differences compared with the effect of roughness on Newtonian fluids. More specifically, the effect of roughness appears to be slightly stronger for HB fluids, in the sense that the bulk Reynolds number, based on the viscosity at the wall, is reduced further due to the increase in viscosity in the troughs of the roughness surface induced by the low shear. At the same time, for the simulated rough surface, the contribution of form drag to the total pressure drop is reduced from 1/4 to about 1/5 due to the persistence of viscous shear in the boundary layer, reducing its shielding effect. As for the friction factor, due to the nonlinearity of the HB constitutive relation, its use with the wall shear rate from the mean wall shear stress underpredicts the minimum viscosity at the wall by up to 18 %. This inevitably leads to uncertainties in the prediction of the friction factor. Finally, it is observed that the rough surface is unable to break the peculiar near-wall flow structure of HB fluids, which consists of long persistent low-speed streaks occupying the entire domain. This means that the small-scale energy is significantly reduced for HB fluids, even in rough channels, with the energy more concentrated in the lower wavenumber range, implying an increase in the slope of the power spectrum to $-7/2$ in the inertial range, as shown by Mitishita et al. (J. Non-Newtonian Fluid Mech., vol. 293, 2021, 104570).

This paper presents numerical simulations of the free fall of homogenous cylinders of length-to-diameter ratios $2$, $3$ and 5 and solid-to-fluid-density ratios $\rho _s/\rho$ going from 0 to 10 in transitional regimes. The path instabilities are shown to be due to two types of transitional states. The well-known fluttering state is a solid mode, characterised by significant oscillations of the cylinder axis due to a strong interaction between the vortex shedding in the wake and the solid degrees of freedom. Weakly oscillating, mostly irregular trajectories, are fluid modes, associated with purely fluid instabilities in the wake. The interplay of solid and fluid modes leads to a varying scenario in which the length-to-diameter and density ratios play an important role. The description is accompanied by the presentation of the identified transitional states in terms of path characteristics and vorticity structure of the wakes and by bifurcation diagrams showing the evolution of asymptotic states with increasing Galileo numbers. There appears to be a strong difference between the behaviour of cylinders of aspect ratio $L/d=3$ and 5. A similar contrast is stated between light cylinders of density ratios $\rho _s/\rho \le 2$ and dense cylinders of density ratios 5 and 10. Finally, the question of the scatter of values of the drag coefficient and of the frequency of oscillations raised in the literature is addressed. It is shown, that in addition to external parameters (Galileo number, density and aspect ratio) the amplitude of oscillations characterising the instability development is to be taken into account to explain this scatter. Fits of the simulation results to simple correlations are proposed. Namely that of the drag coefficient proves to be accurate (better than 1 % of accuracy) but also that of the Strouhal number (a few per cent of accuracy) may be of practical use.

This study presents an approach to investigate the role of eddy viscosity in linearized mean-field analysis of broadband turbulent flows. The procedure is based on spectral proper orthogonal decomposition (SPOD), resolvent analysis and the energy budget of coherent structures and is demonstrated using the example of a turbulent jet. The focus is on the coherent component of the Reynolds stresses, the nonlinear interaction term of the fluctuating velocity component in frequency space, which appears as an unknown in the derivation of the linearized Navier–Stokes equations and which is the quantity modelled by the Boussinesq approach. For the considered jet the coherent Reynolds stresses are found to have a mostly dissipative effect on the energy budget of the dominant coherent structures. Comparison of the energy budgets of SPOD and resolvent modes demonstrates that dissipation caused by nonlinear energy transfer must be explicitly considered within the linear operator to achieve satisfactory results with resolvent analysis. Non-modelled dissipation distorts the energy balance of the resolvent modes and is not, as often assumed, compensated for by the resolvent forcing vector. A comprehensive analysis, considering different predictive and data-driven eddy viscosities, demonstrates that the Boussinesq model is highly suitable for modelling the dissipation caused by nonlinear energy transfer for the considered flow. Suitable eddy viscosities are analysed with regard to their frequency, azimuthal wavenumber and spatial dependence. In conclusion, the energetic considerations reveal that the role of eddy viscosity is to ensure that the energy the structures receive from the mean-field is dissipated.

Energetics of mode-1 internal waves interacting with topographic ridges are investigated using high-resolution two-dimensional simulations at spatial scales of $O(100)$ m that span between classical laboratory-scale ($O(10)$ m) and field-scale simulations ($O(1000\unicode{x2013}10\,000)$ m). This paper focuses on the energetics of wave–topography interaction, with emphasis on systematically examining the partitioning of the incident wave energy as a function of wave forcing and topographic parameters. Partitioning of energy into the transmitted, reflected and dissipated components is quantified as a function of wave Froude number $Fr=U_0/c_{ph}$ ($U_0=$ velocity amplitude of forcing and $c_{ph}=$ internal wave celerity), slope criticality $=\gamma /s$, where $\gamma =$ topographic slope and $s=$ wave characteristic slope, and the ratio of topographic height $h_t$ to water depth $d$. As $Fr$ increases, dense fluid from the base of the stratified water column surges upslope with significant vertical inertia, leading to the formation of internal boluses that plunge over and onto the downstream side of the ridge, resulting in elevated dissipation. Results show that non-hydrostatic contributions to the total energy flux are significant (up to 50 %). Analysis of the energy flux budget shows that transmitted energy flux decreases monotonically as $\gamma /s$ increases for any given $Fr$ and $h_t/d$. At critical slopes ($\gamma /s=1$), the transmitted energy flux scales as a linear function of $h_t/d$, with a mild dependence on $Fr$, a key result that can be useful in energy flux parameterizations. Reflected energy flux exhibits a nonlinear dependence on the ridge height, increasing sharply when $h_t/d > 0.5$. Dissipation is enhanced at critical slopes, with a plateau evident for $\gamma /s \ge 1$ and $h_t/d = 0.5$ for all $Fr$.

Turbulent flow induced by elastorotational instability in viscoelastic Taylor–Couette flow (TCF) with Keplerian rotation is analogous to a turbulent accretion disk destabilized by magnetorotational instability. We examine this novel viscoelastic Keplerian turbulence via direct numerical simulations (DNS) for the shear Reynolds number ($Re$) ranging from $10^2$ to $10^4$. The observed characteristic flow structure consists of penetrating streamwise vortices with axial length scales much smaller than the gap width, distinct from the classic centrifugally induced Taylor vortices, which have axial lengths of the gap width. These intriguing vortices persist for the wide $Re$ range considered and give rise to intriguing scaling behaviour in key flow quantities. Specifically, the characteristic axial length of the penetrating vortices is shown to scale as $Re^{-0.22}$; the angular momentum transport scales as $Re^{0.42}$; the kinetic and elastic boundary-layer thicknesses based on angular velocity and hoop stress near the inner cylinder wall scale as $Re^{-0.48}$ and $Re^{-0.49}$, respectively. This implies that the viscoelastic Keplerian turbulence belongs to the classical turbulent regime of TCF with the Prandtl–Blasius-type boundary layer. Furthermore, we present an analytical relation between the viscous and elastic dissipation rates of kinetic energy and the angular momentum transport and in turn demonstrate its validity using our DNS data. This study has paved the way for future research to explore astrophysics-related Keplerian turbulence and angular momentum transport via the scaling relations of the analogous TCF of dilute polymeric solutions.

Depth-limited overturning wave shape affects water turbulence and sediment suspension. Experiments have shown that wind affects shoaling and overturning wave shape, with uncertain mechanism. Here, we study wind effects (given by the wind Reynolds number) on solitary wave shoaling and overturning with the two-phase direct numerical simulations model Basilisk run in two dimensions on steep bathymetry for fixed wave Reynolds number and Bond number. For all wind, the propagating solitary wave sheds a two-dimensional turbulent air wake and has nearly uniform speed with minimal wave energy changes over the rapidly varying bathymetry. Wave-face slope is influenced by wind, and shoaling wave shape changes are consistent with previous studies. As overturning jet impacts, wind-dependent differences in overturn shape are quantified. The non-dimensional breakpoint location and overturn area have similar wind dependence as previous experience, whereas the overturn aspect ratio has opposite wind dependence. During shoaling, the surface viscous stresses are negligible relative to pressure. Surface tension effects are also small but grow rapidly near overturning. In a wave frame of reference, surface pressure is low in the lee and contributes 2–5 % to the velocity potential rate of change in the surface dynamic boundary condition, which, integrated over time, changes the wave shape. Reasons why the overturn aspect ratio is different than in experiment and why a stronger simulated wind is required are explored. The dramatic wind effects on overturning jet area, and thus to the available overturn potential energy, make concrete the implications of wind-induced changes to wave shape.

Rotor-stator flows are known to exhibit instabilities in the form of circular and spiral rolls. While the spirals are known to emanate from a supercritical Hopf bifurcation, the origin of the circular rolls is still unclear. In the present work we suggest a quantitative scenario for the circular rolls as a response of the system to external forcing. We consider two types of axisymmetric forcing: bulk forcing (based on the resolvent analysis) and boundary forcing using direct numerical simulation. Using the singular value decomposition of the resolvent operator the optimal response is shown to take the form of circular rolls. The linear gain curve shows strong amplification at non-zero frequencies following a pseudo-resonance mechanism. The optimal energy gain is found to scale exponentially with the Reynolds number $Re$ (for $Re$ based on the rotation rate and interdisc spacing $H$). The results for both types of forcing are compared with former experimental works and previous numerical studies. Our findings suggest that the circular rolls observed experimentally are the effect of the high forcing gain together with the roll-like form of the leading response of the linearised operator. For high enough Reynolds number it is possible to delineate between linear and nonlinear responses. For sufficiently strong forcing amplitudes, the nonlinear response is consistent with the self-sustained states found recently for the unforced problem. The onset of such non-trivial dynamics is shown to correspond in state space to a deterministic leaky attractor, as in other subcritical wall-bounded shear flows.

The stability of liquid-film flows is essential in many industrial applications. In the dip-coating process, a liquid film forms over a substrate extracted at a constant speed from a bath. We studied the linear stability of this film considering different thicknesses $\hat {h}$ for four liquids, spanning an extensive range of Kapitza numbers ($Ka$). By solving the Orr–Sommerfeld eigenvalue problem with the Chebyshev–Tau spectral method, we calculated the threshold between growing and decaying perturbations, investigated the instability mechanism, and computed the absolute/convective threshold. The instability mechanism was studied by analysing the perturbations’ vorticity distribution and the kinetic energy balance. It was found that liquids with low $Ka$ (e.g. corn oil, $Ka = 4$) are stable for a smaller range of wavenumbers compared with liquid with high $Ka$ (e.g. liquid zinc, $Ka = 11\,525$). Surface tension has a stabilising and a destabilising effect. For long waves, it curves the vorticity lines near the substrate, reducing the flow under the crests. For short waves, it fosters vorticity production at the interface and creates a region of intense vorticity near the substrate. In addition, we discovered that the surface tension contributes to both the production and dissipation of perturbation's energy depending on the $Ka$ number. Regarding the absolute/convective threshold, we identified a window in the parameter space where unstable waves propagate throughout the entire domain (indicating absolute instability). Perturbations affecting Derjaguin's solution ($\hat {h}=1$) for $Ka<17$ and the Landau–Levich–Derjaguin solution ($\hat {h}=0.945 Re^{1/9}Ka^{-1/6}$), are advected by the flow (indicating convective instability).

By comparing the budget of a data-driven quasi-linear approximation (DQLA) (Holford, Lee & Hwang, J. Fluid Mech., vol. 980, 2024, A12) and direct numerical simulation (DNS) (Lee & Moser, J. Fluid Mech., vol. 860, 2019, pp. 886–938), the energetics of linear models for wall-bounded turbulence are assessed. The DQLA is implemented with the linearised Navier–Stokes equations with a stochastic forcing term and an eddy viscosity diffusion model. The self-consistent nature of the DQLA allows for a global comparison across all wavenumbers to assess the role of the various terms in the linear model in replicating the features present in DNS. Starting from the steady-state second-order statistics of a Fourier mode, a spectral budget equation is derived, connecting Lyapunov-like equations to the transport budget equations obtained from DNS. It is found that the DQLA and DNS are in good qualitative agreement for the streamwise-elongated structures present in DNS, comparing well for production, viscous transport and wall-normal turbulent transport. However, the DQLA does not have an energy-conservative nonlinear term. This results in no dissipation under molecular viscosity, with energy instead being dissipated locally through the eddy viscosity model, which models the energy removal by the nonlinear term at integral length scales. Comparison of the pressure–strain statistics also highlights the absence of the streak instability, with production and forcing mainly being retained in the streamwise and wall-normal components or shifted to the spanwise component. It is demonstrated that the eddy viscosity diffusion term locally enforces a self-similar budget, making the model for the nonlinear term self-consistent with a logarithmic mean profile. Implications and recommendations to improve the current eddy viscosity enhanced linear models are also discussed concerning the comparison with DNS, as well as considerations with regard to pressure statistics to mimic the role of the streak instability through colour of turbulence models.

This DNS study considers transition to turbulence in plane Couette flow (pCf) with a rough stationary wall and a smooth moving wall. The roughness elements are square ribs of height $k=0.2h$ (where $h$ is the half-channel height). Two different pitch separations, $\lambda =2k$ and $10k$, are considered, i.e. d-type and k-type roughness, respectively. The transition in both rough pCf cases takes place through a stage of alternate laminar–turbulent bands aligned in an oblique fashion. However, roughness causes a shift in the transitional Reynolds number ($Re$) range. In the k-type roughness, stable bands are observed in the range $Re \in [300, 325]$, which is a downwards shift from the transitional $Re$ range for the smooth pCf ($Re \in [325,400]$). The d-type roughness, on the other hand, surprisingly shifts the transitional $Re$ range upwards to $Re \in [350,425]$. This peculiar behaviour is associated with the ability of the ribs to shed and regenerate vorticity. Large-scale components extracted using a filtering process relate to the transition bands and flow parallel to the oblique laminar–turbulent boundaries. Counter-rotating vortices are present in the turbulent regions of the flow field, which exist in tandem with the high- and low-velocity streaks. Another interesting observation is the secondary Reynolds shear stresses, $-\overline {v^{\prime }u^{\prime }}$ and $-\overline {w^{\prime }v^{\prime }}$, which are non-zero in the transitional regime, in contrast to the turbulent regime where they are negligible.

We conduct a numerical study on the drag-reduction mechanism of an opposition- controlled turbulent channel flow from the viewpoint of a symbolic dynamics approach. The effect of the virtual wall formed by the opposition control is maximised at the location of the detection plane $y_d^+ \approx 10$. At this wall-normal location, the local link strength of the self-loop of network nodes representing the negative correlation pattern between the streamwise and wall-normal velocity fluctuations is maximised in the uncontrolled flow. In the controlled case, the multiscale complexity–entropy causality plane and the spatial permutation entropy at $y_d^+ \approx 10$ indicate that the drag-reduction effect is attributed to the reduction of the region where streaks actively coalesce and separate and the suppression of the regeneration cycle in the region near the wall.

Extensive three-dimensional boundary-integral simulations are presented for the steady-state, low-Reynolds-number motion of a non-wetting deformable drop in another liquid on an inclined solid wall. The drop remains separated from the wall by a lubricating film. The boundary-integral formulation is based on the half-space Green function. The focus is on the challenging case of small tilt angles $\theta$ combined with high drop-to-medium viscosity ratios $\lambda$, when the drop travels with strong hydrodynamical interaction very close to the wall. Simulations to steady state have required ultrahigh drop surface resolutions (to $3\times 10^5$ boundary elements) achieved through multipole acceleration and combined with novel regularization to fully eliminate near-singular behaviour of the double-layer integrals due to small clearance. Non-dimensional drop speed $U$ is presented for $\theta \geq 7.5^\circ$, $\lambda \leq 300$ and in a broad range of Bond numbers $B$, covering from nearly spherical to strongly pancaked drops. The results are consistent with published experiments on liquid–liquid systems. At small $\theta$ and $\lambda \gg 1$, $U$ is a strong, decreasing function of $B$; the asymptotic regime $U\to 0$ at $B\to 0$ is not observed in the simulated range. For small $B$, the dimpled thin-film geometry is insensitive to $\lambda =1\unicode{x2013}300$. For pancaked drops, the lubrication film is much thicker for $\lambda =1$ than for $\lambda \gg 1$ drops. Approximate thin-film uniformity in the drop motion direction is confirmed for pancaked, but not for $B\ll 1$, drops. Kinematics of drop motion shows that neither perfect tank treading, nor perfect rolling can be approached for liquid–liquid systems in the purely hydrodynamical formulation. The methodology is applicable to other problems and can allow for direct inclusion of short-range colloidal forces in three-dimensional boundary-integral simulations.

We investigate the coupling effects of the two-phase interface, viscosity ratio and density ratio of the dispersed phase to the continuous phase on the flow statistics in two-phase Taylor–Couette turbulence at a system Reynolds number of $6\times 10^3$ and a system Weber number of 10 using interface-resolved three-dimensional direct numerical simulations with the volume-of-fluid method. Our study focuses on four different scenarios: neutral droplets, low-viscosity droplets, light droplets and low-viscosity light droplets. We find that neutral droplets and low-viscosity droplets primarily contribute to drag enhancement through the two-phase interface, whereas light droplets reduce the system's drag by explicitly reducing Reynolds stress due to the density dependence of Reynolds stress. In addition, low-viscosity light droplets contribute to greater drag reduction by further reducing momentum transport near the inner cylinder and implicitly reducing Reynolds stress. While interfacial tension enhances turbulent kinetic energy (TKE) transport, drag enhancement is not strongly correlated with TKE transport for both neutral droplets and low-viscosity droplets. Light droplets primarily reduce the production term by diminishing Reynolds stress, whereas the density contrast between the phases boosts TKE transport near the inner wall. Therefore, the reduction in the dissipation rate is predominantly attributed to decreased turbulence production, causing drag reduction. For low-viscosity light droplets, the production term diminishes further, primarily due to their greater reduction in Reynolds stress, while reduced viscosity weakens the density difference's contribution to TKE transport near the inner cylinder, resulting in a more pronounced reduction in the dissipation rate and consequently stronger drag reduction. Our findings provide new insights into the physics of turbulence modulation by the dispersed phase in two-phase turbulence systems.

Direct numerical simulation (DNS) is performed to explore turbulent Rayleigh–Bénard convection in spherical shells. Our simulations cover six distinct values of radius ratio, $\eta = r_i/r_o = 0.2$, 0.3, 0.4, 0.5, 0.6 and 0.8, under the assumption of a centrally condensed mass with gravity profile $g \sim 1/r^{2}$; where $r_i$, $r_o$ and $r$ denote the inner shell radius, the outer shell radius and the local radial coordinate, respectively. The Prandtl number is kept constant at unity while the Rayleigh number ($Ra$) is varied from $3 \times 10^{3}$ to $5 \times 10^8$. Our primary aim is to analyze how the radius ratio influences the global transport properties and flow physics. To gain insights into the scaling behaviour of the Nusselt number ($Nu$) and the Reynolds number ($Re$) with respect to $Ra$ and $\eta$, we apply the Grossmann–Lohse (GL) theory (Grossmann & Lohse, J. Fluid Mech., vol. 407, 2000, pp. 27–56) to the system. It is observed that the scaling exponents for $Nu$ and $Re$ in relation to $Ra$ are more significant for smaller $\eta$ values, suggesting that the simulations with smaller $\eta$ reach the classical $Nu\sim Ra^{1/3}$ regime at a relatively lower $Ra$. This observation could also imply the systems with smaller $\eta$ might transition to the ultimate regime earlier at a smaller $Ra$. Based on our extensive DNS data, we establish that the thickness of the inner thermal boundary, $\lambda _{\vartheta }^{i}$, follows a scaling relationship of $\lambda _{\vartheta }^{i} \sim \eta ^{1/2}$. This relationship, in turn, leads to a scaling law for $Nu$ in the form of $Nu \sim f(\eta ) Ra^{\gamma }$, where the function $f(\eta )$ is defined as $f(\eta ) = {\eta ^{1/2}}/{(1+\eta ^{4/3})}$, and the exponent $\gamma$ depends on both $Ra$ and $\eta$. Additionally, we characterize and explain the asymmetry in the velocity field by introducing the separate Reynolds numbers for the inner and outer shells. The asymmetry of the kinetic and thermal energy dissipation rates in the inner and outer boundary layers (BLs) is also quantified.

The mechanical behaviour of wet particle aggregates is crucial in many granular processes such as wet granulation and soil degradation. However, the interplay of capillary and viscous forces for aggregate stability and breakage have remained elusive due to the complexity of granular dynamics. We use particle dynamics simulations to analyse the deformation and breakage of wet aggregates colliding with a flat wall. The aggregates are composed of spherical particles and the effect of liquid bonds is modelled through capillary and lubrication forces acting between particles. We perform an extensive parametric study by varying surface tension, impact velocity and liquid viscosity in a broad range of values. We show that when lubrication force is neglected, aggregate breakage is fully controlled by the reduced kinetic energy $\xi$, defined as the ratio of incident kinetic energy to the initial capillary energy. At low values of $\xi$, the aggregate deforms without breakage due to inelastic energy loss induced by rearrangements and loss of capillary bonds, whereas above a critical value of $\xi$ it breaks into smaller aggregates due to the transfer of kinetic energy from aggregate to fragments. In the presence of lubrication forces, the crossover from capillary to viscous regime is controlled by the capillary number, defined as the ratio of viscous dissipation to capillary energy. We find that the critical value of $\xi$ for aggregate breakage in the viscous regime increases as a power law with capillary number while the effective restitution coefficient follows the same trend as in the capillary regime.

We study the effect of geometrical confinement on thermal convection by laboratory experiments and direct numerical simulations using Hele-Shaw geometries (typically the gap-to-height aspect ratio $0.12$) for the Prandtl number $Pr \geq 40$ and the Rayleigh number $Ra \leq 6 \times 10^7$. Under such strong unidirectional confinement, the convective flows are forced to squeeze within the narrow gap and exhibit unique spatiotemporal signatures, which contrast those in unconfined systems. With the increase of $Ra$, we identify that the system experiences five convective regimes that can be classified from two aspects, time dependency and flow dimensionality: (I) quasi-two-dimensional (quasi-2-D) steady flow; (II) quasi-2-D flow with oscillatory corner rolls; (III) three-dimensional (3-D) flow with oscillatory corner rolls; (IV) 3-D steady flow; and (V) 3-D time-dependent motion of plumes around sidewalls. Notably, unsteadiness does not emerge globally, but is localised near the sidewalls as oscillatory corner rolls, resulting in the regime transitions happening in a quasi-steady manner. We confirm that these regime transitions show less dependence on both $Pr$ and the other (wider) horizontal scale of the geometry. Moreover, we find that a recently proposed criterion ‘degree of confinement’ (Noto et al., Proc. Natl Acad. Sci. USA, vol. 121, issue 28, 2024, e2403699121) successfully explains the emergence of 3-D structures, expanding its applicable range to smaller $Ra$. This study deepens the comprehension of the thermal convection emerging in tight geometries, impacting across disciplines, such as Earth and planetary science, and thermal engineering.

Heavy particles suspended in turbulent flow possess inertia and are ejected from violent vortical structures by centrifugal forces. Once piled up along particle paths, this small-scale mechanism leads to an effective large-scale drift. This phenomenon, known as ‘turbophoresis’, causes particles to leave highly turbulent regions and migrate towards calmer regions, explaining why particles transported by non-homogeneous flows tend to concentrate near the minima of turbulent kinetic energy. It is demonstrated here that turbophoretic effects are just as crucial in statistically homogeneous flows. Although the average turbulent activity is uniform, instantaneous spatial fluctuations are responsible for inertial-range inhomogeneities in the particle distribution. Direct numerical simulations are used to probe particle accelerations, specifically how they correlate to local turbulent activity, yielding an effective coarse-grained dynamics that accounts for particle detachment from the fluid and ejection from excited regions through a space- and time-dependent non-Fickian diffusion. This leads to cast fluctuations in particle distributions in terms of a scale-dependent Péclet number ${\textit {Pe}}_\ell$, which measures the importance of turbulent advection compared with inertial turbophoresis at a given scale $\ell$. Multifractal statistics of energy dissipation indicate that $ {\textit {Pe}}_\ell \sim \ell ^\delta /\tau _{p}$ with $\delta \approx 0.84$. Numerical simulations support this behaviour and emphasise the relevance of the turbophoretic Péclet number in characterising how particle distributions, including their radial distribution function, depends on $\ell$. This approach also explains the presence of voids with inertial-range sizes, and the fact that their volumes have a non-trivial distribution with a power-law tail $p(\mathcal {V}) \propto \mathcal {V}^{-\alpha }$, with an exponent $\alpha$ that tends to 2 as ${\textit {Pe}}_\ell \to 0$.

We present herein the derivation of a lubrication-mediated (LM) quasi-potential model for droplet rebounds off deep liquid baths, assuming the presence of a persistent dynamic air layer which acts as a lubricating pressure transfer. We then present numerical simulations of the LM model for axisymmetric rebounds of solid spheres and compare quantitatively to current results in the literature, including experimental data in the low-speed impact regime. In this regime the LM model has the advantage of being far more computationally tractable than direct numerical simulation (DNS) and is also able to provide detailed behaviour within the micro-metric thin lubrication region. The LM system has an interesting mathematical structure, with the lubrication layer providing a free-boundary elliptic problem mediating the drop and bath free-boundary evolutionary equations.

We consider the steady flow of a viscoelastic film over an inclined plane featuring periodic trenches normal to the main flow direction. The trenches have a square cross-section and side length 5–8 times the capillary length. Owing to the orientation of the substrate, the film fails to coat the topographical feature entirely, forming a second gas–liquid interface inside the trench and two three-phase contact lines at the points where the free surface meets the wall of the trench. The volume of entrapped air depends on material and flow parameters and geometric conditions. We develop a computational model and carry out detailed numerical simulations based on the finite element method to investigate this flow. We solve the two-dimensional mass and momentum conservation equations using the exponential Phan-Thien & Tanner constitutive model to account for the rheology of the viscoelastic material. Due to the strong nonlinearity, multiple steady solutions possibly connected by turning points forming hysteresis loops, transcritical bifurcations and isolated solution branches are revealed by pseudo-arc-length continuation. We perform a thorough parametric analysis to investigate the combined effects of elasticity, inertia, capillarity and viscosity on the characteristics of each steady flow configuration. The results of our simulations indicate that fluid inertia and elasticity may or may not promote wetting, while shear thinning and hydrophilicity always promote the wetting of the substrate. Interestingly, there are conditions under which the transition to the inertial regime is not smooth, but a hysteresis loop arises, signifying an abrupt change in the film hydrodynamics. Additionally, we investigate the effect of the geometrical characteristics of the substrate, and our results indicate that there is a unique combination of the geometry and viscoelastic properties that either maximizes or minimizes the wetting lengths.