pith. sign in

arxiv: 2603.00308 · v1 · submitted 2026-02-27 · ⚛️ physics.flu-dyn

Using surface acoustic waves to drive thin film flow over an obstacle

Pith reviewed 2026-05-15 18:21 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords surface acoustic wavesthin film flowacoustic streamingcapillary flowobstacle coatingpiezoelectric actuatorultrasonic propulsion
0
0 comments X

The pith

Surface acoustic waves drive silicone oil films to climb and traverse obstacles.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

Surface acoustic waves at MHz frequencies propel thin films of silicone oil across solid obstacles placed on a piezoelectric substrate. Experiments demonstrate that nanometer-scale wave amplitudes suffice for macroscopic films to climb and cross these obstacles. The flow outcome depends on the interplay of ultrasonic forcing, capillarity, and gravity. A simplified two-dimensional model includes obstacle height directly in the thin-film evolution equation together with a representation of acoustic streaming, and the resulting simulations match experimental observations qualitatively.

Core claim

Nanometer-amplitude surface acoustic waves propagating in the substrate propel macroscopic silicone oil films to climb and traverse solid obstacles. The oil dynamics reveal rich coupling between ultrasonic forcing, capillarity, and gravity that determines coating success. A simplified two-dimensional theoretical model incorporates obstacle geometry directly in the oil thin-film evolution equation and introduces a new representation of acoustic streaming in the presence of substrate height variations. Simulations show qualitative agreement with the experiments.

What carries the argument

Simplified two-dimensional thin-film evolution equation that incorporates substrate height variations directly and includes a new representation of acoustic streaming from surface acoustic waves.

If this is right

  • The balance of ultrasonic forcing, capillarity, and gravity decides whether a film successfully coats the obstacle.
  • Qualitative agreement indicates that the two-dimensional model captures the dominant physical mechanisms.
  • The framework supplies a method to predict coating outcomes from wave parameters and obstacle geometry.
  • Ultrasonic propulsion offers a route to drive thin-film flow over non-flat surfaces without mechanical contact.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same acoustic streaming representation may apply to other wave-driven coating problems on textured substrates.
  • Testing the model on multiple obstacles or three-dimensional topographies could reveal additional flow patterns.
  • Adjusting wave frequency or fluid viscosity might produce controlled deposition profiles for practical coating tasks.

Load-bearing premise

The simplified two-dimensional model with a representation of acoustic streaming accurately captures the essential three-dimensional physics and actual acoustic forcing mechanism.

What would settle it

A quantitative mismatch between predicted film thickness profiles or climbing speeds and experimental data when obstacle height or wave amplitude is varied.

Figures

Figures reproduced from arXiv: 2603.00308 by Avital R. Einhorn, Javier A. Diez, Linda J. Cummings, Lou Kondic, Mark Fasano, Ofer Manor, Yifan Li.

