Interface Fragmentation via Horizontal Vibration: A Pathway to Scalable Monodisperse Emulsification
Pith reviewed 2026-05-19 10:47 UTC · model grok-4.3
The pith
Horizontal vibration of two stratified liquids produces monodisperse micro-emulsions by breaking ordered Faraday waves along the container walls.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
When viscous forces dominate interfacial dynamics, horizontal vibration excites a single line of ordered Faraday waves along each end wall of a rectangular container holding two stably stratified immiscible liquids; the wave tips then break up into a regular array of droplets once the non-dimensional acceleration reaches a threshold that scales as N to the minus one-half times omega star to the three-halves, where N is the kinematic viscosity ratio and omega star is the forcing frequency on the viscous-capillary scale.
What carries the argument
The single line of ordered Faraday waves excited along each end wall under viscous-dominated conditions, whose tips fragment into droplets.
Load-bearing premise
Viscous forces must dominate the interfacial dynamics so that only a single line of ordered Faraday waves forms along the walls.
What would settle it
Measure the critical acceleration at several viscosity ratios and forcing frequencies and test whether the values collapse onto the predicted N to the minus one-half times omega star to the three-halves scaling.
Figures
read the original abstract
We present a scalable method for producing monodisperse micro-scale emulsions in a rectangular container holding two stably stratified layers of immiscible liquids by applying horizontal vibration. This setup enables the excitation of a single line of ordered Faraday waves along each end wall when viscous forces dominate interfacial dynamics. Our experiments and theoretical modelling show that the critical non-dimensional acceleration for the breakup of the wave tips in a regular array of droplets scales as $N^{-1/2} \omega^{*3/2}$, where $N$ is the kinematic viscosity ratio and $\omega^{*}$ is the frequency of forcing on the viscous-capillary scale. The droplet diameter can be easily tuned by varying the forcing parameters, and the number of droplets generated per cycle is proportional to the width of the container.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents a scalable method for monodisperse micro-emulsion production in a rectangular container with two immiscible, stably stratified liquids subjected to horizontal vibration. Under viscous dominance, a single line of ordered Faraday waves forms along each end wall; experiments and modelling indicate that the critical non-dimensional acceleration for tip breakup in the resulting droplet array scales as N^{-1/2} ω^{*3/2}, where N is the kinematic viscosity ratio and ω* is the forcing frequency on the viscous-capillary scale. Droplet diameter is tunable via forcing parameters and the number of droplets per cycle scales with container width.
Significance. If the reported scaling is rigorously derived and experimentally confirmed without fitting, the work supplies a simple, vibration-driven route to uniform emulsions that avoids microfluidic complexity. The combination of a claimed parameter-free scaling with direct visualization of ordered wave breakup constitutes a potentially useful contribution to soft-matter fluid dynamics and emulsification technology.
major comments (2)
- [§3] §3 (theoretical modelling): the central scaling N^{-1/2} ω^{*3/2} is stated to emerge from viscous-capillary balance at the wave tips under horizontal forcing. The derivation must be shown explicitly from the Navier-Stokes equations with the horizontal body force and free-surface boundary conditions; it is not evident how the viscosity ratio N enters at the -1/2 power or how the ω^{*3/2} exponent is obtained without additional assumptions on tip curvature or neglect of inertial corrections.
- [§4] §4 (experimental results) and associated data tables: the abstract asserts that both experiments and modelling support the scaling, yet no quantitative comparison (measured critical accelerations versus predicted values across the tested range of N and ω*) is provided. Without these data or a clear statement that the exponent and prefactor are not obtained by fitting the target dataset, the claim of predictive power remains unverified.
minor comments (2)
- [Figures] Figure captions should explicitly state the range of N and ω* shown and whether the images correspond to the critical acceleration or a supercritical state.
- [Notation] The definition of the viscous-capillary time scale used to non-dimensionalize ω* should be written once in the methods and used consistently; the current notation leaves the precise form of the capillary velocity ambiguous.
Simulated Author's Rebuttal
We thank the referee for their careful review and insightful comments, which have prompted us to clarify and strengthen several aspects of the manuscript. We provide detailed responses to each major comment below.
read point-by-point responses
-
Referee: [§3] §3 (theoretical modelling): the central scaling N^{-1/2} ω^{*3/2} is stated to emerge from viscous-capillary balance at the wave tips under horizontal forcing. The derivation must be shown explicitly from the Navier-Stokes equations with the horizontal body force and free-surface boundary conditions; it is not evident how the viscosity ratio N enters at the -1/2 power or how the ω^{*3/2} exponent is obtained without additional assumptions on tip curvature or neglect of inertial corrections.
Authors: We appreciate the referee pointing out the need for an explicit derivation. In the revised manuscript, we will add a detailed derivation in §3, starting from the incompressible Navier-Stokes equations including the horizontal body force term ρ a cos(ω t) in the x-direction, and applying the kinematic and stress balance conditions at the free surface. The viscosity ratio N = ν_upper / ν_lower enters the scaling through the interfacial shear stress continuity, which determines the effective damping rate scaling as N^{-1/2} when matching the boundary layer solutions. The ω^{*3/2} dependence follows from non-dimensionalizing the capillary pressure and viscous dissipation terms using the viscous-capillary frequency scale, with the wave tip curvature obtained from the linear wave solution amplitude. We include order-of-magnitude analysis showing that inertial terms are smaller by a factor of the inverse Reynolds number in the relevant regime. revision: yes
-
Referee: [§4] §4 (experimental results) and associated data tables: the abstract asserts that both experiments and modelling support the scaling, yet no quantitative comparison (measured critical accelerations versus predicted values across the tested range of N and ω*) is provided. Without these data or a clear statement that the exponent and prefactor are not obtained by fitting the target dataset, the claim of predictive power remains unverified.
Authors: We agree that including a quantitative comparison is essential to substantiate the predictive nature of the scaling. In the revised manuscript, we will add a figure in §4 that plots the experimentally measured critical accelerations against the theoretical prediction N^{-1/2} ω^{*3/2} for the range of viscosity ratios and frequencies tested. We will also add an explicit statement in the text that the scaling and prefactor were obtained from the theoretical model independently of the experimental data, with no fitting performed on the target dataset. The agreement between theory and experiment will be discussed in terms of the observed collapse of the data. revision: yes
Circularity Check
Scaling emerges from viscous-capillary analysis of governing equations without reduction to fitted inputs or self-citation
full rationale
The paper states that experiments and theoretical modelling under viscous dominance yield the reported scaling for critical acceleration. No quoted derivation step reduces the N^{-1/2} ω^{*3/2} form to a fitted parameter, a self-defined quantity, or a load-bearing self-citation. The central claim is presented as following from Navier-Stokes with horizontal body force, viscous-capillary time scale, and free-surface conditions at end-wall waves; this remains an independent scaling argument rather than a tautology or statistical forcing. The derivation chain is therefore self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Viscous forces dominate interfacial dynamics
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
critical non-dimensional acceleration ... scales as N^{-1/2} ω^{*3/2}
-
IndisputableMonolith/Foundation/AlphaCoordinateFixation.leanJ_uniquely_calibrated_via_higher_derivative unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
shear-dominated breakup mechanism ... viscous-capillary scale
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
D. J. McClements, Food Emulsions: Principles, Prac- tices, and Techniques, Second Edition (CRC Press, 2004)
work page 2004
- [2]
-
[3]
T. M. Ho, A. Razzaghi, A. Ramachandran, and K. S. Mikkonen, Adv. Colloid Interface Sci. 299, 102541 (2022)
work page 2022
-
[4]
A. L. Ahmad, A. Kusumastuti, C. J. C. Derek, and B. S. Ooi, J. Chem. Eng. 171, 870 (2011)
work page 2011
-
[5]
A. Kumar, A.and Thakur and P. S. Panesar, Rev. Envi- ron. Sci. Biotechnol. 18, 153 (2019)
work page 2019
- [6]
- [7]
- [8]
-
[9]
S. L. Anna, N. Bontoux, and H. A. Stone, Appl. Phys. Lett. 82, 364 (2003)
work page 2003
-
[10]
A. Evangelio, F. Campo-Cort´ es, and J. M. Gordillo, J. Fluid Mech. 804, 550 (2016)
work page 2016
-
[11]
R. F. Meyer and J. C. Crocker, Phys. Rev. Lett. 102, 194501 (2009)
work page 2009
-
[12]
D. R. Link, S. L. Anna, D. A. Weitz, and H. A. Stone, Phys. Rev. Lett. 92, 054503 (2004)
work page 2004
-
[13]
R. K. Shah, H. C. Shum, A. C. Rowat, D. Lee, J. J. Agresti, A. S. Utada, L. Y. Chu, J. W. Kim, A. Fernandez-Nieves, C. J. Martinez, and D. A. Weitz, Mater. Today 11, 18 (2008)
work page 2008
-
[14]
P. Garstecki, M. J. Fuerstman, H. A. Stone, and G. M. Whitesides, Lab Chip 6, 437 (2006)
work page 2006
-
[15]
M. De Menech, P. Garstecki, F. Jousse, and H. A. Stone, J. Fluid Mech. 595, 141 (2008)
work page 2008
-
[16]
S. L. Anna, Annu. Rev. Fluid Mech. 48, 285 (2016)
work page 2016
-
[17]
M. M. Fryd and T. G. Mason, Annu. Rev. Phys. Chem. 63, 493 (2012)
work page 2012
-
[18]
Q. Chen, N. Singh, K. Schirrmann, Q. Zhou, I. L. Chernyavsky, and A. Juel, Soft Matter 19, 5249 (2023)
work page 2023
-
[19]
C. W. Visser, T. Kamperman, L. P. Karbaat, D. Lohse, and M. Karperien, Sci. Adv. 4, eaao1175 (2018)
work page 2018
-
[20]
Y. Y. Su, D. W. Pan, T. X. Zhang, R. Xie, X. J. Ju, Z. Liu, N. N. Deng, W. Wang, and L. Y. Chu, Sci. Adv. 11, eads1065 (2025)
work page 2025
-
[21]
L. D. Zarzar, V. Sresht, E. M. Sletten, J. A. Kalow, D. Blankschtein, and T. M. Swager, Nature 518, 521 (2015)
work page 2015
-
[22]
E. S. Quintero, A. Evangelio, and J. M. Gordillo, J. Fluid Mech. 846, R3 (2018)
work page 2018
-
[23]
See Supplemental Material at URL-will-be-inserted-by- publisher for movies of regular droplet formation
- [24]
-
[25]
S. V. Jalikop and A. Juel, J. Fluid Mech.640, 131 (2009)
work page 2009
-
[26]
Density and viscosity ( ρu, νu) of silicone oils 10cS (935 kg/m3, 10.3 × 10−6 m2/s), 20cS (950 kg/m3, 21.7 × 10−6 m2/s), 50cS (961 kg /m3, 54.8 × 10−6 m2/s), 100cS (961 kg /m3, 113.7 × 10−6 m2/s), and ( ρl, νl) for HT135 (1752 kg /m3, 1.12 × 10−6 m2/s) and HT270 (1856 kg/m3, 11.7 × 10−6 m2/s)
- [27]
-
[28]
J. Eggers and M. A. Fontelos, in Singularities: Forma- tion, Structure, and Propagation , Vol. 53 (Cambridge University Press, 2015) Chap. 7
work page 2015
-
[29]
J. R. Lister and H. A. Stone, Phys. Fluids 10, 2758 (1998)
work page 1998
-
[30]
C. L. Goodridge, W. T. Shi, H. G. E. Hentschel, and D. P. Lathrop, Phys. Rev. E. 56, 472 (1997)
work page 1997
- [31]
-
[32]
R. J. Lang, J. Acoust. Soc. Am. 34, 6 (1962)
work page 1962
- [33]
- [34]
- [35]
-
[36]
M. Abkarian and H. A. Stone, Phys. Rev. Fluids 5, 102301 (2020). 6 End Matter Horizontal vibration system As illustrated in Fig. A1, the horizontal vibration sys- tem comprises a permanent magnet shaker (LDS-V450), a linear amplifier (LDS, PA1000L), a vibration controller (LDS-COMET USB), a laptop or PC, and a linear air- bearing slide. The container is r...
work page 2020
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.