Capturing the Behavior of Bubbles and Drops Using the Variational Level Set Approach

We reproduce the general behavior of complicated bubble and droplet motions using the variational level set formulation introduced by the authors earlier. Our approach here ignores inertial effects; thus the motion is only correct as an approximation for very viscous problems. However, the steady states are true equilibrium solutions. Inertial forces will be added in future work. The problems include: soap bubbles colliding and merging, drops falling or remaining attached to a (generally irregular) ceiling, and liquid penetrating through a funnel in both two and three dimensions. Each phase is identified with a particular “level set” function. The zero level set of this function is that particular phase boundary. The level set functions all evolve in time through a constrained gradient descent procedure so as to minimize an energy functional. The functions are coupled through physical constraints and through the requirements that different phases do not overlap and vacuum regions do not develop. Both boundary conditions and inequality constraints are cast in terms of (either local or global) equality constraints. The gradient projection method leads to a system of perturbed (by curvature, if surface tension is involved) Hamilton?Jacobi equations coupled through a constraint. The coupling is enforced using the Lagrange multiplier associated with this constraint. The numerical implementation requires much of the modern level set technology; in particular, we achieve a significant speed up by using the fast localization algorithm of H.-K. Zhao, M. Kang, B. Merriman, D. Peng, and S. Osher.

[1]  S. Osher,et al.  Regular Article: A PDE-Based Fast Local Level Set Method , 1999 .

[2]  Stanley Osher,et al.  Regularization of Ill-Posed Problems Via the Level Set Approach , 1998, SIAM J. Appl. Math..

[3]  Steven J. Ruuth.  Efficient Algorithms for Diffusion-Generated Motion by Mean Curvature , 1998 .

[4]  J. Sethian.  METHODS FOR PROPAGATING INTERFACES , 1998 .

[5]  S. Osher,et al.  Regular Article: A PDE-Based Fast Local Level Set Method , 1999 .

[6]  S. Osher,et al.  A Simple Level Set Method for Solving Stefan Problems , 1997, Journal of Computational Physics.

[7]  Steven J. Ruuth.  Eecient Algorithms for Diiusion-generated Motion by Mean Curvature , 1997 .

[8]  T. Chan,et al.  A Variational Level Set Approach to Multiphase Motion , 1996 .

[9]  Stanley Osher,et al.  An Eulerian Approach for Vortex Motion Using a Level Set Regularization Procedure , 1996 .

[10]  S. Osher,et al.  A Level Set Formulation of Eulerian Interface Capturing Methods for Incompressible Fluid Flows , 1996 .

[11]  J. Sethian,et al.  A Fast Level Set Method for Propagating Interfaces , 1995 .

[12]  Stanley Osher,et al.  Subscale Capturing in Numerical Analysis , 1995 .

[13]  J. Hass,et al.  The double bubble conjecture , 1995 .

[14]  S. Osher,et al.  A level set approach for computing solutions to incompressible two-phase flow , 1994 .

[15]  M. Brenner,et al.  A Cascade of Structure in a Drop Falling from a Faucet , 1994, Science.

[16]  S. Osher,et al.  Motion of multiple junctions: a level set approach , 1994 .

[17]  D. Chopp.  Computing Minimal Surfaces via Level Set Curvature Flow , 1993 .

[18]  S. Osher,et al.  High-order essentially nonsocillatory schemes for Hamilton-Jacobi equations , 1990 .

[19]  ShuChi-Wang,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes, II , 1989 .

[20]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[21]  S. Osher,et al.  Algorithms Based on Hamilton-Jacobi Formulations , 1988 .

[22]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[23]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .

[24]  S. Osher,et al.  Uniformly high order accuracy essentially non-oscillatory schemes III , 1987 .

[25]  S. Osher,et al.  Some results on uniformly high-order accurate essentially nonoscillatory schemes , 1986 .

[26]  M. Scott.  An initial value method for the eigenvalue problem for systems of ordinary differential equations , 1973 .

[27]  J. B. Rosen.  The gradient projection method for nonlinear programming: Part II , 1961 .