Steerable drops on heated concentric microgroove arrays

Apr 19, 2023

Nature Communications volume 13, Article number: 3141 (2022) Cite this article

3854 Accesses

8 Citations

1 Altmetric

Metrics details

Guided drop transport is of great importance in various water and thermal management technologies. Unidirectional drop transport on a hot surface has been widely developed, but a bidirectional reversal is still challenging. Here, we report a steerable transport of drop impinging on heated concentric microgroove arrays, on which the directionality of drop transport is dictated by the drop boiling modes. In the transition boiling state, the driving force originated from the Laplace pressure difference rendered by the microgrooves, which enables the drop rebounding towards the center of curvature. While in the film boiling state, a net force towards the opposite side is generated between the grooves and the penetrated liquid, that drives the drop far away from the center of curvature. Our experimental and theoretical results uncover that the lateral displacement is controlled by both the Weber number and off-center distance. These findings strengthen our fundamental understanding of drop impact dynamics at high temperatures and are essential for effective cooling of hot-spot cores and drop sieving.

Rectification of droplet transport on hot surfaces is both of fundamental interest and practical importance in spray cooling1,2, drag reduction3,4,5, and power generation6,7. In recent years, particular attention have been paid to the surface temperature above the so-called Leidenfrost point, where the drop levitates on its own vapor layer. Distinct from ambient conditions where the drop rectification is achieved by harnessing gradients of surface energy and often hindered by contact line pinning, in this Leidenfrost state, the underlying vapor layer of the drop allows the elimination of the contact line pinning and accordingly the associated interfacial friction. Therefore, various asymmetric structures such as macroscale or nanoscale ratchets, have been developed to realize unidirectional drop transport in the Leidenfrost state, resulting from the shear stress generated by the asymmetric vapor ejection8,9,10. However, the underlying vapor layer of the drop that allows negligible interfacial friction also causes a huge heat transfer resistance, which is adverse to effective thermal management.

Spray cooling, as an important cooling technique in many thermal management devices, can release the heat of devices by evaporating drops violently. The transient contact of the impinging drop on the hot surface is accompanied by transient phase change and intensive vibration. Thus, an impinging drop allows random guidance instead of full evaporation in situ, causing unsatisfactory cooling efficiency. Moreover, the temperature distribution tends to be asymmetrical on hot surfaces due to the random droplet transport, which in turn causes a thermal Marangoni effect11,12. Recently, researchers have proposed that directional drop rectification to the preferential region for enhanced heat transfer can be achieved on straight post arrays with a density gradient1 and regularly patterned posts with Janus-mushroom structures13. However, in all of these studies, the drop transport direction is not steerable for a fixed surface design, and the movement of drops are difficult to control quantitatively although the movement directions are explicit. Moreover, the unidirectional drop rectification relies on the elaborated surface structure design and scrupulous control of the specific location where the drop impacts. Therefore, controlled drop vectoring on hot surfaces remains a big challenge.

In this research, we design concentric microgroove arrays to achieve bidirectional drop transport by mediating the surface temperature, which goes beyond the traditional ways to realize unidirectional drop transport by chemical gradient14,15, ratchets8,10, and curvatures16,17. Impacting drop transports towards the center of curvature at temperatures below the Leidenfrost point, while towards the direction far away from the center of curvature at temperatures above the Leidenfrost point, i.e., the drop transport direction is dictated by the synergistic action of surface structure and boiling states. Simple scaling shows that the drop lateral transport distance depends on the Weber number and the off-center distance between the impact position and center of curvature. In addition, we demonstrate that this versatile and robust strategy could be applied in high-efficiency thermal transfer and drop sieving.

