Optimal free-surface pumping by an undulating carpet | Nature Communications

Nature Communications volume 14, Article number: 7735 (2023) Cite this article

3855 Accesses

3 Citations

191 Altmetric

Metrics details

Examples of fluid flows driven by undulating boundaries are found in nature across many different length scales. Even though different driving mechanisms have evolved in distinct environments, they perform essentially the same function: directional transport of liquid. Nature-inspired strategies have been adopted in engineered devices to manipulate and direct flow. Here, we demonstrate how an undulating boundary generates large-scale pumping of a thin liquid near the liquid-air interface. Two dimensional traveling waves on the undulator, a canonical strategy to transport fluid at low Reynolds numbers, surprisingly lead to flow rates that depend non-monotonically on the wave speed. Through an asymptotic analysis of the thin-film equations that account for gravity and surface tension, we predict the observed optimal speed that maximizes pumping. Our findings reveal how proximity to free surfaces, which ensure lower energy dissipation, can be leveraged to achieve directional transport of liquids.

The necessity to manipulate flow and transport liquids is primitive to many biophysical processes such as embryonic growth and development1,2, mucus transport in bronchial tree3,4,5, the motion of food within intestine6,7, animal drinking8,9. Engineered systems also rely on efficient liquid transport such as in heat sinks and exchangers for integrated circuits10,11, micropumps12,13, and lab-on-a-chip devices14. Transporting liquids at small scales requires non-reciprocal motion to overcome the time reversibility of low Reynolds number flows. Deformable boundaries in the form of rhythmic undulation of cilia beds and peristaltic waves are nature’s resolutions to overrule this reversibility and achieve directional liquid transport. While peristaltic pumps have become an integral component of biomedical devices, artificial ciliary metasurfaces that can actuate, pump, and mix flow have been realized only recently15,16,17,18,19.

The design strategy of valveless micropumps essentially relies on a similar working principle as cilia-lined walls; sequential actuation of a channel wall by electrical or magnetic fields creates a traveling wave that drags the liquid along with it20,21. While the primary focus of micropumps has been on the transport of liquids enclosed within a channel, numerous technological applications require handling liquids near fluid–fluid interfaces. In particular, processes such as self-assembly, encapsulation, and emulsification involving micron-sized particles critically rely on the liquid flow near interfaces22,23. Liquid–air interfaces also act as a cradle to Neuston, organisms that inhabit at and below the water surface. For example, the underwater apple snail Pomacea canaliculata exploits the water surface to drive a large-scale surface flow and fetch floating food particles from afar in a process called pedal surface collection24,25,26,27. Whirligig beetles and water striders have developed morphological features to manipulate flow around them for effective locomotion at the interface28,29. Understanding the physics behind these natural phenomena could open up new bio-inspired strategies for flow actuation and sensing at interfaces.

Here, we reveal how a rhythmically deforming solid boundary pumps viscous liquid at the interface, and transports floating objects from distances much larger than its size. Our design, inspired in part by the capability of underwater snails in creating flow through undulations on their flexible foot26,27, produces traveling waves on an undulator. Even though traveling boundaries are a canonical strategy to drive flow within enclosed spaces when placed near a liquid–air interface, the undulator gives rise to non-intuitive observations; pumping does not increase proportionally to the wave speed, and we observe non-monotonicity in the average motion of surface floaters as the wave speed is gradually increased. Detailed measurements of the velocity field in combination with an analysis of the lubrication theory unravel the interfacial hydrodynamics of the problem that emerges from a coupling between capillary, gravity, and viscous forces. We find that the non-monotonic flow is a direct consequence of whether the interface remains flat or conforms to the phase of the undulator. Through the theoretical analysis, we are able to predict the optimal wave speed that maximizes pumping, and this prediction is in excellent agreement with experiments. Finally, we show how pumping near an interface is a less dissipative strategy to transport liquid compared to pumping near a rigid boundary.

A 3D-printed undulator capable of generating traveling waves is attached to the bottom of an acrylic tank. The tank is filled with a viscous liquid (silicone oil or glycerin–water mixture) such that the mean depth of liquid above the undulator (H) remains much smaller the undulator wavelength (λ), i.e., H/λ ≪ 1. The undulator is driven by a servo motor attached to a DC power source. Millimetric styrofoam spheres are sprinkled on the liquid surface and their motion is tracked during the experiment to estimate the large-scale flow of liquid. Additionally, we characterize the flow within the thin film of liquid directly in contact with the undulator by performing 2D particle image velocimetry (PIV) measurements. Our experimental design is essentially a mesoscale realization of the Taylor’s sheet30 placed near a free surface31,32; the crucial difference, however, is that the sheet or undulator is held stationary here, in contrast to free swimming.

Images of the undulator are shown in Fig. 1a, b. The primary component of this design is a helical spine encased with a series of hollow, rectangular links that are interconnected through a thin top surface33 (see SI and Supplementary Movies 1 and 2 for details). The links along with the top surface form an outer shell that transforms the helix rotation into a planar traveling wave of the form, \(\delta \sin [2\pi (x-{V}_{w}t)/\lambda ]\). The pitch and radius of the helix determine the wavelength λ and amplitude δ of the undulations respectively. By modulating the angular frequency of the helix, we are able to vary the wave speed Vw = ωλ from 15 to 120 mm/s (λ is fixed at 50 mm). We perform experiments with undulators of length λ and 2λ, and the results remain invariant of the undulator size. For given Vw, shapes of the undulator surface are shown in Fig. 1c for one period of oscillation.

