A hybrid particle-mesh method for viscous, incompressible, multiphase flows

A hybrid method to simulate unsteady multiphase flows in which a sharp interface separates incompressible fluids of different density and viscosity is described. One phase is represented by moving particles and the other phase is defined on stationary mesh. The flow field is discretized by a conservative finite volume approximation on the stationary mesh, and the interface is automatically captured by the distribution of particles moving through the stationary mesh. The effects of surface tension and wall adhesion are evaluated by the continuum surface force model. The different phases are treated as one fluid with variable material properties. Advection of fluid properties such as density and viscosity is done by following the motion of the particles. The method simplifies the calculation of interface interaction, enables accurate modeling of two- and three-dimensional multiphase flows and does not impose any modeling restrictions on the dynamic evolutions of fluid interfaces having surface tension. Several two-dimensional numerical simulations are given to illustrate the efficiency of the hybrid method.

[1]  G. Tryggvason,et al.  A front-tracking method for viscous, incompressible, multi-fluid flows , 1992 .

[2]  S. Koshizuka,et al.  International Journal for Numerical Methods in Fluids Numerical Analysis of Breaking Waves Using the Moving Particle Semi-implicit Method , 2022 .

[3]  Yoshiaki Oka,et al.  A particle method for elastic and visco-plastic structures and fluid-structure interactions , 2001 .

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

[5]  S. Koshizuka,et al.  Moving-Particle Semi-Implicit Method for Fragmentation of Incompressible Fluid , 1996 .

[6]  Yoshiaki Oka,et al.  Numerical analysis of boiling on high heat-flux and high subcooling condition using MPS-MAFL , 2002 .

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

[8]  D. Jacqmin Regular Article: Calculation of Two-Phase Navier–Stokes Flows Using Phase-Field Modeling , 1999 .

[9]  Yoshiaki Oka,et al.  Direct calculation of bubble growth, departure, and rise in nucleate pool boiling , 2001 .

[10]  Xiaolin Li,et al.  Robust Computational Algorithms for Dynamic Interface Tracking in Three Dimensions , 1999, SIAM J. Sci. Comput..

[11]  J. Monaghan,et al.  Kernel estimates as a basis for general particle methods in hydrodynamics , 1982 .

[12]  Jerzy M. Floryan,et al.  Numerical Methods for Viscous Flows With Moving Boundaries , 1989 .

[13]  Elaine S. Oran,et al.  Surface tension and viscosity with Lagrangian hydrodynamics on a triangular mesh. Memorandum report , 1988 .

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

[15]  Koji Okamoto,et al.  Development of computational techniques for nuclear engineering , 1998 .

[16]  G. Dubner,et al.  Study of the Rayleigh-Taylor instability in Tycho's supernova remnant , 1998 .

[17]  Yoshiaki Oka,et al.  Numerical analysis of fragmentation mechanisms in vapor explosions , 1999 .

[18]  F. Harlow,et al.  Numerical Calculation of Time‐Dependent Viscous Incompressible Flow of Fluid with Free Surface , 1965 .

[19]  Qiang Zhang,et al.  Three-Dimensional Front Tracking , 1998, SIAM J. Sci. Comput..

[20]  H. Fernando,et al.  Evolution of Kelvin-Helmholtz billows in nature and laboratory , 1996 .

[21]  Yoshiaki Oka,et al.  Numerical Analysis of Droplet Breakup Behavior using Particle Method , 2001 .

[22]  Seiichi KOSHIZUKA,et al.  Numerical Analysis of Jet Injection Behavior for Fuel-Coolant Interaction using Particle Method , 2001 .