I. Surface Tension and the Kelvin-Helmholtz Instability
Surface tension poses especially severe difficulties in simulating the evolution of free-surface flows. The Laplace-Young (or Gibbs-Thompson) boundary condition introduces high-order, nonlinear and nonlocal terms, which give strong, time-dependent stability constraints on explicit time-stepping schemes.
The figures below are from simulations of the effect of surface tension mixing of immiscible fluids by the classical Kelvin-Helmholtz instability. The numerical method is based upon a boundary integral formulation of the equations of motion and employs the Small-Scale Decomposition (SSD) to remove high-order time-stepping constraints. This has allowed the evolving flow to be computed with very high accuracy, and over much longer times than had been possible previously. Adaptive spatial grids methods have also been developed and employed within this methodology.
Click a box to see a short mpeg clip.
The Weber number (We) measures the relative importance of inertial to surface tension forces; Here We=200, and there are many unstable scales along the dividing interface. Each simulation begins with a nearly flat interface, and typically ends with the disparate parts of the interface colliding, simultaneously forming corners. This is a "topological singularity", and implies the divergence of velocity gradients. This singularity is driven by the surface tension, and presumeably signals an incipient topological transition in the flow (a bubble forms).
These simulations were performed on a Cray C-90 and a dual processor Silicon Graphics Onyx RE2/R8000. Related simulations use the GMRES iterative method for rapid solution of associated integral equations (using SSD preconditioning), and the Fast Multipole Method (FMM) of Greengard & Rokhlin for rapid velocity evaluations. The animations were made on the SGI Onyx using OpenGL-based software.
The SSD has since been employed and extended by ourselves and other researchers to study many interfacial flow problems driven or mediated by a surface tension or energy. These include studies of the Kelvin-Helmholtz and Rayleigh-Taylor instabilities, water waves, and pattern formation and singularity development in Hele-Shaw flows. In the materials context, studies include pattern formation in directional solidification, the symmetric solidification model, Ostwald ripening with many inclusions (5000-10000), and the coarsening of elastic solids
(ii.) These simulations have revealed that a rich variety of behavior arises from a simple model of fluid mixing. For small Weber numbers, there are no unstable length-scales, and the flow is dispersively dominated. For intermediate Weber numbers, where there are only a few unstable length-scales, the interface forms elongating fingers of fluid that interpenetrate each fluid into the other. At larger Weber numbers, with many unstable scales, the interface rolls-up into a ``Kelvin-Helmholtz'' spiral, with its late evolution terminated by the collision of the interface with itself, forming at that instant bubbles of fluid at the core of the spiral. We study carefully this singular event (a ``topological'' or ``pinching'' singularity) using adaptive grid methods. The numerical techniques and simulational results are discussed in
The Long-Time Evolution of Vortex Sheets with Surface
The Physics of Fluids , Vol. 9, p. 1933, 1997.
by T.Y. Hou (Caltech), J.S. Lowengrub (Minnesota), and M.J. Shelley
Download gzipped postscript PDF
This paper is cited in our being awarded the 1998 Frenkiel
Award by the Division of Fluid Dynamics of the American Physical Society.
See also our review:
Boundary Integral Methods for Multicomponent Materials and Multiphase Fluids,
(w. T. Hou and J. Lowengrub)
Journal of Computational Physics 169, p. 302 (2001).