The actuator, shown in panels a and b, is comprised of a helix rotating inside a blue shell. Rotation of the helix causes an oscillatory motion of the shell forming a traveling wave on the surface. It is placed at a mean depth H below the liquid surface. c Shape of the undulations over a period of oscillation. These shapes are captured by a traveling sine wave of \(\delta \sin [2\pi(x-{V}_{w}t)/\lambda]\). d Trajectories representing motion of styrofoam particles at the interface due to 30 min of continuous oscillation of the undulator in silicone oil (viscosity 0.97 Pa ⋅ s) at a constant Vw. This panel is a top-view image with the actuator position marked at the bottom of the frame. The color coding of dark to light indicates the arrow of time. e Magnified trajectories of particles located straight ahead of the actuator. The filled circles represent the initial positions of the styrofoam particles. f Particle velocity as a function of distance for increasing wave speeds (Vw). Different wave speeds are marked by the color coding. Distances are measured from the edge of the actuator, as shown in panel e. Each of the curves is an average of over 20 trajectories. Particle velocity exhibits a non-monotonic behavior with Vw with maximum velocities measured at intermediate wave speeds. The inset confirms this behavior by showing particle velocity at a fixed location, x = 50 mm for different Vw. Error bars in this plot represent standard deviation in velocity magnitude. The gray line is the prediction from eq. (8).

Figure 1d shows the trajectories of floating styrofoam particles generated by 30 min of continuous oscillations in 1000 cSt silicone oil contained in an acrylic tank of dimensions 61 × 46 cm (Supplementary Movie 3 shows the motion of surface floaters for different Vw). Traveling waves on the actuator move in the downward direction as shown by the direction of Vw in Fig. 1d. The large-scale flow caused by these undulations drags the floaters. The color code on the trajectories represents the arrow of time: blue and yellow colors represent the initial and final positions, respectively. Placing the undulator near a side wall of the tank, we measure the floaters’ motion over a decade in distance. We quantify pumping from the net flow that the undulator creates near the free surface. Unlike closed pipes and channels, this net flow is not driven by an imposed pressure gradient. The shape of the free surface, which is an unknown a priori, determines the local pressure inside the liquid film. The undulator’s open ends and proximity to the tank walls result in flow re-circulation observed in Fig. 1d.

Typical Stokes number, St = ρpRpVw/η remains very small, ≃ 10−2 for the styrofoam floaters (based on wave speed Vw = 100 mm/s, particle radius Rp = 1 mm, particle density ρp ≃ 50 kg/m3 and silicone oil viscosity ηs = 0.97 Pa ⋅ s). Thus we utilize the trajectories of these particles to estimate fluid motion at the interface. The particle tracks in Fig. 1d reveal that some particles are recirculated back in the flow due to the nearby tank wall, which we ignore in the analysis, and only focus on the floaters whose initial positions are straight ahead of the undulator. These trajectories are shown in Fig. 1e with black circles representing the initial positions. For a given Vw, we interpolate 20 trajectories to construct a velocity-distance curve which is shown in Fig. 1f (see SI for details of these measurements). Here, \(| \bar{V}|={({V}_{x}^{2}+{V}_{y})}^{1/2}\) is the magnitude of the velocity at the liquid–air interface, and x is the distance from the edge of the actuator. At the undulator edge, a sharp change in the boundary gives rise to a transition region in the flow field where surface floaters exhibit non-uniform and unsteady motion (cf. Fig. S2b in SI). Thus, we disregard the first 20 mm of data to avoid edge effects. The color code on the curves in Fig. 1f represents the magnitude of Vw. Interestingly, we observe a non-monotonic response in the surface fluid velocity; for any distance ( − 20 ≤ x ≤ 120) speed of fluid parcels initially increases with Vw, but subsequently drops down with a further increase in Vw. Once \(| \bar{V}|\) at a fixed location (x = − 50 mm) is plotted against Vw (inset of Fig. 1f), it becomes apparent that the maximum surface flow is achieved for an intermediate speed, Vw ≃ 80 mm/s. Since the overall flow in the liquid is driven by the hydrodynamics within the thin film of liquid atop the undulator, we focus on quantifying the velocity field and flow rate in this region.

Before we discuss the experimental results further, it is instructive to identify the relevant dimensionless groups that dictate the response of the system. The system has eight-dimensional parameters: three length scales given by film thickness H, amplitude (δ) and wavelength (λ) of the undulator, the velocity scale Vw, gravitational constant g, and three fluid properties set by surface tension (γ), density (ρ), dynamic viscosity (η). These parameters lead to 5 dimensionless groups, namely, ϵ = δ/H, a = H/λ, Reynolds number Re = ρVwλa2/η, Capillary number Ca = ηVw/(γa3), and Bond number Bo = ρgλ2/γ. Here both Re and Ca are defined for the thin-film limit, a ≪ 1. We choose two working liquids, silicone oil (ηs = 0.97 Pa ⋅ s, γs = 0.021 N/m) and glycerin–water mixture (85% glycerin and 15% water by volume), GW (ηGW = 0.133 Pa ⋅ s, γGW = 0.067 N/m). For each of the liquids, the thickness, H (maintaining a ≪ 1), and wave speed, Vw are varied independently. Across all experiments, Re remains lower than 1 (0.01 − 0.45). Thus inertial effects are subdominant and the problem is fully described by ϵ, Ca, and Bo. We vary Ca, the ratio of viscous to capillary forces, over three orders in magnitude, 7 − 6049. The value of Bo, representing the strength of gravitational forces to surface tension, is 1133 and 426 for silicone oil and glycerin–water mixture respectively. For a full list of experimental parameters and dimensionless numbers, see Table S1 in SI. As we will demonstrate in the next sections, Ca/Bo, which represents the ratio of viscous force to gravitational force, turns out to be the key governing parameter.