Figure 1
Figure 1. Figure 1: (a) Top view of the SAW actuator used in the experiment (dimensions: 24.5 mm × 10.8 mm × 0.5 mm), comprising an inter-digital transducer to the left made from metal electrodes fabricated atop the piezoelectric lithium-niobate substrate of the actuator, and a PDMS obstacle placed atop the SAW actuator (to the right). (b) The SAW actuator, supporting a 3D printed PDMS obstacle viewed from above, attached to … view at source ↗
Figure 2
Figure 2. Figure 2: Side view snapshots of a typical experiment (An=1.43 nm), where oil climbs over a ramp-shaped obstacle, indicated by the sketch in (a). During the experiment, (b) initially, the oil film (highlighted using a thin, white curve) moves along the upper surface of the SAW transducer to make contact with the obstacle; the oil film then deforms and (c,d,e) climbs up the obstacle and (f) reaches the peak. Time t i… view at source ↗
Figure 3
Figure 3. Figure 3: Time variation of the climbing height of the film atop the ramp for (a) An = 0.52 nm to 1.69 nm and (b) the same results with a logarithmic time axis, indicating a change in the mechanism that drives the climb above a height of approximately 1 mm. Symbols are experimental data and the connecting curves are guides for the eye. The uncertainty in climb height is 5% based on the spatial resolution of the came… view at source ↗
Figure 4
Figure 4. Figure 4: Maximum climbing height, Hmax, of the film on the ramp, vs SAW normal amplitude displacement (An), for 8 µl of silicone oil volume (blue dots) and for 8+8 µl of oil volume (red squares). The dashed line shows the top of the ramp. The uncertainty in climb height is 5% based on the spatial resolution of the camera. (c) Bump obstacle To better understand how substrate topography affects spreading, we also con… view at source ↗
Figure 5
Figure 5. Figure 5: (a) A sketch of the oil covering bump obstacle and the (b1–e1) side view and (b2–e2) corresponding top view snapshots of a SAW driven silicone oil film climbing over and coating a bump; the oil film outline is marked by a dark, thin line. The snapshot sequence commences at t = 0 (b1, b2) when the oil front reaches the rear of the bump. The oil then climbs over it under the influence of SAW (c1–d1; c2–d2), … view at source ↗
Figure 6
Figure 6. Figure 6: (a) Front position, Xf, measured relative to the rear (left edge) of the bump, and (b) front height, Hf, as functions of time, for a bump of height ho = 0.92 mm. Time is measured from the moment the oil reaches the rear of the bump. The SAW amplitude, An, is varied from 0.78 nm to 1.69 nm. Symbols are experimental data and the connecting curves are guides for the eye. These measurements have a 5% error due… view at source ↗
Figure 7
Figure 7. Figure 7: Time τ to traverse the bump obstacle as a function of An for different bump heights: (a) ho = 0.26 mm, (b) ho = 0.65 mm, and (c) ho = 0.92 mm. The red dots correspond to a drop of volume Vd = 8 µl. The blue dots in (a) are for a drop of volume Vd = 2 µl, while in (c) they stand for the upside down experiment. Climbing time is the time elapsed from when the drop front touches the rear of the bump to when it… view at source ↗
Figure 8
Figure 8. Figure 8: Evolution of the thickness profile as the film climbs the ramp, for A = 4 nm at different times. The black solid line denotes the underlying ramp obstacle. The blue solid line represents the height profile of the liquid film at the current time. The red line indicates the thickness-averaged horizontal velocity component, v, given by Eq. (3.7), dimensionalized by the scale L/T. Full simulation videos availa… view at source ↗
Figure 9
Figure 9. Figure 9: (a) Climbing height of the film on the ramp, Hf, as a function of time from simulations for A = 0.5 nm (magenta) to A = 6 nm (olive color) with ∆A = 0.5 nm increasing according to the black arrow. (b) Steady-state (at large times) height profile of the film on a ramp for different values of the SAW amplitude, A. comparable to capillary forces. As a result, the front climbing height increases slowly in time… view at source ↗
Figure 10
Figure 10. Figure 10: (a) Climbing heights as a function of time for SAW amplitudes A = 1 nm (black), A = 4 nm (orange), and A = 6 nm (green) for ramp slopes of 31.4 ◦ (dashed lines), 42.4 ◦ (solid lines), and 53.9 ◦ (dot-dashed lines); colors shown here are as in [PITH_FULL_IMAGE:figures/full_fig_p015_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: (a) Climbing heights, Hf , of the drop’s front as a function of time on the ramp obstacle for amplitude A = 4 nm calculated for four different drop volumes, Vd: 3 µl (blue), 5 µl (green), and 8 µl (red). (b) Steady-state thickness profile for each drop volume, Vd, for large times (i.e., when there is no flow). the total fluid thickness reaches a fixed cutoff value H = h + s = h ∗ , independent of the tota… view at source ↗
Figure 12
Figure 12. Figure 12: Side view of the bump obstacle used in experiments, with the circular arc employed in simulations overlaid to illustrate the close agreement between the experimental geometry and the model representation [PITH_FULL_IMAGE:figures/full_fig_p016_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: (a)-(f) Evolution of thickness profile (blue line) and thickness-averaged horizontal velocity component (red, dashed line) for A = 4 nm at indicated times as the oil climbs a bump obstacle (black). The gray dashed lines denote the positions of x ∗ r (t), x∗ f (t). Full simulation videos are available upon request for a range of A values. the moment the front reaches it, yet this clearly does not happen: t… view at source ↗
Figure 14
Figure 14. Figure 14: Positions of the front cutoff line x ∗ f (green), rear cutoff line x ∗ r (magenta), the front contact line xf (red, dashed), and the rear contact line xr (blue, dashed) of the oil film for the case of A = 4 nm shown in [PITH_FULL_IMAGE:figures/full_fig_p018_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Maximum height on the ramp obstacle reached by the film front, Hmax, as a function of the SAW amplitude, An (experiments) or A (simulations), as obtained from the experiments (black circles) and simulations (red points connected by a solid line). The shaded gray area indicates the general trend of the experimental data. While it is not clear why this difference arises, some plausible reasons, in addition … view at source ↗
Figure 16
Figure 16. Figure 16: Climbing time, τ , from the experiments (black dots) and simulations (red squares) for selected values of the SAW amplitude, An (experiments) or A (simulations), for four different sized bump obstacles. The black (red) dashed line corresponds to an approximating function for the experimental (numerical) data of the type, τ = aA¯b , where A¯ stands for either A or An. The exponents b for both experimental … view at source ↗
read the original abstract

We study a new paradigm for ultrasonic driven object coating by using a model system where MHz-level surface acoustic waves (SAWs) drive the spreading of a silicone oil film atop topographical obstacles. We use experiments to show that nanometer-amplitude SAWs, propagating in the substrate of a piezoelectric actuator, propel macroscopic oil films to climb and traverse solid obstacles placed on the actuator. The oil dynamics reveal rich coupling between ultrasonic forcing, capillarity, and gravity; the balance of which determines coating success. We formulate a simplified two-dimensional theoretical model that incorporates obstacle geometry directly in the oil thin-film evolution equation, introducing a new representation of acoustic streaming in the presence of substrate height variations. Despite the simplifications inherent in the modeling, simulations show qualitative agreement with the experiments, providing evidence that the model captures the key physics.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The paper demonstrates experimentally that MHz surface acoustic waves (SAWs) propagating in a piezoelectric substrate can drive macroscopic silicone-oil films to climb and traverse solid topographical obstacles. It formulates a simplified 2D lubrication model that incorporates obstacle geometry directly into the thin-film evolution equation via a new effective acoustic-streaming body-force term that accounts for substrate height variations. Simulations of this model are reported to achieve qualitative agreement with the observed climbing and coating dynamics.

Significance. If the central modeling claim holds, the work establishes a new, controllable mechanism for ultrasonic thin-film transport over obstacles, with potential relevance to coating processes and microfluidic actuation. The explicit inclusion of topography in the evolution equation and the qualitative match to experiment constitute a useful step beyond purely flat-substrate SAW-driven flows. The absence of quantitative validation metrics, however, limits the strength of the evidence that the new streaming representation faithfully captures the underlying 3D forcing.

major comments (3)
  1. [Theoretical model (evolution equation)] The new representation of acoustic streaming that incorporates substrate height variations is introduced directly into the 2D evolution equation without a derivation from 3D acoustic boundary-layer analysis or validation against full-wave simulations. Because this term is the key modeling innovation enabling the predicted climbing, its accuracy for finite-width obstacles (where lateral streaming and meniscus effects are expected) remains unverified.
  2. [Results and comparison] The reported agreement between simulations and experiments is characterized only as qualitative. No quantitative metrics—such as measured versus predicted film-height profiles, front speeds, or L2 errors—are provided across multiple obstacle aspect ratios or SAW amplitudes. This omission makes it impossible to determine whether the observed match arises from faithful capture of the forcing physics or from effective adjustment within the simplified 2D framework.
  3. [Model assumptions and limitations] The reduction to a strictly two-dimensional model necessarily omits spanwise variations in streaming and capillary effects around the sides of finite-width obstacles. The manuscript does not discuss the expected magnitude of these 3D corrections or provide any auxiliary 3D simulation or scaling argument to justify their neglect for the geometries studied.
minor comments (2)
  1. [Figures] Figure captions should explicitly state the SAW frequency, amplitude, and oil viscosity used in each panel to allow direct comparison with the model parameters.
  2. [Abstract and Theory] The abstract states that the model 'incorporates obstacle geometry directly'; the precise functional form of this incorporation (e.g., how the height function h(x) enters the streaming term) should be written explicitly in the main text for reproducibility.

Simulated Author's Rebuttal

3 responses · 1 unresolved

We thank the referee for the constructive and detailed comments. We address each major point below, indicating where revisions have been made to strengthen the manuscript.

read point-by-point responses
  1. Referee: The new representation of acoustic streaming that incorporates substrate height variations is introduced directly into the 2D evolution equation without a derivation from 3D acoustic boundary-layer analysis or validation against full-wave simulations. Because this term is the key modeling innovation enabling the predicted climbing, its accuracy for finite-width obstacles (where lateral streaming and meniscus effects are expected) remains unverified.

    Authors: We agree that additional justification for the effective streaming term would improve clarity. We have added a new subsection (Section 3.2) providing a heuristic derivation based on extending the standard flat-substrate acoustic body force to account for local substrate height variations, under the assumption that the acoustic boundary layer thickness remains much smaller than the film height. Scaling arguments are included to support applicability to the obstacle geometries studied. A complete first-principles 3D boundary-layer derivation and full-wave validation lie outside the present scope due to computational demands and are noted as future work. revision: partial

  2. Referee: The reported agreement between simulations and experiments is characterized only as qualitative. No quantitative metrics—such as measured versus predicted film-height profiles, front speeds, or L2 errors—are provided across multiple obstacle aspect ratios or SAW amplitudes. This omission makes it impossible to determine whether the observed match arises from faithful capture of the forcing physics or from effective adjustment within the simplified 2D framework.

    Authors: We acknowledge the benefit of quantitative metrics. The revised manuscript now includes direct comparisons of front propagation speeds versus SAW amplitude and obstacle height, together with L2-norm errors between simulated and experimental film profiles for multiple cases. These are presented in a new Figure 6 and accompanying table, showing agreement within approximately 15% for speeds and profile deviations, providing stronger evidence for the model's fidelity. revision: yes

  3. Referee: The reduction to a strictly two-dimensional model necessarily omits spanwise variations in streaming and capillary effects around the sides of finite-width obstacles. The manuscript does not discuss the expected magnitude of these 3D corrections or provide any auxiliary 3D simulation or scaling argument to justify their neglect for the geometries studied.

    Authors: We have expanded the model limitations paragraph in Section 5 to include scaling estimates of 3D effects. For the obstacle widths employed (several millimeters, much larger than the capillary length and film thickness), spanwise flows are suppressed by geometric confinement, with estimated corrections below 10-15% based on order-of-magnitude analysis of lateral meniscus curvature and streaming. This supports the 2D approximation for the reported experiments, while narrower obstacles would indeed require 3D treatment. revision: yes

standing simulated objections not resolved
  • A complete derivation of the acoustic streaming term from 3D boundary-layer analysis together with validation against full-wave simulations.

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained from lubrication theory

full rationale

The paper formulates a simplified 2D thin-film evolution equation by extending standard lubrication theory to include obstacle geometry and a representation of acoustic streaming. Simulations are then compared qualitatively to experiments. No load-bearing step reduces a claimed prediction or result to a fitted parameter, self-defined quantity, or self-citation chain by construction. The central claim of qualitative agreement rests on independent numerical solution of the stated PDE system rather than tautological equivalence to inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the lubrication approximation for thin films and a novel but simplified effective representation of acoustic streaming; no explicit free parameters or new invented entities are identified from the abstract.

axioms (2)
  • domain assumption The lubrication approximation holds for the macroscopic oil film dynamics under MHz SAW forcing.
    Invoked when formulating the two-dimensional thin-film evolution equation.
  • ad hoc to paper Acoustic streaming effects can be incorporated via an effective forcing term that accounts for substrate height variations.
    This is the new representation introduced in the model.

pith-pipeline@v0.9.0 · 5459 in / 1305 out tokens · 67255 ms · 2026-05-15T18:21:33.123781+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

42 extracted references · 42 canonical work pages

  1. [1]

    2005 Controlled microfluidic interfaces.Nature437, 648–655

    Atencia J, Beebe DJ. 2005 Controlled microfluidic interfaces.Nature437, 648–655

  2. [2]

    2004 Engineering flows in small devices.Ann

    Stone HA, Stroock AD, Ajdari A. 2004 Engineering flows in small devices.Ann. Rev. Fluid Mech.36, 381–411. 24royalsocietypublishing.org/journal/rspa Proc R Soc A 0000000

  3. [3]

    2006 The origins and the future of microfluidics.Nature442, 368–373

    Whitesides GM. 2006 The origins and the future of microfluidics.Nature442, 368–373

  4. [4]

    2001 MEMS-enabled thermal management of high-heat-flux devices EDIFICE: embedded droplet impingement for integrated cooling of electronics.Exp

    Amon CH, Murthy J, Yao SC, Narumanchi S, Wu CF, Hsieh CC. 2001 MEMS-enabled thermal management of high-heat-flux devices EDIFICE: embedded droplet impingement for integrated cooling of electronics.Exp. Thermal Fluid Sci.25, 231–242

  5. [5]

    2006 Direct liquid cooling of high flux micro and nano electronic components.Proc

    Bar-Cohen A, Arik M, Ohadi M. 2006 Direct liquid cooling of high flux micro and nano electronic components.Proc. IEEE94, 1549–1570

  6. [6]

    1974 Evaporation from thin water films on horizontal tubes.Ind

    Fletcher LS, Sernas V, Galowin LS. 1974 Evaporation from thin water films on horizontal tubes.Ind. Eng. Chem. Process Des. Dev.13, 265–269

  7. [7]

    2004 Gravity-driven flow of continuous thin liquid films on non-porous substrates with topography.J

    Gaskell PH, Jimack PK, Sellier M, Thompson HM, Wilson MCT. 2004 Gravity-driven flow of continuous thin liquid films on non-porous substrates with topography.J. Fluid Mech. 509, 253–280

  8. [8]

    1988 Fundamentals of topographic substrate leveling.J

    Stillwagon LE, Larson RG. 1988 Fundamentals of topographic substrate leveling.J. Appl. Phys.63, 5251–5258

  9. [9]

    2003 Surface tension driven flow on a moving curved surface.J

    Howell PD. 2003 Surface tension driven flow on a moving curved surface.J. Eng. Math.45, 283–308

  10. [10]

    1995 Modeling of coating flows on curved surfaces.J

    Schwartz LW, Weidner DE. 1995 Modeling of coating flows on curved surfaces.J. Eng. Math. 29, 91–103

  11. [11]

    1997 Long-scale evolution of thin liquid films.Rev

    Oron A, Davis SH, Bankoff SG. 1997 Long-scale evolution of thin liquid films.Rev. Mod. Phys.69, 931

  12. [12]

    2009 Dynamics and stability of thin liquid films.Rev

    Craster R, Matar O. 2009 Dynamics and stability of thin liquid films.Rev. Mod. Phys.81, 1131

  13. [13]

    1989 Liquid streaming and droplet formation caused by leaky Rayleigh waves.Proc

    Shiokawa S, Matsui Y, Ueda T. 1989 Liquid streaming and droplet formation caused by leaky Rayleigh waves.Proc. IEEE Ultras. Symp.1, 643–646

  14. [14]

    2018 Vibration-driven mass transfer and dynamic wetting.Curr

    Morozov M, Manor O. 2018 Vibration-driven mass transfer and dynamic wetting.Curr. Opin. Colloid Interface Sci.36, 37–45

  15. [15]

    2015 Spreading dynamics of a partially wetting water film atop a MHz substrate vibration.Phys

    Altshuler G, Manor O. 2015 Spreading dynamics of a partially wetting water film atop a MHz substrate vibration.Phys. Fluids27, 102103

  16. [16]

    HoreshA,KhaikinD,KarnilawM,ZigelmanA,ManorO.2019Acoustogravitationalbalance in climbing films.Phys. Rev. Fluids4, 022001(R)

  17. [17]

    1934 On the acoustic radiation pressure on spheres.Proc

    King L. 1934 On the acoustic radiation pressure on spheres.Proc. R. Soc. Lond. A147, 212–240

  18. [18]

    1982 Acoustic radiation pressure produced by a beam of sound.J

    Chu B, Apfel RE. 1982 Acoustic radiation pressure produced by a beam of sound.J. Acoust. Society. Am.72, 1673–1687

  19. [19]

    2000 A general theory of Rayleigh and Langevin radiation pressure.J

    Hasegawa T, Kido T, Iizuka T, Matsuoka C. 2000 A general theory of Rayleigh and Langevin radiation pressure.J. Acoust. Soc. Jap. E21, 146–152

  20. [20]

    1884 On the Circulation of Air Observed in Kundt’s Tubes, and on Some Allied Acoustical Problems.Phil

    Rayleigh L. 1884 On the Circulation of Air Observed in Kundt’s Tubes, and on Some Allied Acoustical Problems.Phil. Trans. Royal Soc. London175, 1–21

  21. [21]

    1978 Acoustic streaming.J

    Lighthill J. 1978 Acoustic streaming.J. Sound Vib.61, 391–418

  22. [22]

    1953 Mass Transport in Water Waves.Philos

    Longuet-Higgins MS. 1953 Mass Transport in Water Waves.Philos. Trans. R. Soc. A.245, 535–581

  23. [23]

    2012 The appearance of boundary layers and drift flows due to high-frequency surface waves.J

    Manor O, Yeo LY, Friend JR. 2012 The appearance of boundary layers and drift flows due to high-frequency surface waves.J. Fluid Mech.707, 482–495

  24. [24]

    2014 Surface acoustic wave microfluidics.Annu

    Yeo L, Friend J. 2014 Surface acoustic wave microfluidics.Annu. Rev. Fluid Mech.46, 379–406

  25. [25]

    2017 An extended Landau Levich model for the dragging of a thin liquid film with a propagating surface acoustic wave.J

    Morozov M, Manor O. 2017 An extended Landau Levich model for the dragging of a thin liquid film with a propagating surface acoustic wave.J. Fluid Mech.810, 307–322

  26. [26]

    1948 Vortices and streams caused by sound waves.Phys

    Eckart C. 1948 Vortices and streams caused by sound waves.Phys. Rev.73, 68–76

  27. [27]

    1952 Acoustic streaming due to attenuated plane waves.J

    Nyborg WL. 1952 Acoustic streaming due to attenuated plane waves.J. Acoust. Soc. Am. 25, 1–8

  28. [28]

    2007 Vibration-induced climbing of drops.Phys

    Brunet P, Eggers J, Deegan RD. 2007 Vibration-induced climbing of drops.Phys. Rev. Lett. 99, 144501

  29. [29]

    2009 Motion of a drop driven by substrate vibrations.Eur

    Brunet P, Eggers J, Deegan RD. 2009 Motion of a drop driven by substrate vibrations.Eur. Phys. J. Special Topics166, 11

  30. [30]

    2025 Modeling the dynamics of an oil drop driven by a surface acoustic wave in the underlying substrate

    Fasano M, Li Y, Diez JA, D’Addesa J, Manor O, Cummings LJ, Kondic L. 2025 Modeling the dynamics of an oil drop driven by a surface acoustic wave in the underlying substrate. J. Fluid Mech.1022, A49

  31. [31]

    Langmuir33, 7352–7363

    ParkJ,KumarS.2017DropletSlidingonanInclinedSubstratewithaTopographicalDefect. Langmuir33, 7352–7363. PMID: 28650652. 25royalsocietypublishing.org/journal/rspa Proc R Soc A 0000000

  32. [32]

    1996Elastic Waves in Solids I

    Royer D, Dieulesaint E. 1996Elastic Waves in Solids I. Springer-Verlag, Berlin/Heidelberg

  33. [33]

    2012 Unique fingering instabilities and soliton-like wave propagation in thin acoustowetting films.Nat

    Rezk AR, Manor O, Friend JR, Yeo LY. 2012 Unique fingering instabilities and soliton-like wave propagation in thin acoustowetting films.Nat. Commun.3, 1167

  34. [34]

    2014 Double flow reversal in thin liquid films driven by megahertz-order surface vibration.Proc

    Rezk A, Manor O, Yeo L, Friend J. 2014 Double flow reversal in thin liquid films driven by megahertz-order surface vibration.Proc. R. Soc. A470, 20130765

  35. [35]

    2002 Computing three-dimensional thin film flows including contact lines.J

    Diez JA, Kondic L. 2002 Computing three-dimensional thin film flows including contact lines.J. Comput. Phys.183, 274

  36. [36]

    1971Acoustic streaming

    Zarembo LK. 1971Acoustic streaming. High-Intensity Ultrasonic Fields. Boston, MA: Springer

  37. [37]

    2021 Unraveling the complex dynamics of acoustofluidics

    Orosco J, Friend J. 2021 Unraveling the complex dynamics of acoustofluidics

  38. [38]

    2023 Theory of acoustic streaming for arbitrary Reynolds number flow.J

    Dubrovski O, Friend J, Manor O. 2023 Theory of acoustic streaming for arbitrary Reynolds number flow.J. Fluid Mech.975, A4

  39. [39]

    2017 An extended Landau-Levich model for the dragging of a thin liquid film with a propagating surface acoustic wave.J

    Morozov M, Manor O. 2017 An extended Landau-Levich model for the dragging of a thin liquid film with a propagating surface acoustic wave.J. Fluid Mech.810, 307

  40. [40]

    2011 Streaming by leaky surface acoustic waves.Proc

    Vanneste J, Bühler O. 2011 Streaming by leaky surface acoustic waves.Proc. R. Soc. A. 467, 1779–1800

  41. [41]

    1970 Propagation of Surface Waves at the Boundary Between a Piezoelectric Crystal and a Fluid Medium.IEEE Trans

    Campbell JJ, Jones WR. 1970 Propagation of Surface Waves at the Boundary Between a Piezoelectric Crystal and a Fluid Medium.IEEE Trans. Son. Ultrason.17, 71–76

  42. [42]

    2003 Instabilities in Gravity Driven Flow of Thin Fluid Films.SIAM Review45, 95

    Kondic L. 2003 Instabilities in Gravity Driven Flow of Thin Fluid Films.SIAM Review45, 95