The disc surface with a diameter d = 2 cm, as shown in the schematic diagram in Fig. 1a and Supplementary Fig. 1, consists of an array of concentric microgrooves with uniform ridge width W = 40 μm, groove width S = 40 μm, and height H = 20 μm, and the surface wettability at room temperature is discussed in Supplementary Note 1. We consider a water drop with radius R = 1.18 mm impinging on the surface with an off-center distance between the impact point and center of curvature r = d/4 under We = 16.9 at different temperatures. The Weber number is defined as We = ρU2R/γ, with ρ and γ being the mass density and surface tension of water, and U the impact velocity. At T = 200 °C, after the drop reaches its maximum spreading diameter, it boils violently until completely evaporates at 102.1 ms at its original location (Supplementary Fig. 3). In this case, the drop is in a contact boiling state that enables a high-efficiency heat transfer. When the temperature increases, i.e., T = 250 °C, a weak boiling is observed. Note that during recoiling, when the width of the drop along the radial direction reaches its minimum, the drop forms a bean shape (i.e., at t = 9.8 ms in Fig. 1b), and finally bounces towards the center of curvature, resulting in a lateral displacement ΔlL = 1.84 mm between the impact point and the point where the drop falls back to the substrate (as shown in Fig. 1b and left column of Supplementary Movie 1). Surprisingly, the impinging drops converge into the center of curvature regardless of the releasing point (Supplementary Movie 2) at this temperature. When the temperature further increases, i.e., T = 350 °C, there is no obvious boiling, and the drop bounces back entirely, which is similar to that on superhydrophobic surfaces18,19. But note that, conversely, the drop bounces far away from the center of curvature with a lateral displacement ΔlR = 1.76 mm, as shown in Fig. 1c and right column of Supplementary Movie 1.

a Schematic diagrams showing bidirectional motions of drops impacting on the heated concentric microgroove arrays at different temperatures. b, c Sequential images of the drop impact on the surfaces at T = 250 and 350 °C, respectively. Here, We = 16.9 and r = d/4. ΔlL and ΔlR denote the lateral displacements along the left and right directions, respectively. Scale bars showing in b, c are 1 mm.

The directional drop transport is obvious when comparing the spreading and retracting processes during the impingement. Fig. 2a, b show the time evolution of normalized drop contact line length ratio K = L(t)/R, where L(t) represents the distance between the left/right contact point and the initial impact point over time. At T = 250 °C, as shown in Fig. 2a, the length ratios of the left and right contact lines are comparable with each other with the left one slightly larger, and finally, the whole drop bounces towards left. Note that, the right contact line recoils faster than the left one after t = 8 ms, due to the sheltering of the left side of the center rim by the front-most liquid. At T = 350 °C, however, the maximum length ratio of the left contact line is remarkably larger than the right one. After t = 4 ms, the left contact line recoils much faster than the right one (Fig. 2b), and finally, the drop bounces to the right.

a, b Variation of contact line length ratio K of an impinging drop on concentric microgroove arrays (r = d/4) on the left (red) and right (black) sides under We = 16.9 at T = 250 and 350 °C, respectively. The insets in a, b represent the side view of maximum spreading. c Lateral displacements Δl of bouncing drops as a function of temperature under various Weber numbers. A positive Δl, i.e., ΔlL, indicates that drops rebound towards the center of curvature, while a negative one, i.e., ΔlR, towards the opposite direction. The insets show the recoiling of impacting drops shortly after they reach the maximum spreading at different boiling states. The error bars of the data in a–c denote the standard deviation of three measurements. Source data are provided as a Source Data file.

We next plot the variation of the lateral displacement Δl of the impinging drops as a function of the substrate temperature T for different We. As shown in Fig. 2c, when T is higher than the boiling point but less than 215 °C, an in situ contact boiling is observed, in which the drop is difficult to take off from the substrate until its complete evaporation. The slightly off-center distance results from measurement error due to the violent boiling of the drop. When the temperature increases, e.g., 215 °C ≤ T ≤ 295 °C, the impacting drop is in a transition boiling state, in which the drop always rebounds towards the center of curvature (Δl > 0). However, when the temperature further increases, e.g., T > 315 °C, the drop is in a film boiling state, in which the drop always rebounds far away from the center of curvature (Δl < 0). Note that, in the temperature range 295 °C ≤ T ≤ 315 °C, there is a transition during which a hybrid boiling exists, and the drop rebound direction is random with a relatively small lateral displacement. The drop would rebound to the left when the transition boiling dominates under a high Weber number, e.g., data points We = 13.7 and 16.9 at T = 300 °C, while to the right when the film boiling prevails under a low Weber number, e.g., data points We = 4.1 and 7.3 at T = 300 °C. This is because the Leidenfrost temperature increases with the Weber number and it is not a fixed value20.