The flow field within the thin film of liquid above the undulator is characterized by performing PIV at five longitudinal planes along the width of the undulator (see “Methods” section and SI section III for details). We find that the velocity field remains invariant of the PIV plane and exhibits low magnitudes of divergence, signifying a dominant 2D flow in the thin film. Figure 2a shows a long-exposure image of illuminated tracer particles, giving a qualitative picture of the flow. The particles essentially oscillate up-down with the actuator, but exhibit net horizontal displacement over a period due to the traveling wave. The presence of an interface is also crucial to the transport mechanism; the interfacial curvature induces a capillary pressure that modifies the local flow field. The coupling between the two deforming boundaries determines the flow within the gap. Snapshots of typical velocity fields for the two liquids are shown in Fig. 2b. The top panel is a silicone oil flow-field with Vw = 23 mm/s and H = 5.7 mm (Ca = 714, Bo = 1133), while the bottom panel represents flow-field of glycerin–water mixture with Vw = 17 mm/s and H = 6.3 mm (Ca = 17, Bo = 426). For both liquids, the instantaneous profile of Vx remains half-parabolic along the depth of the thin film, which signifies a shear-free liquid–air interface. Higher Ca leads to larger deformation of the free surface. Colors in the plot represent the magnitude of the horizontal velocity component, Vx; a portion of the liquid that follows the wave is shown in red, whereas a blue region represents part of the liquid that moves in the opposite direction to the wave. In fact, the velocity vectors at a given location switch directions depending on the phase of the actuator (see Supplementary Movies 4 and 5). Thus, to estimate the net horizontal transport of liquid across a section, we first integrate Vx across thickness in the middle of the undulator (marked by the black dashed line in Fig. 2b, i), which yields an instantaneous flow rate

Here, ha and hf are the positions of the bottom and top boundaries from the reference point, respectively. Figure 3a plots Q as a function of time, measured in silicone oil for three distinct wave speeds. It shows that Q oscillates with the same time period as the undulator (τ = λ/Vw), but there is a net flow of liquid along the traveling wave. Thus a time-averaged flow rate,

gives a measure of liquid transport by the undulator.

a A sketch of the actuator and a long-exposure image of a typical flow-field measurement, showing the motion of the tracer particles in the thin film. The free surface deforms in response to the flow. b Results of PIV for two different capillary numbers, Ca = 714 (top panel), Ca = 17 (bottom panel). In both these panels the bottom boundary is the actuator surface, while the top boundary is the liquid interface. The color coding represents the horizontal component of the velocity field, Vx; red signifies flow along the wave speed (Vw) while blue signifies flow opposite to Vw.

a Instantaneous flow rate in silicone oil over multiple periods of oscillation. The data sets represent increasing Vw, from white to black. These measurements are taken at a cross-section marked by the dashed line in Fig. 2b, i. b Time-averaged flow rate, \(\left\langle Q\right\rangle\) is plotted against the flux scale ϵ2VwH of the problem. The circles represent the experiment in silicone oil with bigger markers denoting larger height, H: 4.3 mm (red), 5.7 mm (green), 6.8 mm (blue), 8 mm (orange). The squares represent gylcerin-water (GW) experiments with H = 6.3 mm. Error bars are based on standard deviation. The dashed line is the theoretical prediction, given by \(\left\langle Q\right\rangle=3{\epsilon }^{2}{V}_{w}H/2\).

Figure 3b gives a comprehensive picture of the flow rate measured across all the experiments. \(\left\langle Q\right\rangle\) is plotted against the characteristic flow rate VwH. The geometric prefactor of ϵ2 is a direct consequence of the thin geometry of the flow34. Two interesting observations are in order; (i) regardless of the fluid properties, the flow rates at first increase linearly with ϵ2VwH, (ii) all the data sets other than the GW exhibit a non-monotonic behavior with flow rates reaching maximum values at intermediate ϵ2VwH. Note that Re in GW experiments tends to reach \({{{{{{{\mathscr{O}}}}}}}}(1)\) at the higher speeds. Thus, it is expected to behave differently from the silicone oil experiments where Re ≤ 10−1. Additional data presented in Fig. S5 of SI for experiments in GW (performed at higher ϵ2VwH than what is shown in Fig. 3b) and pure water show that inertial effects lead to a continuous increase of \(\left\langle Q\right\rangle\) with wave speeds, and confirm that the non-monotonic flow rates are a feature of the low Re regime. Importantly, the non-monotonic surface flow observed in Fig. 1f seems to be directly correlated to the flow within the thin film above the undulator. We would like to point out that these measurements remain invariant of the undulator size as shown in the SI (cf. Fig. S4) where we compare the time-averaged flow rates measured in single and double-wave undulators. In the next section, we develop a theoretical model to explain how the geometrical and material parameters combine to give the optimal wave speed that maximizes the flow rate.

