In this paper, drop impact onto a sphere is numerically investigated at moderate Reynolds and Weber numbers. It is naturally expected that the aspect ratio of the sphere to the drop, $\unicode[STIX]{x1D706}_{r}$ , would make a big difference to drop spreading and retraction on the sphere, compared with drop impact onto a flat substrate. To quantitatively assess the effect of $\unicode[STIX]{x1D706}_{r}$ , a diffuse-interface immersed-boundary method is adopted after being validated against experiments. With the help of numerical simulations, we identify the key regimes in the spreading and retraction, analyse the results by scaling laws, and quantitatively evaluate the effect of $\unicode[STIX]{x1D706}_{r}$ on the impact dynamics. We find that the thickness of the liquid film spreading on the sphere can be well approximated by $h_{L,\infty }(1+3/4\unicode[STIX]{x1D706}_{r}^{-3/2})$ , where $h_{L,\infty }$ represents the film thickness of drop impact on a flat substrate. At the early stage of spreading, the temporal variation of the wetted area is independent of $\unicode[STIX]{x1D706}_{r}$ when the time is rescaled by the thickness of the liquid film. Drops are observed to retract on the sphere at a roughly constant speed, and the predictions of theoretical analysis are in good agreement with numerical results.

We numerically study the impact of a compound drop on a hydrophobic substrate using a ternary-fluid diffuse-interface method, aiming to understand how the presence of the inner droplet affects the spreading dynamics and maximal spreading of the compound drop. First, it is interesting to see that the numerical results for an impacting pure drop agree well with the universal rescaling of maximal spreading ratio proposed by Lee et al. (J. Fluid Mech., vol. 786, 2016, R4). Second, two flow regimes have been identified for an impacting compound drop: namely jammed spreading and joint rim formation. The maximal spreading ratio of the compound drop is found to depend on the volume fraction of the inner droplet $\unicode[STIX]{x1D6FC}$ , the surface tension ratio $\unicode[STIX]{x1D6FE}$ , the Weber number and the flow regime. Moreover, we propose a universal rescaling of maximal spreading ratio for compound drops, by integrating the one for pure drops with a corrected Weber number that takes $\unicode[STIX]{x1D6FC}$ , $\unicode[STIX]{x1D6FE}$ and the flow regime into account. The predictions of the universal rescaling are in good agreement with the numerical results for impacting compound drops.

We investigate the cavity formation during the impact of spheres and cylinders into a liquid pool by using a combination of experiments, simulations and theoretical analysis, with particular interest in contact-line pinning and its relation with the subsequent cavity evolution. The flows are simulated by a Navier–Stokes diffuse-interface solver that allows for moving contact lines. On the basis of agreement on experimentally measured quantities such as the position of the pinned contact line and the interface shape, we investigate flow details that are not accessible experimentally, identify the interface regions in the cavity formation and examine the geometric effects of impact objects. We connect wettability, inertia, geometry of the impact object, interface bending and contact-line position with the contact-line pinning by analysing the force balance at a pinned meniscus, and the result compares favourably with those from simulations and experiments. In addition to adjusting the interface bending, the object geometry also has a significant effect on the magnitude of low pressure in the liquid and the occurrence of flow separation. As a result, it is easier for an object with sharp edges to generate a cavity than a smooth object. A theoretical model based on the Rayleigh–Besant equation is developed to provide a quantitative description of the radial expansion of the cavity after the pinning of the contact line. The accuracy of the solution is greatly affected by the geometrical information on the interface connected to the pinned meniscus, showing the dependence of the global cavity dynamics on the local flows around the pinned contact line. Vertical ripple propagation on the cavity wall is found to follow the dispersion relation for the perturbation evolution on a hollow jet.

We numerically investigate the mechanism leading to the entrapment of spheres at the gas–liquid interface after impact. Upon impact onto a liquid pool, a hydrophobic sphere is seen to follow one of the three regimes identified in the experiment (Lee & Kim, Langmuir, vol. 24, 2008, pp. 142–145): sinking, bouncing or being entrapped at the interface. It is important to understand the role of wettability in this process of flow–structure interaction with dynamic wetting, and in particular, to what extent the wettability can determine whether the sphere is entrapped at the interface. For this purpose, a diffuse-interface immersed boundary method is adopted in the numerical simulations. We expand the parameter space considered previously, provide the phase diagrams and identify the key phenomena in the impact dynamics. Then, we propose the scaling models to interpret the critical conditions for the occurrence of sphere entrapment, accounting for the wettability of the sphere. The models are shown to provide a good correlation among the impact inertia of the drop, the surface tension, the wettability and the density ratio of the sphere to the liquid.

We report direct numerical simulation results that clearly elucidate the mechanism that leads to curvature dependence of drag enhancement (DE) in viscoelastic turbulent Taylor–Couette flow. Change in the angular momentum transport and its inherent link to transitions in vortical flow structures have been explored to depict the influence of the curvature of the flow geometry on DE. Specifically, it has been demonstrated that a transition in vortical structures with increasing radius ratio leads to weakening and elimination of the small-scale Görtler vortices and development and better organization (occupying the entire gap) of large-scale Taylor vortices as also evinced by the patterns of angular momentum current. The commensurate change in DE and its underlying mechanism are examined by contributions of convective flux and polymeric stress to the angular momentum current. The present finding paves the way for capturing highly localized elastic turbulence structures in direct numerical simulation by increasing geometry curvature in traditional turbulent curvilinear flows.