Further, we investigate the influence of off-center distance r and Weber number We on the lateral displacement Δl at different temperatures. As shown in Fig. 3a, b, in the transition boiling state at T = 250 °C, it is obvious that ΔlL decreases with r for a fixed We, but increase with We for a specific r. Similar tendencies are observed in the boiling state at T = 350 °C, as shown in Fig. 3c, d. These facts suggest that, despite remarkably different temperatures of the substrates and rebounding directions of the drops, ΔlL and ΔlR are closely related to r and We.

a, b Variation of ΔlL at T = 250 °C as a function of r under various We and as a function of We for various r, respectively. c, d Variation of ΔlR at T = 350 °C as a function of r under various We and as a function of We for various r, respectively. The error bars of the data in a–d denote the standard deviation of three measurements. Source data are provided as a Source Data file.

We now explore the underlying mechanisms accounting for the directional rebounding in the transition boiling at T = 250 °C. To address the question, we first consider the drop deformation in the spreading and receding processes, which plays an important role in the interaction between the drop and the substrate. As shown in Fig. 1b, violent boiling occurs when the drop contacts the substrate, which may render a partial wetting of the contact region. In addition, due to a smaller contact angle hysteresis along the circumferential direction compared with the radial direction of the substrate, the drop spreads larger along the grooves and evolves into an elliptical shape (i.e., 3.6 ms) rather than a perfect circle at its maximum extension. In this case, more liquid flows to the two ends along the circumferential direction, leaving a relatively thinner layer on the left and right sides. Since the recoiling drop dewets at the Taylor–Culick velocity which is inversely proportional to one-half of the thickness of the liquid layer21,22, the left and right sides recoil faster than the upper and lower ends, and finally, a bean-shaped drop is formed at the end of the retraction, i.e., at 9.8 ms (Supplementary Discussion 1). This behavior is remarkably different from that in previous research, in which an impacting drop recoils in a circular shape23 or evaporates immediately20. Note that the arc length of the bean-shaped drop at t = 9.8 ms is close to its maximum spreading width Dmax at t = 3.6 ms (Fig. 1b), and therefore we have 2rα ≈ Dmax, where α is the corresponding central half-angle, as shown in the inset of Fig. 4a. Following the scaling24 of Dmax/D ~ We1/4, we obtain

where D = 2 R is the initial drop diameter. Indeed, this is observed for different values of r and We, where the experimental data is in good agreement with Eq. 1 with a prefactor 1.25 based on the best fit (Fig. 4a).

a Relationships between 2αr/D and We. b The lateral displacement ΔlL as a function of We and the initial off-center distance r. Here, R = 1.18 mm, and the red solid lines in a, b are respectively the best fits of Eqs. 1 and 3 to the experimental data. c The top panel sketching the morphology of the drop during recoiling from the top view. The red and blue arrows represent the forces that the grooves exert on the liquid. The bottom panel sketching the wetting state of the sections a-a. The left and right parts of the drop are marked by the light red and light blue shadows, respectively. d The lateral displacement ΔlR as a function of We and r. The solid line is the best fit of Eq. 5. The error bars of the data in a, b, d denote the standard deviation of three measurements. Source data are provided as a Source Data file.