We consider the two-dimensional geometry depicted in Fig. 4a for the theoretical model. An infinite train of periodic undulations of the form \({h}_{a}=\delta \sin 2\pi (x-{V}_{w}t)/\lambda\) propagates on the actuator located at a mean depth of H from the free surface. We analyze the flow in the thin-film limit, such that a = H/λ ≪ 1. A key aspect of the problem is that the shape of the interface, hf is unknown along with the flow-field. The explicit time dependence in this problem is a direct manifestation of the traveling wave on the boundary. Thus in a coordinate system (X,  Z) moving with the wave, the flow becomes steady. A simple Galilean transformation relates these coordinates to the laboratory coordinates (x,  z): X = x − Vwt, and Z = z. Thus, we first solve the problem in the wave frame and then transform the solution to the lab frame. Leaving the details of the derivation in Materials & Methods, here we present the key results.

a The thin-film geometry with relevant quantities. We consider an infinite train of traveling undulations of amplitude δ and wavelength λ moving at a speed of Vw. The coordinate frame (X, Z) travels with the undulations. The red curve represents the bottom boundary in motion. A liquid layer of mean thickness H reside on top the deformable bottom boundary. Shape of the free surface is given by hf, while the bottom surface is given by ha. b Numerical solution of (3) and (4), plotted in terms of the time-averaged flow rate \(\left\langle \bar{Q}\right\rangle=\bar{q}+1\) as a function of Capillary number (Ca), for different Bond numbers (Bo). The two large Bond numbers correspond to the experimental values. c The rescaled experimental data of Fig. 3b are in excellent agreement with the theoretical prediction of (7), plotted as the solid black line. Error bars are based on standard deviation. The small Vw (Ca/Bo ≪ 1) limit is given by \(\left\langle \bar{Q}\right\rangle=3{\epsilon }^{2}/2\), while the large Vw (Ca/Bo ≫ 1) limit is given by \(\left\langle \bar{Q}\right\rangle=2{\pi }^{2}{\epsilon }^{2}{({{{{{{{\rm{Ca}}}}}}}}/{{{{{{{\rm{Bo}}}}}}}})}^{-2}/3\).

In the thin-film limit, the separation of vertical and horizontal scales leads to a predominantly horizontal flow-field, and both mass and momentum conservation equations are integrated across the film thickness to reach an ordinary differential equation involving the free-surface shape, hf, and volume flow rate q. Introducing dimensionless variables \(\bar{X}=X/\lambda,{\bar{h}}_{f}={h}_{f}/H,{\bar{h}}_{a}={h}_{a}/H\), and \(\bar{q}=q/{V}_{w}H\) we get

where both \(\bar{q}\) and \({\bar{h}}_{f}\) are unknowns, and \({\bar{h}}_{a}=\epsilon \sin 2\pi \bar{X}\) is known. We close the problem by imposing the following constraint on \({\bar{h}}_{f}\):

which states that the mean film thickness over one wavelength does not change due to deformation. Along with periodic boundary conditions, equations (3) and (4) form a set of nonlinear coupled equations whose solutions depend on the three parameters, Ca, Bo, and ϵ. For chosen Bo and ϵ, these equations are solved by a shooting method for a wide range of Ca. To be able to compare the numerical results with the experimental data of Fig. 3, we transform the results to the lab frame using the relation \(\bar{Q}=\bar{q}+({\bar{h}}_{f}(\bar{x},\bar{t})-{\bar{h}}_{a}(\bar{x},\bar{t}))\). Owing to the periodic nature of \({\bar{h}}_{f}\) and \({\bar{h}}_{a}\), the time-averaged flow rate simplifies to \(\left\langle \bar{Q}\right\rangle=\bar{q}+1\). Figure 4b shows the numerical solution of \(\left\langle \bar{Q}\right\rangle\) as a function of Ca for ϵ = 0.3 and different Bo. All curves exhibit the same qualitative behavior; at low Ca, the scaled flow rate reaches a constant value as \(\left\langle Q\right\rangle \sim {V}_{w}H\), which is analogous to what we observe in Fig. 3b. At large Ca, however, we recover a decreasing flow rate as \(\left\langle Q\right\rangle \sim {({V}_{w}H)}^{-\alpha }\) with α > 0. The transition between the two regimes scales with the Bo. Thus the thin-film equation captures the qualitative behavior found in the experiments. Interestingly, the non-monotonic response is also found in the Bo = 0 limit (light blue curve in Fig. 4b), broadening the significance of our results to smaller scales where capillary forces dominate over gravity.

For Bo ≫ 1, the third-order term in Eq. (3) can be neglected, which simplifies the governing equation to

Indeed Bo values in experiments are large (433 and 1132) justifying the above simplification. Furthermore, we assume that the amplitude of the wave, δ is much smaller than H, ϵ ≪ 1. Interestingly Ca/Bo, the single parameter dictating the solution of eq. (5), does not contain surface tension. This ratio is reciprocal to the Galileo number which plays a crucial role in the stability of thin films driven by gravity35. Here we look for asymptotic solutions of the form, \({\bar{h}}_{f}=1+\epsilon {\bar{h}}_{f1}+{\epsilon }^{2}{\bar{h}}_{f2}+{{{{{{{\mathscr{O}}}}}}}}({\epsilon }^{3})\) and \(\bar{q}={q}_{0}+\epsilon {\bar{q}}_{1}+{\epsilon }^{2}{\bar{q}}_{2}+{{{{{{{\mathscr{O}}}}}}}}({\epsilon }^{3})\)36. We insert these expansions in eqns. (5) and (4), and solve the equations in orders of ϵ. Leaving the solution of \({\bar{h}}_{f}\) in the SI, here we present the solution of that becomes

