Nonlinear oscillations of inviscid free drops

Tadeusz Patzek*, R. E. Benner, O. A. Basaran, L. E. Scriven

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

53 Scopus citations


Inviscid oscillations of free liquid drops are analyzed by solving Bernoulli's equation for the free surface shape and Laplace's equation for the velocity potential field. The means are: (a) Galerkin's weighted residual method which converts the governing equations into a large system of nonlinear, time-dependent ordinary differential equations; (b) an implicit predictor-corrector method for time integration which automatically adjusts time steps; and (c) Newton's method which solves the large system of nonlinear algebraic equations that results from time discretization. Results presented include sequences of drop shapes, pressure distributions, particle paths, and evolution with time of kinetic and surface energies. Accuracy is attested by virtual constancy of drop volume and total energy and smallness of mass and momentum fluxes across drop surfaces. Dynamic response to small amplitude disturbances agrees with linear theory. Large-amplitude oscillations are compared to the predictions by the marker- and-cell method and second-order perturbation theory. Mode interactions and frequency shifts are analyzed by Fourier power spectra and lend further insight into the nature of the oscillations.

Original languageEnglish (US)
Pages (from-to)489-515
Number of pages27
JournalJournal of Computational Physics
Issue number2
StatePublished - Jan 1 1991

ASJC Scopus subject areas

  • Numerical Analysis
  • Modeling and Simulation
  • Physics and Astronomy (miscellaneous)
  • Physics and Astronomy(all)
  • Computer Science Applications
  • Computational Mathematics
  • Applied Mathematics


Dive into the research topics of 'Nonlinear oscillations of inviscid free drops'. Together they form a unique fingerprint.

Cite this