The driving force that propels the drop to the curvature center of the substrate could be deduced based on the asymmetric profiles of the drop. From the top view in the left column of Supplementary Movie 1, the asymmetry of the left and right profiles of the drop is particularly obvious in the recoiling state, resulting in the different curvatures on the left and right sides of the torus. Consequently, a Laplace pressure difference δP is rendered, which not only accounts for the asymmetric deformation of the drop but also the lateral motion of the entire drop. An exact deduction of δP in the whole recoiling stage is challenging due to the complex and changeable morphology of the drop. However, δP could be perceived based on the bean-shaped drop at 9.8 ms in Fig. 1b. In other words, δP in the recoiling stage must have a strength that is able to force the circular-shape drop into a bean-shaped one, and δP exerted on the bean-shaped drop could be considered as the characteristic value of which in the recoiling stage. In this regard, we make analyses based on the bean-shaped drop where there are two principal curvature radii, i.e., azimuthal radii r1 = r − e/2 and r2 = r + e/2 closed to the left (point A) and right (point B) sides of the liquid-vapor menisci (Supplementary Fig. 5), respectively, while meridional radii e/2 for both of the left and right sides (points A and B), where e is denoted as the width of the bean-shaped drop (inset of Fig. 4b). Based on the Young–Laplace equation25, a pressure difference between points A and B is obtained, i.e., δP = PB – PA ≈ γ(1/r1 + 1/r2) > 0, which renders an inner pressure difference that propels the drop to move to the left side. Since e ≪ r, δP ≈ γ(1/r1 + 1/r2) ≈ 2γ/r ~ γ/r is expected (Supplementary Discussion 2). Considering δP works along the horizontal direction, the force exerting on the drop along the horizontal direction is FL = δP·A, where A ~ Dmax2 characterizes the contact area between the solid and liquid interfaces in the recoiling stage rather than the specific moment of the triggered bean shape. Finally, we obtain

Thus, the acceleration along the lateral direction that propels the drop to the center of the curvature is characterized as aL = FL/m, with m = 4πρR3/3 being the mass of the drop. Accordingly, the corresponding lateral speed could be estimated as VL = aLτ0, with τ0 ~ (ρR3/γ)1/2 being the classic solid-liquid contact time23 during the drop impingement. Considering a free falling of the drop after detaching from the substrate with an upward speed \({U}_{\perp }\), the horizontal flying time span is \(\varDelta t \sim {U}_{\perp }/g\), with g being the acceleration of gravity. Here, \({U}_{\perp }\) is the characteristic rebounding velocity26 and could be estimated as \({U}_{\perp }={[\gamma /(\rho R)]}^{1/2}\). Finally, we have \(\varDelta {l}_{{{\mathrm{L}}}} \sim {V}_{{{\mathrm{L}}}}\varDelta t \sim ({F}_{{{\mathrm{L}}}}{\tau }_{0}/m)({U}_{\perp }/g)\) and a combination of Eq. 2 leads to

In Eq. 3, lc = [γ/(ρg)]1/2 is defined as the capillary length, which is about 2.73 mm for water in the ambient laboratory environment. As shown in Fig. 4b, the lateral distance ΔlL estimated by Eq. 3 is very well consistent with the experimental data, with a prefactor 0.32 based on the best fit.

We now explore the underlying mechanism accounting for the lateral displacement of the drop in the film boiling state. Different from the transition boiling state, the temperature of the film boiling state is higher, and a thin layer of air film is generated underneath the drop that avoids direct contact between the liquid and solid substrate. Moreover, the existence of the air film significantly reduces the adhesion between the drop and the grooves, which endows a circular shape of the drop during the whole impingement (Fig. 1c). Here, we hypothesize that the force propelling the drop to move laterally would be a result of the interaction between the liquid and the substrate. To address this point, we first check the morphology of the liquid in the contact region. Upon impact, a balance reaches between the dynamic pressure ρU2/2 and the capillary pressure γ/R*, resulting from the dynamics of the drop and the liquid-vapor meniscus generated in the grooves, respectively. Here, R* ~ S2/ε represents the curvature radius of the liquid-vapor meniscus in the grooves27, with ε being the penetration depth. The above analysis leads to an estimate of the penetration depth, i.e., ε ~ S2We/R (see details in Supplementary Discussion 3). The parameters used in our experiments (i.e., S = 40 μm, R ~ 1 mm, and We ~ 10) give rise to ε ~ 10 μm, whose order is comparable to the depth of the groove (ε ~ H) and suggests a penetration of the liquid into the grooves, and therefore, the recoiling liquid will exert forces on the inner wall of the grooves during the drop retraction process. In turn, the inner wall of the grooves will exert equivalent forces on the liquid in the opposite direction, as represented by the red and blue arrows in Fig. 4c. To give a quantitative analysis, we use μ to denote the aforementioned force per length (N/m) in a single groove, i.e., μ ~ ρU//2H, with U// ~ [γ/(ρh)]1/2 and h = RWe−1/2 being the characteristic lateral velocity and thickness of the recoiling liquid, respectively28. After integration over the whole contact region, the resultant acting force could be obtained (see details in Supplementary Discussion 4). For simplicity, we consider the drop in two parts, i.e., the left and right parts which are marked by the red and blue colors, respectively. On the left side of the drop (red shadow in Fig. 4c), the acting force on the drop from the grooves could be estimated as Fleft ≈ 2μ(R2 − R3/r)/(W + S). Similarly, the force on the right (light blue shadow in Fig. 4c) is Fright ≈ 2μ(R2 + R3/r)/(W + S). Moreover, the net force FR = Fright − Fleft exerted on the drop is obtained as