Thus the time-averaged flow rate in the lab frame is given by

This is the key result of the theoretical model. It demonstrates that the flow rate is quadratic in the amplitude of the traveling wave, which is why we incorporated ϵ2 in the horizontal scale of Fig. 3b. Importantly, eq. (7) captures the non-monotonic behavior of the experiments. Once the data in Fig. 3b are rescaled, all the different cases collapse onto a master curve which is in excellent, overall agreement with the black solid line representing eq. (7), as shown in Fig. 4c. The difference between theory and experiment at large Ca/Bo may be attributed to large deformations of the interface where the lubrication approximation of δ/H ≪ 1 breaks down.

The physical picture behind the non-monotonic nature of the flow rate becomes clear once the free-surface shapes are found. For a given Bo with a low Ca, the liquid–air interface behaves as an infinitely taut membrane with minimal deformations. Thus, a liquid parcel moves primarily in the horizontal direction, and the flow rate is given purely by the kinematics (Vw, H, δ). Indeed for Ca/Bo ≪ 1 eq. (7) simplifies to give \(\left\langle \bar{Q}\right\rangle /{\epsilon }^{2}=3/2\). In the dimensional form, this relation explains the increase in the flow rate with the wave speed, \(\left\langle Q\right\rangle=3{\delta }^{2}{V}_{w}/2H\). Thus, the flow rate is independent of the liquid properties which we have noted in Fig. 3b. As Ca increases, the interface starts to deform up and down by conforming with the undulating actuator. In this limit, the translational velocity of tracer particles decreases thereby lowering the flow rate. Indeed, in the limit of Ca/Bo ≫ 1, we find a decreasing flow rate given by \(\left\langle \bar{Q}\right\rangle=2{\pi }^{2}{\epsilon }^{2}{(Ca/Bo)}^{-2}/3\). These two asymptotic limits are shown as dashed lines in Fig. 4c. The flow rate attains a maximum at the intersection of these two lines where Ca/Bo = 2π/3. In the dimensional form, this particular value of Ca/Bo gives the optimal wave speed at which flow rate peaks,

The optimal wave speed emerges from a competition between hydrostatic pressure (~ρgH) and lubrication pressure (~ηVwλ/H2); surface tension drops out in the above expression. Eq. (8) gives the optimal speed at which the undulator maximizes pumping. Now we are in a position to examine whether eq. (8) captures the peak surface velocities observed in Fig. 1f. Plugging in the density (ρ = 970 kg/m3), viscosity (ηs = 0.97 Pa ⋅ s), H = 5.7 mm, λ = 50 mm, we find \({V}_{w}^{{{{{{{{\rm{(max)}}}}}}}}}=76.1\) mm/s, which matches very well with the observation (shown as the gray line in Fig. 1f inset). Note that the blue and orange data sets in Fig. 4c exhibit optimal flow rates at lower Ca/Bo than the above asymptotic prediction. This difference could be attributed to higher H/λ where the thin-film equation can not capture the exact flow-field.

The flow rate achieved by this mechanism comes at the expense of the power needed to drive the undulator. This power expenditure equals the viscous dissipation within the flow. To this end, we estimate the efficacy of the mechanism by comparing the output, \(\left\langle \bar{Q}\right\rangle\) to the input, viscous dissipation \(\bar{{{{{{{{\mathscr{E}}}}}}}}}\) (see “Methods” for a derivation of \(\bar{{{{{{{{\mathscr{E}}}}}}}}}\)). To demonstrate the benefit of having an interface on the pumping capability of this mechanism, we compare it with the flux and dissipation for a rigid top boundary. These results are shown in Fig. 5. The data points represent dimensionless flux plotted against dissipation for a wide range of Ca/Bo. The ϵ ≪ 1 asymptotic result, shown as the black dashed line in Fig. 5, captures these results perfectly, giving the following algebraic relation between the two

The importance of the free surface becomes apparent when the above result is compared to the scenario of the thin film bounded by a rigid, solid wall on top. As shown in the SI, for a rigid top boundary, both the flow rate and dissipation is given purely by the ratio of δ/H. We find that the flow dissipates 4 times more energy to achieve the same amount of flow,

This is plotted as the solid black line in Fig. 5. Thus it is clear that the liquid–air interface facilitates pumping by promoting horizontal transport of fluid parcels at a lower power consumption.

The dimensionless flow rate is plotted against the dimensionless dissipation for a wide range of Ca/Bo values. The data points represent numerical results, which are obtained for a fixed amplitude-to-depth ratio of ϵ = 0.3. The dashed line is the asymptotic prediction of eq. (9). The solid line is the result of a top rigid boundary and represents eq. (10).

