An improved unified solver for compressible and incompressible fluids involving free surfaces. II. Multi-time-step integration and applications

Abstract An improved numerical solver for the unified solution of compressible and incompressible fluids involving interfaces is proposed. The present method is based on the CIP–CUP (Cubic Interpolated Propagation/Combined, Unified Procedure) method, which is a pressure-based semi-implicit solver for the Euler equations of fluid flows. In Part I of this series of articles [M. Ida, Comput. Phys. Commun. 132 (2000) 44], we proposed an improved scheme for the convection terms in the equations, which allowed us discontinuous descriptions of the density interface by replacing the cubic interpolation function used in the CIP scheme with a quadratic extrapolation function only around the interface. In this paper, as Part II of this series, the multi-time-step integration technique is adapted to the CIP–CUP integration. Because the CIP–CUP treats different-nature components in the fluid equations separately, the adaptation of the technique is straightforward. This modification allows us flexible determinations of the time interval, which results in an efficient and accurate integration. Furthermore, some additional discussion on our methods is presented. Finally, the application results to composite flow problems such as compressible and incompressible Kelvin–Helmholtz instabilities and the dynamics of two acoustically coupled deformable bubbles in a viscous liquid are provided.

[1]  De-Kang Mao.  A treatment of discontinuities in shock-capturing finite difference methods , 1991 .

[2]  O Louisnard,et al.  High bubble concentrations produced by ultrasounds in binary mixtures. , 2001, Ultrasonics sonochemistry.

[3]  R. Fedkiw,et al.  A numerical method for two-phase flow consisting of separate compressible and incompressible regions , 2000 .

[4]  Takashi Yabe,et al.  Unified Numerical Procedure for Compressible and Incompressible Fluid , 1991 .

[5]  A. A. Amsden,et al.  A numerical fluid dynamics calculation method for all flow speeds , 1971 .

[6]  Junichiro Makino,et al.  A Modified Aarseth Code for GRAPE and Vector Processors , 1991 .

[7]  Claus-Dieter Ohl,et al.  Experimental and Theoretical Bubble Dynamics , 2007 .

[8]  J. Scott,et al.  Singular perturbation theory applied to the collective oscillation of gas bubbles in a liquid , 1981, Journal of Fluid Mechanics.

[9]  R. Klein.  Semi-implicit extension of a Godunov-type scheme based on low Mach number asymptotics , 1995 .

[10]  Masato Ida.  An improved unified solver for compressible and incompressible fluids involving free surfaces. Part I. Convection , 2000 .

[11]  Lawrence A. Crum,et al.  Bjerknes forces on bubbles in a stationary sound field , 1975 .

[12]  R. Gleim.  Legal notice , 1989 .

[13]  Feng Xiao,et al.  Description of Complex and Sharp Interface during Shock Wave Interaction with Liquid Drop , 1993 .

[14]  Mark E. Tuckerman,et al.  Reversible multiple time scale molecular dynamics , 1992 .

[15]  John R. Blake,et al.  Interaction of two cavitation bubbles with a rigid boundary , 1993, Journal of Fluid Mechanics.

[16]  Yabe,et al.  Two- and three-dimensional behavior of Rayleigh-Taylor and Kelvin-Helmholtz instabilities. , 1991, Physical review. A, Atomic, molecular, and optical physics.

[17]  Andrea Prosperetti,et al.  Bubble phenomena in sound fields: part one , 1984 .

[18]  Qiang Zhang,et al.  A numerical study of bubble interactions in Rayleigh–Taylor instability for compressible fluids , 1990 .

[19]  R. D. Jackson,et al.  Heat Transfer 1 , 1965 .

[20]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[21]  A. Shima.  The Natural Frequencies of Two Spherical Bubbles Oscillating in Water , 1971 .

[22]  P Charrier,et al.  On front-tracking methods applied to hyperbolic systems of nonlinear conservation laws , 1986 .

[23]  Takashi Yabe,et al.  A universal solver for hyperbolic equations by cubic-polynomial interpolation I. One-dimensional solver , 1991 .

[24]  D. Youngs,et al.  Three-dimensional numerical simulation of turbulent mixing by Rayleigh-Taylor instability , 1991 .