which points to the right and enables the drop to bounce far away from the center of the curvature. Considering \(\varDelta {l}_{{{\mathrm{R}}}} \sim {V}_{{{\mathrm{R}}}}\varDelta t \sim ({F}_{{{\mathrm{R}}}}{\tau }_{0}/m)({U}_{\perp }/g)\) as analyzed before, we finally obtain

Equation (5) suggests that the lateral displacement ΔlR of the drop impacting the concentric microgroove arrays at the film boiling state is also controlled by both the Weber number and the initial off-center distance, which is well confirmed by the data in Fig. 4d. The prefactor of Eq. 5 obtained from the best fit to the experimental data is 1.15. In addition, the obtained forces Fleft and Fright in the above give us fundamental insights to understand the drop spreading and recoiling dynamics in Fig. 2b. The above result Fleft < Fright indicates a smaller resistance exerted on the liquid, that the drop spreads faster and the contact length is larger to the left than that to the right. On the contrary, the left contact line moves faster than the right one during the recoiling stage.

To further validate the universal applicability of the steerable drop transport method by employing concentric microgroove arrays, as well as the proposed models, systematic experiments were conducted to elaborate on the effect of structural parameters on the drop lateral transport distance (see details in Supplementary Discussion 6). Specifically, substrates with different ridge width W (20 and 40 μm), groove width S (10, 40, and 80 μm), and ridge height H (10, 20, 50 μm) were fabricated, and labeled as W20S10H20 (i.e., W = 20 μm, S = 10 μm, H = 20 μm), W20S40H20, W20S80H20, W40S10H20, W40S40H20, W40S80H20, W40S40H10, and W40S40H50, respectively (Supplementary Fig. 9). It is obvious that the proposed approach is well validated and the lateral displacement of impacting drops on different substrates follows the same scaling laws, e.g., ΔlL ~ lc2r−1We1/2 in the transition boiling state and ΔlR ~ lc2H(W + S)−1r−1We1/2 in the film boiling, respectively, as shown in Supplementary Fig. 10.

The ability to deposit drops to obtain desired transport directions in a highly repeatable manner is desired for thermal dissipation with high efficiency, and can also be utilized to perform practical tasks such as drop sieving for fluidic-based applications29,30. The advantages of our work will largely benefit these concepts, which will be illustrated as follows. It is well-known that the boiling state of the liquid as a function of the substrate temperature is related to the liquid property31,32, and Fig. 5a shows a phase diagram describing the bidirectional transport of various liquids (Table S1) at different temperatures. First, for a specific value of the substrate temperature, we could use a liquid to guarantee that the impinging drop is always in the transition boiling state. As shown in Fig. 5b, when successive water drops are deposited, they can be converged in the center of curvature (Supplementary Movie 3), which is promising to achieve a high-efficient heat transfer at the center of the substrate, especially in cooling hot spots. Second, as shown in Supplementary Fig. 12a and Supplementary Movie 4, by leveraging the opposite rebounding direction of two kinds of droplets with different liquid properties (for example, ethanol and n-hexane) at a certain temperature, e.g., T = 160 °C, the synchronously deposited drops rebound to opposite directions to realize drop sieving (Fig. 5c and Supplementary Movie 5). These intriguing transport behaviors on our surfaces are general, and more results are given in Supplementary Fig. 12b–d and Supplementary Movies 6–8.