In summary, we have demonstrated the pumping capability of a sub-surface undulating carpet; the traveling wave triggers a large-scale flow beyond its body size. A direct observation of the liquid motion above the undulator in combination with a quantitative analysis of the thin-film equations, yields the optimal speed at which this device transports the maximum amount of liquid for given geometric and fluid properties. This optimal wave speed scales inversely with the wavelength of the undulations and linearly with the cube of the film thickness. It is interesting to note that the key governing parameter, Ca/Bo can be interpreted as a ratio of two velocities - wave speed (Vw) to a characteristic relaxation or leveling speed (Vr = ρgH3/ηλ) at which surface undulations flatten out. This leveling process is dominated by gravity since the scale of undulations (~λ) is much larger than the capillary length (\(\sqrt{\gamma /\rho g}\)). Thus for Vr ≫ Vw, the undulator essentially works against a relaxed, flat interface, and liquid parcels primarily exhibit horizontal displacement over a period. In the other limit of Vr ≪ Vw, the free surface tends to beat in phase with the traveling boundary amplifying the vertical displacement, and subsequently reducing the net transport.

Our study demonstrates that the large-scale surface flow is a direct manifestation of the thin-film hydrodynamics above the undulator by showing how the optimal pumping speed captures the peak velocities in surface floaters. However, a quantitative analysis connecting the above two aspects of the flow field is necessary to exploit the full potential of this mechanism. Additionally, in the unexplored inertial regime, we expect the mechanism to showcase interesting dynamics due to the coupling between surface waves and finite-size particles37. In particular, for the transport of heavy or sedimentary particles, the Basset history force becomes important. This history force is proportional to the particle’s size and density, especially for large or dense particles in an unsteady flow. Thus future studies on driven heavy particles and floating rafts could exhibit intricate, new dynamics due to the Basset force. We believe that this work opens up new pathways for self-assembly and patterning at the mesoscale38, 39, bio-inspired strategies for remote sensing and actuation within liquids40,41, and control of interfacial flows using active boundaries42,43.

The models are designed in Fusion 360 (Autodesk). The helix is 3D printed in a Formlab Form 2 SLA printer by photo-crosslinking resin, whereas the outer shell comprising the top surface and rectangular links is printed in an Ultimaker S5 (Ultimaker Ltd.) using a blue TPE (thermoplastic elastomer). Due to the relative flexibility of TPE, the outer shell conforms to the helix. The helix is connected to a mini servo motor which is driven by DC power supply. All other parts (Base, Undulator holders, etc.) are printed using PLA (Polylactic acid) filaments on an Ultimaker S5 (Ultimaker Ltd.) printer.

We perform particle image velocimetry measurements on the thin liquid layer above the undulator. The viscous liquid is seeded with 10 μm glass microspheres (LaVision). A 520 nm 1W laser sheet (Laserland) illuminates a longitudinal plane in the middle of the undulator. Images are recorded by a Photron Fastcam SAZ camera at 500 frames per second. Image cross-correlation is performed in the open-source PIVlab44 to construct the velocity field.

The separation of scales (H ≪ λ) in the thin-film geometry leads to a set of reduced momentum equations and a flow field that is predominantly horizontal. Thus integration of the X-momentum equation with no-slip boundary condition on the undulator (Z = ha) and no shear stress condition at the free surface (Z = hf) results in

Similarly, we integrate the Z-momentum equation and apply the Young-Laplace equation at the free surface, which yields the following expression for the pressure p,

Integration of the continuity equation gives the volume flow rate, \(q=\int\nolimits_{{h}_{a}}^{{h}_{f}}{v}_{X}{{{{{{{\rm{d}}}}}}}}Z\), (per unit depth in this two-dimensional case). Plugging eqs. (11) and (12) into the expression of flow rate gives the following ODE,

This equation relates the yet unknown constant q and the unknown free-surface shape hf. We close the problem by imposing the following additional constraint on hf:

which states that the mean film thickness over one wavelength remains the same as that of the unperturbed interface, H. In dimensionless form, the above set of equations take the form of eqs. (3) and (4) of the main text.

For a direct comparison with experiments, we transform the flow rate, q back to the lab frame which is an explicit function of time, Q(x, t) = q + Vw(hf(x, t) − ha(x, t)). We seek for periodic free-surface shapes, such that hf = H + periodic terms. The time-averaged flow rate thus simplifies to

where the integration is performed at a fixed spatial location. In dimensionless form, this equation becomes, \(\left\langle \bar{Q}\right\rangle=\bar{q}+1\), as mentioned in the main text.

To estimate the efficacy of the pumping mechanism, we compare the output, and flow rate to the energy dissipation within the flow, which is given by,

Using the expression for the velocity field of eq. (11), we integrate on Z to find that the free-surface shape hf fully determines the amount of dissipation in the flow. In dimensionless form, dissipation becomes

where \(\bar{{{{{{{{\mathscr{E}}}}}}}}}={{{{{{{\mathscr{E}}}}}}}}H/(\eta \lambda {V}_{w}^{2})\).

Source data are provided in this paper. Additionally, all data and Mathematica scripts are available on the Open Science Framework (DOI 10.17605/OSF.IO/ERZ79). Source data are provided in this paper.

Okada, Y., Takeda, S., Tanaka, Y., Belmonte, J.-C. I. & Hirokawa, N. Mechanism of nodal flow: a conserved symmetry breaking event in left-right axis determination. Cell 121, 633–644 (2005).

Article CAS PubMed Google Scholar

