Convergence of a non-stiff boundary integral method for interfacial flows with surface tension
by Héctor D. Ceniceros and Thomas Y. Hou PDF
Math. Comp. 67 (1998), 137-182 Request permission
Boundary integral methods to simulate interfacial flows are very sensitive to numerical instabilities. In addition, surface tension introduces nonlinear terms with high order spatial derivatives into the interface dynamics. This makes the spatial discretization even more difficult and, at the same time, imposes a severe time step constraint for stable explicit time integration methods. A proof of the convergence of a reformulated boundary integral method for two-density fluid interfaces with surface tension is presented. The method is based on a scheme introduced by Hou, Lowengrub and Shelley [ J. Comp. Phys. 114 (1994), pp. 312–338] to remove the high order stability constraint or stiffness. Some numerical filtering is applied carefully at certain places in the discretization to guarantee stability. The key of the proof is to identify the most singular terms of the method and to show, through energy estimates, that these terms balance one another. The analysis is at a time continuous-space discrete level but a fully discrete case for a simple Hele-Shaw interface is also studied. The time discrete analysis shows that the high order stiffness is removed and also provides an estimate of how the CFL constraint depends on the curvature and regularity of the solution. The robustness of the method is illustrated with several numerical examples. A numerical simulation of an unstably stratified two-density interfacial flow shows the roll-up of the interface; the computations proceed up to a time where the interface is about to pinch off and trapped bubbles of fluid are formed. The method remains stable even in the full nonlinear regime of motion. Another application of the method shows the process of drop formation in a falling single fluid.
Uri M. Ascher, Steven J. Ruuth, and Brian T. R. Wetton, Implicit-explicit methods for time-dependent partial differential equations, SIAM J. Numer. Anal. 32 (1995), no. 3, 797–823. MR 1335656, DOI 10.1137/0732037
G. Baker, Generalized vortex methods for free-surface flows, Waves on Fluid Interfaces (R. Meyer, ed.), Univ. Wisc. Press, 1983, pp. 53–81.
G. R. Baker and M. J. Shelley, On the connection between thin vortex layers and vortex sheets, J. Fluid Mech. 215 (1990), 161–194. MR 1061499, DOI 10.1017/S0022112090002609
Gregory Baker and André Nachbin, Stable methods for vortex sheet motion in presence of surface tension, Submitted to J. Comp. Phys.
Gregory R. Baker, Daniel I. Meiron, and Steven A. Orszag, Vortex simulations of the Rayleigh-Taylor instability, Phys. Fluids 23 (1980), 1485–1490.
Gregory R. Baker, Daniel I. Meiron, and Steven A. Orszag, Generalized vortex methods for free-surface flow problems, J. Fluid Mech. 123 (1982), 477–501. MR 687014, DOI 10.1017/S0022112082003164
J. T. Beale, Thomas Y. Hou, and John S. Lowengrub, Convergence of a boundary integral method for water waves, SIAM J. Num. Anal. 33 (1996), 1797–1843.
—, The stability of two-fluid flows with surface tension, part 1: Growth rates for the linear motion far from equilibrium; part 2: Convergence of suitably modified vortex methods, In preparation.
J. Thomas Beale, Thomas Y. Hou, and John S. Lowengrub, Growth rates for the linearized motion of fluid interfaces away from equilibrium, Comm. Pure Appl. Math. 46 (1993), no. 9, 1269–1301. MR 1231428, DOI 10.1002/cpa.3160460903
—, On the well-posedness of two fluid interfacial flows with surface tension, Singularities in Fluids, Plasmas and Optics (London) (R. C. Caflisch and G. C. Papanicolaou, eds.), Kluwer Academic, 1993.
J. T. Beale, T. Y. Hou, J. S. Lowengrub, and M. J. Shelley, Spatial and temporal stability issues for interfacial flows with surface tension, Math. Comput. Modelling 20 (1994), no. 10-11, 1–27. Theory and numerical methods for initial-boundary value problems. MR 1306284, DOI 10.1016/0895-7177(94)90167-8
J. B. Bell and D. L. Marcus, A second-order projection method for variable-density flows, J. Comp. Phys. 101 (1992), 334–348.
J. U. Brackbill, D. B. Kothe, and C. Zemach, A continuum method for modeling surface tension, J. Comput. Phys. 100 (1992), no. 2, 335–354. MR 1167749, DOI 10.1016/0021-9991(92)90240-Y
H. D. Ceniceros and T. Y. Hou, Numerical studies of surface tension effects in interfacial flows, In preparation.
Héctor D. Ceniceros, Convergence of a reformulated boundary integral method for two fluid interfaces with surface tension, Ph.D. thesis, New York University, May 1995.
Y. C. Chang, T. Y. Hou, B. Merriman, and S. Osher, A level set formulation of Eulerian interface capturing methods for incompressible fluid flows, J. Comput. Phys. 124 (1996), no. 2, 449–464. MR 1383769, DOI 10.1006/jcph.1996.0072
W. S. Dai and M. J. Shelley, A numerical study of the effect of surface tension and noise on an expanding Hele-Shaw bubble, Phys. Fluids A 5 (1993), no. 9, 2131–2146.
Bart J. Daly, Numerical study of the effect of surface tension on interface instability, Phys. Fluids 12 (1969), no. 7, 1340–1354.
J. W. Dold, An efficient surface-integral algorithm applied to unsteady gravity waves, J. Comput. Phys. 103 (1992), no. 1, 90–115. MR 1188090, DOI 10.1016/0021-9991(92)90327-U
P. G. Drazin and William Hill Reid, Hydrodynamic stability, Cambridge Monographs on Mechanics and Applied Mathematics, Cambridge University Press, Cambridge-New York, 1981. MR 604359
Raymond E. Goldstein and Dean M. Petrich, The Korteweg-de Vries hierarchy as dynamics of closed curves in the plane, Phys. Rev. Lett. 67 (1991), no. 23, 3203–3206. MR 1135964, DOI 10.1103/PhysRevLett.67.3203
Günther Hämmerlin and Karl-Heinz Hoffmann, Numerical mathematics, Undergraduate Texts in Mathematics, Springer-Verlag, New York, 1991. Translated from the German by Larry Schumaker; Readings in Mathematics. MR 1088482, DOI 10.1007/978-1-4612-4442-4
Thomas Y. Hou, John S. Lowengrub, and Michael J. Shelley, The roll-up and self-intersection of vortex sheets under surface tension, Preprint. Courant Institute.
Thomas Y. Hou, John S. Lowengrub, and Michael J. Shelley, Removing the stiffness from interfacial flows with surface tension, J. Comput. Phys. 114 (1994), no. 2, 312–338. MR 1294935, DOI 10.1006/jcph.1994.1170
D.A. Kessler, J. Koplik, and H. Levine, Geometrical models of interface evolution. II, Phys. Rev. A 30 (1984), 3161–3174.
Robert Krasny, Desingularization of periodic vortex sheet roll-up, J. Comp. Phys. 65 (1986), 292–313.
Robert Krasny, A study of singularity formation in a vortex sheet by the point-vortex approximation, J. Fluid Mech. 167 (1986), 65–93. MR 851670, DOI 10.1017/S0022112086002732
L. D. Landau and E. M. Lifshitz, Fluid mechanics, Course of Theoretical Physics, Vol. 6, Pergamon Press, London-Paris-Frankfurt; Addison-Wesley Publishing Company, Inc., Reading, Mass., 1959. Translated from the Russian by J. B. Sykes and W. H. Reid. MR 0108121
J. S. Langer, Instabilities and pattern formation in crystal growth, Rev. Modern Phys. 52 (1980), 1–28.
M. S. Longuet-Higgins and E. D. Cokelet, The deformation of steep surface waves on water. I. A numerical method of computation, Proc. Roy. Soc. London Ser. A 350 (1976), no. 1660, 1–26. MR 411355, DOI 10.1098/rspa.1976.0092
A. I. Markushevich, Theory of functions of a complex variable. Vol. I, II, III, Second English edition, Chelsea Publishing Co., New York, 1977. Translated and edited by Richard A. Silverman. MR 0444912
E. Meiburg and G. M. Homsy, Nonlinear unstable viscous fingers in Hele-Shaw flows .2. Numerical simulation, Phys. Fluids 31 (1988), no. 3, 429–439.
Stanley Osher and James A. Sethian, Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations, J. Comput. Phys. 79 (1988), no. 1, 12–49. MR 965860, DOI 10.1016/0021-9991(88)90002-2
C. Pozrikidis, Boundary integral and singularity methods for linearized viscous flow, Cambridge Texts in Applied Mathematics, Cambridge University Press, Cambridge, 1992. MR 1156495, DOI 10.1017/CBO9780511624124
D. I. Pullin, Numerical studies of surface-tension effects in nonlinear Kelvin-Helmholtz and Rayleigh-Taylor instability, J. Fluid Mech. 119 (1982), 507–532.
R.H. Rangel and W.A. Sirignano, Nonlinear growth of Kelvin-Helmholtz instability: Effect of surface tension and density ratio, Phys. Fluids 31 (1988), no. 7, 1845–1855.
A. J. Roberts, A stable and accurate numerical method to calculate the motion of a sharp interface between fluids, IMA J. Appl. Math. 31 (1983), 13–35.
M. J. Shelley, A study of singularity formation in vortex-sheet motion by a spectrally accurate vortex method, J. Fluid Mech. 244 (1992), 493–526. MR 1192255, DOI 10.1017/S0022112092003161
Avram Sidi and Moshe Israeli, Quadrature methods for periodic singular and weakly singular Fredholm integral equations, J. Sci. Comput. 3 (1988), no. 2, 201–231. MR 981760, DOI 10.1007/BF01061258
John Strain, A boundary integral approach to unstable solidification, J. Comput. Phys. 85 (1989), no. 2, 342–389. MR 1029193, DOI 10.1016/0021-9991(89)90155-1
Gilbert Strang, Accurate partial difference methods. II. Non-linear problems, Numer. Math. 6 (1964), 37–46. MR 166942, DOI 10.1007/BF01386051
M. Sussman, P. Smereka, and S. Osher, A level set approach for computing solutions to incompressible 2-phase flow, J. Comp. Phys. 114(1) (1994), 146–159.
Eitan Tadmor, Stability analysis of finite difference, pseudospectral and Fourier-Galerkin approximations for time-dependent problems, SIAM Rev. 29 (1987), no. 4, 525–555. MR 917694, DOI 10.1137/1029110
Grétar Tryggvason, Numerical simulations of the Rayleigh-Taylor instability, J. Comp. Phys. 75 (1988), 253–282.
Retrieve articles in Mathematics of Computation with MSC (1991): 65M12, 76B15
Retrieve articles in all journals with MSC (1991): 65M12, 76B15
Héctor D. Ceniceros
Affiliation: Department of Applied Mathematics, California Institute of Technology, Pasadena, California 91125
Address at time of publication: Centro de Investigación en Computación Instituto Politécnico Nacional Col. Lindavista, Mexico City, Mexico 07300.
Email: hdc@jsbach.cic.ipn.mx
Thomas Y. Hou
Email: hou@ama.caltech.edu
Received by editor(s): December 7, 1995
Received by editor(s) in revised form: June 5, 1996
Additional Notes: The first author was partially supported by the Office of Naval Research under Grant N00014-94-1-0310.
The second author was partially supported by the Office of Naval Research under Grant N00014-94-1-0310 and the National Science Foundation under grant DMS-9407030.
Journal: Math. Comp. 67 (1998), 137-182
MSC (1991): Primary 65M12, 76B15
MathSciNet review: 1443116