Abstract
Droplet rebound is ubiquitous on super-repellent surfaces. Conversion between kinetic and surface energies suggests that rebound suppression is unachievable due to negligible energy dissipation. Here, we present an effective approach to suppressing rebounds by incorporating bubbles into droplets, even in super-repellent states. This suppression arises from the counteractive capillary effects within bubble-encapsulated hollow droplets. The capillary flows induced by the deformed inner-bubble surface counterbalance those driven by the outer-droplet surface, resulting in a reduction of the effective take-off momentum. We propose a double-spring system with reduced effective elasticity for hollow droplets, wherein the competing springs offer distinct behavior from the classical single-spring model employed for single-phase droplets. Through experimental, analytical, and numerical validations, we establish a comprehensive and unified understanding of droplet rebound, by which the behavior of single-phase droplets represents the exceptional case of zero bubble volume and can be encompassed within this overarching framework.
Similar content being viewed by others
Introduction
Liquid droplets bounce deftly when they impact super-repellent surfaces, exhibiting remarkable elasticity owing to the efficient conversion between kinetic and surface energies1,2,3,4,5,6. Previous attempts to prevent droplet rebound rely on augmenting energy dissipation, for example, by using sticky surfaces7,8 for enhanced liquid-solid adhesion, designing compound drops9,10,11 to form a self-lubricating liquid layer for droplet arrest, elevating the temperature of water12 to produce adhesive droplets through reduced surface tension and vapor condensation, and introducing droplet additives such as polymers13,14 to increase the non-Newtonian elongational viscosity, surfactants15,16 to decrease the surface tension, and oppositely charged polyelectrolytes17 to produce hydrophilic surface defects. However, wetting transition occurs in these systems and unduly compromises liquid repellency.
With negligible energy dissipation, a bouncing droplet is analogous to a spring2,18,19. The solid-liquid contact time tc scales with the natural oscillation period as, tc ~ (m/γ)1/2 ~ (ρR3/γ)1/2, where m, γ, ρ, and R represent the mass, surface tension, volumetric density, and radius of the droplet, respectively. The contact time of a millimetric droplet impacting super-repellent surfaces is of the order of several milliseconds. In contrast, non-rebound droplets imply an infinitely large contact time tc, which requires either infinite mass or vanishing surface tension in the mass–spring model. Heavy droplets cannot bounce off solid surfaces when R is larger than the capillary length, lc ~ (γ/ρg)1/2 with g the gravitational acceleration, because gravity dominates capillary effects18. However, the non-rebound behavior of capillary-dominant droplets (R < lc) on super-repellent surfaces remains a challenge, as it is inherently impossible to achieve a surface tension of zero for any liquid.
Hollow droplets, which consist of a gas bubble surrounded by a liquid shell, possess a distinct core-shell structure that imparts them with unique physical and chemical properties compared to single-phase droplets. This distinctive morphology is prevalent in various practical processes20,21,22,23, including raindrops colliding with the Earth’s surface20, deposition of hollow spherical particles during thermal spray coating24, bubble-bursting in the preparation of nanoemulsions25, and aerosol transfer from sea26. Prior research has contributed to our understanding of the dynamic behavior exhibited by hollow droplets, such as shape oscillations27, droplet spreading21,28,29,30,31,32, and counter-jet formation21,28,29,30,31,32 during impact. However, previous studies have predominantly focused on the interactions of hollow droplets with hydrophilic or wetting solid surfaces23,28,30,31,32, as well as liquid pools33. In contrast, there has been limited exploration regarding the impact of hollow droplets on nonwetting surfaces31.
Here we present evidence of the suppressed rebound of hollow droplets upon impact with super-repellent surfaces, a phenomenon that cannot be achieved with single-phase droplets. The ability to suppress droplet rebound presents exciting prospects for the advancement of droplet-based shock absorbers, which can facilitate non-sticking liquid deposition, offering practical implications in various fields such as spray cooling34,35, self-cleaning36,37, inkjet printing38,39, agricultural spraying13,14, liquid transport40,41, and fire extinguishing42,43.
Results
Non-rebound of hollow droplets
To distinguish the bouncing dynamics of a single-phase droplet (SD) and hollow droplet (HD), we released them from the same height H0 with identical impacting velocities V0 = (2gH0)1/2 and kinetic energies. To ensure a fair comparison, the mass and liquid material of the two droplets were the same. Consequently, both droplets had the same liquid volume and characteristic radius Rh = (R03 – Rb3)1/3, where R0 and Rb denote the apparent radius of the droplet (Fig. 1a, b) and radius of the encapsulated bubble, respectively (Fig. 1b). The volume fraction of the bubble was determined as Φ = Rb3/R03. Notably, SDs can be considered a special instance of HDs with a completely filled core (Φ = 0) and, consequently, the two radii are identical, Rh = R0. HDs were generated and controlled in a co-flow microfluidic device (Supplementary Fig. 1), in which air and water were used as the core and shell fluids, respectively. Sodium dodecyl benzene sulfonate (SDBS) was added to water as the surfactant and the concentration was adjusted to stabilize HDs, by which the surface tension was in the range of 35–53 mN m−1 (see details in “Methods”). To ensure a negligible effect of gravity, we kept Rh < lc in experiments.
Figure 1c, d show the impact of an SD and HD on a superhydrophobic surface (Supplementary Fig. 2), respectively. Unlike the SD, which bounces off the surface (Fig. 1c), the HD wobbles and promptly rests on the surface after collision (Fig. 1d and Supplementary Movie 1). Compared with the SD, the HD exhibits a considerably smaller oscillating magnitude in the position of its centroid and the aspect ratio of its shape (Supplementary Fig. 3). The difference in the dynamics can be further illustrated by considering the temporal evolution of the droplet’s bottom height Hb. Figure 1e shows that the bottom height of the HD remains zero after impact, whereas that of the SD increases to be higher than zero to form multiple positive peaks owing to the rebound. Considering that non-rebound HDs remain in contact with the solid surface after impingement, we can assume that the contact time would tend towards infinity when compared to the bouncing case, neglecting droplet evaporation. In experiments, a transition from rebound to non-rebound behavior occurred at a moderate value of the bubble volume fraction (Φ ≈ 0.4), as indicated by a sharp increase in contact time (Fig. 1f). The HD-enabled rebound suppression is ubiquitous, irrespective of the droplets and super-repellent surfaces, such as air-in-hexadecane HDs on superamphiphobic surfaces, air-in-water HDs on superhydrophobic surfaces, air-in-SDBS HDs on horizontal and inclined surfaces in Leidenfrost regime, air-in-SDBS HDs impacting vertical superamphiphobic surfaces, and air-in-hexadecane HDs impacting solid surfaces under-liquid (Supplementary Movie 2).
Suppressed rebound of hollow droplets
The release height of the droplets was increased to characterize the bouncing dynamics at higher Weber numbers, We = ρV02R0/γ, where ρ represents the density of the liquid (see in “Methods—Definition of Weber number” for discussion). Upon increasing the kinetic energy of impact, both SDs and HDs were observed to rebound (Supplementary Fig. 4 and Supplementary Movie 3). We compared the contact times (tc) of SDs and HDs and made several notable observations. Firstly, tc was found to be independent of the impacting velocity (V0) for both SDs and HDs (Supplementary Fig. 5a). Secondly, tc exhibited an increase with the characteristic radius Rh in the form of tc ~ Rh3/2 (Supplementary Fig. 5b). As such, for both bouncing SDs and HDs, the contact time scales according to the inertial-capillary timescale, τ ~ (m/γ)1/2 ~ (ρRh3/γ)1/2 with m ~ ρRh3 (Supplementary Fig. 5c), consistent with previous studies2,19. Astonishingly, bouncing HDs displayed a shorter contact time compared to SDs by approximately 25% (Supplementary Fig. 5).
Despite the shorter contact time, we demonstrate that the rebound of HDs is still suppressed in comparison to SDs, as evidenced by their distinct dynamic behavior during impact. With the internal bubble, the HD exhibits a flying-saucer shape at its maximal deformation (Fig. 2a), instead of the pancake shape commonly observed for the SD. After the bounce, the HD is considerably less vibrant than the SD. The dynamics of the HD differ from those of the SD in three aspects: (i) The spreading factor β = Rmax/R0 (where Rmax is the maximum spreading radius, Fig. 2a) is smaller and decreases with Φ (Supplementary Fig. 6a), (ii) the retraction velocity Vret (defined as the time derivative of the contact radius Rcontact, as shown in Fig. 2a) is higher (Fig. 2b) and increases with Φ (Supplementary Fig. 6b), and (iii) the restitution coefficient ε = Vreb/V0 (where Vreb is the rebound velocity) is lower and decreases with Φ (Supplementary Fig. 6c). These results suggest that the rebound is markedly suppressed by encapsulating a bubble in droplets. The reduced contact time of bouncing HDs is thus attributed to the smaller spreading factor and faster retraction velocity.
To clarify the influence of the internal bubble on the droplet dynamics, we performed a scaling analysis considering Φ and We. We assumed that the bubble consists of a hemisphere with a maximum radius of Rb-max and a deformed spherical cap with the same maximum radius and a height of h0 (Fig. 2a). To conserve the volume of the bubble, we had the relationship \({R}_{{{{{{\rm{b}}}}}}-{{\max }}}^{3}+{R}_{{{{{{\rm{b}}}}}}-{{\max }}}^{2}{h}_{0} \sim {R}_{{{{{{\rm{b}}}}}}}^{3}\), which leads to \({\left({R}_{{{{{{\rm{b}}}}}}-{{\max }}}/{R}_{{{{{{\rm{b}}}}}}}\right)}^{2} \sim {R}_{{{{{{\rm{b}}}}}}}/\left({R}_{{{{{{\rm{b}}}}}}-{{\max }}}+{h}_{0}\right)\). In the case of a slightly deformed bubble, we can approximate \({R}_{{{{{{\rm{b}}}}}}-{{\max }}}\,\approx \,{R}_{{{{{{\rm{b}}}}}}}\), giving us \({R}_{{{{{{\rm{b}}}}}}-{{\max }}} \sim {R}_{{{{{{\rm{b}}}}}}}{\left[{R}_{{{{{{\rm{b}}}}}}}/\left({R}_{{{{{{\rm{b}}}}}}}+{h}_{0}\right)\right]}^{1/2}\). The height h0 can be determined as \({h}_{0} \sim {R}_{0}{{We}}^{-1/2}\) by balancing capillarity and gravity18, \(\gamma /{h}_{0} \sim \rho a{h}_{0}\), with a reinforced gravitation acceleration of \(a={V}_{0}^{2}/{R}_{0}\). Substituting Rb with R0Φ1/3, we can estimate \({R}_{{{{{{\rm{b}}}}}}-{{\max }}} \sim {R}_{0}{\varPhi }^{1/3}{(1+{\varPhi }^{-1/3}{{We}}^{-1/2})}^{-1/2}\). To model the maximally deformed HD composed of a hemispherical cap (with radius Rb-max) and a pancake-like lamella (with radius Rmax and thickness h0), volume conservation of the HD yields \({R}_{{{{{{\rm{b}}}}}}-{{\max }}}^{3}+{R}_{{{\max }}}^{2}{h}_{0} \sim {R}_{0}^{3}\) (Fig. 2a). Using the estimated values of Rb-max and h0 from the previous analysis, we find that Rmax scales as \({R}_{{{\max }}} \sim {R}_{0}{{We}}^{1/4}{[1-\varPhi {(1+{\varPhi }^{-1/3}{{We}}^{-1/2})}^{-3/2}]}^{1/2}\). The spreading factor can then be determined based on these estimations:
where C1 is a pre-factor. When Φ = 0, Eq. (1) reduces to β ~ We1/4, identical to a previous prediction of β for SDs18,44. Overall, our model provides a unified description of how β increases with We and decreases with Φ, consistent with the experimental results (Supplementary Fig. 6a). Figure 2c shows the experimental results against the theoretical prediction of β. The comparison validates Eq. (1), and the linear fitting indicates that C1 = 0.73.
After maximal spreading, the three-phase contact line of the HD retracts at a constant velocity Vret (Fig. 2b) that is determined by the balance between the capillary and inertial forces, \(\gamma /{h}_{0} \sim \rho {V}_{{{{{{\rm{ret}}}}}}}^{2}\). Considering that both the HD and SD have the same mass and volume of the liquid phase, the lamellar thickness (h0) of the HD is smaller compared to that of the SD. This is due to the liquid being pushed away by the central bubble, resulting in the spread of the liquid into a lamella with a larger maximum radius (see Fig. 2a). To obtain a more accurate estimation of h0 that highlights the differences between HD and SD by taking the bubble into consideration, we need to re-evaluate the determination of h0 by considering the force balance of \(\gamma /{h}_{0} \sim \rho a{h}_{0}\). The reinforced gravitation acceleration is reconsidered as \(a \sim {V}_{0}^{2}/({R}_{0}-{R}_{{{{{{\rm{b}}}}}}})={V}_{0}^{2}/{R}_{0}(1-{\varPhi }^{1/3})\). Consequently, h0 is updated into \({h}_{0} \sim {R}_{0}{{We}}^{-1/2}{(1-{\varPhi }^{1/3})}^{1/2}\), which is a function of and decreases with both We and Φ. This validates the assertion that the lamellar thickness of the HD is smaller than that of the SD. As a result, the driving capillary force (~γ/h0) is larger, leading to an increased retraction velocity (Vret) for the HD. Using the value of h0 estimated above, we determine Vret to be,
The Vret estimated by the proposed model is consistent with the experimental results shown in Fig. 2d, which yields a pre-factor of C2 = 0.41. For SDs with Φ = 0, Eq. (2) reduces to Vret ~ (γ/ρR0)1/2We1/4, consistent with the previous studies5,44.
The characteristic take-off velocity of a droplet scales as18 Vc ~ (γ/ρR0)1/2. However, the bubble acts as a baffle that obstructs the upward liquid flow in the HD. Assuming that the bubble-induced resisting momentum scaling as (Rb-max/Rmax)1/2mVc, the momentum balance can be obtained as \(m{V}_{{{{{{\rm{reb}}}}}}}=(1-{R}_{{{{{{\rm{b}}}}}}-{{\max }}}^{2}/{R}_{{{\max }}\,}^{2})m{V}_{{{{{{\rm{c}}}}}}}\) for the HD. By substituting the estimated values of Rb-max and Rmax into the analysis, we can derive the rebound velocity Vreb as a function of We and Φ, which is \({V}_{{{{{{\rm{reb}}}}}}} \sim [1-{{We}}^{-1/2}{\varPhi }^{2/3}{(A-\varPhi {A}^{-1/2})}^{-1}]{V}_{{{{{{\rm{c}}}}}}}\) with \(A=1+{\varPhi }^{-1/3}{{We}}^{-1/2}\). By defining ε = Vreb/V0, the restitution coefficient can be predicted as,
Equation (3) suggests that ε decreases with both We and Φ, consistent with the experimental results (Supplementary Fig. 6c). When Φ = 0, Eq. (3) reduces to ε ~ We–1/2, which represents a conventional scaling for the SD. Figure 2e shows that the results obtained experimentally and theoretically agree with each other, with a pre-factor of C3 = 0.94.
Mechanism for rebound suppression
Subsequently, we discuss the mechanism underlying rebound suppression. A superficial conjecture may attribute the suppressed rebound to the energy dissipation (Ed) in compound fluid systems45,46. Numerical simulations have shown that the energy dissipation ratio of HDs is even smaller than that of SDs during impact. Supplementary Fig. 7a, b demonstrate that only approximately 4% of the total energy (E0) is lost during the spreading and recoiling of a rebounded HD at an impact velocity of V0 = 0.270 m s−1, while this energy loss proportion is reduced to around 1% when the HD does not rebound at V0 = 0.135 m s−1. It is noteworthy that the SD bounces despite having a higher energy dissipation ratio compared to the HD at both impact velocities. Thus, the suppression of rebound cannot be attributed to energy dissipation.
With negligible energy dissipation, the sum of kinetic energy (Ek), surface energy (Es), and potential energy (Ep) remains nearly conserved, and they can be transformed among one another (Fig. 3). However, HDs and SDs follow different energy conversion pathways. For bouncing SDs, almost all the kinetic energy is converted to surface energy during spreading (Fig. 3a). In contrast, the non-rebound HDs maintain a high and nearly constant value of kinetic energy after impact (Fig. 3b). This disparity is also evident at higher impact velocities, where both HDs and SDs rebound (Fig. 3c, d). Numerical simulations reveal that only approximately 30% of the kinetic energy is converted to surface energy during the spreading stage of bouncing HDs (Fig. 3d). The suppression of HD rebound with such a low conversion efficiency from kinetic energy to surface energy calls for further exploration.
The rebound suppression can potentially be understood by the distribution of the flow field (Fig. 4 and Supplementary Movie 4). In the SD, the liquid presents a unidirectional, upward internal flow responsible for the take-off. In contrast, in the HD, the distribution of the flow velocity is inhomogeneous during droplet retraction. The encapsulated bubble always hovers in the upper part of the HD owing to buoyancy. The retraction-induced upward flow is blocked by the bubble, and the direction of the flow velocity may change from upward to downward. As both the inner and outer gas-liquid surfaces deform, the direction of the flow velocity alternately shifts between upward and downward. The velocity oscillations experienced by HDs significantly attenuate the net momentum, resulting in either a negative net momentum (Supplementary Fig. 7c) for non-rebound or a low positive value (Supplementary Fig. 7d) that lifts the HD with a suppressed height. This scenario is distinct from that of SDs, where higher rebound is associated with larger net momentum (Fig. 4).
Two capillary effects compete during the HD retraction. The outer-droplet surface retracts to push the liquid upward, whereas the dimple-shaped bottom of the inner-bubble surface produces a Laplace pressure (ΔP) that prevents the liquid from upwelling (Fig. 5a). The dimple-shaped surface is reminiscent of the shape of counter-jet reported in previous studies21,30,31,32. However, in our case, the upward growth of the dimple is inhibited due to the low We, preventing the bubble from breaking up (Figs. 4 and 5a). Considering these aspects, the recoiling HD can be modeled as a double-spring system in which a mass (representing the liquid) is sandwiched between a lower spring (representing the outer-droplet surface) with stiffness k and length lo and an upper spring (representing the inner-bubble surface) with stiffness k and length li (Fig. 5b). The two springs counteract each other: the elastic energy released by one spring is partly absorbed by the other. The evolutions of the inner surface area (Si) and outer surface area (So) exhibit opposite trends when the HD retracts (Fig. 5c, d), which validates the results of the previous analysis. The decrease in So (or Si) is compensated by the increase in Si (or So), which accounts for the increase in the total surface area (Si + So) during HD retraction, in contrast to the SD case (Supplementary Fig. 7e, f).
The double-spring model is characterized by a reduced effective elasticity. Unlike in the case of SDs, where released surface energy is directly converted into kinetic energy, impacting HDs exhibit an extra energy conversion between the inner and outer surfaces. In this process, the surface energy released from one surface is absorbed by the other, resulting in a reduced conversion efficiency between surface energy and kinetic energy. With the two surfaces analogous to the two competing springs in the double-spring system, we can assume that k ~ γ, lo ~ R0, and li ~ Rb. We propose a qualitative estimation of the apparent surface tension (γa), which accounts for the overall counteractive capillary effect as the difference in capillarity between the two surfaces, R0γa ~ (R0γ-Rbγ). Scaling analysis suggests that γa ~ (R0 − Rb)γ/R0 = (1 − Φ1/3)γ. This apparent surface tension is lower than the actual liquid surface tension (γ) and decreases with the bubble volume fraction (Φ).
Based on this reasoning, the apparent surface tension decreases asymptotically to zero when the bubble volume fraction approaches unity. This result suggests that high-Φ HDs can behave as droplets with nearly zero surface tension when impacting super-repellent surfaces. The rebound suppression can be strengthened by increasing the number of bubbles in the HD (Supplementary Movie 5), in which case the two-spring model extends to a multiple-spring model (Supplementary Fig. 8) with enhanced counteractive capillary effects.
Phase diagram
According to the energy-conversion argument, the effects of the inner-bubble surface are twofold: (i) The dimple-like surface deformation absorbs part of the kinetic energy, which attenuates the upward motion of the liquid, and (ii) the recovery of the inner-surface deformation converts the surface energy back into the kinetic energy of the liquid flowing downward. Despite the high total kinetic energy of the HD, only a small fraction is responsible for droplet rebound. We thus speculate that droplet rebound will be completely suppressed when a critical amount of the impacting kinetic energy (~\({R}_{{{{{{\rm{h}}}}}}}^{3}\rho {V}^{2}\)) is absorbed by the inner-bubble surface (~\({R}_{b}^{2}\gamma\)). This hypothesis implies an energy-balance scaling of \(({R}_{0}^{3}-{R}_{b}^{3})\rho {V}^{2} \sim {R}_{b}^{2}\gamma\). Therefore, a critical Weber number pertaining to the transition between the rebound and non-rebound conditions can be found:
with C4 ≈ 1/7 determined empirically from the experimental results. The phase diagram in Fig. 6 shows that our model can accurately predict the bouncing dynamics of HDs.
Discussion
Overall, the results of this study highlight that the droplet structure can be simply modified to suppress the rebound, without altering the droplet composition or surface properties. Notably, bubbles are spontaneously entrained during droplet impact3. The influence of these entrained bubbles on the bouncing dynamics may bring awareness of their new scientific interest. The enhanced deposition of hollow droplets on super-repellent surfaces can potentially be leveraged, for example, for increasing the self-cleaning efficiency via the enlarged contact area of the sliding/rolling droplets (Supplementary Fig. 9 and Supplementary Movie 6) and increasing the efficiency of spray cooling owing to the prolonged contact time. Similar phenomena can also be observed in macroscopic systems. For example, the rebound of a compound balloon (containing water and an air-filled balloon) is conspicuously suppressed (Supplementary Fig. 10 and Supplementary Movie 7) compared with that of a single-phase balloon (filled with only water). This finding can extend the applications of shock-absorbing hollow systems to a wider horizon, for instance, to develop throw-type fire-extinguishing balls43 that must be precisely positioned at the target fire to release the interior content without bouncing.
Methods
Hollow droplet generation
A co-flow microfluidic device was used to generate hollow droplets, in which air and liquid were used as the inner and outer phase fluids, respectively (Supplementary Fig. 1). The inner diameter of the inner capillary ranged from 0.06 mm to 0.1 mm, and that of the outer capillary ranged from 0.25 mm and 0.9 mm. Two syringe pumps (Longer Pump, LSP01-1A) were used to drive the two working fluids and control the flow rates. To maintain the dripping mode for the generation of monodisperse hollow droplets, the flow rate of the outer phase liquid was set in a range between 2 mL h−1 and 8 mL h−1. The volume fraction of the encapsulated bubble was adjusted by controlling the inner airflow in a range between 1 mL h−1 and 5 mL h−1.
Water, with sodium dodecyl benzene sulfonate (SDBS, Macklin) as the surfactant, was used as the outer fluid for hollow droplet impact on superhydrophobic, superamphiphobic, and Leidenfrost-state surfaces. In addition to pure water, aqueous solutions with three different SDBS concentrations were used in experiments: impingement of 0.003 wt% SDBS (surface tension of 53 mN m−1) droplets on superhydrophobic surfaces; impingement of 0.1 wt% SDBS (surface tension of 43 mN m−1) droplets on smooth aluminum plates in the Leidenfrost regime; impingement of 0.2 wt% SDBS (surface tension of 35 mN m−1) droplets on superamphiphobic surfaces. Besides, n-hexadecane (98%, Aladdin) was used as the outer liquid phase to produce air-in-oil HDs for impact experiments. The selection of various liquids and solid surfaces allowed us to demonstrate the universality of the phenomenon of rebound suppression by HDs across different liquid compositions and solid surface properties.
Fabrication of superhydrophobic surfaces
The superhydrophobic coating suspension was made by dispersing fumed silica nanoparticles (Aerosil R805, Evonik) in ethanol at a concentration of 2 wt%. A glass substrate was cleaned with deionized water first and then ethanol in an ultrasonic bath. After drying, the substrate was coated with the superhydrophobic solution twice by dip-coating. Then the coated surface was baked at 60 °C for 10 min in a drying oven.
Fabrication of superamphiphobic surfaces
The superamphiphobic coating was fabricated by the following steps. First, 15 mL ammonia solution (28–30%, Aladdin) was added to 45 mL ethanol (absolute, Anaqua), which was stirred at an agitation speed of 400 RPM and heated to a temperature of 60°C in a water bath. Then, 2 mL tetraethyl orthosilicate (>99%, Macklin) was added into the solution drop by drop. When the solution became opaque, 0.1 mL 1H,1H,2H,2H-perfluorodecyltriethoxysilane (>98.0%, TCI) was added to the mixture. After 24 h of reaction, a superamphiphobic solution was obtained, which contained fluorinated silica nanoparticles. A glass substrate was dip-coated with the solution and then incubated at 80 °C for 1 h in a drying oven. The coating and drying processes were repeated two to three times until a uniform superamphiphobic nanoparticle coating was obtained on the glass surface.
Droplet impact
Hollow droplets impact the solid super-repellent surface directly after their generation. The height of the droplet impact was tunable by a sliding rail. The impact process was recorded using a high-speed camera (Fastcam Mini, Photron) at a frame rate ranging from 5000 frames per second (fps) to 6520 fps. The shape variation of droplets during impingement was analyzed using ImageJ (National Institutes of Health).
Leidenfrost experiment
A smooth, flat, and well-polished aluminum plate was placed on a copper heater equipped with a temperature controller. The real-time surface temperature of the aluminum plate was measured using a K-type thermocouple and monitored using a recorder. In Leidenfrost experiments, the temperature of the aluminum plate was controlled at 210 ± 5 °C. A tilting stage was used to maintain the tilt angle at 10° for droplets impacting inclined surfaces.
Self-cleaning
A superamphiphobic surface was covered with a layer of Fe3O4 nanoparticles (Zhonghangzhongmai Metal Material Co. Ltd) as the dust. The surface was tilted at an angle of 10° by a tilting stage. A hollow droplet impacted the tilted surface and absorbed the dust particles when rolling down the surface.
Balloon impact
The compound core-shell balloon was made by containing water as the shell and a hollow balloon (filled with air) as the core. In the control group, a single-phase balloon was filled with only water. For a fair comparison, the two balloon systems contained the same mass of water (500 g) and were released from the same height to impact the ground.
Definition of Weber number
In the study, the Weber number (We) was defined using the density of the liquid phase (ρ) rather than the apparent density (ρa) of the HD. The mass of the HD was calculated by \(m=4\pi {R}_{0}^{3}{\rho }_{a}/3=4\pi {R}_{{{{{{\rm{h}}}}}}}^{3}\rho /3\) where \({R}_{{{{{{\rm{h}}}}}}}^{3}={R}_{0}^{3}-{R}_{{{{{{\rm{b}}}}}}}^{3}\). We then have \({\rho }_{a}={R}_{{{{{{\rm{h}}}}}}}^{3}\rho /{R}_{0}^{3}=({R}_{0}^{3}-{R}_{{{{{{\rm{b}}}}}}}^{3})\rho /{R}_{0}^{3}=(1-{R}_{{{{{{\rm{b}}}}}}}^{3}/{R}_{0}^{3})\rho=\left(1-\varPhi \right)\rho\). Based on this, the ρa-based Weber number is defined as \({{We}}_{{{{{{\rm{a}}}}}}}={\rho }_{a}{V}_{0}^{2}{R}_{0}/\gamma={We}\left(1-\varPhi \right)\) where \({We}=\rho {V}_{0}^{2}{R}_{0}/\gamma\). It is important to note that Wea depends on Φ while We does not. To investigate the bouncing behavior of HDs systematically and explicitly, it is crucial to consider independent groups of non-dimensional parameters (such as We and Φ). With this consideration, the ρ-based We was used instead of ρa-based Wea in this study, and this choice does not affect the accuracy of the models presented in Eqs. (1)–(4), because Wea can be readily determined using the values of We and Φ.
Numerical simulation
For the numerical simulation of the droplet impact, we used the open-source code Basilisk (basilisk.fr)47. The code uses a second-order in space, time-splitting projection method for incompressible two-phase flows, and a geometric Volume of Fluid (VoF) method to accurately track the liquid-gas interface.
Basilisk discretizes the domain using an Octree grid (Quadtree grid in 2D), which can dynamically refine and coarsen its resolution during computation. This adaptive mesh refinement (AMR) feature allows for the finest grid resolution around the droplet interface while keeping relatively lower resolutions elsewhere to reduce the computational cost.
In our case, the simulation was enforced to be axis symmetry, so simulating a 2D cross-section of the droplet is adequate. Given a square numerical domain of the side length L, a cell refined by n times has the side-length of L/2n. With L = 5.74 mm and the maximum level of the refinement n = 13, the solver can resolve spatial scales as small as 0.7 microns. This fine spatial resolution is crucial to avoid numerical breakups, and Basilisk’s AMR capability makes it possible at a reasonable computational cost.
Data availability
The data that support the findings of this study are available from the corresponding authors. Source data are provided with this paper.
Code availability
The code that supports the findings of this study is available from the corresponding authors.
References
Schutzius, T. M. et al. Spontaneous droplet trampolining on rigid superhydrophobic surfaces. Nature 527, 82–85 (2015).
Richard, D., Clanet, C. & Quéré, D. Contact time of a bouncing drop. Nature 417, 811 (2002).
Josserand, C. & Thoroddsen, S. T. Drop impact on a solid surface. Annu. Rev. Fluid Mech. 48, 365–391 (2016).
Yarin, A. L. Drop impact dynamics: splashing, spreading, receding, bouncing…. Annu. Rev. Fluid Mech. 38, 159–192 (2006).
Zhu, P., Chen, C., Nandakumar, K. & Wang, L. Nonspecular reflection of droplets. Small 17, 2006695 (2021).
Vasileiou, T., Gerber, J., Prautzsch, J., Schutzius, T. M. & Poulikakos, D. Superhydrophobicity enhancement through substrate flexibility. Proc. Natl Acad. Sci. USA 113, 13307–13312 (2016).
Zhu, P., Kong, T., Tang, X. & Wang, L. Well-defined porous membranes for robust omniphobic surfaces via microfluidic emulsion templating. Nat. Commun. 8, 15823 (2017).
Bhushan, B. & Her, E. K. Fabrication of superhydrophobic surfaces with high and low adhesion inspired from rose petal. Langmuir 26, 8207–8217 (2010).
Damak, M., Ruiter, J. D., Panat, S. & Varanasi, K. K. Dynamics of an impacting emulsion droplet. Sci. Adv. 8, eabl7160 (2022).
Han, X. et al. Slippery damper of an overlay for arresting and manipulating droplets on nonwetting surfaces. Nat. Commun. 12, 3154 (2021).
Blanken, N., Saleem, M. S., Thoraval, M.-J. & Antonini, C. Impact of compound drops: a perspective. Curr. Opin. Colloid Interface Sci. 51, 101389 (2021).
Zhu, P., Chen, R. & Wang, L. Topography-directed hot-water super-repellent surfaces. Adv. Sci. 6, 1900798 (2019).
Bergeron, V., Bonn, D., Martin, J. Y. & Vovelle, L. Controlling droplet deposition with polymer additives. Nature 405, 772–775 (2000).
Bartolo, D., Boudaoud, A., Narcy, G. & Bonn, D. Dynamics of non-Newtonian droplets. Phys. Rev. Lett. 99, 174502 (2007).
Song, M. et al. Controlling liquid splash on superhydrophobic surfaces by a vesicle surfactant. Sci. Adv. 3, e1602188 (2017).
Gaskin, R. E., Steele, K. D. & Forster, W. A. Characterising plant surfaces for spray adhesion and retention. N. Z. Plant Prot. 58, 179–183 (2005).
Damak, M., Mahmoudi, S. R., Hyder, M. N. & Varanasi, K. K. Enhancing droplet deposition through in-situ precipitation. Nat. Commun. 7, 12560 (2016).
Biance, A.-L., Chevy, F., Clanet, C., Lagubeau, G. & Quéré, D. On the elasticity of an inertial liquid shock. J. Fluid Mech. 554, 47–66 (2006).
Bird, J. C., Dhiman, R., Kwon, H.-M. & Varanasi, K. K. Reducing the contact time of a bouncing drop. Nature 503, 385–388 (2013).
Villermaux, E. & Bossa, B. Single-drop fragmentation determines size distribution of raindrops. Nat. Phys. 5, 697–702 (2009).
Gulyaev, I. P. & Solonenko, O. P. Hollow droplets impacting onto a solid surface. Exp. Fluids 54, 1432 (2012).
Lohse, D. Inside a micro-reactor. Nature 418, 381–383 (2002).
Li, D., Duan, X., Zheng, Z. & Liu, Y. Dynamics and heat transfer of a hollow droplet impact on a wetted solid surface. Int. J. Heat Mass Transf. 122, 1014–1023 (2018).
Solonenko, O. P., Gulyaev, I. P. & Smirnov, A. V. Plasma processing and deposition of powdered metal oxides consisting of hollow spherical particles. Tech. Phys. Lett. 34, 1050–1052 (2008).
Feng, J. et al. Nanoemulsions obtained via bubble-bursting at a compound interface. Nat. Phys. 10, 606–612 (2014).
Tripathi, M. K., Sahu, K. C. & Govindarajan, R. Dynamics of an initially spherical bubble rising in quiescent liquid. Nat. Commun. 6, 6268 (2015).
Balla, M., Tripathi, M. K. & Sahu, K. C. A numerical study of a hollow water droplet falling in air. Theor. Comput. Fluid Dyn. 34, 133–144 (2020).
Wei, Y. & Thoraval, M.-J. Maximum spreading of an impacting air-in-liquid compound drop. Phys. Fluids 33, 061703 (2021).
Nasiri, M., Amini, G., Moreau, C. & Dolatabadi, A. Flattening of a hollow droplet impacting a solid surface. J. Fluid Mech. 962, A1 (2023).
Gulyaev, I. P., Solonenko, O. P., Gulyaev, P. Y. & Smirnov, A. V. Hydrodynamic features of the impact of a hollow spherical drop on a flat surface. Tech. Phys. Lett. 35, 885–888 (2009).
Nasiri, M., Amini, G., Moreau, C. & Dolatabadi, A. Hollow droplet impact on a solid surface. Int. J. Multiphase Flow 143, 103740 (2021).
Naidu, D. P. & Dash, S. Impact dynamics of air-in-liquid compound droplets. Phys. Fluids 34, 073604 (2022).
Zhu, S., Kherbeche, A., Feng, Y. & Thoraval, M.-J. Impact of an air-in-liquid compound drop onto a liquid surface. Phys. Fluids 32, 041705 (2020).
Jiang, M. et al. Inhibiting the Leidenfrost effect above 1,000 °C for sustained thermal cooling. Nature 601, 568–572 (2022).
Shiri, S. & Bird, J. C. Heat exchange between a bouncing drop and a superhydrophobic substrate. Proc. Natl Acad. Sci. USA 114, 6930–6935 (2017).
Blossey, R. Self-cleaning surfaces — virtual realities. Nat. Mater. 2, 301–306 (2003).
Lu, Y. et al. Robust self-cleaning surfaces that function when exposed to either air or oil. Science 347, 1132–1135 (2015).
Minas, C., Carnelli, D., Tervoort, E. & Studart, A. R. 3D printing of emulsions and foams into hierarchical porous ceramics. Adv. Mater. 28, 9993–9999 (2016).
Visser, C. W., Amato, D. N., Mueller, J. & Lewis, J. A. Architected polymer foams via direct bubble writing. Adv. Mater. 31, 1904668 (2019).
Feng, S. L. et al. Three-dimensional capillary ratchet-induced liquid directional steering. Science 373, 1344–1348 (2021).
Prakash, M., Quéré, D. & Bush, J. W. Surface tension transport of prey by feeding shorebirds: the capillary ratchet. Science 320, 931–934 (2008).
Ratzer, A. History and development of foam as a fire extinguishing medium. Ind. Eng. Chem. 48, 2013–2016 (1956).
Held, J. M. New technique for fire fighting-large scale open fires. US patent US11/546,272 (2008).
Clanet, C., Béguin, C., Richard, D. & Quéré, D. Maximal deformation of an impacting drop. J. Fluid Mech. 517, 199–208 (2004).
Sauret, A., Boulogne, F., Cappello, J., Dressaire, E. & Stone, H. A. Damping of liquid sloshing by foams. Phys. Fluids 27, 022103 (2015).
Chiu, S.-L. & Lin, T.-H. Experiment on the dynamics of a compound drop impinging on a hot surface. Phys. Fluids 17, 122103 (2005).
Popinet, S. A quadtree-adaptive multigrid solver for the Serre–Green–Naghdi equations. J. Comput. Phys. 302, 336–358 (2015).
Acknowledgements
The financial support from the Research Grants Council of Hong Kong (Project No. CityU21213621), Shenzhen Science and Technology Program (Project No. JCYJ20220530140812028), and City University of Hong Kong (Project Nos. 9610502 and 7005936) is gratefully acknowledged.
Author information
Authors and Affiliations
Contributions
P.Z. conceived the research. P.Z. and Y.Z. designed the experiments. W.Z. fabricated the superamphiphobic coating. S.W. and Y.Z. fabricated the microfluidic device for droplet generation. Y.Z. performed the experiments. P.Z. performed the theoretical analysis. C.Z. performed the numerical simulation. P.Z., Y.Z., and C.Z. analyzed the data. P.Z. and Y.Z. wrote the manuscript. P.Z. supervised the research. All authors commented on the paper.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. A peer review file is available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Source data
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Zhou, Y., Zhang, C., Zhao, W. et al. Suppression of hollow droplet rebound on super-repellent surfaces. Nat Commun 14, 5386 (2023). https://doi.org/10.1038/s41467-023-40941-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-023-40941-3