a Phase diagram for different transport directions of various liquids at different temperatures. "Left" represents the drop rebounding towards the center of curvature, while "Right" represents the drop rebounding far away from the center of curvature. b Sequential images showing the convergence of successive water drops on the center of curvature at T = 250 °C and We = 10.5. c Comparison between opposite rebounding directions of ethanol and n-hexane drops showing the ability of drop sieving at T = 160 °C, We = 10.5, and r = d/4. The ethanol drop bounces towards the center of curvature, while the n-hexane drop bounces far away from the center of curvature. Scale bars in b, c are 2 mm.

In summary, we have reported a strategy to steer drops on hot concentric microgrooved substrates. In the transition boiling, the impinging drop forms a bean shape during the retraction process, which generates a Laplace pressure difference, driving the drop towards the center of curvature. In the film boiling state, however, the impacting drop rebounds towards the opposite direction, i.e., far away from the center of curvature, owning to the interaction force between the penetrated liquid and the grooves. This work breaks through our understanding that drops on hot textured surfaces could only be transported in one direction, such as the Leidenfrost drops. Recently, Jiang et al. report a remarkable phenomenon that the Leidenfrost effect could be inhibited up to a record-high temperature33, e.g., 1100 °C, which is rationalized by constructing the termed "structured thermal armor" to decouple the wetting phenomena from the vapor dynamics. Here, we propose an alternative approach that emphasizes the richness of thermal management. This steerable and highly repeatable manner to realize directional transportation of drops, as well as detailed fundamental understandings of the physics that we have revealed, provides a paradigm to manipulate drop transportation on hot surfaces available in various applications, including spray cooling and drop sieving, among many others.

Still, some open questions remain. For example, the contact state due to boiling modes between the liquid and the substrate during the impingement plays a crucial role in dictating the drop rebound dynamics, however, direct observation of the contact state on our structured surface challenges the current technology20, but which is essential for a scientific breakthrough. In the present experiment, hydrophilic substrates are employed. Impingement on substrates with different wettabilities under even higher temperatures and Weber number remain to be further explored, which could be of considerable interest for practical applications. Moreover, the vertical component of the pressure resulting from the thin vapor film may affect the overall drop bouncing dynamics, and further investigation is still needed. Therefore, a thorough clarification of these influences accounting for the steerable drop transport would deserve a dedicated study in the future.

The concentric microgroove arrays were fabricated using the standard photolithography technique on a silicon wafer. Specifically, a SiO2 layer with a thickness of 2 μm was first deposited on the silicon wafer with a thickness of 500 μm at high temperature. The photoresist was uniformly coated on the SiO2 layer. Then the uncovered SiO2 was etched by Reactive Ion Etching (RIE). Deep RIE was used to further etch the silicon substrate to form the micropattern. A series of substrates were fabricated with different ridge width W (20 and 40 μm), groove width S (10, 40, and 80 μm), and ridge height H (10, 20, 50 μm), and are labeled as W20S10H20 (i.e., W = 20 μm, S = 10 μm, H = 20 μm), W20S40H20, W20S80H20, W40S10H20, W40S40H20, W40S80H20, W40S40H10, and W40S40H50, respectively, as shown in Supplementary Figs. 1, 9.

The sample was heated on a hot plate, in which a thermocouple was used to monitor the surface temperature. Deionized water droplets of radius R = 1.18 mm were released from a syringe pump equipped with a fine needle at a flow rate of 2 μL/s. The height between the droplets and the substrates was adjusted to change the impact velocity U from 0.23 to 1.02 m/s, corresponding to the Weber number (We) of 0.8 < We < 16.9, where We = ρU2R/γ is defined as the Weber number with the liquid density ρ = 998 kg/m3 and the liquid surface tension γ = 72 mN/m. The drop dynamic behavior is filmed simultaneously from the side and top views using two synchronous high-speed cameras (Photron SA5) at a frame rate of 10,000 fps.

No statistical method was used to predetermine sample size. No data were excluded from the analyses. The experiments were not randomized. The Investigators were not blinded to allocation during experiments and outcome assessment.

