REMOVING THE STIFFNESS FROM INTERFACIAL FLOW WITH SURFACE-TENSION

被引:441
作者
HOU, TY
LOWENGRUB, JS
SHELLEY, MJ
机构
[1] UNIV MINNESOTA,SCH MATH,MINNEAPOLIS,MN 55455
[2] NYU,COURANT INST MATH SCI,NEW YORK,NY 10012
关键词
D O I
10.1006/jcph.1994.1170
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
A new formulation and new methods are presented for computing the motion of fluid interfaces with surface tension in two-dimensional, irrotational, and incompressible fluids. Through the Laplace-Young condition at the interface, surface tension introduces high-order terms, both nonlinear and nonlocal, into the dynamics. This leads to severe stability constraints for explicit time integration methods and makes the application of implicit methods difficult. This new formulation has all the nice properties for time integration methods that are associated with having a linear, constant coefficient, highest order term. That is, using this formulation, we give implicit time integration methods that have no high order time step stability constraint associated with surface tension and are explicit in Fourier space. The approach is based on a boundary integral formulation and applies more generally, even to problems beyond the fluid mechanical context. Here they are applied to computing with high resolution the motion of interfaces in Hele-Shaw flows and the motion of free surfaces in inviscid flows governed by the Euler equations. One Hele-Shaw computation shows the behavior of an expanding gas bubble over long-time as the interface undergoes successive tip-splittings and finger competition. A second computation shows the formation of a very ramified interface through the interaction of surface tension with an unstable density stratification. In Euler flows, the computation of a vortex sheet shows its roll-up through the Kelvin-Helmholtz instability. This motion culminates in the late time self-intersection of the interface, creating trapped bubbles of fluid. This is, we believe, a type of singularity formation previously unobserved for such flows in 2D. Finally, computations of falling plumes in an unstably stratified Boussinesq fluid show a very similar behavior. (C) 1994 Academic Press, Inc.
引用
收藏
页码:312 / 338
页数:27
相关论文
共 51 条
[1]   SCALING BEHAVIOR IN ANISOTROPIC HELE-SHAW FLOW [J].
ALMGREN, R ;
DAI, WS ;
HAKIM, V .
PHYSICAL REVIEW LETTERS, 1993, 71 (21) :3461-3464
[2]  
Ames WF, 1977, NUMERICAL METHODS PA, V2nd
[3]   A METHOD OF LOCAL CORRECTIONS FOR COMPUTING THE VELOCITY-FIELD DUE TO A DISTRIBUTION OF VORTEX BLOBS [J].
ANDERSON, CR .
JOURNAL OF COMPUTATIONAL PHYSICS, 1986, 62 (01) :111-123
[4]  
BAKER G, PREPRINT
[5]   GENERALIZED VORTEX METHODS FOR FREE-SURFACE FLOW PROBLEMS [J].
BAKER, GR ;
MEIRON, DI ;
ORSZAG, SA .
JOURNAL OF FLUID MECHANICS, 1982, 123 (OCT) :477-501
[6]   ON THE CONNECTION BETWEEN THIN VORTEX LAYERS AND VORTEX SHEETS [J].
BAKER, GR ;
SHELLEY, MJ .
JOURNAL OF FLUID MECHANICS, 1990, 215 :161-194
[7]   BOUNDARY INTEGRAL TECHNIQUES FOR MULTICONNECTED DOMAINS [J].
BAKER, GR ;
SHELLEY, MJ .
JOURNAL OF COMPUTATIONAL PHYSICS, 1986, 64 (01) :112-132
[8]  
BAKER GR, 1983, WAVES FLUID INTERFAC
[9]  
BEALE JT, UNPUB
[10]  
BEALE JT, IN PRESS COMMUN PURE