[25]  Feng Xiao,et al.  Effect of EOS on Break-Up of Shoemaker-Levy 9 Entering Jovian Atmosphere , 1994 .

[26]  Edward D Harder,et al.  Efficient multiple time step method for use with Ewald and particle mesh Ewald for large biomolecular systems , 2001 .

[27]  Ali Nadim,et al.  Coupled pulsation and translation of two gas bubbles in a liquid , 2001, Journal of Fluid Mechanics.

[28]  Takashi Yabe,et al.  The compact CIP (cubic-interpolated pseudo-particle) method as a general hyperbolic solver , 1991 .

[29]  Suhas V. Patankar,et al.  A NEW FINITE-DIFFERENCE SCHEME FOR PARABOLIC DIFFERENTIAL EQUATIONS , 1978 .

[30]  Feng Xiao,et al.  Regular Article: A Computational Model for Suspended Large Rigid Bodies in 3D Unsteady Viscous Flows , 1999 .

[31]  A. Doinikov,et al.  Translational motion of two interacting bubbles in a strong acoustic field. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[32]  M. Rudman.  A Volume-Tracking Method for Incompressible Multifluid Flows with Large Density Variations , 1998 .

[33]  Feng Xiao,et al.  Constructing oscillation preventing scheme for advection equation by rational function , 1996 .

[34]  D. Sharp.  An overview of Rayleigh-Taylor instability☆ , 1984 .

[35]  C. Feuillade,et al.  Acoustically coupled gas bubbles in fluids: time-domain phenomena. , 2001, The Journal of the Acoustical Society of America.

[36]  Takashi Yabe,et al.  The unified simulation for incompressible and compressible flow by the predictor-corrector scheme based on the CIP method , 1999 .

[37]  O. C. Zienkiewicz,et al.  A general explicit or semi-explicit algorithm for compressible and incompressible flows , 1992 .

[38]  A. Gosman,et al.  Solution of the implicitly discretised reacting flow equations by operator-splitting , 1986 .

[39]  Winston C. Chao.  Formulation of an Explicit-Multiple-Time-Step Time Integration Method for Use in a Global Primitive Equation Grid Model , 1982 .

[40]  Masato Ida,et al.  An Eulerian Scheme for Direct Numerical Simulation of Multibubble Dynamics in an Acoustic Field , 2001 .

[41]  J. Brackbill,et al.  A continuum method for modeling surface tension , 1992 .

[42]  S. J. Aarseth,et al.  Dynamical Evolution of Clusters Of Galaxies, II , 1963 .

[43]  Masato Ida.  A conservative semi-Lagrangian method for oscillation-free computation of advection processes , 2001 .

[44]  J. Haile,et al.  A multiple time-step method for molecular dynamics simulations of fluids of chain molecules , 1984 .

[45]  S. Osher,et al.  An improved level set method for incompressible two-phase flows , 1998 .

[46]  Marilyn Schneider,et al.  LARGE AND SMALL SCALE STRUCTURE IN RAYLEIGH-TAYLOR MIXING , 1998 .

[47]  M. Strasberg.  The Pulsation Frequency of Nonspherical Gas Bubbles in Liquids , 1953 .

[48]  Zhen Ye,et al.  ACOUSTIC LOCALIZATION IN BUBBLY LIQUID MEDIA , 1998 .

[49]  K. Nightingale,et al.  A preliminary evaluation of the effects of primary and secondary radiation forces on acoustic contrast agents , 1997, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[50]  M. Rabaud,et al.  Gap size effects for the Kelvin-Helmholtz instability in a Hele-Shaw cell. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[51]  C. W. Hirt,et al.  Volume of fluid (VOF) method for the dynamics of free boundaries , 1981 .

[52]  S. K. Zhdanov.  Nonlinear theory of Kelvin-Helmholtz instability , 1995 .

[53]  Rupert Klein,et al.  Regular Article: Extension of Finite Volume Compressible Flow Solvers to Multi-dimensional, Variable Density Zero Mach Number Flows , 1999 .

[54]  T. Yabe,et al.  The constrained interpolation profile method for multiphase analysis , 2001 .

[55]  S. Osher,et al.  Computing interface motion in compressible gas dynamics , 1992 .

[56]  Feng Xiao.  A Class of Single-Cell High-Order Semi-Lagrangian Advection Schemes , 2000 .