The data that support the findings of this study are available from the corresponding authors upon request. Source data are provided with this paper.

Li, J. et al. Directional transport of high-temperature Janus droplets mediated by structural topography. Nat. Phys. 12, 606–612 (2016).

Article ADS CAS Google Scholar

Wendelstorf, J., Spitzer, K. H. & Wendelstorf, R. Spray water cooling heat transfer at high temperatures and liquid mass fluxes. Int. J. Heat. Mass. Transf. 51, 4902–4910 (2008).

Article CAS Google Scholar

Vakarelski, I. U., Marston, J. O., Chan, D. Y. C. & Thoroddsen, S. T. Drag reduction by Leidenfrost vapor layers. Phys. Rev. Lett. 106, 214501 (2011).

Article ADS Google Scholar

Quéré, D. Leidenfrost dynamics. Annu. Rev. Fluid Mech. 45, 197–215 (2013).

Article ADS MathSciNet Google Scholar

Vakarelski, I. U., Patankar, N. A., Marston, J. O., Chan, D. Y. C. & Thoroddsen, S. T. Stabilization of Leidenfrost vapour layer by textured superhydrophobic surfaces. Nature 489, 274–277 (2012).

Article ADS CAS Google Scholar

Wells, G. G., Ledesma-Aguilar, R., McHale, G. & Sefiane, K. A sublimation heat engine. Nat. Commun. 6, 6390 (2015).

Article ADS CAS Google Scholar

Agrawal, P. et al. Leidenfrost heat engine: Sustained rotation of levitating rotors on turbine-inspired substrates. Appl. Energ. 240, 399–408 (2019).

Article Google Scholar

Linke, H. et al. Self-propelled Leidenfrost droplets. Phys. Rev. Lett. 96, 154502 (2006).

Article ADS CAS Google Scholar

Dupeux, G. et al. Viscous mechanism for Leidenfrost propulsion on a ratchet. Europhys. Lett. 96, 58001 (2011).

Article ADS Google Scholar

Lagubeau, G., Le Merrer, M., Clanet, C. & Quéré, D. Leidenfrost on a ratchet. Nat. Phys. 7, 395–398 (2011).

Article CAS Google Scholar

Thimbleby, H. The Leidenfrost phenomenon. Phys. Educ. 23, 300–303 (1989).

Article ADS Google Scholar

de Ruiter, J., Soto, D. & Varanasi, K. K. Self-peeling of impacting droplets. Nat. Phys. 14, 35–39 (2018).

Article Google Scholar

Liu, M. et al. Inhibiting random droplet motion on hot surfaces by engineering symmetry-breaking Janus-mushroom structure. Adv. Mater. 32, 1907999 (2020).

Article CAS Google Scholar

Chaudhury, M. K. & Whitesides, G. M. How to make water run uphill. Science 256, 1539–1541 (1992).

Article ADS CAS Google Scholar

Ichimura, K., Oh, S. K. & Nakagawa, M. Light-driven motion of liquids on a photoresponsive surface. Science 288, 1624–1626 (2000).

Article ADS CAS Google Scholar

Lv, C. et al. Substrate curvature gradient drives rapid droplet motion. Phys. Rev. Lett. 113, 026101 (2014).

Article ADS Google Scholar

Ju, J. et al. A multi-structural and multi-functional integrated fog collection system in cactus. Nat. Commun. 3, 1247 (2012).

Article ADS Google Scholar

Liu, Y. et al. Pancake bouncing on superhydrophobic surfaces. Nat. Phys. 10, 515–519 (2014).

Article CAS Google Scholar

Zhan, H. et al. Horizontal motion of a superhydrophobic substrate affects the drop bouncing dynamics. Phys. Rev. Lett. 126, 234503 (2021).

Article ADS CAS Google Scholar

Tran, T., Staat, H. J. J., Prosperetti, A., Sun, C. & Lohse, D. Drop impact on superheated surfaces. Phys. Rev. Lett. 108, 036101 (2012).

Article ADS Google Scholar