Cartwright, J. H., Piro, O. & Tuval, I. Fluid-dynamical basis of the embryonic development of left-right asymmetry in vertebrates. Proc. Natl Acad. Sci. USA 101, 7234–7239 (2004).

Article ADS CAS PubMed PubMed Central Google Scholar

Sleigh, M. A., Blake, J. R. & Liron, N. The propulsion of mucus by cilia. Am. Rev. Respir. Dis. 137, 726–741 (1988).

Article CAS PubMed Google Scholar

Blake, J. On the movement of mucus in the lung. J. Biomech. 8, 179–190 (1975).

Article CAS PubMed Google Scholar

Bustamante-Marin, X. M. & Ostrowski, L. E. Cilia and mucociliary clearance. Cold Spring Harbor Perspect. Biol. 9, a028241 (2017).

Article Google Scholar

Burns, J. & Parkes, T. Peristaltic motion. J. Fluid Mech. 29, 731–743 (1967).

Article ADS Google Scholar

Agbesi, R. J. A. & Chevalier, N. R. Flow and mixing induced by single, colinear, and colliding contractile waves in the intestine. Phys. Rev. Fluids 7, 043101 (2022).

Article ADS Google Scholar

Reis, P. M., Jung, S., Aristoff, J. M. & Stocker, R. How cats lap: water uptake by Felis catus. Science 330, 1231–1234 (2010).

Article ADS CAS PubMed Google Scholar

Gart, S., Socha, J. J., Vlachos, P. P. & Jung, S. Dogs lap using acceleration-driven open pumping. Proc. Natl Acad. Sci. USA 112, 15798–15802 (2015).

Article ADS CAS PubMed PubMed Central Google Scholar

Tuckerman, D. B. & Pease, R. F. W. High-performance heat sinking for VLSI. IEEE Electr. Device Lett. 2, 126–129 (1981).

Article ADS Google Scholar

Das, S. K., Choi, S. U. & Patel, H. E. Heat transfer in nanofluids-a review. Heat Transf. Eng. 27, 3–19 (2006).

Article ADS CAS Google Scholar

Laser, D. J. & Santiago, J. G. A review of micropumps. J. Micromech. Microeng. 14(apr), R35–R64 (2004).

Article Google Scholar

Iverson, B. D. & Garimella, S. V. Recent advances in microscale pumping technologies: a review and evaluation. Microfluid. Nanofluidics 5, 145–174 (2008).

Article CAS Google Scholar

Kirby, B. J. Micro-and Nanoscale Fluid Mechanics: Transport in Microfluidic Devices. (Cambridge University Press, 2010).

Shields, A. et al. Biomimetic cilia arrays generate simultaneous pumping and mixing regimes. Proc. Natl Acad. Sci. USA 107, 15670–15675 (2010).

Article ADS CAS PubMed PubMed Central Google Scholar

Milana, E. et al. Metachronal patterns in artificial cilia for low Reynolds number fluid propulsion. Sci. Adv. 6, eabd2508 (2020).

Article ADS CAS PubMed PubMed Central Google Scholar

Wang, W. et al. Cilia metasurfaces for electronically programmable microfluidic manipulation. Nature 605, 681–686 (2022).

Article ADS CAS PubMed Google Scholar

Gu, H. et al. Magnetic cilia carpets with programmable metachronal waves. Nat. Commun. 11, 1–10 (2020).

Article ADS Google Scholar

Hanasoge, S., Hesketh, P. J. & Alexeev, A. Microfluidic pumping using artificial magnetic cilia. Microsyst. Nanoeng. 4, 1–9 (2018).

Article ADS CAS Google Scholar

Liu, G. and Zhang, W. in Traveling-Wave Micropumps. 1017–1035. (Springer Singapore, 2018).

Ogawa, J., Kanno, I., Kotera, H., Wasa, K. & Suzuki, T. Development of liquid pumping devices using vibrating microchannel walls. Sens. Actuator A Phys. 152, 211–218 (2009).

Article CAS Google Scholar

Chatzigiannakis, E., Jaensson, N. & Vermant, J. Thin liquid films: where hydrodynamics, capillarity, surface stresses and intermolecular forces meet. Curr. Opin. Colloid Interface Sci. 53, 101441 (2021).

Article CAS Google Scholar

Langevin, D. Influence of interfacial rheology on foam and emulsion properties. Adv. Colloid Interface Sci. 88, 209–222 (2000).

Article CAS PubMed Google Scholar

Saveanu, L. & Martín, P. R. Pedal surface collecting as an alternative feeding mechanism of the invasive apple snail Pomacea canaliculata (caenogastropoda: Ampullariidae). J. Molluscan Stud. 79, 11–18 (2013).

Article Google Scholar

Saveanu, L. & Martín, P. R. Neuston: a relevant trophic resource for apple snails? Limnologica 52, 75–82 (2015).

Article Google Scholar

Joo, S., Jung, S., Lee, S., Cowie, R. H. & Takagi, D. Freshwater snail feeding: lubrication-based particle collection on the water surface. J. R. Soc. Interface 17, 20200139 (2020).

Article PubMed PubMed Central Google Scholar

Huang, Z., Yang, H., Xu, K., Wu, J. & Zhang, J. Collecting differently sized particles on water surface by maneuvering pedal waves on the foot of the water snail Pomacea canaliculata. Soft Matter 18, 7850–7858 (2022).

