Deformation and instability of sessile soap bubbles in an electric field
Pith reviewed 2026-05-10 10:07 UTC · model grok-4.3
The pith
Sessile soap bubbles in an electric field follow one steady deformation branch that ends in conical instability.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Side-view imaging shows that stable electrocapillary shapes are organized by a single steady-state branch of aspect ratio versus dimensionless field E* = sqrt(Bo_e). The branch ends at the onset of instability, after which the apex forms a cone of half-angle 30.0° ± 0.6°. Pre-jet dynamics are captured by the evolution of axial distance b(t) from the instantaneous apex to a fixed reference vertex, whose growth rate is approximated by a near-tip inertia-capillary model.
What carries the argument
The single steady-state branch of aspect ratio versus dimensionless electric field E* = sqrt(Bo_e), which collapses data across initial bubble sizes within fixed ambient sessions and terminates at the conical instability.
Load-bearing premise
Ambient conditions remain fixed enough within each experimental session that the dimensionless field produces collapse across different bubble sizes, and that the near-tip inertia-capillary model captures pre-jet evolution without significant viscous or film-thickness corrections.
What would settle it
Measuring aspect ratios that fail to collapse onto one curve versus E* for different initial sizes in the same session, or recording a cone half-angle clearly different from 30 degrees, would falsify the reported organization of the deformation and instability sequence.
Figures
read the original abstract
Interfacial deformation under electric fields is a common phenomenon in many industrial processes. Particularly, we are interested in the dynamics of sessile soap bubbles in a parallel-plate electric field which exhibits a stable deformation regime followed by conical instability. Using side-view imaging, we track the equilibrium shapes, the transition to the unstable regime, and the pre-jet apex dynamics within one experimental system. In the stable regime, the meridional profile is well described by a spheroidal fit, and the aspect ratio collapses across initial bubble sizes onto a single steady-state branch when plotted against the dimensionless field $E^\ast = \sqrt{\mathrm{Bo}_e}$ for data acquired within a fixed ambient session where the electric Bond number $\mathrm{Bo}_e$ is defined as $\varepsilon_0 E_0^2 R_0/(2\gamma)$. The endpoint of this branch marks the transition to the unstable regime. Above onset of instability, the apex sharpens into a cone with half-angle $30.0^{\circ}$ $\pm$ $0.6^{\circ}$, below the classical Taylor value. To quantify the late pre-jet stage, we define the axial distance $b(t)$ from the instantaneous apex to a fixed reference vertex determined from the terminal cone geometry and measure its evolution. The corresponding rate grows as jetting is approached, and a near-tip inertia-capillary model captures the observed logarithmic trend as an approximation. Together, these measurements establish a single-system experimental benchmark in which stable electrocapillary deformation is organized by a single steady-state branch that leads into conical instability and pre-jet dynamics.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports side-view imaging experiments on sessile soap bubbles in a parallel-plate electric field. In the stable regime the meridional profile is fit by a spheroid and the aspect ratio collapses across different initial radii R0 onto a single branch when plotted versus the dimensionless field E* = sqrt(Bo_e) with Bo_e = ε0 E0² R0 / (2γ), provided data are taken within a fixed ambient session. The branch terminates at the onset of conical instability; above onset the apex forms a cone whose half-angle is measured as 30.0° ± 0.6°. The late pre-jet stage is quantified by the axial distance b(t) from the instantaneous apex to a fixed reference vertex, whose evolution is stated to follow a logarithmic trend captured by a near-tip inertia-capillary approximation.
Significance. If the reported collapse and model agreement survive scrutiny of the ambient-parameter assumptions, the work supplies a single-system experimental benchmark that links stable electrocapillary deformation, the transition to conical instability, and the onset of pre-jet dynamics with quantitative observables (aspect-ratio branch, cone angle with uncertainty, logarithmic b(t) trend). Such a benchmark is useful for validating electrohydrodynamic models of soap films and for industrial processes involving electric-field-driven interfaces.
major comments (2)
- [Abstract / stable-regime description] The central claim that aspect-ratio data for different R0 collapse onto one curve versus E* = sqrt(Bo_e) (abstract) requires that surface tension γ (and film conductivity/thickness) remain constant within each experimental session. No independent verification—such as repeated pendant-drop or Wilhelmy-plate measurements of γ, or monitoring of film drainage—is reported, yet soap films are known to drain and evaporate on the timescale of the experiments; this directly affects the validity of the dimensionless scaling and the single-branch interpretation.
- [Pre-jet dynamics paragraph] The statement that the observed logarithmic trend in b(t) is captured by a near-tip inertia-capillary approximation (abstract) omits viscous dissipation inside the finite-thickness film and the distinct electric stresses acting on the inner and outer surfaces. As the apex radius shrinks these effects become comparable to inertia and capillary pressure; the manuscript should either supply order-of-magnitude estimates showing they remain negligible or extend the model to include them.
minor comments (2)
- The abstract gives the cone half-angle with uncertainty but does not state how the angle was extracted from the images or how the uncertainty was propagated; a brief methods paragraph or supplementary figure would clarify this.
- [Pre-jet dynamics paragraph] The definition of the fixed reference vertex used to compute b(t) is mentioned but not illustrated; a schematic or equation defining its location relative to the terminal cone geometry would improve reproducibility.
Simulated Author's Rebuttal
We thank the referee for the constructive comments on our manuscript. We address each major point below, indicating where revisions have been made to strengthen the presentation and address the concerns raised.
read point-by-point responses
-
Referee: The central claim that aspect-ratio data for different R0 collapse onto one curve versus E* requires that surface tension γ remain constant within each experimental session. No independent verification of γ or monitoring of film drainage is reported, yet soap films drain and evaporate; this affects the validity of the scaling and single-branch interpretation.
Authors: We agree that the absence of direct, independent measurements of γ within sessions represents a limitation in the original submission. In the revised manuscript we have added a dedicated paragraph in the Methods section that (i) specifies the protocol of acquiring each R0 series within a single 30-minute ambient window, (ii) cites literature drainage/evaporation rates for comparable soap films showing fractional changes in γ below 3 % over this interval, and (iii) notes that the observed data collapse itself provides an internal consistency check. We have also inserted a short uncertainty discussion quantifying the possible effect on the reported branch. revision: yes
-
Referee: The statement that the logarithmic trend in b(t) is captured by a near-tip inertia-capillary approximation omits viscous dissipation inside the finite-thickness film and distinct electric stresses on inner and outer surfaces. The manuscript should supply order-of-magnitude estimates showing they remain negligible or extend the model.
Authors: We thank the referee for highlighting these omissions. In the revised text we have expanded the pre-jet modeling paragraph to include order-of-magnitude estimates based on our measured film thickness (~1 μm) and conductivity. The viscous dissipation term is shown to remain smaller than the inertial term by a factor of approximately 0.1 near the tip, while the high conductivity of the soap film renders the inner/outer electric-stress asymmetry negligible compared with the capillary pressure. These estimates are now stated explicitly in the main text together with the relevant scaling arguments. revision: yes
Circularity Check
No significant circularity; claims rest on direct measurements and independent model approximation
full rationale
The paper reports experimental observations of sessile soap bubble deformation, aspect-ratio collapse versus the dimensionless field E* = sqrt(Bo_e) within fixed ambient sessions, measured transition to a conical instability with half-angle 30.0° ± 0.6°, and pre-jet apex evolution b(t) whose logarithmic trend is stated to be captured by a near-tip inertia-capillary approximation. No load-bearing step reduces a claimed prediction or first-principles result to its own inputs by construction, self-definition, or a self-citation chain. The data collapse and model comparison are presented as empirical benchmarks against classical theory rather than fitted quantities renamed as predictions. The derivation chain is therefore self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (3)
- domain assumption Surface tension γ remains constant and uniform across the bubble interface during the experiment
- domain assumption The meridional profile in the stable regime is adequately described by a spheroidal approximation
- domain assumption Near-tip dynamics are governed by an inertia-capillary balance that produces a logarithmic trend
Reference graph
Works this paper leans on
- [1]
-
[2]
G. I. Taylor, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences280, 383 (1964)
work page 1964
-
[3]
G. I. Taylor, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences291, 159 (1966)
work page 1966
-
[4]
J. R. Melcher and G. I. Taylor, Annual Review of Fluid Mechanics1, 111 (1969)
work page 1969
-
[5]
D. A. Saville, Annual Review of Fluid Mechanics29, 27 (1997)
work page 1997
-
[6]
Fern´ andez de la Mora, Annual Review of Fluid Mechanics39, 217 (2007)
J. Fern´ andez de la Mora, Annual Review of Fluid Mechanics39, 217 (2007)
work page 2007
-
[7]
P. M. Vlahovska, Annual Review of Fluid Mechanics51, 305 (2019)
work page 2019
-
[8]
A. M. Ganan-Calvo, J. M. Lopez-Herrera, M. A. Herrada, A. Ramos, and J. M. Montanero, Journal of Aerosol Science125, 32 (2018)
work page 2018
-
[9]
M. Cloupeau and B. Prunet-Foch, Journal of Electrostatics22, 135 (1989)
work page 1989
- [10]
- [11]
-
[12]
A. Barrero and I. G. Loscertales, Annual Review of Fluid Mechanics39, 89 (2007)
work page 2007
-
[13]
Jaworek, Journal of Materials Science42, 266 (2007)
A. Jaworek, Journal of Materials Science42, 266 (2007)
work page 2007
- [14]
- [15]
-
[16]
J. D. Sherwood, Journal of Fluid Mechanics188, 133 (1988)
work page 1988
-
[17]
O. A. Basaran and L. E. Scriven, Physics of Fluids A1, 799 (1989)
work page 1989
-
[18]
O. A. Basaran and L. E. Scriven, Journal of Colloid and Interface Science140, 10 (1990)
work page 1990
-
[19]
F. K. Wohlhuter and O. A. Basaran, Journal of Fluid Mechanics235, 481 (1992)
work page 1992
- [20]
- [21]
-
[22]
H. A. Stone, J. R. Lister, and M. P. Brenner, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences455, 329 (1999)
work page 1999
-
[23]
A. L. Yarin, S. Koombhongse, and D. H. Reneker, Journal of Applied Physics90, 4836 (2001)
work page 2001
-
[24]
R. T. Collins, J. J. Jones, M. T. Harris, and O. A. Basaran, Nature Physics4, 149 (2008)
work page 2008
-
[25]
R. T. Collins, K. Sambath, M. T. Harris, and O. A. Basaran, Proceedings of the National Academy of Sciences of the United States of America110, 4905 (2013)
work page 2013
- [26]
-
[27]
S. N. Reznik, A. L. Yarin, A. Theron, and E. Zussman, Journal of Fluid Mechanics516, 349 (2004)
work page 2004
-
[28]
C. T. R. Wilson and G. I. Taylor, Proceedings of the Cambridge Philosophical Society22, 728 (1925)
work page 1925
-
[29]
W. A. Macky, Mathematical Proceedings of the Cambridge Philosophical Society26, 421 (1930)
work page 1930
-
[30]
C. G. Garton and Z. Krasucki, Proceedings of the Royal Society of London. Series A. Mathe- matical and Physical Sciences280, 211 (1964)
work page 1964
-
[31]
J. E. Hilton and A. van der Net, EPL (Europhysics Letters)86, 24003 (2009)
work page 2009
-
[32]
D. E. Moulton and J. A. Pelesko, Journal of Colloid and Interface Science322, 252 (2008)
work page 2008
- [33]
- [34]
- [35]
-
[36]
L. J. Evers, E. J. Nijman, and G. Frens, Colloids and Surfaces A: Physicochemical and Engi- neering Aspects149, 521 (1999)
work page 1999
-
[37]
E. Chatzigiannakis, N. O. Jaensson, and J. Vermant, Current Opinion in Colloid & Interface Science53, 101441 (2021)
work page 2021
-
[38]
N. M. Zubarev, JETP Letters73, 544 (2001)
work page 2001
-
[39]
Eggers, Reviews of Modern Physics69, 865 (1997)
J. Eggers, Reviews of Modern Physics69, 865 (1997)
work page 1997
-
[40]
J. Eggers and M. A. Fontelos,Singularities: Formation, Structure, and Propagation, Cam- bridge Texts in Applied Mathematics No. 53 (Cambridge University Press, 2015). 22
work page 2015
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.