Taylor, G. The dynamics of thin sheets of fluid. III. Disintegration of fluid sheets. Proc. R. Soc. Lond. A 253, 313–321 (1959).

Article ADS MathSciNet Google Scholar

Culick, F. E. C. Comments on a ruptured soap film. J. Appl. Phys. 31, 1128–1129 (1960).

Article ADS Google Scholar

Richard, D., Clanet, C. & Quéré, D. Contact time of a bouncing drop. Nature 417, 811–811 (2002).

Article ADS CAS Google Scholar

Clanet, C., Beguin, C., Richard, D. & Quéré, D. Maximal deformation of an impacting drop. J. Fluid Mech. 517, 199–208 (2004).

Article ADS Google Scholar

Quéré, D. Wetting and roughness. Ann. Rev. Mater. Res. 38, 71–99 (2008).

Article ADS Google Scholar

Reyssat, M., Pardo, F. & Quéré, D. Drops onto gradients of texture. Europhys. Lett. 87, 36603 (2009).

Article Google Scholar

Soto, D., Lagubeau, G., Clanet, C. & Quéré, D. Surfing on a herringbone. Phys. Rev. Fluids 1, 013902 (2016).

Article ADS Google Scholar

Mouterde, T. et al. Two recipes for repelling hot water. Nat. Commun. 10, 1410 (2019).

Article ADS Google Scholar

Feng, S. et al. Three-dimensional capillary ratchet-induced liquid directional steering. Science 373, 1344–1348 (2021).

Article ADS CAS Google Scholar

Mazaltarim, A. J., Bowen, J. J., Taylor, J. M. & Morin, S. A. Dynamic manipulation of droplets using mechanically tunable microtextured chemical gradients. Nat. Commun. 12, 3114 (2021).

Article ADS CAS Google Scholar

Gottfried, B. S., Lee, C. J. & Bell, K. J. The leidenfrost phenomenon: film boiling of liquid droplets on a flat plate. Int. J. Heat. Mass. Transf. 9, 1167–1187 (1966).

Article CAS Google Scholar

Biance, A. L., Clanet, C. & Quéré, D. Leidenfrost drops. Phys. Fluids 15, 1632–1637 (2003).

Article ADS CAS Google Scholar

Jiang, M. et al. Inhibiting the Leidenfrost effect above 1,000 °C for sustained thermal cooling Leidenfrost drops. Nature 601, 568–572 (2022).

Article ADS CAS Google Scholar

Download references

The work is supported by the National Natural Science Foundation of China (NSFC, No. 52075071, Y.L., 11872227 and 11921002, C.L. (Cunjing Lv)), Fundamental Research Funds for the Central Universities (DUT18RC3048, Y.L.), and Liao Ning Revitalization Talents Program (XLYC1807092, Y.L.).

Key Laboratory for Precision and Non-Traditional Machining Technology of Ministry of Education, Dalian University of Technology, Dalian, 116024, P. R. China

Cong Liu, Chenguang Lu, Zichao Yuan & Yahua Liu

Department of Engineering Mechanics and Center for Nano and Micro Mechanics, AML, Tsinghua University, Beijing, 100084, P. R. China

Cunjing Lv

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

Y.L., C.L. (Cong Liu), and C.L. (Cunjing Lv) conceived the research. Y.L. and C.L. (Cunjing Lv) supervised the research. C.L. (Cong Liu), C.L. (Chenguang Lu), and Z.Y. carried out the experiment. C.L. (Cunjing Lv) built the models. All authors analyzed the data and wrote the paper.

Correspondence to Cunjing Lv or Yahua Liu.

The authors declare no competing interests.

Nature Communications thanks Rajneesh Bhardwaj, Li Jia, and Zuankai Wang for their contribution to the peer review of this work.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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/.

Reprints and Permissions

Liu, C., Lu, C., Yuan, Z. et al. Steerable drops on heated concentric microgroove arrays. Nat Commun 13, 3141 (2022). https://doi.org/10.1038/s41467-022-30837-z

Download citation

Received: 02 November 2021

Accepted: 20 May 2022

Published: 06 June 2022

DOI: https://doi.org/10.1038/s41467-022-30837-z

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.