Article ADS CAS PubMed Google Scholar

Voise, J. & Casas, J. The management of fluid and wave resistances by whirligig beetles. J. R. Soc. Interface 7, 343–352 (2010).

Article PubMed Google Scholar

Bush, J. W. & Hu, D. L. Walking on water: biolocomotion at the interface. Annu. Rev. Fluid Mech. 38, 339–369 (2006).

Article ADS MathSciNet MATH Google Scholar

Taylor, G. I. Analysis of the swimming of microscopic organisms. Proc. R. Soc. London. Ser. A. Math. Phys. Sci. 209, 447–461 (1951).

ADS MathSciNet MATH Google Scholar

Shaik, V. A. & Ardekani, A. M. Swimming sheet near a plane surfactant-laden interface. Phys. Rev. E 99, 033101 (2019).

Article ADS CAS PubMed Google Scholar

Dias, M. A. & Powers, T. R. Swimming near deformable membranes at low Reynolds number. Phys. Fluids 25, 101901 (2013).

Article ADS Google Scholar

Zarrouk, D. et al. Single actuator wave-like robot (saw): design, modeling, and experiments. Bioinspir. Biomim. 11, 046004 (2016).

Article ADS PubMed Google Scholar

Oron, A., Davis, S. H. & Bankoff, S. G. Long-scale evolution of thin liquid films. Rev. Modern Phys. 69, 931 (1997).

Article ADS CAS Google Scholar

Craster, R. V. & Matar, O. K. Dynamics and stability of thin liquid films. Rev. Modern Phys. 81, 1131 (2009).

Article ADS CAS Google Scholar

Lee, S., Bush, J. W., Hosoi, A. & Lauga, E. Crawling beneath the free surface: water snail locomotion. Phys. Fluids 20, 082106 (2008).

Article ADS MATH Google Scholar

Punzmann, H., Francois, N., Xia, H., Falkovich, G. & Shats, M. Generation and reversal of surface flows by propagating waves. Nat. Phys. 10, 658–663 (2014).

Article CAS Google Scholar

Zamanian, B. et al. Interface-directed self-assembly of cell-laden microgels. Small 6, 937–944 (2010).

Article CAS PubMed PubMed Central Google Scholar

Snezhko, A., Belkin, M., Aranson, I. & Kwok, W.-K. Self-assembled magnetic surface swimmers. Phys. Rev. Lett. 102, 118103 (2009).

Article ADS CAS PubMed Google Scholar

Santhanakrishnan, A., Dollinger, M., Hamlet, C. L., Colin, S. P. & Miller, L. A. Flow structure and transport characteristics of feeding and exchange currents generated by upside-down cassiopea jellyfish. J. Exp. Bio.215, 2369–2381 (2012).

Article Google Scholar

Ryu, S., Pepper, R. E., Nagai, M. & France, D. C. Vorticella: a protozoan for bio-inspired engineering. Micromachines 8, 4 (2016).

Article PubMed Central Google Scholar

Manna, R. K., Laskar, A., Shklyaev, O. E. & Balazs, A. C. Harnessing the power of chemically active sheets in solution. Nat. Rev. Phys. 4, 125–137 (2022).

Article Google Scholar

Laskar, A., Shklyaev, O. E. & Balazs, A. C. Designing self-propelled, chemically active sheets: Wrappers, flappers, and creepers. Sci. Adv. 4, eaav1745 (2018).

Article ADS CAS PubMed PubMed Central Google Scholar

Thielicke, W. & Stamhuis, E. Pivlab–towards user-friendly, affordable and accurate digital particle image velocimetry in Matlab. J. Open Res. Softw. 2, 1 (2014).

Article Google Scholar

Download references

We thank Yohan Sequeira and Sarah MacGregor for their initial contributions. C.R., D.T., S.L., and S.J. acknowledge the support of NSF through grant no CMMI-2042740. A.P. acknowledges startup funding from Syracuse University.

Mechanical & Aerospace Engineering Department and BioInspired Syracuse, Syracuse University, Syracuse, NY, 13244, USA

Anupam Pandey

Department of Mechanical Engineering, University of Minnesota, Minneapolis, MN, 55455, USA

Zih-Yin Chen & Sungyon Lee

Department of Biological & Environmental Engineering, Cornell University, Ithaca, NY, 14853, USA

Jisoo Yuk, Chris Roh & Sunghwan Jung

Sibley School of Mechanical & Aerospace Engineering, Cornell University, Ithaca, NY, 14853, USA

Yuming Sun

Department of Mathematics, University of Hawaii at Manoa, Honolulu, HI, 96822, USA

Daisuke Takagi

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

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

A.P., S.J. conceived the idea. A.P., J.Y., Y.S., C.R. and S.J. designed and performed experiments, and analyzed data. Z.C., D.T. and S.L. developed the theoretical and numerical models. All authors wrote the paper collectively.

Correspondence to Anupam Pandey or Sunghwan Jung.

The authors declare no competing interests.

Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. A peer review file is available.

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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

Pandey, A., Chen, ZY., Yuk, J. et al. Optimal free-surface pumping by an undulating carpet. Nat Commun 14, 7735 (2023). https://doi.org/10.1038/s41467-023-43059-8

Download citation

Received: 23 February 2023

Accepted: 30 October 2023

Published: 25 November 2023

DOI: https://doi.org/10.1038/s41467-023-43059-8

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